Splicing factor proline and glutamine rich intron retention, reduced expression and aggregate formation are pathological features of amyotrophic lateral sclerosis.

AIM: Splicing factor proline and glutamine rich (SFPQ) is an RNA-DNA binding protein that is dysregulated in Alzheimer's disease and frontotemporal dementia. Dysregulation of SFPQ, specifically increased intron retention and nuclear depletion, has been linked to several genetic subtypes of amyotrophic lateral sclerosis (ALS), suggesting that SFPQ pathology may be a common feature of this heterogeneous disease. Our study aimed to investigate this hypothesis by providing the first comprehensive assessment of SFPQ pathology in large ALS case-control cohorts. METHODS: We examined SFPQ at the RNA, protein and DNA levels. SFPQ RNA expression and intron retention were examined using RNA-sequencing and quantitative PCR. SFPQ protein expression was assessed by immunoblotting and immunofluorescent staining. At the DNA level, SFPQ was examined for genetic variation novel to ALS patients. RESULTS: At the RNA level, retention of SFPQ intron nine was significantly increased in ALS patients' motor cortex. In addition, SFPQ RNA expression was significantly reduced in the central nervous system, but not blood, of patients. At the protein level, neither nuclear depletion nor reduced expression of SFPQ was found to be a consistent feature of spinal motor neurons. However, SFPQ-positive ubiquitinated protein aggregates were observed in patients' spinal motor neurons. At the DNA level, our genetic screen identified two novel and two rare SFPQ sequence variants not previously reported in the literature. CONCLUSIONS: Our findings confirm dysregulation of SFPQ as a pathological feature of the central nervous system of ALS patients and indicate that investigation of the functional consequences of this pathology will provide insight into ALS biology.

Genome-wide microRNA profiling in brain and blood samples in a mouse model of epileptogenesis.

OBJECTIVES: This study profiled circulating and hippocampal microRNAs (miRNAs) to identify alterations associated with the risk of epileptogenesis in a mouse temporal lobe epilepsy model. METHODS: Next-generation sequencing was performed to examine the changes in miRNA expression 24 h after pilocarpine-induced status epilepticus (SE) in C57BL/6NCrl mice using both blood and hippocampus samples. Differentially expressed miRNAs were identified from SE animals and matched controls that failed to develop SE after receiving equal doses of pilocarpine (NS animals). Blood and brain miRNA profiles were then compared to identify circulating miRNA alterations reflecting the changes in the brain. RESULTS: We identified 3 miRNAs that were significantly up-regulated and 4 miRNAs that were significantly down-regulated in the blood of SE animals compared with NS animals. When hippocampal miRNAs of SE animals and NS animals were compared, 5 miRNAs were up-regulated and 4 were down-regulated. Of these, miR-434-3p and miR-133a-3p were observed to have greatest changes in both blood and brain of SE animals. SIGNIFICANCE: This study extends current knowledge of changes in miRNAs associated with epileptogenesis by profiling miRNAs in SE and NS animals in an experimental temporal lobe epilepsy model. The study was designed to allow non-specific changes due to the activation of muscarinic cholinergic receptors in peripheral organs by pilocarpine to be ruled out. Significantly altered circulating miRNAs that reflect changes in the brain during epileptogenesis after SE have the potential to be developed as prognostic biomarkers for epileptogenesis.

Monozygotic twins and triplets discordant for amyotrophic lateral sclerosis display differential methylation and gene expression.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterised by the loss of upper and lower motor neurons. ALS exhibits high phenotypic variability including age and site of onset, and disease duration. To uncover epigenetic and transcriptomic factors that may modify an ALS phenotype, we used a cohort of Australian monozygotic twins (n = 3 pairs) and triplets (n = 1 set) that are discordant for ALS and represent sporadic ALS and the two most common types of familial ALS, linked to C9orf72 and SOD1. Illumina Infinium HumanMethylation450K BeadChip, EpiTYPER and RNA-Seq analyses in these ALS-discordant twins/triplets and control twins (n = 2 pairs), implicated genes with consistent longitudinal differential DNA methylation and/or gene expression. Two identified genes, RAD9B and C8orf46, showed significant differential methylation in an extended cohort of >1000 ALS cases and controls. Combined longitudinal methylation-transcription analysis within a single twin set implicated CCNF, DPP6, RAMP3, and CCS, which have been previously associated with ALS. Longitudinal transcriptome data showed an 8-fold enrichment of immune function genes and under-representation of transcription and protein modification genes in ALS. Examination of these changes in a large Australian sporadic ALS cohort suggest a broader role in ALS. Furthermore, we observe that increased methylation age is a signature of ALS in older patients.

Bioorthogonal Metabolic Labeling of Nascent RNA in Neurons Improves the Sensitivity of Transcriptome-Wide Profiling.

Transcriptome-wide expression profiling of neurons has provided important insights into the underlying molecular mechanisms and gene expression patterns that transpire during learning and memory formation. However, there is a paucity of tools for profiling stimulus-induced RNA within specific neuronal cell populations. A bioorthogonal method to chemically label nascent (i.e., newly transcribed) RNA in a cell-type-specific and temporally controlled manner, which is also amenable to bioconjugation via click chemistry, was recently developed and optimized within conventional immortalized cell lines. However, its value within a more fragile and complicated cellular system such as neurons, as well as for transcriptome-wide expression profiling, has yet to be demonstrated. Here, we report the visualization and sequencing of activity-dependent nascent RNA derived from neurons using this labeling method. This work has important implications for improving transcriptome-wide expression profiling and visualization of nascent RNA in neurons, which has the potential to provide valuable insights into the mechanisms underlying neural plasticity, learning, and memory.

Whole-exome sequencing in amyotrophic lateral sclerosis suggests NEK1 is a risk gene in Chinese.

BACKGROUND: Amyotrophic lateral sclerosis (ALS) is a progressive neurological disease characterised by the degeneration of motor neurons, which are responsible for voluntary movement. There remains limited understanding of disease aetiology, with median survival of ALS of three years and no effective treatment. Identifying genes that contribute to ALS susceptibility is an important step towards understanding aetiology. The vast majority of published human genetic studies, including for ALS, have used samples of European ancestry. The importance of trans-ethnic studies in human genetic studies is widely recognised, yet a dearth of studies of non-European ancestries remains. Here, we report analyses of novel whole-exome sequencing (WES) data from Chinese ALS and control individuals. METHODS: WES data were generated for 610 ALS cases and 460 controls drawn from Chinese populations. We assessed evidence for an excess of rare damaging mutations at the gene level and the gene set level, considering only singleton variants filtered to have allele frequency less than 5 × 10-5 in reference databases. To meta-analyse our results with a published study of European ancestry, we used a Cochran-Mantel-Haenszel test to compare gene-level variant counts in cases vs controls. RESULTS: No gene passed the genome-wide significance threshold with ALS in Chinese samples alone. Combining rare variant counts in Chinese with those from the largest WES study of European ancestry resulted in three genes surpassing genome-wide significance: TBK1 (p = 8.3 × 10-12), SOD1 (p = 8.9 × 10-9) and NEK1 (p = 1.1 × 10-9). In the Chinese data alone, SOD1 and NEK1 were nominally significantly associated with ALS (p = 0.04 and p = 7 × 10-3, respectively) and the case/control frequencies of rare coding variants in these genes were similar in Chinese and Europeans (SOD1: 1.5%/0.2% vs 0.9%/0.1%, NEK1 1.8%/0.4% vs 1.9%/0.8%). This was also true for TBK1 (1.2%/0.2% vs 1.4%/0.4%), but the association with ALS in Chinese was not significant (p = 0.14). CONCLUSIONS: While SOD1 is already recognised as an ALS-associated gene in Chinese, we provide novel evidence for association of NEK1 with ALS in Chinese, reporting variants in these genes not previously found in Europeans.

Whole exome sequencing and DNA methylation analysis in a clinical amyotrophic lateral sclerosis cohort.

BACKGROUND: Gene discovery has provided remarkable biological insights into amyotrophic lateral sclerosis (ALS). One challenge for clinical application of genetic testing is critical evaluation of the significance of reported variants. METHODS: We use whole exome sequencing (WES) to develop a clinically relevant approach to identify a subset of ALS patients harboring likely pathogenic mutations. In parallel, we assess if DNA methylation can be used to screen for pathogenicity of novel variants since a methylation signature has been shown to associate with the pathogenic C9orf72 expansion, but has not been explored for other ALS mutations. Australian patients identified with ALS-relevant variants were cross-checked with population databases and case reports to critically assess whether they were "likely causal," "uncertain significance," or "unlikely causal." RESULTS: Published ALS variants were identified in >10% of patients; however, in only 3% of patients (4/120) could these be confidently considered pathogenic (in SOD1 and TARDBP). We found no evidence for a differential DNA methylation signature in these mutation carriers. CONCLUSIONS: The use of WES in a typical ALS clinic demonstrates a critical approach to variant assessment with the capability to combine cohorts to enhance the largely unknown genetic basis of ALS.

Small non-coding RNA expression from anterior cingulate cortex in schizophrenia shows sex specific regulation.

MicroRNAs (miRNAs) are known to regulate the expression of genes that are important for brain development and function, but the roles of other classes of small non-coding RNAs (sncRNAs) are less well understood. Additionally, although miRNA expression studies have been conducted in post-mortem brain samples from schizophrenia (SCZ) patients, other classes of sncRNAs are yet to be investigated in SCZ. We profiled the expression of miRNAs, piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs) and small nuclear RNAs (snRNAs) in SCZ by applying small RNA sequencing (RNA-Seq) to sncRNA isolated from post-mortem anterior cingulate cortex (ACC) of SCZ-affected individuals (n=22) and matched controls (n=22). We identified about one-third of annotated miRNAs, one-quarter of snoRNAs and a small proportion of piRNAs and snRNAs. No sncRNAs were significantly differentially expressed between SCZ and controls, but there was evidence for an interaction between disease status and sex on the expression level of a number of miRNAs and snoRNAs. Many of these transcripts exhibited differential expression between male and female cases, and/or between female cases and controls, suggesting sex based dysregulation in ACC of SCZ. These findings require replication in an independent sample, but our study provides further insights into the potential involvement of sncRNAs in brain function and SCZ.

Purification of neural precursor cells reveals the presence of distinct, stimulus-specific subpopulations of quiescent precursors in the adult mouse hippocampus.

The activity of neural precursor cells in the adult hippocampus is regulated by various stimuli; however, whether these stimuli regulate the same or different precursor populations remains unknown. Here, we developed a novel cell-sorting protocol that allows the purification to homogeneity of neurosphere-forming neural precursors from the adult mouse hippocampus and examined the responsiveness of individual precursors to various stimuli using a clonal assay. We show that within the Hes5-GFP(+)/Nestin-GFP(+)/EGFR(+) cell population, which comprises the majority of neurosphere-forming precursors, there are two distinct subpopulations of quiescent precursor cells, one directly activated by high-KCl depolarization, and the other activated by norepinephrine (NE). We then demonstrate that these two populations are differentially distributed along the septotemporal axis of the hippocampus, and show that the NE-responsive precursors are selectively regulated by GABA, whereas the KCl-responsive precursors are selectively modulated by corticosterone. Finally, based on RNAseq analysis by deep sequencing, we show that the progeny generated by activating NE-responsive versus KCl-responsive quiescent precursors are molecularly different. These results demonstrate that the adult hippocampus contains phenotypically similar but stimulus-specific populations of quiescent precursors, which may give rise to neural progeny with different functional capacity.

A comparative study of techniques for differential expression analysis on RNA-Seq data.

Recent advances in next-generation sequencing technology allow high-throughput cDNA sequencing (RNA-Seq) to be widely applied in transcriptomic studies, in particular for detecting differentially expressed genes between groups. Many software packages have been developed for the identification of differentially expressed genes (DEGs) between treatment groups based on RNA-Seq data. However, there is a lack of consensus on how to approach an optimal study design and choice of suitable software for the analysis. In this comparative study we evaluate the performance of three of the most frequently used software tools: Cufflinks-Cuffdiff2, DESeq and edgeR. A number of important parameters of RNA-Seq technology were taken into consideration, including the number of replicates, sequencing depth, and balanced vs. unbalanced sequencing depth within and between groups. We benchmarked results relative to sets of DEGs identified through either quantitative RT-PCR or microarray. We observed that edgeR performs slightly better than DESeq and Cuffdiff2 in terms of the ability to uncover true positives. Overall, DESeq or taking the intersection of DEGs from two or more tools is recommended if the number of false positives is a major concern in the study. In other circumstances, edgeR is slightly preferable for differential expression analysis at the expense of potentially introducing more false positives.

Cube - an online tool for comparison and contrasting of protein sequences.

When comparing sequences of similar proteins, two kinds of questions can be asked, and the related two kinds of inference made. First, one may ask to what degree they are similar, and then, how they differ. In the first case one may tentatively conclude that the conserved elements common to all sequences are of central and common importance to the protein's function. In the latter case the regions of specialization may be discriminative of the function or binding partners across subfamilies of related proteins. Experimental efforts - mutagenesis or pharmacological intervention - can then be pointed in either direction, depending on the context of the study. Cube simplifies this process for users that already have their favorite sets of sequences, and helps them collate the information by visualization of the conservation and specialization scores on the sequence and on the structure, and by spreadsheet tabulation. All information can be visualized on the spot, or downloaded for reference and later inspection. SERVER HOMEPAGE: http://eopsf.org/cube.

Global transcriptome profiles of Camellia sinensis during cold acclimation.

BACKGROUND: Tea is the most popular non-alcoholic health beverage in the world. The tea plant (Camellia sinensis (L.) O. Kuntze) needs to undergo a cold acclimation process to enhance its freezing tolerance in winter. Changes that occur at the molecular level in response to low temperatures are poorly understood in tea plants. To elucidate the molecular mechanisms of cold acclimation, we employed RNA-Seq and digital gene expression (DGE) technologies to the study of genome-wide expression profiles during cold acclimation in tea plants. RESULTS: Using the Illumina sequencing platform, we obtained approximately 57.35 million RNA-Seq reads. These reads were assembled into 216,831 transcripts, with an average length of 356 bp and an N50 of 529 bp. In total, 1,770 differentially expressed transcripts were identified, of which 1,168 were up-regulated and 602 down-regulated. These include a group of cold sensor or signal transduction genes, cold-responsive transcription factor genes, plasma membrane stabilization related genes, osmosensing-responsive genes, and detoxification enzyme genes. DGE and quantitative RT-PCR analysis further confirmed the results from RNA-Seq analysis. Pathway analysis indicated that the "carbohydrate metabolism pathway" and the "calcium signaling pathway" might play a vital role in tea plants' responses to cold stress. CONCLUSIONS: Our study presents a global survey of transcriptome profiles of tea plants in response to low, non-freezing temperatures and yields insights into the molecular mechanisms of tea plants during the cold acclimation process. It could also serve as a valuable resource for relevant research on cold-tolerance and help to explore the cold-related genes in improving the understanding of low-temperature tolerance and plant-environment interactions.

Cube-DB: detection of functional divergence in human protein families.

Cube-DB is a database of pre-evaluated results for detection of functional divergence in human/vertebrate protein families. The analysis is organized around the nomenclature associated with the human proteins, but based on all currently available vertebrate genomes. Using full genomes enables us, through a mutual-best-hit strategy, to construct comparable taxonomical samples for all paralogues under consideration. Functional specialization is scored on the residue level according to two models of behavior after divergence: heterotachy and homotachy. In the first case, the positions on the protein sequence are scored highly if they are conserved in the reference group of orthologs, and overlap poorly with the residue type choice in the paralogs groups (such positions will also be termed functional determinants). The second model additionally requires conservation within each group of paralogs (functional discriminants). The scoring functions are phylogeny independent, but sensitive to the residue type similarity. The results are presented as a table of per-residue scores, and mapped onto related structure (when available) via browser-embedded visualization tool. They can also be downloaded as a spreadsheet table, and sessions for two additional molecular visualization tools. The database interface is available at http://epsf.bmad.bii.a-star.edu.sg/cube/db/html/home.html.

Determinants, discriminants, conserved residues--a heuristic approach to detection of functional divergence in protein families.

In this work, belonging to the field of comparative analysis of protein sequences, we focus on detection of functional specialization on the residue level. As the input, we take a set of sequences divided into groups of orthologues, each group known to be responsible for a different function. This provides two independent pieces of information: within group conservation and overlap in amino acid type across groups. We build our discussion around the set of scoring functions that keep the two separated and the source of the signal easy to trace back to its source.We propose a heuristic description of functional divergence that includes residue type exchangeability, both in the conservation and in the overlap measure, and does not make any assumptions on the rate of evolution in the groups other than the one under consideration. Residue types acceptable at a certain position within an orthologous group are described as a distribution which evolves in time, starting from a single ancestral type, and is subject to constraints that can be inferred only indirectly. To estimate the strength of the constraints, we compare the observed degrees of conservation and overlap with those expected in the hypothetical case of a freely evolving distribution.Our description matches the experiment well, but we also conclude that any attempt to capture the evolutionary behavior of specificity determining residues in terms of a scalar function will be tentative, because no single model can cover the variety of evolutionary behavior such residues exhibit. Especially, models expecting the same type of evolutionary behavior across functionally divergent groups tend to miss a portion of information otherwise retrievable by the conservation and overlap measures they use.

WLS-dependent secretion of WNT3A requires Ser209 acylation and vacuolar acidification.

Wnt proteins are secreted post-translationally modified proteins that signal locally to regulate development and proliferation. The production of bioactive Wnts requires a number of dedicated factors in the secreting cell whose coordinated functions are not fully understood. A screen for small molecules identified inhibitors of vacuolar acidification as potent inhibitors of Wnt secretion. Inhibition of the V-ATPase or disruption of vacuolar pH gradients by diverse drugs potently inhibited Wnt/ß-catenin signaling both in cultured human cells and in vivo, and impaired Wnt-regulated convergent extension movements in Xenopus embryos. WNT secretion requires its binding to the carrier protein wntless (WLS); we find that WLS is ER-resident in human cells and WNT3A binding to WLS requires PORCN-dependent lipid modification of WNT3A at serine 209. Inhibition of vacuolar acidification results in accumulation of the WNT3A-WLS complex both in cells and at the plasma membrane. Modeling predictions suggest that WLS has a lipid-binding ß-barrel that is similar to the lipocalin-family fold. We propose that WLS binds Wnts in part through a lipid-binding domain, and that vacuolar acidification is required to release palmitoylated WNT3A from WLS in secretory vesicles, possibly to facilitate transfer of WNT3A to a soluble carrier protein.

deconSTRUCT: general purpose protein database search on the substructure level.

deconSTRUCT webserver offers an interface to a protein database search engine, usable for a general purpose detection of similar protein (sub)structures. Initially, it deconstructs the query structure into its secondary structure elements (SSEs) and reassembles the match to the target by requiring a (tunable) degree of similarity in the direction and sequential order of SSEs. Hierarchical organization and judicious use of the information about protein structure enables deconSTRUCT to achieve the sensitivity and specificity of the established search engines at orders of magnitude increased speed, without tying up irretrievably the substructure information in the form of a hash. In a post-processing step, a match on the level of the backbone atoms is constructed. The results presented to the user consist of the list of the matched SSEs, the transformation matrix for rigid superposition of the structures and several ways of visualization, both downloadable and implemented as a web-browser plug-in. The server is available at http://epsf.bmad.bii.a-star.edu.sg/struct_server.html.

Reduced representation of protein structure: implications on efficiency and scope of detection of structural similarity.

BACKGROUND: Computational comparison of two protein structures is the starting point of many methods that build on existing knowledge, such as structure modeling (including modeling of protein complexes and conformational changes), molecular replacement, or annotation by structural similarity. In a commonly used strategy, significant effort is invested in matching two sets of atoms. In a complementary approach, a global descriptor is assigned to the overall structure, thus losing track of the substructures within. RESULTS: Using a small set of geometric features, we define a reduced representation of protein structure, together with an optimizing function for matching two representations, to provide a pre-filtering stage in a database search. We show that, in a straightforward implementation, the representation performs well in terms of resolution in the space of protein structures, and its ability to make new predictions. CONCLUSIONS: Perhaps unexpectedly, a substantial discriminating power already exists at the level of main features of protein structure, such as directions of secondary structural elements, possibly constrained by their sequential order. This can be used toward efficient comparison of protein (sub)structures, allowing for various degrees of conformational flexibility within the compared pair, which in turn can be used for modeling by homology of protein structure and dynamics.

In silico characterization of immunogenic epitopes presented by HLA-Cw*0401.

BACKGROUND: HLA-C locus products are poorly understood in part due to their low expression at the cell surface. Recent data indicate that these molecules serve as major restriction elements for human immunodeficiency virus type 1 (HIV-1) cytotoxic T lymphocyte (CTL) epitopes. We report here a structure-based technique for the prediction of peptides binding to Cw*0401. The models were rigorously trained, tested and validated using experimentally verified Cw*0401 binding and non-binding peptides obtained from biochemical studies. A new scoring scheme facilitates the identification of immunological hot spots within antigens, based on the sum of predicted binding energies of the top four binders within a window of 30 amino acids. RESULTS: High predictivity is achieved when tested on the training (r(2) = 0.88, s = 3.56 kJ/mol, q(2) = 0.84, s(press) = 5.18 kJ/mol) and test (A(ROC) = 0.93) datasets. Characterization of the predicted Cw*0401 binding sequences indicate that amino acids at key anchor positions share common physico-chemical properties which correlate well with existing experimental studies. CONCLUSION: The analysis of predicted Cw*0401-binding peptides showed that anchor residues may not be restrictive and the Cw*0401 binding pockets may possibly accommodate a wide variety of peptides with common physico-chemical properties. The potential Cw*0401-specific T-cell epitope repertoires for HIV-1 p24(gag) and gp160(gag) glycoproteins are well distributed throughout both glycoproteins, with thirteen and nine immunological hot spots for HIV-1 p24(gag) and gp160(gag) glycoproteins respectively. These findings provide new insights into HLA-C peptide selectivity, indicating that pre-selection of candidate HLA-C peptides may occur at the TAP level, prior to peptide loading in the endoplasmic reticulum.

AllerTool: a web server for predicting allergenicity and allergic cross-reactivity in proteins.

UNLABELLED: Assessment of potential allergenicity and patterns of cross-reactivity is necessary whenever novel proteins are introduced into human food chain. Current bioinformatic methods in allergology focus mainly on the prediction of allergenic proteins, with no information on cross-reactivity patterns among known allergens. In this study, we present AllerTool, a web server with essential tools for the assessment of predicted as well as published cross-reactivity patterns of allergens. The analysis tools include graphical representation of allergen cross-reactivity information; a local sequence comparison tool that displays information of known cross-reactive allergens; a sequence similarity search tool for assessment of cross-reactivity in accordance to FAO/WHO Codex alimentarius guidelines; and a method based on support vector machine (SVM). A 10-fold cross-validation results showed that the area under the receiver operating curve (A(ROC)) of SVM models is 0.90 with 86.00% sensitivity (SE) at specificity (SP) of 86.00%. AVAILABILITY: AllerTool is freely available at http://research.i2r.a-star.edu.sg/AllerTool/.

ALLERDB database and integrated bioinformatic tools for assessment of allergenicity and allergic cross-reactivity.

Databases and computational tools are increasingly important in the study of allergies, particularly in the assessment of allergenicity and allergic cross-reactivity. ALLERDB database contains sequences of allergens and information on reported cross-reactivity between allergens. It focuses on analysis of allergenicity and allergic cross-reactivity of clinically relevant protein allergens. The official IUIS allergen data were extracted from the IUIS Allergen Nomenclature Sub-Committee website, and their sequence information from the public databases, and reference publications. The analysis tools assist allergen data analysis and retrieval, and include keyword searching, BLAST, prediction of allergenicity, modification of BLAST that displays cross-reactive allergens, and graphics representation of cross-reactivity data. ALLERDB is new brand of allergen databases with a rich set of tools for sequence comparison, pattern identification, and visualization of results. It is accessible at http://research.i2r.a-star.edu.sg/Templar/DB/Allergen.