Effects of spaced k-mers on alignment-free genotyping.

MOTIVATION: Alignment-free, k-mer based genotyping methods are a fast alternative to alignment-based methods and are particularly well suited for genotyping larger cohorts. The sensitivity of algorithms, that work with k-mers, can be increased by using spaced seeds, however, the application of spaced seeds in k-mer based genotyping methods has not been researched yet. RESULTS: We add a spaced seeds functionality to the genotyping software PanGenie and use it to calculate genotypes. This significantly improves sensitivity and F-score when genotyping SNPs, indels, and structural variants on reads with low (5×) and high (30×) coverage. Improvements are greater than what could be achieved by just increasing the length of contiguous k-mers. Effect sizes are particularly large for low coverage data. If applications implement effective algorithms for hashing of spaced k-mers, spaced k-mers have the potential to become an useful technique in k-mer based genotyping. AVAILABILITY AND IMPLEMENTATION: The source code of our proposed tool MaskedPanGenie is openly available on https://github.com/hhaentze/MaskedPangenie.

Long non-coding RNA 2310069B03Rik functions as a suppressor of Ucp1 expression under prolonged cold exposure in murine beige adipocytes.

Specific conditions, such as exposure to cold, can induce the production of brown-like adipocytes in white adipose tissue. These adipocytes express high levels of uncoupling protein 1 (UCP1) and energy expended by generating heat. Thus, these are a potential target for the prevention or treatment of obesity. The present study involved a comprehensive analysis of the adipose tissue to understand the relationship between long non-coding RNA (lncRNA) 2310069B03Rik and UCP1. Cold exposure increased both lncRNA 2310069B03Rik and Ucp1 expression in inguinal white adipose tissue (iWAT). However, overexpression of lncRNA 2310069B03Rik suppressed the Ucp1 mRNA expression and the promoter activity of UCP1 in the iWAT primary adipocytes. In addition, compared to the early induction of Ucp1 expression by cold stimulation, the induction of lncRNA 2310069B03Rik expression was later. These results suggest that lncRNA 2310069B03Rik functions as a suppression factor of Ucp1 expression.

Sequence homology in eukaryotes (SHOE): interactive visual tool for promoter analysis.

BACKGROUND: Microarray and DNA-sequencing based technologies continue to produce enormous amounts of data on gene expression. This data has great potential to illuminate our understanding of biology and medicine, but the data alone is of limited value without computational tools to allow human investigators to visualize and interpret it in the context of their problem of interest. RESULTS: We created a web server called SHOE that provides an interactive, visual presentation of the available evidence of transcriptional regulation and gene co-expression to facilitate its exploration and interpretation. SHOE predicts the likely transcription factor binding sites in orthologous promoters of humans, mice, and rats using the combined information of 1) transcription factor binding preferences (position-specific scoring matrix (PSSM) libraries such as Transfac32, Jaspar, HOCOMOCO, ChIP-seq, SELEX, PBM, and iPS-reprogramming factor), 2) evolutionary conservation of putative binding sites in orthologous promoters, and 3) co-expression tendencies of gene pairs based on 1,714 normal human cells selected from the Gene Expression Omnibus Database. CONCLUSION: SHOE enables users to explore potential interactions between transcription factors and target genes via multiple data views, discover transcription factor binding motifs on top of gene co-expression, and visualize genes as a network of gene and transcription factors on its native gadget GeneViz, the CellDesigner pathway analyzer, and the Reactome database to search the pathways involved. As we demonstrate here when using the CREB1 and Nf-κB datasets, SHOE can reliably identify experimentally verified interactions and predict plausible novel ones, yielding new biological insights into the gene regulatory mechanisms involved. SHOE comes with a manual describing how to run it on a local PC or via the Garuda platform ( www.garuda-alliance.org ), where it joins other popular gadgets such as the CellDesigner pathway analyzer and the Reactome database, as part of analysis workflows to meet the growing needs of molecular biologists and medical researchers. SHOE is available from the following URL http://ec2-54-150-223-65.ap-northeast-1.compute.amazonaws.com A video demonstration of SHOE can be found here: https://www.youtube.com/watch?v=qARinNb9NtE.

A bioinformatics potpourri.

The 16th International Conference on Bioinformatics (InCoB) was held at Tsinghua University, Shenzhen from September 20 to 22, 2017. The annual conference of the Asia-Pacific Bioinformatics Network featured six keynotes, two invited talks, a panel discussion on big data driven bioinformatics and precision medicine, and 66 oral presentations of accepted research articles or posters. Fifty-seven articles comprising a topic assortment of algorithms, biomolecular networks, cancer and disease informatics, drug-target interactions and drug efficacy, gene regulation and expression, imaging, immunoinformatics, metagenomics, next generation sequencing for genomics and transcriptomics, ontologies, post-translational modification, and structural bioinformatics are the subject of this editorial for the InCoB2017 supplement issues in BMC Genomics, BMC Bioinformatics, BMC Systems Biology and BMC Medical Genomics. New Delhi will be the location of InCoB2018, scheduled for September 26-28, 2018.

MitoFates: improved prediction of mitochondrial targeting sequences and their cleavage sites.

Mitochondria provide numerous essential functions for cells and their dysfunction leads to a variety of diseases. Thus, obtaining a complete mitochondrial proteome should be a crucial step toward understanding the roles of mitochondria. Many mitochondrial proteins have been identified experimentally but a complete list is not yet available. To fill this gap, methods to computationally predict mitochondrial proteins from amino acid sequence have been developed and are widely used, but unfortunately, their accuracy is far from perfect. Here we describe MitoFates, an improved prediction method for cleavable N-terminal mitochondrial targeting signals (presequences) and their cleavage sites. MitoFates introduces novel sequence features including positively charged amphiphilicity, presequence motifs, and position weight matrices modeling the presequence cleavage sites. These features are combined with classical ones such as amino acid composition and physico-chemical properties as input to a standard support vector machine classifier. On independent test data, MitoFates attains better performance than existing predictors in both detection of presequences and in predicting their cleavage sites. We used MitoFates to look for undiscovered mitochondrial proteins from 42,217 human proteins (including isoforms such as alternative splicing or translation initiation variants). MitoFates predicts 1167 genes to have at least one isoform with a presequence. Five-hundred and eighty of these genes were not annotated as mitochondrial in either UniProt or Gene Ontology. Interestingly, these include candidate regulators of parkin translocation to damaged mitochondria, and also many genes with known disease mutations, suggesting that careful investigation of MitoFates predictions may be helpful in elucidating the role of mitochondria in health and disease. MitoFates is open source with a convenient web server publicly available.

A bioinformatician's guide to the forefront of suffix array construction algorithms.

The suffix array and its variants are text-indexing data structures that have become indispensable in the field of bioinformatics. With the uninitiated in mind, we provide an accessible exposition of the SA-IS algorithm, which is the state of the art in suffix array construction. We also describe DisLex, a technique that allows standard suffix array construction algorithms to create modified suffix arrays designed to enable a simple form of inexact matching needed to support 'spaced seeds' and 'subset seeds' used in many biological applications.

GIW and InCoB are advancing bioinformatics in the Asia-Pacific.

GIW/InCoB2015 the joint 26th International Conference on Genome Informatics (GIW) and 14th International Conference on Bioinformatics (InCoB) held in Tokyo, September 9-11, 2015 was attended by over 200 delegates. Fifty-one out of 89 oral presentations were based on research articles accepted for publication in four BMC journal supplements and three other journals. Sixteen articles in this supplement and six articles in the BMC Systems Biology GIW/InCoB2015 Supplement are covered by this introduction. The topics range from genome informatics, protein structure informatics, image analysis to biological networks and biomarker discovery.

GIW and InCoB, two premier bioinformatics conferences in Asia with a combined 40 years of history.

Knowledge discovery in bioinformatics thrives on joint and inclusive efforts of stakeholders. Similarly, knowledge dissemination is expected to be more effective and scalable through joint efforts. Therefore, the International Conference on Bioinformatics (InCoB) and the International Conference on Genome Informatics (GIW) were organized as a joint conference for the first time in 13 years of coexistence. The Asia-Pacific Bioinformatics Network (APBioNet) and the Japanese Society for Bioinformatics (JSBi) collaborated to host GIW/InCoB2015 in Tokyo, September 9-11, 2015. The joint endeavour yielded 51 research articles published in seven journals, 78 poster and 89 oral presentations, showcasing bioinformatics research in the Asia-Pacific region. Encouraged by the results and reduced organizational overheads, APBioNet will collaborate with other bioinformatics societies in organizing co-located bioinformatics research and training meetings in the future. InCoB2016 will be hosted in Singapore, September 21-23, 2016.

Multimodal functional and structural neuroimaging investigation of major depressive disorder following treatment with duloxetine.

BACKGROUND: Longitudinal neuroimaging studies of major depressive disorder (MDD) have most commonly assessed the effects of antidepressants from the serotonin reuptake inhibitor class and usually reporting a single measure. Multimodal neuroimaging assessments were acquired from MDD patients during an acute depressive episode with serial measures during a 12-week treatment with the serotonin-norepinephrine reuptake inhibitor (SNRI) duloxetine. METHODS: Participants were medication-free MDD patients (n = 32; mean age 40.2 years) in an acute depressive episode and healthy controls matched for age, gender, and IQ (n = 25; mean age 38.8 years). MDD patients received treatment with duloxetine 60 mg daily for 12 weeks with an optional dose increase to 120 mg daily after 8 weeks. All participants had serial imaging at weeks 0, 1, 8, and 12 on a 3 Tesla magnetic resonance imaging (MRI) scanner. Neuroimaging tasks included emotional facial processing, negative attentional bias (emotional Stroop), resting state functional MRI and structural MRI. RESULTS: A significant group by time interaction was identified in the anterior default mode network in which MDD patients showed increased connectivity with treatment, while there were no significant changes in healthy participants. In the emotional Stroop task, increased posterior cingulate activation in MDD patients normalized following treatment. No significant group by time effects were observed for happy or sad facial processing, including in amygdala responsiveness, or in regional cerebral volumes. Reduced baseline resting state connectivity within the orbitofrontal component of the default mode network was predictive of clinical response. An early increase in hippocampal volume was predictive of clinical response. CONCLUSIONS: Baseline resting state functional connectivity was predictive of subsequent clinical response. Complementary effects of treatment were observed from the functional neuroimaging correlates of affective facial expressions, negative attentional bias, and resting state. No significant effects were observed in affective facial processing, while the interaction effect in negative attentional bias and individual group effects in resting state connectivity could be related to the SNRI class of antidepressant medication. The specificity of the observed effects to SNRI pharmacological treatments requires further investigation. TRIAL REGISTRATION: Registered at clinicaltrials.gov ( NCT01051466 ).

Identification of cargo proteins specific for the nucleocytoplasmic transport carrier transportin by combination of an in vitro transport system and stable isotope labeling by amino acids in cell culture (SILAC)-based quantitative proteomics.

The human importin-ß family consists of 21 nucleocytoplasmic transport carrier proteins that carry proteins and RNAs across the nuclear envelope through nuclear pores in specific directions. These transport carriers are responsible for the nucleocytoplasmic transport of thousands of proteins, but the cargo allocation of each carrier, which is necessary information if one wishes to understand the physiological context of transport, is poorly characterized. To address this issue, we developed a high-throughput method to identify the cargoes of transport carriers by applying stable isotope labeling by amino acids in cell culture to construct an in vitro transport system. Our method can be outlined in three steps. (1) Cells are cultured in a medium containing a stable isotope. (2) The cell membranes of the labeled cells are permeabilized, and proteins extracted from unlabeled cells are transported into the nuclei of the permeabilized cells. In this step, the reaction system is first depleted of all importin-ß family carriers and then supplemented with a particular importin-ß family carrier of interest. (3) Proteins in the nuclei are extracted and analyzed quantitatively via LC-MS/MS. As an important test case, we used this method to identify cargo proteins of transportin, a representative member of the importin-ß family. As expected, the identified candidate cargo proteins included previously reported transportin cargoes as well as new potential cargoes, which we corroborated via in vitro binding assays. The identified cargoes are predominately RNA-interacting proteins, affirming that cargoes allotted to the same carrier share functional characteristics. Finally, we found that the transportin cargoes possessed at least two classes of signal sequences: the well characterized PY-nuclear localization signals specific for transportin, and Lys/Arg-rich segments capable of binding to both transportin and importin-ß. Thus, our method will be useful for linking a carrier to features shared among its cargoes and to specific nuclear localization signals.

ValidNESs: a database of validated leucine-rich nuclear export signals.

ValidNESs (http://validness.ym.edu.tw/) is a new database for experimentally validated leucine-rich nuclear export signal (NES)-containing proteins. The therapeutic potential of the chromosomal region maintenance 1 (CRM1)-mediated nuclear export pathway and disease relevance of its cargo proteins has gained recognition in recent years. Unfortunately, only about one-third of known CRM1 cargo proteins are accessible in a single database since the last compilation in 2003. CRM1 cargo proteins are often recognized by a classical NES (leucine-rich NES), but this signal is notoriously difficult to predict from sequence alone. Fortunately, a recently developed prediction method, NESsential, is able to identify good candidates in some cases, enabling valuable hints to be gained by in silico prediction, but until now it has not been available through a web interface. We present ValidNESs, an integrated, up-to-date database holding 221 NES-containing proteins, combined with a web interface to prediction by NESsential.

Contesting heteronormativity: the fight for lesbian, gay, bisexual and transgender recognition in India and Vietnam.

Recent public debates about sexuality in India and Vietnam have brought the rights of lesbian, gay, bisexual and transgender people sharply into focus. Drawing on legal documents, secondary sources and ethnographic fieldwork conducted in the urban centres of Delhi and Hanoi, this article shows how the efforts of civil society organisations dedicated to the fight for lesbian, gay, bisexual and transgender rights have had different consequences in these two Asian contexts. The paper considers how these organisations navigated government regulations about their formation and activities, as well as the funding priorities of national and international agencies. The HIV epidemic has had devastating consequences for gay men and other men who have sex with men, and has been highly stigmatising. As a sad irony, the epidemic has provided at the same time a strategic entry point for organisations to struggle for lesbian, gay, bisexual and transgender recognition. This paper examines how the fight for lesbian, gay, bisexual and transgender recognition has been doubly framed through health-based and rights-based approaches and how the struggle for recognition has positioned lesbian, gay, bisexual and transgender people in India and Vietnam differently.

Plus ça change - evolutionary sequence divergence predicts protein subcellular localization signals.

BACKGROUND: Protein subcellular localization is a central problem in understanding cell biology and has been the focus of intense research. In order to predict localization from amino acid sequence a myriad of features have been tried: including amino acid composition, sequence similarity, the presence of certain motifs or domains, and many others. Surprisingly, sequence conservation of sorting motifs has not yet been employed, despite its extensive use for tasks such as the prediction of transcription factor binding sites. RESULTS: Here, we flip the problem around, and present a proof of concept for the idea that the lack of sequence conservation can be a novel feature for localization prediction. We show that for yeast, mammal and plant datasets, evolutionary sequence divergence alone has significant power to identify sequences with N-terminal sorting sequences. Moreover sequence divergence is nearly as effective when computed on automatically defined ortholog sets as on hand curated ones. Unfortunately, sequence divergence did not necessarily increase classification performance when combined with some traditional sequence features such as amino acid composition. However a post-hoc analysis of the proteins in which sequence divergence changes the prediction yielded some proteins with atypical (i.e. not MPP-cleaved) matrix targeting signals as well as a few misannotations. CONCLUSION: We report the results of the first quantitative study of the effectiveness of evolutionary sequence divergence as a feature for protein subcellular localization prediction. We show that divergence is indeed useful for prediction, but it is not trivial to improve overall accuracy simply by adding this feature to classical sequence features. Nevertheless we argue that sequence divergence is a promising feature and show anecdotal examples in which it succeeds where other features fail.

Adaptive seeds tame genomic sequence comparison.

The main way of analyzing biological sequences is by comparing and aligning them to each other. It remains difficult, however, to compare modern multi-billionbase DNA data sets. The difficulty is caused by the nonuniform (oligo)nucleotide composition of these sequences, rather than their size per se. To solve this problem, we modified the standard seed-and-extend approach (e.g., BLAST) to use adaptive seeds. Adaptive seeds are matches that are chosen based on their rareness, instead of using fixed-length matches. This method guarantees that the number of matches, and thus the running time, increases linearly, instead of quadratically, with sequence length. LAST, our open source implementation of adaptive seeds, enables fast and sensitive comparison of large sequences with arbitrarily nonuniform composition.

RecountDB: a database of mapped and count corrected transcribed sequences.

The field of gene expression analysis continues to benefit from next-generation sequencing generated data, which enables transcripts to be measured with unmatched accuracy and resolution. But the high-throughput reads from these technologies also contain many errors, which can compromise the ability to accurately detect and quantify rare transcripts. Fortunately, techniques exist to ameliorate the affects of sequencer error. We present RecountDB, a secondary database derived from primary data in NCBI's short read archive. RecountDB holds sequence counts from RNA-seq and 5' capped transcription start site experiments, corrected and mapped to the relevant genome. Via a searchable and browseable interface users can obtain corrected data in formats useful for transcriptomic analysis. The database is currently populated with 2265 entries from 45 organisms and continuously growing. RecountDB is publicly available at: http://recountdb.cbrc.jp.

Mammalian NUMT insertion is non-random.

It is well known that remnants of partial or whole copies of mitochondrial DNA, known as Nuclear MiTochondrial sequences (NUMTs), are found in nuclear genomes. Since whole genome sequences have become available, many bioinformatics studies have identified putative NUMTs and from those attempted to infer the factors involved in NUMT creation. These studies conclude that NUMTs represent randomly chosen regions of the mitochondrial genome. There is less consensus regarding the nuclear insertion sites of NUMTs - previous studies have discussed the possible role of retrotransposons, but some recent ones have reported no correlation or even anti-correlation between NUMT sites and retrotransposons. These studies have generally defined NUMT sites using BLAST with default parameters. We analyze a redefined set of human NUMTs, computed with a carefully considered protocol. We discover that the inferred insertion points of NUMTs have a strong tendency to have high-predicted DNA curvature, occur in experimentally defined open chromatin regions and often occur immediately adjacent to A + T oligomers. We also show clear evidence that their flanking regions are indeed rich in retrotransposons. Finally we show that parts of the mitochondrial genome D-loop are under-represented as a source of NUMTs in primate evolution.

'I thought I was the only one': the misrecognition of LGBT youth in contemporary Vietnam.

While recent LGBT rights demonstrations and discussions about same-sex marriage have thrust the issue of homosexuality into the spotlight, it was not long ago that the issue of homosexuality was notable by its absence in Vietnam. Drawing on semi-structured interviews with young gay, lesbian and bisexual people in Vietnam's capital city Hanoi, this paper considers the increasing visibility of homosexuality through the theoretical lens of recognition, and illustrates the heterosexist misrecognition that LGBT young people have been subjected to in legislation, the media, their families, and through the education system. Drawing on the narratives of LGBT young people, the paper highlights the potentially negative impact such misrecognition may have on psychological and social wellbeing.

Paving the future: finding suitable ISMB venues.

The International Society for Computational Biology, ISCB, organizes the largest event in the field of computational biology and bioinformatics, namely the annual international conference on Intelligent Systems for Molecular Biology, the ISMB. This year at ISMB 2012 in Long Beach, ISCB celebrated the 20th anniversary of its flagship meeting. ISCB is a young, lean and efficient society that aspires to make a significant impact with only limited resources. Many constraints make the choice of venues for ISMB a tough challenge. Here, we describe those challenges and invite the contribution of ideas for solutions.