ALBERT

Beschreibung: Motivation: Mapping of high-throughput sequencing data and other bulk sequence comparison applications have motivated a search for high-efficiency sequence alignment algorithms. The bit-parallel approach represents individual cells in an alignment scoring matrix as bits in computer words and emulates the calculation of scores by a series of logic operations composed of AND, OR, XOR, complement, shift and addition. Bit-parallelism has been successfully applied to the longest common subsequence (LCS) and edit-distance problems, producing fast algorithms in practice. Results: We have developed BitPAl, a bit-parallel algorithm for general, integer-scoring global alignment. Integer-scoring schemes assign integer weights for match, mismatch and insertion/deletion. The BitPAl method uses structural properties in the relationship between adjacent scores in the scoring matrix to construct classes of efficient algorithms, each designed for a particular set of weights. In timed tests, we show that BitPAl runs 7–25 times faster than a standard iterative algorithm. Availability and implementation: Source code is freely available for download at http://lobstah.bu.edu/BitPAl/BitPAl.html . BitPAl is implemented in C and runs on all major operating systems. Contact : jloving@bu.edu or yhernand@bu.edu or gbenson@bu.edu Supplementary information : Supplementary data are available at Bioinformatics online.

Beschreibung: : Next-generation sequencing (NGS) has a large potential in HIV diagnostics, and genotypic prediction models have been developed and successfully tested in the recent years. However, albeit being highly accurate, these computational models lack computational efficiency to reach their full potential. In this study, we demonstrate the use of graphics processing units (GPUs) in combination with a computational prediction model for HIV tropism. Our new model named gCUP, parallelized and optimized for GPU, is highly accurate and can classify 〉175 000 sequences per second on an NVIDIA GeForce GTX 460. The computational efficiency of our new model is the next step to enable NGS technologies to reach clinical significance in HIV diagnostics. Moreover, our approach is not limited to HIV tropism prediction, but can also be easily adapted to other settings, e.g. drug resistance prediction. Availability and implementation: The source code can be downloaded at http://www.heiderlab.de Contact: d.heider@wz-straubing.de

Beschreibung: : We present a new method to incrementally construct the FM-index for both short and long sequence reads, up to the size of a genome. It is the first algorithm that can build the index while implicitly sorting the sequences in the reverse (complement) lexicographical order without a separate sorting step. The implementation is among the fastest for indexing short reads and the only one that practically works for reads of averaged kilobases in length. Availability and implementation: https://github.com/lh3/ropebwt2 Contact: hengli@broadinstitute.org

Beschreibung: : AliView is an alignment viewer and editor designed to meet the requirements of next-generation sequencing era phylogenetic datasets. AliView handles alignments of unlimited size in the formats most commonly used, i.e. FASTA, Phylip, Nexus, Clustal and MSF. The intuitive graphical interface makes it easy to inspect, sort, delete, merge and realign sequences as part of the manual filtering process of large datasets. AliView also works as an easy-to-use alignment editor for small as well as large datasets. Availability and implementation: AliView is released as open-source software under the GNU General Public License, version 3.0 (GPLv3), and is available at GitHub ( www.github.com/AliView ). The program is cross-platform and extensively tested on Linux, Mac OS X and Windows systems. Downloads and help are available at http://ormbunkar.se/aliview Contact: anders.larsson@ebc.uu.se Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: The ability to accurately read the order of nucleotides in DNA and RNA is fundamental for modern biology. Errors in next-generation sequencing can lead to many artifacts, from erroneous genome assemblies to mistaken inferences about RNA editing. Uneven coverage in datasets also contributes to false corrections. Result: We introduce Trowel, a massively parallelized and highly efficient error correction module for Illumina read data. Trowel both corrects erroneous base calls and boosts base qualities based on the k -mer spectrum. With high-quality k -mers and relevant base information, Trowel achieves high accuracy for different short read sequencing applications.The latency in the data path has been significantly reduced because of efficient data access and data structures. In performance evaluations, Trowel was highly competitive with other tools regardless of coverage, genome size read length and fragment size. Availability and implementation: Trowel is written in C++ and is provided under the General Public License v3.0 (GPLv3). It is available at http://trowel-ec.sourceforge.net . Contact: euncheon.lim@tue.mpg.de or weigel@tue.mpg.de Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: : The application of protein–protein docking in large-scale interactome analysis is a major challenge in structural bioinformatics and requires huge computing resources. In this work, we present MEGADOCK 4.0, an FFT-based docking software that makes extensive use of recent heterogeneous supercomputers and shows powerful, scalable performance of 〉97% strong scaling. Availability and Implementation: MEGADOCK 4.0 is written in C++ with OpenMPI and NVIDIA CUDA 5.0 (or later) and is freely available to all academic and non-profit users at: http://www.bi.cs.titech.ac.jp/megadock . Contact: akiyama@cs.titech.ac.jp Supplementary information: Supplementary data are available at Bioinformatics online

Beschreibung: Motivation: The identification of active transcriptional regulatory elements is crucial to understand regulatory networks driving cellular processes such as cell development and the onset of diseases. It has recently been shown that chromatin structure information, such as DNase I hypersensitivity (DHS) or histone modifications, significantly improves cell-specific predictions of transcription factor binding sites. However, no method has so far successfully combined both DHS and histone modification data to perform active binding site prediction. Results: We propose here a method based on hidden Markov models to integrate DHS and histone modifications occupancy for the detection of open chromatin regions and active binding sites. We have created a framework that includes treatment of genomic signals, model training and genome-wide application. In a comparative analysis, our method obtained a good trade-off between sensitivity versus specificity and superior area under the curve statistics than competing methods. Moreover, our technique does not require further training or sequence information to generate binding location predictions. Therefore, the method can be easily applied on new cell types and allow flexible downstream analysis such as de novo motif finding. Availability and implementation: Our framework is available as part of the Regulatory Genomics Toolbox. The software information and all benchmarking data are available at http://costalab.org/wp/dh-hmm . Contact: ivan.costa@rwth-aachen.de or eduardo.gusmao@rwth-aachen.de Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: A proper target or marker is essential in any diagnosis (e.g. an infection or cancer). An ideal diagnostic target should be both conserved in and unique to the pathogen. Currently, these targets can only be identified manually, which is time-consuming and usually error-prone. Because of the increasingly frequent occurrences of emerging epidemics and multidrug-resistant ‘superbugs’, a rapid diagnostic target identification process is needed. Results: A new method that can identify uniquely conserved regions (UCRs) as candidate diagnostic targets for a selected group of organisms solely from their genomic sequences has been developed and successfully tested. Using a sequence-indexing algorithm to identify UCRs and a k -mer integer-mapping model for computational efficiency, this method has successfully identified UCRs within the bacteria domain for 15 test groups, including pathogenic, probiotic, commensal and extremophilic bacterial species or strains. Based on the identified UCRs, new diagnostic primer sets were designed, and their specificity and efficiency were tested by polymerase chain reaction amplifications from both pure isolates and samples containing mixed cultures. Availability and implementation: The UCRs identified for the 15 bacterial species are now freely available at http://ucr.synblex.com . The source code of the programs used in this study is accessible at http://ucr.synblex.com/bacterialIdSourceCode.d.zip Contact: yazhousun@synblex.com Supplementary Information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: A popular method for classification of protein domain movements apportions them into two main types: those with a ‘hinge’ mechanism and those with a ‘shear’ mechanism. The intuitive assignment of domain movements to these classes has limited the number of domain movements that can be classified in this way. Furthermore, whether intended or not, the term ‘shear’ is often interpreted to mean a relative translation of the domains. Results: Numbers of occurrences of four different types of residue contact changes between domains were optimally combined by logistic regression using the training set of domain movements intuitively classified as hinge and shear to produce a predictor for hinge and shear. This predictor was applied to give a 10-fold increase in the number of examples over the number previously available with a high degree of precision. It is shown that overall a relative translation of domains is rare, and that there is no difference between hinge and shear mechanisms in this respect. However, the shear set contains significantly more examples of domains having a relative twisting movement than the hinge set. The angle of rotation is also shown to be a good discriminator between the two mechanisms. Availability and implementation: Results are free to browse at http://www.cmp.uea.ac.uk/dyndom/interface/ . Contact: sjh@cmp.uea.ac.uk . Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Recent studies on human disease have revealed that aberrant interaction between proteins probably underlies a substantial number of human genetic diseases. This suggests a need to investigate disease inheritance mode using interaction, and based on which to refresh our conceptual understanding of a series of properties regarding inheritance mode of human disease. Results: We observed a strong correlation between the number of protein interactions and the likelihood of a gene causing any dominant diseases or multiple dominant diseases, whereas no correlation was observed between protein interaction and the likelihood of a gene causing recessive diseases. We found that dominant diseases are more likely to be associated with disruption of important interactions. These suggest inheritance mode should be understood using protein interaction. We therefore reviewed the previous studies and refined an interaction model of inheritance mode, and then confirmed that this model is largely reasonable using new evidences. With these findings, we found that the inheritance mode of human genetic diseases can be predicted using protein interaction. By integrating the systems biology perspectives with the classical disease genetics paradigm, our study provides some new insights into genotype–phenotype correlations. Contact: haodapeng@ems.hrbmu.edu.cn or biofomeng@hotmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: : Recently, several high profile studies collected cell viability data from panels of cancer cell lines treated with many drugs applied at different concentrations. Such drug sensitivity data for cancer cell lines provide suggestive treatments for different types and subtypes of cancer. Visualization of these datasets can reveal patterns that may not be obvious by examining the data without such efforts. Here we introduce Drug/Cell-line Browser (DCB), an online interactive HTML5 data visualization tool for interacting with three of the recently published datasets of cancer cell lines/drug-viability studies. DCB uses clustering and canvas visualization of the drugs and the cell lines, as well as a bar graph that summarizes drug effectiveness for the tissue of origin or the cancer subtypes for single or multiple drugs. DCB can help in understanding drug response patterns and prioritizing drug/cancer cell line interactions by tissue of origin or cancer subtype. Availability and implementation: DCB is an open source Web-based tool that is freely available at: http://www.maayanlab.net/LINCS/DCB Contact: avi.maayan@mssm.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Clinical and neuropathological similarities between dementia with Lewy bodies (DLB), Parkinson's and Alzheimer's diseases (PD and AD, respectively) suggest that these disorders may share etiology. To test this hypothesis, we have performed an association study of 54 genomic regions, previously implicated in PD or AD, in a large cohort of DLB cases and controls. The cohort comprised 788 DLB cases and 2624 controls. To minimize the issue of potential misdiagnosis, we have also performed the analysis including only neuropathologically proven DLB cases (667 cases). The results show that the APOE is a strong genetic risk factor for DLB, confirming previous findings, and that the SNCA and SCARB2 loci are also associated after a study-wise Bonferroni correction, although these have a different association profile than the associations reported for the same loci in PD. We have previously shown that the p.N370S variant in GBA is associated with DLB, which, together with the findings at the SCARB2 locus, suggests a role for lysosomal dysfunction in this disease. These results indicate that DLB has a unique genetic risk profile when compared with the two most common neurodegenerative diseases and that the lysosome may play an important role in the etiology of this disorder. We make all these data available.

Beschreibung: Mitochondrial DNA mutations at MT-ATP6 gene are relatively common in individuals suffering from striatal necrosis syndromes. These patients usually do not show apparent histochemical and/or biochemical signs of oxidative phosphorylation dysfunction. Because of this, MT-ATP6 is not typically analyzed in many other mitochondrial disorders that have not been previously associated to mutations in this gene. To correct this bias, we have performed a screening of the MT-ATP6 gene in a large collection of patients suspected of suffering different mitochondrial DNA (mtDNA) disorders. In three cases, biochemical, molecular-genetics and other analyses in patient tissues and cybrids were also carried out. We found three new pathologic mutations. Two of them in patients showing phenotypes that have not been commonly associated to mutations in the MT-ATP6 gene. These results remark the importance of sequencing the MT-ATP6 gene in patients with striatal necrosis syndromes, but also within other mitochondrial pathologies. This gene should be sequenced at least in all those patients suspected of suffering an mtDNA disorder disclosing normal results for histochemical and biochemical analyses of respiratory chain.

Beschreibung: Immunoglobulin-like domain containing receptor 1 ( ILDR1 ) is a poorly characterized gene that was first identified in lymphoma cells. Recently, ILDR1 has been found to be responsible for autosomal recessive hearing impairment DFNB42. Patients with ILDR1 mutations cause bilateral non-progressive moderate-to-profound sensorineural hearing impairment. However, the etiology and mechanism of ILDR1 -related hearing loss remains to be elucidated. In order to uncover the pathology of DFNB42 deafness, we used the morpholino injection technique to establish an ildr1b -morphant zebrafish model. Ildr1b -morphant zebrafish displayed defective hearing and imbalanced swimming, and developmental delays were seen in the semicircular canals of the inner ear. The gene expression profile and real-time PCR revealed down-regulation of atp1b2b (encoding Na + /K + transporting, beta 2b polypeptide) in ildr1b -morphant zebrafish. We found that injection of atp1b2b mRNA into ildr1b -knockdown zebrafish could rescue the phenotype of developmental delay of the semicircular canals. Moreover, ildr1b -morphant zebrafish had reduced numbers of lateral line neuromasts due to the disruption of lateral line primordium migration. In situ hybridization showed the involvement of attenuated FGF signaling and the chemokine receptor 4b ( cxcr4b ) and chemokine receptor 7b ( cxcr7b ) in posterior lateral line primordium of ildr1b -morphant zebrafish. We concluded that Ildr1b is crucial for the development of the inner ear and the lateral line system. This study provides the first evidence for the mechanism of Ildr1b on hearing in vivo and sheds light on the pathology of DFNB42.

Beschreibung: Huntington's disease (HD) is an inherited neurodegenerative disorder caused by abnormal expansion of CAG repeats in the gene encoding huntingtin. Mutant huntingtin undergoes proteolytic processing and its N-terminal fragment containing polyglutamine repeat accumulates as inclusion not only in nucleus but also in cytoplasm and neuronal processes. Here, we demonstrate that removal of ubiquitin ligase Ube3a selectively from HD mice brain resulted in accelerated disease phenotype and shorter lifespan in comparison with HD mice. The deficiency of Ube3a in HD mice brain also caused significant increase in global aggregates load, and these aggregates were less ubiquitinated when compared with age-matched HD mice. These Ube3a -maternal deficient HD mice also showed drastic reduction of DARPP-32, a dopamine-regulated phoshphoprotein in their striatum. These results emphasize the crucial role of Ube3a in the progression of HD and its immense potential as therapeutic target.

Beschreibung: Parent-of-origin-specific expression at imprinted genes is regulated by allele-specific DNA methylation at imprinting control regions (ICRs). This mechanism of gene regulation, where one element controls allelic expression of multiple genes, is not fully understood. Furthermore, the mechanism of gene dysregulation through ICR epimutations, such as loss or gain of DNA methylation, remains a mystery. We have used genetic mouse models to dissect ICR-mediated genetic and epigenetic regulation of imprinted gene expression. The H19/insulin-like growth factor 2 (Igf2) ICR has a multifunctional role including insulation, activation and repression. Microdeletions at the human H19/IGF2 ICR (IC1) are proposed to be responsible for IC1 epimutations associated with imprinting disorders such as Beckwith–Wiedemann syndrome (BWS). Here, we have generated and characterized a mouse model that mimics BWS microdeletions to define the role of the deleted sequence in establishing and maintaining epigenetic marks and imprinted expression at the H19/IGF2 locus. These mice carry a 1.3 kb deletion at the H19/Igf2 ICR [2,3] removing two of four CCCTC-binding factor (CTCF) sites and the intervening sequence, ~75% of the ICR. Surprisingly, the 2,3 deletion does not perturb DNA methylation at the ICR; however, it does disrupt imprinted expression. While repressive functions of the ICR are compromised by the deletion regardless of tissue type, insulator function is only disrupted in tissues of mesodermal origin where a significant amount of CTCF is poly(ADP-ribosyl)ated. These findings suggest that insulator activity of the H19/Igf2 ICR varies by cell type and may depend on cell-specific enhancers as well as posttranslational modifications of the insulator protein CTCF.

Beschreibung: Cerebral cavernous malformation (CCM) is a disease of vascular malformations known to be caused by mutations in one of three genes: CCM1 , CCM2 or CCM3 . Despite several studies, the mechanism of CCM lesion onset remains unclear. Using a Ccm1 knockout mouse model, we studied the morphogenesis of early lesion formation in the retina in order to provide insight into potential mechanisms. We demonstrate that lesions develop in a stereotypic location and pattern, preceded by endothelial hypersprouting as confirmed in a zebrafish model of disease. The vascular defects seen with loss of Ccm1 suggest a defect in endothelial flow response. Taken together, these results suggest new mechanisms of early CCM disease pathogenesis and provide a framework for further study.

Beschreibung: Genomic imprinting is the epigenetic process that results in monoallelic expression of genes depending on parental origin. These genes are known to be critical for placental development and fetal growth in mammals. Aberrant epigenetic profiles at imprinted loci, such as DNA methylation defects, are surprisingly rare in pregnancies with compromised fetal growth, while variations in transcriptional output from the expressed alleles of imprinted genes are more commonly reported in pregnancies complicated with intrauterine growth restriction (IUGR). To determine if PLAGL1 and HYMAI , two imprinted transcripts deregulated in Transient Neonatal Diabetes Mellitus, are involved in non-syndromic IUGR we compared the expression and DNA methylation levels in a large cohort of placental biopsies from IUGR and uneventful pregnancies. This revealed that despite appropriate maternal methylation at the shared PLAGL1 / HYMAI promoter, there was a loss of correlation between PLAGL1 and HYMAI expression in IUGR. This incongruity was due to higher HYMAI expression in IUGR gestations, coupled with PLAGL1 down-regulation in placentas from IUGR girls, but not boys. The PLAGL1 protein is a zinc-finger transcription factor that has been shown to be a master coordinator of a genetic growth network in mice. We observe PLAGL1 binding to the H19 / IGF2 shared enhancers in placentae, with significant correlations between PLAGL1 levels with H19 and IGF2 expression levels. In addition, PLAGL1 binding and expression also correlate with expression levels of metabolic regulator genes SLC2A4 , TCF4 and PPAR1 . Our results strongly suggest that fetal growth can be influenced by altered expression of the PLAGL1 gene network in human placenta.

Beschreibung: MicroRNAs (miRNAs) have emerged as a class of small, endogenous, regulatory RNAs that exhibit the ability to epigenetically modulate the translation of mRNAs into proteins. This feature enables them to control cell phenotypes and, consequently, modify cell function in a disease context. The role of inflammatory miRNAs in Alzheimer's disease (AD) and their ability to modulate glia responses are now beginning to be explored. In this study, we propose to disclose the functional role of miR-155, one of the most well studied immune-related miRNAs in AD-associated neuroinflammatory events, employing the 3xTg AD animal model. A strong upregulation of miR-155 levels was observed in the brain of 12-month-old 3xTg AD animals. This event occurred simultaneously with an increase of microglia and astrocyte activation, and before the appearance of extracellular Aβ aggregates, suggesting that less complex Aβ species, such as Aβ oligomers may contribute to early neuroinflammation. In addition, we investigated the contribution of miR-155 and the c-Jun transcription factor to the molecular mechanisms that underlie Aβ-mediated activation of glial cells. Our results suggest early miR-155 and c-Jun upregulation in the 3xTg AD mice, as well as in Aβ-activated microglia and astrocytes, thus contributing to the production of inflammatory mediators such as IL-6 and IFN-β. This effect is associated with a miR-155-dependent decrease of suppressor of cytokine signaling 1. Furthermore, since c-Jun silencing decreases the levels of miR-155 in Aβ-activated microglia and astrocytes, we propose that miR-155 targeting can constitute an interesting and promising approach to control neuroinflammation in AD.

Beschreibung: Huntington's disease (HD) is a neurodegenerative disorder caused by a CAG trinucleotide repeat expansion in the huntingtin ( HTT ) gene. Disease pathogenesis derives, at least in part, from the long polyglutamine tract encoded by mutant HTT . Therefore, considerable effort has been dedicated to the development of therapeutic strategies that significantly reduce the expression of the mutant HTT protein. Antisense oligonucleotides (ASOs) targeted to the CAG repeat region of HTT transcripts have been of particular interest due to their potential capacity to discriminate between normal and mutant HTT transcripts. Here, we focus on phosphorodiamidate morpholino oligomers (PMOs), ASOs that are especially stable, highly soluble and non-toxic. We designed three PMOs to selectively target expanded CAG repeat tracts (CTG22, CTG25 and CTG28), and two PMOs to selectively target sequences flanking the HTT CAG repeat (HTTex1a and HTTex1b). In HD patient–derived fibroblasts with expanded alleles containing 44, 77 or 109 CAG repeats, HTTex1a and HTTex1b were effective in suppressing the expression of mutant and non-mutant transcripts. CTGn PMOs also suppressed HTT expression, with the extent of suppression and the specificity for mutant transcripts dependent on the length of the targeted CAG repeat and on the CTG repeat length and concentration of the PMO. PMO CTG25 reduced HTT-induced cytotoxicity in vitro and suppressed mutant HTT expression in vivo in the N171-82Q transgenic mouse model. Finally, CTG28 reduced mutant HTT expression and improved the phenotype of Hdh Q7/Q150 knock-in HD mice. These data demonstrate the potential of PMOs as an approach to suppressing the expression of mutant HTT.

Beschreibung: Reduced expression of SMN protein causes spinal muscular atrophy (SMA), a neurodegenerative disorder leading to motor neuron dysfunction and loss. However, the molecular mechanisms by which SMN regulates neuronal dysfunction are not fully understood. Here, we report that reduced SMN protein level alters miRNA expression and distribution in neurons. In particular, miR-183 levels are increased in neurites of SMN-deficient neurons. We demonstrate that miR-183 regulates translation of mTor via direct binding to its 3' UTR. Interestingly, local axonal translation of mTor is reduced in SMN-deficient neurons, and this can be recovered by miR-183 inhibition. Finally, inhibition of miR-183 expression in the spinal cord of an SMA mouse model prolongs survival and improves motor function of Smn -mutant mice. Together, these observations suggest that axonal miRNAs and the mTOR pathway are previously unidentified molecular mechanisms contributing to SMA pathology.

Beschreibung: Microphthalmia-associated transcription factor ( MITF ) is a master regulator of pigmented cell survival and differentiation with direct transcriptional links to cell cycle, apoptosis and pigmentation. In mouse, Mitf is expressed early and uniformly in optic vesicle (OV) cells as they evaginate from the developing neural tube, and null Mitf mutations result in microphthalmia and pigmentation defects. However, homozygous mutations in MITF have not been identified in humans; therefore, little is known about its role in human retinogenesis. We used a human embryonic stem cell (hESC) model that recapitulates numerous aspects of retinal development, including OV specification and formation of retinal pigment epithelium (RPE) and neural retina progenitor cells (NRPCs), to investigate the earliest roles of MITF. During hESC differentiation toward a retinal lineage, a subset of MITF isoforms was expressed in a sequence and tissue distribution similar to that observed in mice. In addition, we found that promoters for the MITF-A , -D and -H isoforms were directly targeted by Visual Systems Homeobox 2 (VSX2), a transcription factor involved in patterning the OV toward a NRPC fate. We then manipulated MITF RNA and protein levels at early developmental stages and observed decreased expression of eye field transcription factors, reduced early OV cell proliferation and disrupted RPE maturation. This work provides a foundation for investigating MITF and other highly complex, multi-purposed transcription factors in a dynamic human developmental model system.

Beschreibung: The p.N478D missense mutation in human mitochondrial poly(A) polymerase (mtPAP) has previously been implicated in a form of spastic ataxia with optic atrophy. In this study, we have investigated fibroblast cell lines established from family members. The homozygous mutation resulted in the loss of polyadenylation of all mitochondrial transcripts assessed; however, oligoadenylation was retained. Interestingly, this had differential effects on transcript stability that were dependent on the particular species of transcript. These changes were accompanied by a severe loss of oxidative phosphorylation complexes I and IV, and perturbation of de novo mitochondrial protein synthesis. Decreases in transcript polyadenylation and in respiratory chain complexes were effectively rescued by overexpression of wild-type mtPAP. Both mutated and wild-type mtPAP localized to the mitochondrial RNA-processing granules thereby eliminating mislocalization as a cause of defective polyadenylation. In vitro polyadenylation assays revealed severely compromised activity by the mutated protein, which generated only short oligo(A) extensions on RNA substrates, irrespective of RNA secondary structure. The addition of LRPPRC/SLIRP, a mitochondrial RNA-binding complex, enhanced activity of the wild-type mtPAP resulting in increased overall tail length. The LRPPRC/SLIRP effect although present was less marked with mutated mtPAP, independent of RNA secondary structure. We conclude that (i) the polymerase activity of mtPAP can be modulated by the presence of LRPPRC/SLIRP, (ii) N478D mtPAP mutation decreases polymerase activity and (iii) the alteration in poly(A) length is sufficient to cause dysregulation of post-transcriptional expression and the pathogenic lack of respiratory chain complexes.

Beschreibung: During postnatal development, neuronal activity controls the remodeling of initially imprecise neuronal connections through the regulation of gene expression. MeCP2 binds to methylated DNA and modulates gene expression during neuronal development and MECP2 mutation causes the autistic disorder Rett syndrome. To investigate a role for MeCP2 in neuronal circuit refinement and to identify activity-dependent MeCP2 transcription regulations, we leveraged the precise organization and accessibility of olfactory sensory axons to manipulation of neuronal activity through odorant exposure in vivo . We demonstrate that olfactory sensory axons failed to develop complete convergence when Mecp2 is deficient in olfactory sensory neurons (OSNs) in an otherwise wild-type animal. Furthermore, we demonstrate that expression of selected adhesion genes was elevated in Mecp2 -deficient glomeruli, while acute odor stimulation in control mice resulted in significantly reduced MeCP2 binding to these gene loci, correlating with increased expression. Thus, MeCP2 is required for both circuitry refinement and activity-dependent transcriptional responses in OSNs.

Beschreibung: Simultaneous generation of neural cells and that of the nutrient-supplying vasculature during brain development is called neurovascular coupling. We report on a transgenic mouse with impaired transforming growth factor β (TGFβ)-signalling in forebrain-derived neural cells using a Foxg1-cre knock-in to drive the conditional knock-out of the Tgfbr2 . Although the expression of FOXG1 is assigned to neural progenitors and neurons of the telencephalon, Foxg1 cre/+ ; Tgfbr2 flox/flox (Tgfbr2-cKO) mutants displayed intracerebral haemorrhage. Blood vessels exhibited an atypical, clustered appearance were less in number and displayed reduced branching. Vascular endothelial growth factor (VEGF) A, insulin-like growth factor (IGF) 1, IGF2, TGFβ, inhibitor of DNA binding (ID) 1, thrombospondin (THBS) 2, and a disintegrin and metalloproteinase with thrombospondin motifs (ADAMTS) 1 were altered in either expression levels or tissue distribution. Accordingly, human umbilical vein endothelial cells (HUVEC) displayed branching defects after stimulation with conditioned medium (CM) that was derived from primary neural cultures of the ventral and dorsal telencephalon of Tgfbr2-cKO. Supplementing CM of Tgfbr2-cKO with VEGFA rescued these defects, but application of TGFβ aggravated them. HUVEC showed reduced migration towards CM of mutants compared with controls. Supplementing the CM with growth factors VEGFA, fibroblast growth factor (FGF) 2 and IGF1 partially restored HUVEC migration. In contrast, TGFβ supplementation further impaired migration of HUVEC. We observed differences along the dorso-ventral axis of the telencephalon with regard to the impact of these factors on the phenotype. Together these data establish a TGFBR2-dependent molecular crosstalk between neural and endothelial cells during brain vessel development. These findings will be useful to further elucidate neurovascular interaction in general and to understand pathologies of the blood vessel system such as intracerebral haemorrhages, hereditary haemorrhagic telangiectasia, Alzheimeŕs disease, cerebral amyloid angiopathy or tumour biology.

Beschreibung: Pneumoconiosis is the most serious occupational disease in China and its leading cause is occupational silica exposure. Pneumoconiosis takes several years to develop depending on the exposure level of silica. However, individual variation in the susceptibility to pneumoconiosis has been observed among the subjects with similar exposure. We conducted a genome-wide screening with 710 999 single nucleotide polymorphisms (SNPs) in a cohort of 400 coal workers (202 cases and 198 exposed controls) for pneumoconiosis susceptible loci. Seven promising variants were evaluated in an independent cohort of 568 coal workers (323 cases and 245 exposed controls), followed by a second replication on 463 iron ore workers (167 cases and 296 exposed controls). By pooling all of the genome-wide association studies and replication stages together, we found a genome-wide significant ( P 〈 5.0 x 10 –8 ) association for rs73329476 ( P = 1.74 x 10 –8 , OR = 2.17, 95% CI = 1.66–2.85) and two additional replicated associations for rs4320486 ( P 〈 0.05) and rs117626015 ( P 〈 0.05) with combined P -values of 4.29 x 10 –6 and 5.05 x 10 –6 , respectively. In addition, the risk allele T of rs73329476 was significantly associated with lower mRNA expression levels of carboxypeptidase M ( CPM ) in total cellular RNA from whole blood of 156 healthy individuals ( P = 0.0252). The identified pneumoconiosis susceptibility loci may provide new insights into the pathogenesis of pneumoconiosis, and may also have some clinical utility for risk prediction for pneumoconiosis and high-risk population screening for workers with occupational silica exposure.

Beschreibung: Hearing function is known to be heritable, but few significant and reproducible associations of genetic variants have been identified to date in the adult population. In this study, genome-wide association results of hearing function from the G-EAR consortium and TwinsUK were used for meta-analysis. Hearing ability in eight population samples of Northern and Southern European ancestry ( n = 4591) and the Silk Road ( n = 348) was measured using pure-tone audiometry and summarized using principal component (PC) analysis. Genome-wide association analyses for PC1–3 were conducted separately in each sample assuming an additive model adjusted for age, sex and relatedness of subjects. Meta-analysis was performed using 2.3 million single-nucleotide polymorphisms (SNPs) tested against each of the three PCs of hearing ability in 4939 individuals. A single SNP lying in intron 6 of the salt-inducible kinase 3 ( SIK3 ) gene was found to be associated with hearing PC2 ( P = 3.7 x 10 –8 ) and further supported by whole-genome sequence in a subset. To determine the relevance of this gene in the ear, expression of the Sik3 protein was studied in mouse cochlea of different ages. Sik3 was expressed in murine hair cells during early development and in cells of the spiral ganglion during early development and adulthood. Our results suggest a developmental role of Sik3 in hearing and may be required for the maintenance of adult auditory function.

Beschreibung: Genome instability, epigenetic remodelling and structural chromosomal rearrangements are hallmarks of cancer. However, the coordinated epigenetic effects of constitutional chromosomal rearrangements that disrupt genes associated with congenital neurodevelopmental diseases are poorly understood. To understand the genetic–epigenetic interplay at breakpoints of chromosomal translocations disrupting CG-rich loci, we quantified epigenetic modifications at DLGAP4 ( SAPAP4), a key post-synaptic density 95 (PSD95) associated gene, truncated by the chromosome translocation t(8;20)(p12;q11.23), co-segregating with cerebellar ataxia in a five-generation family. We report significant epigenetic remodelling of the DLGAP4 locus triggered by the t(8;20)(p12;q11.23) translocation and leading to dysregulation of DLGAP4 expression in affected carriers. Disruption of DLGAP4 results in monoallelic hypermethylation of the truncated DLGAP4 promoter CpG island. This induced hypermethylation is maintained in somatic cells of carriers across several generations in a t(8;20) dependent-manner however, is erased in the germ cells of the translocation carriers. Subsequently, chromatin remodelling of the locus-perturbed monoallelic expression of DLGAP4 mRNAs and non-coding RNAs in haploid cells having the translocation. Our results provide new mechanistic insight into the way a balanced chromosomal rearrangement associated with a neurodevelopmental disorder perturbs allele-specific epigenetic mechanisms at breakpoints leading to the deregulation of the truncated locus.

Beschreibung: We conducted blinded psychiatric assessments of 26 Amish subjects (52 ± 11 years) from four families with prevalent bipolar spectrum disorder, identified 10 potentially pathogenic alleles by exome sequencing, tested association of these alleles with clinical diagnoses in the larger Amish Study of Major Affective Disorder (ASMAD) cohort, and studied mutant potassium channels in neurons. Fourteen of 26 Amish had bipolar spectrum disorder. The only candidate allele shared among them was rs78247304, a non-synonymous variant of KCNH7 (c.1181G〉A, p.Arg394His). KCNH7 c.1181G〉A and nine other potentially pathogenic variants were subsequently tested within the ASMAD cohort, which consisted of 340 subjects grouped into controls subjects and affected subjects from overlapping clinical categories (bipolar 1 disorder, bipolar spectrum disorder and any major affective disorder). KCNH7 c.1181G〉A had the highest enrichment among individuals with bipolar spectrum disorder ( 2 = 7.3) and the strongest family-based association with bipolar 1 ( P = 0.021), bipolar spectrum ( P = 0.031) and any major affective disorder ( P = 0.016). In vitro, the p.Arg394His substitution allowed normal expression, trafficking, assembly and localization of HERG3/Kv11.3 channels, but altered the steady-state voltage dependence and kinetics of activation in neuronal cells. Although our genome-wide statistical results do not alone prove association, cumulative evidence from multiple independent sources (parallel genome-wide study cohorts, pharmacological studies of HERG-type potassium channels, electrophysiological data) implicates neuronal HERG3/Kv11.3 potassium channels in the pathophysiology of bipolar spectrum disorder. Such a finding, if corroborated by future studies, has implications for mental health services among the Amish, as well as development of drugs that specifically target HERG3/Kv11.3.

Beschreibung: Complex III (cytochrome bc 1 ) is a protein complex of the mitochondrial inner membrane that transfers electrons from ubiquinol to cytochrome c . Its assembly requires the coordinated expression of mitochondrial-encoded cytochrome b and nuclear-encoded subunits and assembly factors. Complex III deficiency is a severe multisystem disorder caused by mutations in subunit genes or assembly factors. Sequence-profile-based orthology predicts C11orf83 , hereafter named UQCC3 , to be the ortholog of the fungal complex III assembly factor CBP4 . We describe a homozygous c.59T〉A missense mutation in UQCC3 from a consanguineous patient diagnosed with isolated complex III deficiency, displaying lactic acidosis, hypoglycemia, hypotonia and delayed development without dysmorphic features. Patient fibroblasts have reduced complex III activity and lower levels of the holocomplex and its subunits than controls. They have no detectable UQCC3 protein and have lower levels of cytochrome b protein. Furthermore, in patient cells, cytochrome b is absent from a high-molecular-weight complex III. UQCC3 is reduced in cells depleted for the complex III assembly factors UQCC1 and UQCC2. Conversely, absence of UQCC3 in patient cells does not affect UQCC1 and UQCC2. This suggests that UQCC3 functions in the complex III assembly pathway downstream of UQCC1 and UQCC2 and is consistent with what is known about the function of Cbp4 and of the fungal orthologs of UQCC1 and UQCC2, Cbp3 and Cbp6. We conclude that UQCC3 functions in complex III assembly and that the c.59T〉A mutation has a causal role in complex III deficiency.

Beschreibung: Mutations in the photoreceptor-specific gene peripherin-2 ( PRPH-2 , also known as retinal degeneration slow/RDS) cause incurable retinal degeneration with a high degree of phenotypic variability. Patient phenotypes range from retinitis pigmentosa to various forms of macular and pattern dystrophy. Macular and pattern dystrophy in particular are associated with complex, poorly understood disease mechanisms, as severe vision loss is often associated both with defects in the photoreceptors, as well as the choroid and retinal pigment epithelium (RPE). Since there is currently no satisfactory model to study pattern dystrophy disease mechanisms, we generated a knockin mouse model expressing an RDS pattern dystrophy mutation, Y141C. Y141C mice exhibited clinical signs similar to those in patients including late-onset fundus abnormalities characteristic of RPE and choroidal defects and electroretinogram defects. Ultrastructural examination indicated that disc formation was initiated by the Y141C protein, but proper sizing and alignment of discs required wild-type RDS. The biochemical mechanism underlying these abnormalities was tied to defects in the normal process of RDS oligomerization which is required for proper RDS function. Y141C-RDS formed strikingly abnormal disulfide-linked complexes which were localized to the outer segment (OS) where they impaired the formation of proper OS structure. These data support a model of pattern dystrophy wherein a primary molecular defect occurring in all photoreceptors leads to secondary sequellae in adjacent tissues, an outcome which leads to macular vision loss. An understanding of the role of RDS in the interplay between these tissues significantly enhances our understanding of RDS-associated pathobiology and our ability to design rational treatment strategies.

Beschreibung: Motivation : Structural variation is common in human and cancer genomes. High-throughput DNA sequencing has enabled genome-scale surveys of structural variation. However, the short reads produced by these technologies limit the study of complex variants, particularly those involving repetitive regions. Recent ‘third-generation’ sequencing technologies provide single-molecule templates and longer sequencing reads, but at the cost of higher per-nucleotide error rates. Results : We present MultiBreak-SV, an algorithm to detect structural variants (SVs) from single molecule sequencing data, paired read sequencing data, or a combination of sequencing data from different platforms. We demonstrate that combining low-coverage third-generation data from Pacific Biosciences (PacBio) with high-coverage paired read data is advantageous on simulated chromosomes. We apply MultiBreak-SV to PacBio data from four human fosmids and show that it detects known SVs with high sensitivity and specificity. Finally, we perform a whole-genome analysis on PacBio data from a complete hydatidiform mole cell line and predict 1002 high-probability SVs, over half of which are confirmed by an Illumina-based assembly. Availability and implementation : MultiBreak-SV is available at http://compbio.cs.brown.edu/software/ . Contact : annaritz@vt.edu or braphael@cs.brown.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Insertions play an important role in genome evolution. However, such variants are difficult to detect from short-read sequencing data, especially when they exceed the paired-end insert size. Many approaches have been proposed to call short insertion variants based on paired-end mapping. However, there remains a lack of practical methods to detect and assemble long variants. Results: We propose here an original method, called M ind T he G ap , for the integrated detection and assembly of insertion variants from re-sequencing data. Importantly, it is designed to call insertions of any size, whether they are novel or duplicated, homozygous or heterozygous in the donor genome. M ind T he G ap uses an efficient k -mer-based method to detect insertion sites in a reference genome, and subsequently assemble them from the donor reads. M ind T he G ap showed high recall and precision on simulated datasets of various genome complexities. When applied to real Caenorhabditis elegans and human NA12878 datasets, M ind T he G ap detected and correctly assembled insertions 〉1 kb, using at most 14 GB of memory. Availability and implementation: http://mindthegap.genouest.org Contact: guillaume.rizk@inria.fr or claire.lemaitre@inria.fr

Beschreibung: Motivation: Most tumor samples are a heterogeneous mixture of cells, including admixture by normal (non-cancerous) cells and subpopulations of cancerous cells with different complements of somatic aberrations. This intra-tumor heterogeneity complicates the analysis of somatic aberrations in DNA sequencing data from tumor samples. Results: We describe an algorithm called THetA2 that infers the composition of a tumor sample—including not only tumor purity but also the number and content of tumor subpopulations—directly from both whole-genome (WGS) and whole-exome (WXS) high-throughput DNA sequencing data. This algorithm builds on our earlier Tumor Heterogeneity Analysis (THetA) algorithm in several important directions. These include improved ability to analyze highly rearranged genomes using a variety of data types: both WGS sequencing (including low ~7 x coverage) and WXS sequencing. We apply our improved THetA2 algorithm to WGS (including low-pass) and WXS sequence data from 18 samples from The Cancer Genome Atlas (TCGA). We find that the improved algorithm is substantially faster and identifies numerous tumor samples containing subclonal populations in the TCGA data, including in one highly rearranged sample for which other tumor purity estimation algorithms were unable to estimate tumor purity. Availability and implementation: An implementation of THetA2 is available at http://compbio.cs.brown.edu/software Contact: layla@cs.brown.edu or braphael@brown.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: High-throughput sequencing (HTS) technologies have made low-cost sequencing of large numbers of samples commonplace. An explosion in the type, not just number, of sequencing experiments has also taken place including genome re-sequencing, population-scale variation detection, whole transcriptome sequencing and genome-wide analysis of protein-bound nucleic acids. Results: We present Artemis as a tool for integrated visualization and computational analysis of different types of HTS datasets in the context of a reference genome and its corresponding annotation. Availability: Artemis is freely available (under a GPL licence) for download (for MacOSX, UNIX and Windows) at the Wellcome Trust Sanger Institute websites: http://www.sanger.ac.uk/resources/software/artemis/ . Contact: artemis@sanger.ac.uk ; tjc@sanger.ac.uk

Beschreibung: : microRibonucleic acid (miRNAs) are small regulatory molecules that act by mRNA degradation or via translational repression. Although many miRNAs are ubiquitously expressed, a small subset have differential expression patterns that may give rise to tissue-specific complexes. Motivation: This work studies gene targeting patterns amongst miRNAs with differential expression profiles, and links this to control and regulation of protein complexes. Results: We find that, when a pair of miRNAs are not expressed in the same tissues, there is a higher tendency for them to target the direct partners of the same hub proteins. At the same time, they also avoid targeting the same set of hub-spokes. Moreover, the complexes corresponding to these hub-spokes tend to be specific and nonoverlapping. This suggests that the effect of miRNAs on the formation of complexes is specific. Contact: wongls@comp.nus.edu.sg Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Small interfering RNAs (siRNAs) are produced from much longer sequences of double-stranded RNA precursors through cleavage by Dicer or a Dicer-like protein. These small RNAs play a key role in genetic and epigenetic regulation; however, a full understanding of the mechanisms by which they operate depends on the characterization of the precursors from which they are derived. Results: High-throughput sequencing of small RNA populations allows the locations of the double-stranded RNA precursors to be inferred. We have developed methods to analyse small RNA sequencing data from multiple biological sources, taking into account replicate information, to identify robust sets of siRNA precursors. Our methods show good performance on both a set of small RNA sequencing data in Arabidopsis thaliana and simulated datasets. Availability: Our methods are available as the Bioconductor ( www.bioconductor.org ) package segmentSeq (version 1.5.6 and above). Contact: tjh48@cam.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Intrinsically disordered regions are key for the function of numerous proteins, and the scant available experimental annotations suggest the existence of different disorder flavors. While efficient predictions are required to annotate entire genomes, most existing methods require sequence profiles for disorder prediction, making them cumbersome for high-throughput applications. Results: In this work, we present an ensemble of protein disorder predictors called ESpritz. These are based on bidirectional recursive neural networks and trained on three different flavors of disorder, including a novel NMR flexibility predictor. ESpritz can produce fast and accurate sequence-only predictions, annotating entire genomes in the order of hours on a single processor core. Alternatively, a slower but slightly more accurate ESpritz variant using sequence profiles can be used for applications requiring maximum performance. Two levels of prediction confidence allow either to maximize reasonable disorder detection or to limit expected false positives to 5%. ESpritz performs consistently well on the recent CASP9 data, reaching a S w measure of 54.82 and area under the receiver operator curve of 0.856. The fast predictor is four orders of magnitude faster and remains better than most publicly available CASP9 methods, making it ideal for genomic scale predictions. Conclusions: ESpritz predicts three flavors of disorder at two distinct false positive rates, either with a fast or slower and slightly more accurate approach. Given its state-of-the-art performance, it can be especially useful for high-throughput applications. Availability: Both a web server for high-throughput analysis and a Linux executable version of ESpritz are available from: http://protein.bio.unipd.it/espritz/ Contact: silvio.tosatto@unipd.it Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Clustering protein structures is an important task in structural bioinformatics. De novo structure prediction, for example, often involves a clustering step for finding the best prediction. Other applications include assigning proteins to fold families and analyzing molecular dynamics trajectories. Results: We present Pleiades, a novel approach to clustering protein structures with a rigorous mathematical underpinning. The method approximates clustering based on the root mean square deviation by first mapping structures to Gauss integral vectors—which were introduced by Røgen and co-workers—and subsequently performing K-means clustering. Conclusions: Compared to current methods, Pleiades dramatically improves on the time needed to perform clustering, and can cluster a significantly larger number of structures, while providing state-of-the-art results. The number of low energy structures generated in a typical folding study, which is in the order of 50 000 structures, can be clustered within seconds to minutes. Contact: thamelry@binf.ku.dk ; harder@binf.ku.dk Supplementary Information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Transmembrane β barrel proteins (TMBs) are found in the outer membrane of Gram-negative bacteria, chloroplast and mitochondria. They play a major role in the translocation machinery, pore formation, membrane anchoring and ion exchange. TMBs are also promising targets for antimicrobial drugs and vaccines. Given the difficulty in membrane protein structure determination, computational methods to identify TMBs and predict the topology of TMBs are important. Results: Here, we present BOCTOPUS; an improved method for the topology prediction of TMBs by employing a combination of support vector machines (SVMs) and Hidden Markov Models (HMMs). The SVMs and HMMs account for local and global residue preferences, respectively. Based on a 10-fold cross-validation test, BOCTOPUS performs better than all existing methods, reaching a Q3 accuracy of 87%. Further, BOCTOPUS predicted the correct number of strands for 83% proteins in the dataset. BOCTOPUS might also help in reliable identification of TMBs by using it as an additional filter to methods specialized in this task. Availability: BOCTOPUS is freely available as a web server at: http://boctopus.cbr.su.se/ . The datasets used for training and evaluations are also available from this site. Contact: arne@bioinfo.se Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: High-dimensional data such as microarrays have created new challenges to traditional statistical methods. One such example is on class prediction with high-dimension, low-sample size data. Due to the small sample size, the sample mean estimates are usually unreliable. As a consequence, the performance of the class prediction methods using the sample mean may also be unsatisfactory. To obtain more accurate estimation of parameters some statistical methods, such as regularizations through shrinkage, are often desired. Results: In this article, we investigate the family of shrinkage estimators for the mean value under the quadratic loss function. The optimal shrinkage parameter is proposed under the scenario when the sample size is fixed and the dimension is large. We then construct a shrinkage-based diagonal discriminant rule by replacing the sample mean by the proposed shrinkage mean. Finally, we demonstrate via simulation studies and real data analysis that the proposed shrinkage-based rule outperforms its original competitor in a wide range of settings. Contact: tongt@hkbu.edu.hk

Beschreibung: Motivation: The advent of high-throughput sequencing technologies is revolutionizing our ability in discovering and genotyping DNA copy number variants (CNVs). Read count-based approaches are able to detect CNV regions with an unprecedented resolution. Although this computational strategy has been recently introduced in literature, much work has been already done for the preparation, normalization and analysis of this kind of data. Results: Here we face the many aspects that cover the detection of CNVs by using read count approach. We first study the characteristics and systematic biases of read count distributions, focusing on the normalization methods designed for removing these biases. Subsequently, we compare the algorithms designed to detect the boundaries of CNVs and we investigate the ability of read count data to predict the exact number of DNA copy. Finally, we review the tools publicly available for analysing read count data. To better understand the state of the art of read count approaches, we compare the performance of the three most widely used sequencing technologies (Illumina Genome Analyzer, Roche 454 and Life Technologies SOLiD) in all the analyses that we perform. Contact: albertomagi@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: We investigate and quantify the generalizability of the white blood cell (WBC) transcriptome to the general, multiorgan transcriptome. We use data from the NCBI's Gene Expression Omnibus (GEO) public repository to define two datasets for comparison, WBC and OO (Other Organ) sets. Results: Comprehensive pair-wise correlation and expression level profiles are calculated for both datasets (with sizes of 81 and 1463, respectively). We have used mapping and ranking across the Gene Ontology (GO) categories to quantify similarity between the two sets. GO mappings of the most correlated and highly expressed genes from the two datasets tightly match, with the notable exceptions of components of the ribosome, cell adhesion and immune response. That is, 10 877 or 48.8% of all measured genes do not change 〉10% of rank range between WBC and OO; only 878 (3.9%) change rank 〉50%. Two trans -tissue gene lists are defined, the most changing and the least changing genes in expression rank. We also provide a general, quantitative measure of the probability of expression rank and correlation profile in the OO system given the expression rank and correlation profile in the WBC dataset. Contact: vvaltchinov@partners.org Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: The understanding of the molecular sources for diseases like cancer can be significantly improved by computational models. Recently, Boolean networks have become very popular for modeling signaling and regulatory networks. However, such models rely on a set of Boolean functions that are in general not known. Unfortunately, while detailed information on the molecular interactions becomes available in large scale through electronic databases, the information on the Boolean functions does not become available simultaneously and has to be included manually into the models, if at all known. Results: We propose a new Boolean approach which can directly utilize the mechanistic network information available through modern databases. The Boolean function is implicitly defined by the reaction mechanisms. Special care has been taken for the treatment of kinetic features like inhibition. The method has been applied to a signaling model combining the Wnt and MAPK pathway. Availability: A sample C++ implementation of the proposed method is available for Linux and compatible systems through http://code.google.com/p/libscopes/wiki/Paper2011 Contact: handorf@physik.hu-berlin.de Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Multiple sequence alignment (MSA) is a core method in bioinformatics. The accuracy of such alignments may influence the success of downstream analyses such as phylogenetic inference, protein structure prediction, and functional prediction. The importance of MSA has lead to the proliferation of MSA methods, with different objective functions and heuristics to search for the optimal MSA. Different methods of inferring MSAs produce different results in all but the most trivial cases. By measuring the differences between inferred alignments, we may be able to develop an understanding of how these differences (i) relate to the objective functions and heuristics used in MSA methods, and (ii) affect downstream analyses. Results: We introduce four metrics to compare MSAs, which include the position in a sequence where a gap occurs or the location on a phylogenetic tree where an insertion or deletion (indel) event occurs. We use both real and synthetic data to explore the information given by these metrics and demonstrate how the different metrics in combination can yield more information about MSA methods and the differences between them. Availability: MetAl is a free software implementation of these metrics in Haskell. Source and binaries for Windows, Linux and Mac OS X are available from http://kumiho.smith.man.ac.uk/whelan/software/metal/ . Contact: simon.whelan@manchester.ac.uk

Beschreibung: Motivation: Peptide detection is a crucial step in mass spectrometry (MS) based proteomics. Most existing algorithms are based upon greedy isotope template matching and thus may be prone to error propagation and ineffective to detect overlapping peptides. In addition, existing algorithms usually work at different charge states separately, isolating useful information that can be drawn from other charge states, which may lead to poor detection of low abundance peptides. Results: BPDA2d models spectra as a mixture of candidate peptide signals and systematically evaluates all possible combinations of possible peptide candidates to interpret the given spectra. For each candidate, BPDA2d takes into account its elution profile, charge state distribution and isotope pattern, and it combines all evidence to infer the candidate's signal and existence probability. By piecing all evidence together—especially by deriving information across charge states—low abundance peptides can be better identified and peptide detection rates can be improved. Instead of local template matching, BPDA2d performs global optimization for all candidates and systematically optimizes their signals. Since BPDA2d looks for the optimal among all possible interpretations of the given spectra, it has the capability in handling complex spectra where features overlap. BPDA2d estimates the posterior existence probability of detected peptides, which can be directly used for probability-based evaluation in subsequent processing steps. Our experiments indicate that BPDA2d outperforms state-of-the-art detection methods on both simulated data and real liquid chromatography–mass spectrometry data, according to sensitivity and detection accuracy. Availability: The BPDA2d software package is available at http://gsp.tamu.edu/Publications/supplementary/sun11a/ Contact: Michelle.Zhang@utsa.edu ; edward@ece.tamu.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: The continued progress in developing technological platforms, availability of many published experimental datasets, as well as different statistical methods to analyze those data have allowed approaching the same research question using various methods simultaneously. To get the best out of all these alternatives, we need to integrate their results in an unbiased manner. Prioritized gene lists are a common result presentation method in genomic data analysis applications. Thus, the rank aggregation methods can become a useful and general solution for the integration task. Results: Standard rank aggregation methods are often ill-suited for biological settings where the gene lists are inherently noisy. As a remedy, we propose a novel robust rank aggregation (RRA) method. Our method detects genes that are ranked consistently better than expected under null hypothesis of uncorrelated inputs and assigns a significance score for each gene. The underlying probabilistic model makes the algorithm parameter free and robust to outliers, noise and errors. Significance scores also provide a rigorous way to keep only the statistically relevant genes in the final list. These properties make our approach robust and compelling for many settings. Availability: All the methods are implemented as a GNU R package R obust R ank A ggreg , freely available at the Comprehensive R Archive Network http://cran.r-project.org/ . Contact: vilo@ut.ee Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: : CLARE is a computational method designed to reveal sequence encryption of tissue-specific regulatory elements. Starting with a set of regulatory elements known to be active in a particular tissue/process, it learns the sequence code of the input set and builds a predictive model from features specific to those elements. The resulting model can then be applied to user-supplied genomic regions to identify novel candidate regulatory elements. CLARE's model also provides a detailed analysis of transcription factors that most likely bind to the elements, making it an invaluable tool for understanding mechanisms of tissue-specific gene regulation. Availability: CLARE is freely accessible at http://clare.dcode.org/ . Contact: taherl@ncbi.nlm.nih.gov ; ovcharen@nih.gov Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: We present a pipeline for the pre-processing, quality assessment, read distribution and methylation estimation for methylated DNA immunoprecipitation (MeDIP)-sequence datasets. This is the first MeDIP-seq-specific analytic pipeline that starts at the output of the sequencers. This pipeline will reduce the data analysis load on staff and allows the easy and straightforward analysis of sequencing data for DNA methylation. The pipeline integrates customized scripting and several existing tools, which can deal with both paired and single end data. Availability: The package and extensive documentation, and comparison to public data is available at http://life.tongji.edu.cn/meqa/ Contact: jhuang@cephb.fr

Beschreibung: Motivation: A plethora of bioinformatics analysis has led to the discovery of numerous gene sets, which can be interpreted as discrete measurements emitted from latent signaling pathways. Their potential to infer signaling pathway structures, however, has not been sufficiently exploited. Existing methods accommodating discrete data do not explicitly consider signal cascading mechanisms that characterize a signaling pathway. Novel computational methods are thus needed to fully utilize gene sets and broaden the scope from focusing only on pairwise interactions to the more general cascading events in the inference of signaling pathway structures. Results: We propose a gene set based simulated annealing (SA) algorithm for the reconstruction of signaling pathway structures. A signaling pathway structure is a directed graph containing up to a few hundred nodes and many overlapping signal cascades, where each cascade represents a chain of molecular interactions from the cell surface to the nucleus. Gene sets in our context refer to discrete sets of genes participating in signal cascades, the basic building blocks of a signaling pathway, with no prior information about gene orderings in the cascades. From a compendium of gene sets related to a pathway, SA aims to search for signal cascades that characterize the optimal signaling pathway structure. In the search process, the extent of overlap among signal cascades is used to measure the optimality of a structure. Throughout, we treat gene sets as random samples from a first-order Markov chain model. We evaluated the performance of SA in three case studies. In the first study conducted on 83 KEGG pathways, SA demonstrated a significantly better performance than Bayesian network methods. Since both SA and Bayesian network methods accommodate discrete data, use a ‘search and score’ network learning strategy and output a directed network, they can be compared in terms of performance and computational time. In the second study, we compared SA and Bayesian network methods using four benchmark datasets from DREAM. In our final study, we showcased two context-specific signaling pathways activated in breast cancer. Availibility: Source codes are available from http://dl.dropbox.com/u/16000775/sa_sc.zip Contact: dzhu@wayne.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: : We provide a Bioconductor package with quality assessment, processing and visualization tools for high-throughput sequencing data, with emphasis in ChIP-seq and RNA-seq studies. It includes detection of outliers and biases, inefficient immuno-precipitation and overamplification artifacts, de novo identification of read-rich genomic regions and visualization of the location and coverage of genomic region lists. Availability: www.bioconductor.org Contact: david.rossell@irbbarcelona.org Supplementary information: Supplementary data available at Bioinformatics online.

Beschreibung: Motivation: We study a stochastic method for approximating the set of local minima in partial RNA folding landscapes associated with a bounded-distance neighbourhood of folding conformations. The conformations are limited to RNA secondary structures without pseudoknots. The method aims at exploring partial energy landscapes p L induced by folding simulations and their underlying neighbourhood relations. It combines an approximation of the number of local optima devised by Garnier and Kallel (2002) with a run-time estimation for identifying sets of local optima established by Reeves and Eremeev (2004). Results: The method is tested on nine sequences of length between 50 nt and 400 nt, which allows us to compare the results with data generated by RNAsubopt and subsequent barrier tree calculations. On the nine sequences, the method captures on average 92% of local minima with settings designed for a target of 95%. The run-time of the heuristic can be estimated by O ( n 2 D ln), where n is the sequence length, is the number of local minima in the partial landscape p L under consideration and D is the maximum number of steepest descent steps in attraction basins associated with p L . Contact: a.albrecht@qub.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: RNA-seq is a powerful technology for the study of transcriptome profiles that uses deep-sequencing technologies. Moreover, it may be used for cellular phenotyping and help establishing the etiology of diseases characterized by abnormal splicing patterns. In RNA-Seq, the exact nature of splicing events is buried in the reads that span exon–exon boundaries. The accurate and efficient mapping of these reads to the reference genome is a major challenge. Results: We developed PASSion, a pattern growth algorithm-based pipeline for splice site detection in paired-end RNA-Seq reads. Comparing the performance of PASSion to three existing RNA-Seq analysis pipelines, TopHat, MapSplice and HMMSplicer, revealed that PASSion is competitive with these packages. Moreover, the performance of PASSion is not affected by read length and coverage. It performs better than the other three approaches when detecting junctions in highly abundant transcripts. PASSion has the ability to detect junctions that do not have known splicing motifs, which cannot be found by the other tools. Of the two public RNA-Seq datasets, PASSion predicted ~ 137 000 and 173 000 splicing events, of which on average 82 are known junctions annotated in the Ensembl transcript database and 18% are novel. In addition, our package can discover differential and shared splicing patterns among multiple samples. Availability: The code and utilities can be freely downloaded from https://trac.nbic.nl/passion and ftp://ftp.sanger.ac.uk/pub/zn1/passion Contact: y.zhang@lumc.nl ; k.ye@lumc.nl Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: The completion of 168 genome sequences from a single population of Drosophila melanogaster provides a global view of genomic variation and an understanding of the evolutionary forces shaping the patterns of DNA polymorphism and divergence along the genome. Results: We present the ‘Population Drosophila Browser’ (PopDrowser), a new genome browser specially designed for the automatic analysis and representation of genetic variation across the D. melanogaster genome sequence. PopDrowser allows estimating and visualizing the values of a number of DNA polymorphism and divergence summary statistics, linkage disequilibrium parameters and several neutrality tests. PopDrowser also allows performing custom analyses on-the-fly using user-selected parameters. Availability: PopDrowser is freely available from http://PopDrowser.uab.cat . Contact: miquel.ramia@uab.cat

Beschreibung: Motivation: Probabilistic approaches for inferring transcription factor binding sites (TFBSs) and regulatory motifs from DNA sequences have been developed for over two decades. Previous work has shown that prediction accuracy can be significantly improved by incorporating features such as the competition of multiple transcription factors (TFs) for binding to nearby sites, the tendency of TFBSs for co-regulated TFs to cluster and form cis -regulatory modules and explicit evolutionary modeling of conservation of TFBSs across orthologous sequences. However, currently available tools only incorporate some of these features, and significant methodological hurdles hampered their synthesis into a single consistent probabilistic framework. Results: We present MotEvo, a integrated suite of Bayesian probabilistic methods for the prediction of TFBSs and inference of regulatory motifs from multiple alignments of phylogenetically related DNA sequences, which incorporates all features just mentioned. In addition, MotEvo incorporates a novel model for detecting unknown functional elements that are under evolutionary constraint, and a new robust model for treating gain and loss of TFBSs along a phylogeny. Rigorous benchmarking tests on ChIP-seq datasets show that MotEvo's novel features significantly improve the accuracy of TFBS prediction, motif inference and enhancer prediction. Availability: Source code, a user manual and files with several example applications are available at www.swissregulon.unibas.ch . Contact: erik.vannimwegen@unibas.ch Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: : We present LaTcOm, a new web tool, which offers several alternative methods for ‘rare codon cluster’ (RCC) identification from a single and simple graphical user interface. In the current version, three RCC detection schemes are implemented: the recently described %MinMax algorithm and a simplified sliding window approach, along with a novel modification of a linear-time algorithm for the detection of maximally scoring subsequences tailored to the RCC detection problem. Among a number of user tunable parameters, several codon-based scales relevant for RCC detection are available, including tRNA abundance values from Escherichia coli and several codon usage tables from a selection of genomes. Furthermore, useful scale transformations may be performed upon user request (e.g. linear, sigmoid). Users may choose to visualize RCC positions within the submitted sequences either with graphical representations or in textual form for further processing. Availability: LaTcOm is freely available online at the URL http://troodos.biol.ucy.ac.cy/latcom.html . Contact: vprobon@ucy.ac.cy Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: : Many existing databases annotate experimentally characterized single nucleotide polymorphisms (SNPs). Each non-synonymous SNP (nsSNP) changes one amino acid in the gene product (single amino acid substitution;SAAS). This change can either affect protein function or be neutral in that respect. Most polymorphisms lack experimental annotation of their functional impact. Here, we introduce SNPdbe—SNP database of effects, with predictions of computationally annotated functional impacts of SNPs. Database entries represent nsSNPs in dbSNP and 1000 Genomes collection, as well as variants from UniProt and PMD. SAASs come from 〉2600 organisms; ‘human’ being the most prevalent. The impact of each SAAS on protein function is predicted using the SNAP and SIFT algorithms and augmented with experimentally derived function/structure information and disease associations from PMD, OMIM and UniProt. SNPdbe is consistently updated and easily augmented with new sources of information. The database is available as an MySQL dump and via a web front end that allows searches with any combination of organism names, sequences and mutation IDs. Availability: http://www.rostlab.org/services/snpdbe Contact: schaefer@rostlab.org ; snpdbe@rostlab.org

Beschreibung: : We have implemented in a single package all the features required for extracting, visualizing and manipulating fully conserved positions as well as those with a family-dependent conservation pattern in multiple sequence alignments. The program allows, among other things, to run different methods for extracting these positions, combine the results and visualize them in protein 3D structures and sequence spaces. Availability and implementation: JDet is a multiplatform application written in Java. It is freely available, including the source code, at http://csbg.cnb.csic.es/JDet . The package includes two of our recently developed programs for detecting functional positions in protein alignments ( Xdet and S3Det ), and support for other methods can be added as plug-ins. A help file and a guided tutorial for JDet are also available. Contact: pazos@cnb.csic.es

Beschreibung: : VarSifter is a graphical software tool for desktop computers that allows investigators of varying computational skills to easily and quickly sort, filter, and sift through sequence variation data. A variety of filters and a custom query framework allow filtering based on any combination of sample and annotation information. By simplifying visualization and analyses of exome-scale sequence variation data, this program will help bring the power and promise of massively-parallel DNA sequencing to a broader group of researchers. Availability and Implementation: VarSifter is written in Java, and is freely available in source and binary versions, along with a User Guide, at http://research.nhgri.nih.gov/software/VarSifter/ . Contact: mullikin@mail.nih.gov Supplementary Information: Additional figures and methods available online at the journal's website.

Beschreibung: Insect cells are widely used for recombinant glycoprotein production, but they cannot provide the glycosylation patterns required for some biotechnological applications. This problem has been addressed by genetically engineering insect cells to express mammalian genes encoding various glycoprotein glycan processing functions. However, for various reasons, the impact of a mammalian cytosine-5'-monophospho (CMP)-sialic acid transporter has not yet been examined. Thus, we transformed Spodoptera frugiperda (Sf9) cells with six mammalian genes to generate a new cell line, SfSWT-4, that can produce sialylated glycoproteins when cultured with the sialic acid precursor, N -acetylmannosamine. We then super-transformed SfSWT-4 with a human CMP-sialic acid transporter (hCSAT) gene to isolate a daughter cell line, SfSWT-6, which expressed the hCSAT gene in addition to the other mammalian glycogenes. SfSWT-6 cells had higher levels of cell surface sialylation and also supported higher levels of recombinant glycoprotein sialylation, particularly when cultured with low concentrations of N -acetylmannosamine. Thus, hCSAT expression has an impact on glycoprotein sialylation, can reduce the cost of recombinant glycoprotein production and therefore should be included in ongoing efforts to glycoengineer the baculovirus-insect cell system. The results of this study also contributed new insights into the endogenous mechanism and potential mechanisms of CMP-sialic acid accumulation in the Golgi apparatus of lepidopteran insect cells.

Beschreibung: Chondroitin sulfate (CS) chains regulate the development of the central nervous system in vertebrates and are linear polysaccharides consisting of variously sulfated repeating disaccharides, [–4GlcUAβ1-3GalNAcβ1–] n , where GlcUA and GalNAc represent d -glucuronic acid and N -acetyl- d -galactosamine, respectively. CS chains containing D-disaccharide units [GlcUA(2- O -sulfate)-GalNAc(6- O -sulfate)] are involved in the development of cerebellar Purkinje cells and neurite outgrowth-promoting activity through interaction with a neurotrophic factor, pleiotrophin, resulting in the regulation of signaling. In this study, to obtain further structural information on the CS chains containing d -disaccharide units involved in brain development, oligosaccharides containing D-units were isolated from a shark fin cartilage. Seven novel hexasaccharide sequences, O-D-D, A-D-D, C-D-D, E-A-D, D-D-C, E-D-D and A-B-D, in addition to three previously reported sequences, C-A-D, C-D-C and A-D-A, were isolated from a CS preparation of shark fin cartilage after exhaustive digestion with chondroitinase AC-I, which cannot act on the galactosaminidic linkages bound to D-units. The symbol stands for a 4,5-unsaturated bond of uronic acids, whereas A, B, C, D, E and O represent [GlcUA-GalNAc(4- O -sulfate)], [GlcUA(2- O -sulfate)-GalNAc(4- O -sulfate)], [GlcUA-GalNAc(6- O -sulfate)], [GlcUA(2- O -sulfate)-GalNAc(6- O -sulfate)], [GlcUA-GalNAc(4- O -, 6- O -sulfate)] and [GlcUA-GalNAc], respectively. In binding studies using an anti-CS monoclonal antibody, MO-225, the epitopes of which are involved in cerebellar development in mammals, novel epitope structures, A-D-A, A-D-D and A-B-D, were revealed. Hexasaccharides containing two consecutive D-units or a B-unit will be useful for the structural and functional analyses of CS chains particularly in the neuroglycobiological fields.

Beschreibung: Protein O -fucosyltransferase 1 (Pofut1) and protein O -fucosyltransferase 2 (Pofut2) add O -linked fucose at distinct consensus sequences in properly folded epidermal growth factor (EGF)-like repeats and thrombospondin type-1 (TSR) repeats, respectively. Glycan chain elongation past O -fucose can occur to yield a tetrasaccharide on EGF repeats and a disaccharide on TSRs. Elimination of Pofut1 in mice causes embryonic lethality with Notch-like phenotypes demonstrating that O -fucosylation of Notch is essential for its function. Similarly, elimination of Pofut2 results in an early embryonic lethal phenotype in mice, although the molecular mechanism for the lethality is unknown. The recent development of sugar analogs has revolutionized the study of glycans by providing a convenient method for labeling and tracking glycosylation. In order to study O -fucosylation, we took advantage of the recently developed reporter, 6-alkynyl fucose. Using the Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC), or "click" reaction, azido-biotin allows tagging and detection of 6AF-modified proteins. Here we examine whether proteins containing EGF repeats or TSRs with O -fucose consensus sequences are specifically modified with 6AF in cell culture. Using mass spectrometry (MS), we demonstrate that 6AF is efficiently incorporated onto the appropriate consensus sequences on EGF repeats and TSRs. Furthermore, the elongation of the O -fucose monosaccharide on EGF repeats and TSRs is not hampered when 6AF is used. These results show that 6AF is efficiently utilized in a truly bioorthogonal manner by Pofut1, Pofut2 and the enzymes that elongate O -fucose, providing evidence that 6AF is a significant new tool in the study of protein O -fucosylation.

Beschreibung: Xanthan is a polysaccharide secreted by Xanthomonas campestris that contains pentameric repeat units. The biosynthesis of xanthan involves an operon composed of 12 genes ( gumB to gumM ). In this study, we analyzed the proteins encoded by gumB and gumC . Membrane fractionation showed that GumB was mainly associated with the outer membrane, whereas GumC was an inner membrane protein. By in silico analysis and specific globomycin inhibition, GumB was characterized as a lipoprotein. By reporter enzyme assays, GumC was shown to contain two transmembrane segments flanking a large periplasmic domain. We confirmed that gumB and gumC mutant strains uncoupled the synthesis of the lipid-linked repeat unit from the polymerization process. We studied the effects of gumB and gumC gene amplification on the production, composition and viscosity of xanthan. Overexpression of GumB, GumC or GumB and GumC simultaneously did not affect the total amount or the chemical composition of the polymer. GumB overexpression did not affect xanthan viscosity; however, a moderate increase in xanthan viscosity was achieved when GumC protein levels were increased 5-fold. Partial degradation of GumC was observed when only that protein was overexpressed; but co-expression of GumB and GumC diminished GumC degradation and resulted in higher xanthan viscosity than individual GumB or GumC overexpression. Compared with xanthan from the wild-type strain, longer polymer chains from the strain that simultaneously overexpressed GumB and GumC were observed by atomic force microscopy. Our results suggest that GumB–GumC protein levels modulate xanthan chain length, which results in altered polymer viscosity.

Beschreibung: Bifidobacterium bifidum is one of the most frequently found bifidobacteria in the intestines of newborn infants. We previously reported that B. bifidum possesses unique metabolic pathways for O -linked glycans on gastrointestinal mucin (Yoshida E, Sakurama H, Kiyohara M, Nakajima M, Kitaoka M, Ashida H, Hirose J, Katayama T, Yamamoto K, Kumagai H. 2012. Bifidobacterium longum subsp. infantis uses two different β-galactosidases for selectively degrading type-1 and type-2 human milk oligosaccharides. Glycobiology . 22:361–368). The nonreducing termini of O -linked glycans on mucin are frequently covered with histo-blood group antigens. Here, we identified a gene agabb from B. bifidum JCM 1254, which encodes glycoside hydrolase (GH) family 110 α-galactosidase. AgaBb is a 1289-amino acid polypeptide containing an N-terminal signal sequence, a GH110 domain, a carbohydrate-binding module (CBM) 51 domain, a bacterial Ig-like (Big) 2 domain and a C-terminal transmembrane region, in this order. The recombinant enzyme expressed in Escherichia coli hydrolyzed α1,3-linked Gal in branched blood group B antigen [Galα1-3(Fucα1-2)Galβ1-R], but not in a linear xenotransplantation antigen (Galα1-3Galβ1-R). The enzyme also acted on group B human salivary mucin and erythrocytes. We also revealed that CBM51 specifically bound blood group B antigen using both isothermal titration calorimetry and a solid-phase binding assay, and it enhanced the affinity of the enzyme toward substrates with multivalent B antigens. We suggest that this enzyme plays an important role in degrading B antigens to acquire nutrients from mucin oligosaccharides in the gastrointestinal tracts.

Beschreibung: Alg3 of Saccharomyces cerevisiae catalyzes the mannosyl transfer from Man-P-Dol to Man 5 GlcNAc 2 -PP-Dol resulting in the formation of Man 6 GlcNAc 2 -PP-Dol, which is then further processed to the final precursor oligosaccharide Glc 3 Man 9 GlcNAc 2 for N-glycosylation of proteins. Here, we identified the alg3 gene of the mushroom-forming fungus Schizophyllum commune by homology search. Its function was confirmed by the complementation of the alg3 strain of S. cerevisiae . Inactivation of alg3 in S. commune resulted in the production of predominantly Man 3 GlcNAc 2 protein-linked N -glycans. No impact on growth nor a developmental phenotype due to the deletion was observed. This provides a first step toward engineering of a homogeneous, humanized N-glycosylation pattern for the production of therapeutic glycoproteins in mushrooms.

Beschreibung: We previously demonstrated that Siglec-15, a member of the Siglec family of glycan-recognition proteins, is expressed on a subset of macrophages and preferentially recognizes the sialyl-Tn (sTn) antigen, a tumor-associated glycan structure. In this study, we report on the biological significance of the Siglec-15-mediated interaction between monocytes/macrophages and cancer cells. Siglec-15 is expressed on tumor-associated macrophages (TAMs) in various human tumor tissues. We further demonstrated that its expression is substantially elevated in macrophage colony-stimulating factor-induced M2-like macrophages, which produced more transforming growth factor-β (TGF-β) in response to sTn-positive cells than to negative cells. We designed a co-culture model of THP-1 (human monocytic leukemia) cells and H157 (human lung carcinoma) cells mimicking the interaction between monocytes/macrophages and cancer cells that recapitulated the enhanced TGF-β production in Siglec-15 expressing THP-1 cells by the cellular interaction with sTn expressing H157 cells. The enhanced TGF-β production required a direct interaction between the two cell lines through sialic acids. Siglec-15 associates with adaptor protein DNAX activation protein of 12 kDa (DAP12) at the binding determinant Lys 274 in the transmembrane domain and transduces a signal to spleen tyrosine kinase (Syk). The enhanced TGF-β secretion was significantly attenuated by Syk inhibitor treatment of THP-1 cells or by substitution of the Siglec-15 Lys 274 to Ala, which disrupts the molecular interaction between Siglec15 and DAP12. These findings indicate that Siglec-15 recognizes the tumoral sTn antigen and transduces a signal for enhanced TGF-β secretion in TAMs and further suggest that Siglec-15 on macrophages may contribute to tumor progression by the TGF-β-mediated modulation of intratumoral microenvironments.

Beschreibung: Frontonasal dysplasia (FND) refers to a class of midline facial malformations caused by abnormal development of the facial primordia. The term encompasses a spectrum of severities but characteristic features include combinations of ocular hypertelorism, malformations of the nose and forehead and clefting of the facial midline. Several recent studies have drawn attention to the importance of Alx homeobox transcription factors during craniofacial development. Most notably, loss of Alx1 has devastating consequences resulting in severe orofacial clefting and extreme microphthalmia. In contrast, mutations of Alx3 or Alx4 cause milder forms of FND. Whilst Alx1 , Alx3 and Alx4 are all known to be expressed in the facial mesenchyme of vertebrate embryos, little is known about the function of these proteins during development. Here, we report the establishment of a zebrafish model of Alx -related FND. Morpholino knock-down of zebrafish alx1 expression causes a profound craniofacial phenotype including loss of the facial cartilages and defective ocular development. We demonstrate for the first time that Alx1 plays a crucial role in regulating the migration of cranial neural crest (CNC) cells into the frontonasal primordia. Abnormal neural crest migration is coincident with aberrant expression of foxd3 and sox10 , two genes previously suggested to play key roles during neural crest development, including migration, differentiation and the maintenance of progenitor cells. This novel function is specific to Alx1, and likely explains the marked clinical severity of Alx1 mutation within the spectrum of Alx -related FND.

Beschreibung: Activating somatic and germline mutations of closely related RAS genes (H, K, N) have been found in various types of cancer and in patients with developmental disorders, respectively. The involvement of the RAS signalling pathways in developmental disorders has recently emerged as one of the most important drivers in RAS research. In the present study, we investigated the biochemical and cell biological properties of two novel missense KRAS mutations (Y71H and K147E). Both mutations affect residues that are highly conserved within the RAS family. KRAS Y71H showed no clear differences to KRAS wt , except for an increased binding affinity for its major effector, the RAF1 kinase. Consistent with this finding, even though we detected similar levels of active KRAS Y71H when compared with wild-type protein, we observed an increased activation of MEK1/2, irrespective of the stimulation conditions. In contrast, KRAS K147E exhibited a tremendous increase in nucleotide dissociation generating a self-activating RAS protein that can act independently of upstream signals. As a consequence, levels of active KRAS K147E were strongly increased regardless of serum stimulation and similar to the oncogenic KRAS G12V . In spite of this, KRAS K147E downstream signalling did not reach the level triggered by oncogenic KRAS G12V , especially because KRAS K147E was downregulated by RASGAP and moreover exhibited a 2-fold lower affinity for RAF kinase. Here, our findings clearly emphasize that individual RAS mutations, despite being associated with comparable phenotypes of developmental disorders in patients, can cause remarkably diverse biochemical effects with a common outcome, namely a rather moderate gain-of-function.

Beschreibung: Spinocerebellar ataxia type 2 (SCA2) is an autosomal dominantly inherited disorder, which is caused by a pathological expansion of a polyglutamine (polyQ) tract in the coding region of the ATXN2 gene. Like other ataxias, SCA2 most overtly affects Purkinje cells (PCs) in the cerebellum. Using a transgenic mouse model expressing a full-length ATXN2 Q127 -complementary DNA under control of the Pcp2 promoter (a PC-specific promoter), we examined the time course of behavioral, morphologic, biochemical and physiological changes with particular attention to PC firing in the cerebellar slice. Although motor performance began to deteriorate at 8 weeks of age, reductions in PC number were not seen until after 12 weeks. Decreases in the PC firing frequency first showed at 6 weeks and paralleled deterioration of motor performance with progression of disease. Transcription changes in several PC-specific genes such as Calb1 and Pcp2 mirrored the time course of changes in PC physiology with calbindin-28 K changes showing the first small, but significant decreases at 4 weeks. These results emphasize that in this model of SCA2, physiological and behavioral phenotypes precede morphological changes by several weeks and provide a rationale for future studies examining the effects of restoration of firing frequency on motor function and prevention of future loss of PCs.

Beschreibung: Balancing selection has maintained human leukocyte antigen (HLA) allele diversity, but it is unclear whether this selection is symmetric (all heterozygotes are comparable and all homozygotes are comparable in terms of fitness) or asymmetric (distinct heterozygote genotypes display greater fitness than others). We tested the hypothesis that HLA is under asymmetric balancing selection in populations by estimating allelic branch lengths from genetic sequence data encoding peptide-binding domains. Significant deviations indicated changes in the ratio of terminal to internal branch lengths. Such deviations could arise even if no individual alleles present a strikingly altered branch length (e.g. if there is an overall distortion, with all or many terminal branches being longer than expected). DQ and DP loci were also analyzed as haplotypes. Using allele frequencies for 419 distinct populations in 10 geographical regions, we examined population differentiation in alleles within and between regions, and the relationship between allelic branch length and frequency. The strongest evidence for asymmetrical balancing selection was observed for HLA-DRB1 , HLA-B and HLA-DPA1 , with significant deviation ( P ≤ 1.1 x 10 –4 ) in about half of the populations. There were significant results at all loci except HLA-DQB1 / DQA1 . We observed moderate genetic variation within and between geographic regions, similar to the rest of the genome. Branch length was not correlated with allele frequency. In conclusion, sequence data suggest that balancing selection in HLA is asymmetric (some heterozygotes enjoy greater fitness than others). Because HLA polymorphism is crucial for pathogen resistance, this may manifest as a frequency-dependent selection with fluctuation in the fitness of specific heterozygotes over time.

Beschreibung: Birt–Hogg–Dubé syndrome (BHD) is a human cancer disorder caused by mutations in the tumor suppressor gene Folliculin ( FLCN ) with unknown biological functions. Here, we show that the Drosophila homolog of FLCN, dFLCN (a.k.a. dBHD ) localizes to the nucleolus and physically interacts with the 19S proteasomal ATPase, Rpt4, a nucleolar resident and known regulator of rRNA transcription. Downregulation of dFLCN resulted in an increase in nucleolar volume and upregulation of rRNA synthesis, whereas dFLCN overexpression reduced rRNA transcription and counteracted the effects of Rpt4 on rRNA production by preventing the association of Rpt4 with the rDNA locus. We further show that human FLCN exhibited evolutionarily conserved function and that Rpt4 knockdown inhibits the growth of FLCN-deficient human renal cancer cells in mouse xenografts. Our study suggests that FLCN functions as a tumor suppressor by negatively regulating rRNA synthesis.

Beschreibung: KitL, via its receptor cKit, supports primordial germ cell (PGC) growth, survival, migration and reprogramming to pluripotent embryonic germ cells (EGCs). However, the signaling downstream of KitL and its regulation in PGCs remain unclear. A constitutively activating mutation, cKit V558 , causes gain-of-function phenotypes in mast cells and intestines, and gastrointestinal stromal tumors (GISTs) when heterozygous. Unexpectedly, we find that PGC growth is not significantly affected in cKit V558 heterozygotes, whereas in homozygotes, increased apoptosis and inefficient migration lead to the depletion of PGCs. Through genetic studies, we reveal that this oncogenic cKit allele exhibits loss-of-function behavior in PGCs distinct from that in GIST development. Examination of downstream signaling in GISTs from cKit V558/+ mice confirmed hyperphosphorylation of AKT and ERK, but both remain unperturbed in cKit V558/+ PGCs and EGCs. In contrast, we find reduced activation of ERK1/2 and JNK1 in cKit V558 homozygous PGCs and EGCs. Inhibiting JNK, though not ERK1/2, increased apoptosis of wild-type PGCs, but did not further affect the already elevated apoptosis of cKit V558 / V558 PGCs. These results demonstrate a cell-context-dependent response to the cKit V558 mutation. We propose that AKT overload protection and JNK-mediated survival comprise PGC-specific mechanisms for regulating cKit signaling.

Beschreibung: TDP-43 is an evolutionarily conserved RNA-binding protein currently under intense investigation for its involvement in the molecular pathogenesis of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 is normally localized in the nucleus, but translocated to the cytoplasm in diseased neurons. The endogenous functions of TDP-43 in the nervous system remain poorly understood. Here, we show that the loss of Drosophila TDP-43 (dTDP-43) results in an increased production of sensory bristles and sensory organ precursor (SOP) cells on the notum of some but not all flies. The location of ectopic SOPs varies among mutant flies. The penetrance of this novel phenotype is dependent on the gender and sensitive to environmental influences. A similar SOP phenotype was also observed on the wing and in the embryos. Overexpression of dTDP-43 causes both loss and ectopic production of SOPs. Ectopic expression of ALS-associated mutant human TDP-43 (hTDP-43 M337V and hTDP-43 Q331K ) produces a less severe SOP phenotype than hTDP-43 WT , indicating a partial loss of function of mutant hTDP-43. In dTDP-43 mutants, miR-9a expression is significantly reduced. Genetic interaction studies further support the notion that dTDP-43 acts through miR-9a to control the precision of SOP specification. These findings reveal a novel role for endogenous TDP-43 in neuronal specification and suggest that the FTD/ALS-associated RNA-binding protein TDP-43 functions to ensure the robustness of genetic control programs.

Beschreibung: Mutations in leucine-rich repeat kinase 2 (LRRK2) are the most frequent known cause of late-onset Parkinson's disease (PD). To explore the therapeutic potential of small molecules targeting the LRRK2 kinase domain, we characterized two LRRK2 kinase inhibitors, TTT-3002 and LRRK2-IN1, for their effects against LRRK2 activity in vitro and in Caenorhabditis elegans models of LRRK2-linked neurodegeneration. TTT-3002 and LRRK2-IN1 potently inhibited in vitro kinase activity of LRRK2 wild-type and mutant proteins, attenuated phosphorylation of cellular LRRK2 and rescued neurotoxicity of mutant LRRK2 in transfected cells. To establish whether LRRK2 kinase inhibitors can mitigate pathogenesis caused by different mutations including G2019S and R1441C located within and outside of the LRRK2 kinase domain, respectively, we evaluated effects of TTT-3002 and LRRK2-IN1 against R1441C- and G2019S-induced neurodegeneration in C. elegans models. TTT-3002 and LRRK2-IN1 rescued the behavioral deficit characteristic of dopaminergic impairment in transgenic C. elegans expressing human R1441C- and G2019S-LRRK2. The inhibitors displayed nanomolar to low micromolar rescue potency when administered either pre-symptomatically or post-symptomatically, indicating both prevention and reversal of the dopaminergic deficit. The same treatments also led to long-lasting prevention and rescue of neurodegeneration. In contrast, TTT-3002 and LRRK2-IN1 were ineffective against the neurodegenerative phenotype in transgenic worms carrying the inhibitor-resistant A2016T mutation of LRRK2, suggesting that they elicit neuroprotective effects in vivo by targeting LRRK2 specifically. Our findings indicate that the LRRK2 kinase activity is critical for neurodegeneration caused by R1441C and G2019S mutations, suggesting that kinase inhibition of LRRK2 may represent a promising therapeutic strategy for PD.

Beschreibung: Mitochondrial DNA (mtDNA) mutations leading to the disruption of respiratory complex I (CI) have been shown to exhibit anti-tumorigenic effects, at variance with those impairing only the function but not the assembly of the complex, which appear to contribute positively to cancer development. Owing to the challenges in the analysis of the multi-copy mitochondrial genome, it is yet to be determined whether tumour-associated mtDNA lesions occur as somatic modifying factors or as germ-line predisposing elements. Here we investigated the whole mitochondrial genome sequence of 20 pituitary adenomas with oncocytic phenotype and identified pathogenic and/or novel mtDNA mutations in 60% of the cases. Using highly sensitive techniques, namely fluorescent PCR and allele-specific locked nucleic acid quantitative PCR, we identified the most likely somatic nature of these mutations in our sample set, since none of the mutations was detected in the corresponding blood tissue of the patients analysed. Furthermore, we have subjected a series of 48 pituitary adenomas to a high-resolution array comparative genomic hybridization analysis, which revealed that CI disruptive mutations, and the oncocytic phenotype, significantly correlate with low number of chromosomal aberrations in the nuclear genome. We conclude that CI disruptive mutations in pituitary adenomas are somatic modifiers of tumorigenesis most likely contributing not only to the development of oncocytic change, but also to a less aggressive tumour phenotype, as indicated by a stable karyotype.

Beschreibung: Functional loss of SMN1 causes proximal spinal muscular atrophy (SMA), the most common genetic condition accounting for infant lethality. Hence, the hypomorphic copy gene SMN2 is the only resource of functional SMN protein in SMA patients and influences SMA severity in a dose-dependent manner. Consequently, current therapeutic approaches focus on SMN2 . Histone deacetylase inhibitors (HDACi), such as the short chain fatty acid VPA (valproic acid), ameliorate the SMA phenotype by activating the SMN2 expression. By analyzing blood SMN2 expression in 16 VPA-treated SMA patients, about one-third of individuals were identified as positive responders presenting increased SMN2 transcript levels. In 66% of enrolled patients, a concordant response was detected in the respective fibroblasts. Most importantly, by taking the detour of reprograming SMA patients' fibroblasts, we showed that the VPA response was maintained even in GABAergic neurons derived from induced pluripotent stem cells (iPS) cells. Differential expression microarray analysis revealed a complete lack of response to VPA in non-responders, which was associated with an increased expression of the fatty acid translocase CD36. The pivotal role of CD36 as the cause of non-responsiveness was proven in various in vitro approaches. Most importantly, knockdown of CD36 in SMA fibroblasts converted non- into pos-responders. In summary, the concordant response from blood to the central nervous system (CNS) to VPA may allow selection of pos-responders prior to therapy. Increased CD36 expression accounts for VPA non-responsiveness. These findings may be essential not only for SMA but also for other diseases such as epilepsy or migraine frequently treated with VPA.

Beschreibung: Rett syndrome (RTT) is a neurodevelopmental disorder caused primarily by mutations of the X-linked MECP2 gene. Although the loss of MeCP2 function affects many neural systems, impairments of catecholaminergic function have been hypothesized to underlie several of the cardinal behavioral deficits of RTT patients and Mecp2-deficient mice. Although recent Mecp2 reactivation studies indicate that RTT may be a reversible condition, it remains unclear whether specifically preserving Mecp2 function within a specific system will be sufficient to convey beneficial effects. Here, we test whether the selective preservation of Mecp2 within catecholaminergic cells will improve the phenotype of Mecp2-deficient mice. Our results show that this targeted preservation of Mecp2 significantly improves the lifespan, phenotypic severity and cortical epileptiform discharge activity of both male and female Mecp2-deficient mice. Further, we found that the catecholaminergic preservation of Mecp2 also improves the ambulatory rate, rearing activity, motor coordination, anxiety and nest-building performances of Mecp2-deficient mice of each gender. Interestingly, our results also revealed a gender-specific improvement, as specific cortical and hippocampal electroencephalographic abnormalities were significantly improved in male, but not female, rescue mice. Collectively, these results support the role of the catecholaminergic system in the pathogenesis of RTT and provide proof-of-principle that restoring MeCP2 function within this specific system could represent a treatment strategy for RTT.

Beschreibung: Mutations in COL4A1 have been identified in families with hereditary small vessel disease of the brain presumably due to a dominant-negative mechanism. Here, we report on two novel mutations in COL4A1 in two families with porencephaly, intracerebral hemorrhage and severe white matter disease caused by haploinsufficiency. Two families with various clinical presentations of cerebral microangiopathy and autosomal dominant inheritance were examined. Clinical, neuroradiological and genetic investigations were performed. Electron microscopy of the skin was also performed. In one of the families, sequence analysis revealed a one base deletion, c.2085del, leading to a frameshift and a premature stopcodon, p.(Gly696fs). In the other family, a splice site mutation was identified, c.2194-1G〉A, which most likely leads to skipping of an exon with a frameshift and premature termination as a result. In fibroblasts of affected individuals from both the families, nonsense-mediated decay (NMD) of the mutant COL4A1 messenger RNAs (mRNAs) and a clear reduction of COL4A1 protein expression were demonstrated, indicating haploinsufficiency of COL4A1. Moreover, thickening of the capillary basement membrane in the skin was documented, similar to reports in patients with COL4A1 missense mutations. These findings suggest haploinsufficiency, a different mechanism from the commonly assumed dominant-negative effect, for COL4A1 mutations as a cause of (antenatal) intracerebral hemorrhage and white matter disease.

Beschreibung: Motivation: Given the current costs of next-generation sequencing, large studies carry out low-coverage sequencing followed by application of methods that leverage linkage disequilibrium to infer genotypes. We propose a novel method that assumes study samples are sequenced at low coverage and genotyped on a genome-wide microarray, as in the 1000 Genomes Project (1KGP). We assume polymorphic sites have been detected from the sequencing data and that genotype likelihoods are available at these sites. We also assume that the microarray genotypes have been phased to construct a haplotype scaffold. We then phase each polymorphic site using an MCMC algorithm that iteratively updates the unobserved alleles based on the genotype likelihoods at that site and local haplotype information. We use a multivariate normal model to capture both allele frequency and linkage disequilibrium information around each site. When sequencing data are available from trios, Mendelian transmission constraints are easily accommodated into the updates. The method is highly parallelizable, as it analyses one position at a time. Results: We illustrate the performance of the method compared with other methods using data from Phase 1 of the 1KGP in terms of genotype accuracy, phasing accuracy and downstream imputation performance. We show that the haplotype panel we infer in African samples, which was based on a trio-phased scaffold, increases downstream imputation accuracy for rare variants (R2 increases by 〉0.05 for minor allele frequency 〈1%), and this will translate into a boost in power to detect associations. These results highlight the value of incorporating microarray genotypes when calling variants from next-generation sequence data. Availability: The method (called MVNcall) is implemented in a C++ program and is available from http://www.stats.ox.ac.uk/~marchini/#software . Contact: marchini@stats.ox.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Accurate alignment of high-throughput RNA-seq data is a challenging and yet unsolved problem because of the non-contiguous transcript structure, relatively short read lengths and constantly increasing throughput of the sequencing technologies. Currently available RNA-seq aligners suffer from high mapping error rates, low mapping speed, read length limitation and mapping biases. Results: To align our large (〉80 billon reads) ENCODE Transcriptome RNA-seq dataset, we developed the Spliced Transcripts Alignment to a Reference (STAR) software based on a previously undescribed RNA-seq alignment algorithm that uses sequential maximum mappable seed search in uncompressed suffix arrays followed by seed clustering and stitching procedure. STAR outperforms other aligners by a factor of 〉50 in mapping speed, aligning to the human genome 550 million 2 x 76 bp paired-end reads per hour on a modest 12-core server, while at the same time improving alignment sensitivity and precision. In addition to unbiased de novo detection of canonical junctions, STAR can discover non-canonical splices and chimeric (fusion) transcripts, and is also capable of mapping full-length RNA sequences. Using Roche 454 sequencing of reverse transcription polymerase chain reaction amplicons, we experimentally validated 1960 novel intergenic splice junctions with an 80–90% success rate, corroborating the high precision of the STAR mapping strategy. Availability and implementation: STAR is implemented as a standalone C++ code. STAR is free open source software distributed under GPLv3 license and can be downloaded from http://code.google.com/p/rna-star/ . Contact: dobin@cshl.edu .

Beschreibung: Motivation: Local motifs are patterns of DNA or protein sequences that occur within a sequence interval relative to a biologically defined anchor or landmark. Current protein motif discovery methods do not adequately consider such constraints to identify biologically significant motifs that are only weakly over-represented but spatially confined. Using negatives, i.e. sequences known to not contain a local motif, can further increase the specificity of their discovery. Results: This article introduces the method DLocalMotif that makes use of positional information and negative data for local motif discovery in protein sequences. DLocalMotif combines three scoring functions, measuring degrees of motif over-representation, entropy and spatial confinement, specifically designed to discriminatively exploit the availability of negative data. The method is shown to outperform current methods that use only a subset of these motif characteristics. We apply the method to several biological datasets. The analysis of peroxisomal targeting signals uncovers several novel motifs that occur immediately upstream of the dominant peroxisomal targeting signal-1 signal. The analysis of proline-tyrosine nuclear localization signals uncovers multiple novel motifs that overlap with C2H2 zinc finger domains. We also evaluate the method on classical nuclear localization signals and endoplasmic reticulum retention signals and find that DLocalMotif successfully recovers biologically relevant sequence properties. Availability: http://bioinf.scmb.uq.edu.au/dlocalmotif/ Contact: m.boden@uq.edu.au Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Much evidence show that over-expression of epidermal growth factor receptor (EGFR) plays an important role in regulating carcinogenesis. Genetic variations in 3' untranslated region (3'UTR) of gene have been reported to affect gene expression by interfering with microRNAs (miRNAs), which are thought to function as either tumour suppressors or oncogenes by binding to their target mRNA. In this study, we investigated the association between the EGFR 3'UTR 774T〉C polymorphism and bladder cancer risk. We used the TaqMan technology to genotype this genetic variant in a hospital-based case–control study of 908 bladder cancer patients and 1239 controls in a Chinese population. We found that the 774CC genotype was associated with a statistically significantly increased risk of bladder cancer [adjusted odds ratio = 1.29, 95% confidence interval = 1.05–1.58], compared with the 774TT/TC genotype, and this increased risk was more pronounced among subgroups of age 〉 65 years, non-smokers and patients’ tumour invasive stage. Furthermore, luciferase assays in T24 cell showed that EGFR 3'UTR 774 T to C substitution could increase the expression of EGFR, which was consistent with the association study finding. Additionally, we also provide evidence that 774T〉C polymorphism increasing EGFR expression was not regulated by hsa-miR-214 binding. These findings suggested that EGFR 3'UTR 774T〉C polymorphism may contribute to susceptibility to bladder cancer.

Beschreibung: Motivation: Protein–protein interaction (PPI) plays an important role in understanding gene functions, and many computational PPI prediction methods have been proposed in recent years. Despite the extensive efforts, PPI prediction still has much room to improve. Sequence-based co-evolution methods include the substitution rate method and the mirror tree method, which compare sequence substitution rates and topological similarity of phylogenetic trees, respectively. Although they have been used to predict PPI in species with small genomes like Escherichia coli , such methods have not been tested in large scale proteome like Homo sapiens . Result: In this study, we propose a novel sequence-based co-evolution method, co-evolutionary divergence (CD), for human PPI prediction. Built on the basic assumption that protein pairs with similar substitution rates are likely to interact with each other, the CD method converts the evolutionary information from 14 species of vertebrates into likelihood ratios and combined them together to infer PPI. We showed that the CD method outperformed the mirror tree method in three independent human PPI datasets by a large margin. With the arrival of more species genome information generated by next generation sequencing, the performance of the CD method can be further improved. Availability: Source code and support are available at http://mib.stat.sinica.edu.tw/LAP/tmp/CD.rar . Contact: syuan@stat.sinica.edu.tw Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: : In higher eukaryotes, the identification of translation initiation sites (TISs) has been focused on finding these signals in cDNA or mRNA sequences. Using Arabidopsis thaliana ( A.t. ) information, we developed a prediction tool for signals within genomic sequences of plants that correspond to TISs. Our tool requires only genome sequence, not expressed sequences. Its sensitivity/specificity is for A.t. (90.75%/92.2%), for Vitis vinifera (66.8%/94.4%) and for Populus trichocarpa (81.6%/94.4%), which suggests that our tool can be used in annotation of different plant genomes. We provide a list of features used in our model. Further study of these features may improve our understanding of mechanisms of the translation initiation. Availability and implementation: Our tool is implemented as an artificial neural network. It is available as a web-based tool and, together with the source code, the list of features, and data used for model development, is accessible at http://cbrc.kaust.edu.sa/dts . Contact: vladimir.bajic@kaust.edu.sa Supplementary information : Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: Pathway or gene set analysis has been widely applied to genomic data. Many current pathway testing methods use univariate test statistics calculated from individual genomic markers, which ignores the correlations and interactions between candidate markers. Random forests-based pathway analysis is a promising approach for incorporating complex correlation and interaction patterns, but one limitation of previous approaches is that pathways have been considered separately, thus pathway cross-talk information was not considered. Results: In this article, we develop a new pathway hunting algorithm for survival outcomes using random survival forests, which prioritize important pathways by accounting for gene correlation and genomic interactions. We show that the proposed method performs favourably compared with five popular pathway testing methods using both synthetic and real data. We find that the proposed methodology provides an efficient and powerful pathway modelling framework for high-dimensional genomic data. Availability: The R code for the analysis used in this article is available upon request. Contact: xi.steven.chen@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Motivation: PacBio sequencers produce two types of characteristic reads (continuous long reads: long and high error rate and circular consensus sequencing: short and low error rate), both of which could be useful for de novo assembly of genomes. Currently, there is no available simulator that targets the specific generation of PacBio libraries. Results: Our analysis of 13 PacBio datasets showed characteristic features of PacBio reads (e.g. the read length of PacBio reads follows a log-normal distribution). We have developed a read simulator, PBSIM, that captures these features using either a model-based or sampling-based method. Using PBSIM, we conducted several hybrid error correction and assembly tests for PacBio reads, suggesting that a continuous long reads coverage depth of at least 15 in combination with a circular consensus sequencing coverage depth of at least 30 achieved extensive assembly results. Availability: PBSIM is freely available from the web under the GNU GPL v2 license ( http://code.google.com/p/pbsim/ ). Contact: mhamada@k.u-tokyo.ac.jp Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: : Drugster is a fully interactive pipeline designed to break the command line barrier and introduce a new user-friendly environment to perform drug design, lead and structure optimization experiments through an efficient combination of the PDB2PQR, Ligbuilder, Gromacs and Dock suites. Our platform features a novel workflow that guides the user through each logical step of the iterative 3D structural optimization setup and drug design process, by providing a seamless interface to all incorporated packages. Availability: Drugster can be freely downloaded via our dedicated server system at http://www.bioacademy.gr/bioinformatics/drugster/ . Contact: dvlachakis@bioacademy.gr .

Beschreibung: XiP (eXtensible integrative Pipeline) is a flexible, editable and modular environment with a user-friendly interface that does not require previous advanced programming skills to run, construct and edit workflows. XiP allows the construction of workflows by linking components written in both R and Java, the analysis of high-throughput data in grid engine systems and also the development of customized pipelines that can be encapsulated in a package and distributed. XiP already comes with several ready-to-use pipeline flows for the most common genomic and transcriptomic analysis and ~300 computational components. Availability: XiP is open source, freely available under the Lesser General Public License (LGPL) and can be downloaded from http://xip.hgc.jp . Contact: nagasaki@megabank.tohoku.ac.jp

Beschreibung: Existing repositories for experimental datasets typically capture snapshots of data acquired using a single experimental technique and often require manual population and continual curation. We present a storage system for heterogeneous research data that performs dynamic automated indexing to provide powerful search, discovery and collaboration features without the restrictions of a structured repository. ADAM is able to index many commonly used file formats generated by laboratory assays and therefore offers specific advantages to the experimental biology community. However, it is not domain specific and can promote sharing and re-use of working data across scientific disciplines. Availability and implementation: ADAM is implemented using Java and supported on Linux. It is open source under the GNU General Public License v3.0. Installation instructions, binary code, a demo system and virtual machine image and are available at http://www.imperial.ac.uk/bioinfsupport/resources/software/adam . Contact: m.woodbridge@imperial.ac.uk

Beschreibung: : Drug versus Disease (DvD) provides a pipeline, available through R or Cytoscape, for the comparison of drug and disease gene expression profiles from public microarray repositories. Negatively correlated profiles can be used to generate hypotheses of drug-repurposing, whereas positively correlated profiles may be used to infer side effects of drugs. DvD allows users to compare drug and disease signatures with dynamic access to databases Array Express, Gene Expression Omnibus and data from the Connectivity Map. Availability and implementation: R package (submitted to Bioconductor) under GPL 3 and Cytoscape plug-in freely available for download at www.ebi.ac.uk/saezrodriguez/DVD/ . Contact: saezrodriguez@ebi.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Beschreibung: Angiogenesis, a formation of neovessels, is regulated by the local balance between angiogenesis stimulators and inhibitors. A number of such endogenous regulators of angiogenesis have been found in the body. Recently, vasohibin-1 (VASH1) was isolated as a negative feedback regulator of angiogenesis produced by endothelial cells (ECs) and subsequently vasohibin-2 (VASH2) as a homologue of VASH1. It was then explored that VASH1 is expressed in ECs to terminate angiogenesis, whereas VASH2 is expressed in cells other than ECs to promote angiogenesis in the mouse model of angiogenesis. This review will focus on the vasohibin family members, which are novel regulators of angiogenesis.

Beschreibung: Tertiary dentin is deposited inside teeth after various stimuli and serves as a major defensive wall to preserve pulp cells. However, the molecular mechanisms of the activation of quiescent odontoblasts, immature pulp cells and tertiary dentin formation are still unclear. Therefore, we performed a comprehensive gene expression analysis of pulp cells after cavity preparation of 9-week-old rat molars to clarify the critical molecules in tertiary dentinogenesis. As a result, mRNA expression of various molecules was up- or down-regulated. Notably, several members of the matrix metalloprotease family and their endogenous inhibitors were up-regulated after cavity preparation. In situ hybridization showed that tissue inhibitor of metalloprotease 1 ( Timp1 ) was widely and continuously distributed in the pulp beneath the cavity in vivo . We also observed accumulation of β-catenin in the pulp cells beneath the cavity by fluorescence immunohistochemistry. Furthermore, Timp1 transcription was repressed by a dominant-negative TCF4 in immature undifferentiated mesenchymal cells, but not altered in mature odontoblast-like cells. These results indicate that cavity preparation may activate the Wnt/β-catenin pathway and the Wnt/β-catenin pathway and Timp1 may be correlatively involved in pulp repair. Timp1 might play crucial roles in reactivation of immature pulp cells for tertiary dentinogenesis.