Locality-sensitive hashing for protein classification

Determination of sequence similarity is a central issue in computational biology, a problem addressed primarily through BLAST, an alignment based heuristic which has underpinned much of the analysis and annotation of the genomic era. Despite their success, alignment-based approaches scale poorly with increasing data set size, and are not robust under structural sequence rearrangements. Successive waves of innovation in sequencing technologies – so-called Next Generation Sequencing (NGS) approaches – have led to an explosion in data availability, challenging existing methods and motivating novel approaches to sequence representation and similarity scoring, including adaptation of existing methods from other domains such as information retrieval. In this work, we investigate locality-sensitive hashing of sequences through binary document signatures, applying the method to a bacterial protein classification task. Here, the goal is to predict the gene family to which a given query protein belongs. Experiments carried out on a pair of small but biologically realistic datasets (the full protein repertoires of families of Chlamydia and Staphylococcus aureus genomes respectively) show that a measure of similarity obtained by locality sensitive hashing gives highly accurate results while offering a number of avenues which will lead to substantial performance improvements over BLAST..

National working group meeting on ALK diagnostics in lung cancer

The global landscape of molecular testing is rapidly changing, with the recent publication of the International Association for the Study of Lung Cancer (IASLC)/College of American Pathologists (CAP) guidelines and the ALK Atlas. The IASLC/CAP guidelines recommend that tumors from patients with non-small cell lung cancer (NSCLC) be tested for ALK rearrangements in addition to epidermal growth factor receptor (EGFR) mutations. The spur for this recommendation is the availability of novel therapies that target these rearrangements. This article is based on coverage of a Pfizer-sponsored National Working Group Meeting on ALK Diagnostics in Lung Cancer, held around the 15th World Lung Cancer Conference, in Sydney on October 31, 2013. It is based on the presentations given by the authors at the meeting and the discussion that ensued. The content for this article was discussed and agreed on by the authors.

Double-strand breaks and the concept of short- and long-term epigenetic memory

Double-strand breaks represent an extremely cytolethal form of DNA damage and thus pose a serious threat to the preservation of genetic and epigenetic information. Though it is well-known that double-strand breaks such as those generated by ionising radiation are among the principal causative factors behind mutations, chromosomal aberrations, genetic instability and carcino-genesis, significantly less is known about the epigenetic consequences of double-strand break formation and repair for carcinogenesis. Double-strand break repair is a highly coordinated process that requires the unravelling of the compacted chromatin structure to facilitate repair machinery access and then restoration of the original undamaged chromatin state. Recent experimental findings have pointed to a potential mechanism for double-strand break-induced epigenetic silencing. This review will discuss some of the key epigenetic regulatory processes involved in double-strand break (DSB) repair and how incomplete or incorrect restoration of chromatin structure can leave a DSB-induced epigenetic memory of damage with potentially pathological repercussions

J-Circos : an interactive Circos plotter

Circos plots are graphical outputs that display three dimensional chromosomal interactions and fusion transcripts. However, the Circos plot tool is not an interactive visualization tool, but rather a figure generator. For example, it does not enable data to be added dynamically, nor does it provide information for specific data points interactively. Recently, an R-based Circos tool (RCircos) has been developed to integrate Circos to R, but similarly, Rcircos can only be used to generate plots. Thus, we have developed a Circos plot tool (J-Circos) that is an interactive visualization tool that can plot Circos figures, as well as being able to dynamically add data to the figure, and providing information for specific data points using mouse hover display and zoom in/out functions. J-Circos uses the Java computer language to enable it to be used on most operating systems (Windows, MacOS, Linux). Users can input data into JCircos using flat data formats, as well as from the GUI. J-Circos will enable biologists to better study more complex chromosomal interactions and fusion transcripts that are otherwise difficult to visualize from next-generation sequencing data.

Inferring phylogenies of evolving sequences without multiple sequence alignment

Alignment-free methods, in which shared properties of sub-sequences (e.g. identity or match length) are extracted and used to compute a distance matrix, have recently been explored for phylogenetic inference. However, the scalability and robustness of these methods to key evolutionary processes remain to be investigated. Here, using simulated sequence sets of various sizes in both nucleotides and amino acids, we systematically assess the accuracy of phylogenetic inference using an alignment-free approach, based on D2 statistics, under different evolutionary scenarios. We find that compared to a multiple sequence alignment approach, D2 methods are more robust against among-site rate heterogeneity, compositional biases, genetic rearrangements and insertions/deletions, but are more sensitive to recent sequence divergence and sequence truncation. Across diverse empirical datasets, the alignment-free methods perform well for sequences sharing low divergence, at greater computation speed. Our findings provide strong evidence for the scalability and the potential use of alignment-free methods in large-scale phylogenomics.

Epigenetic modelling

The field of epigenetics looks at changes in the chromosomal structure that affect gene expression without altering DNA sequence. A large-scale modelling project to better understand these mechanisms is gaining momentum. Early advances in genetics led to the all-genetic paradigm: phenotype (an organism's characteristics/behaviour) is determined by genotype (its genetic make-up). This was later amended and expressed by the well-known formula P = G + E, encompassing the notion that the visible characteristics of a living organism (the phenotype, P) is a combination of hereditary genetic factors (the genotype, G) and environmental factors (E). However, this method fails to explain why in diseases such as schizophrenia we still observe differences between identical twins. Furthermore, the identification of environmental factors (such as smoking and air quality for lung cancer) is relatively rare. The formula also fails to explain cell differentiation from a single fertilized cell. In the wake of early work by Waddington, more recent results have emphasized that the expression of the genotype can be altered without any change in the DNA sequence. This phenomenon has been tagged as epigenetics. To form the chromosome, DNA strands roll over nucleosomes, which are a cluster of nine proteins (histones), as detailed in Figure 1. Epigenetic mechanisms involve inherited alterations in these two structures, eg through attachment of a functional group to the amino acids (methyl, acetyl and phosphate). These 'stable alterations' arise during development and cell proliferation and persist through cell division. While information within the genetic material is not changed, instructions for its assembly and interpretation may be. Modelling this new paradigm, P = G + E + EpiG, is the object of our study.

WT1 interacts with MAD2 and regulates mitotic checkpoint function

Tumour suppressors safeguard the fidelity of the mitotic checkpoint by transcriptional regulation of genes that encode components of the mitotic checkpoint complex (MCC). Here we report a new role for the tumour suppressor and transcription factor, WT1, in the mitotic checkpoint. We show that WT1 regulates the MCC by directly interacting with the spindle assembly checkpoint protein, MAD2. WT1 colocalizes with MAD2 during mitosis and preferentially binds to the functionally active, closed-conformer, C-MAD2. Furthermore, WT1 associates with the MCC containing MAD2, BUBR1 and CDC20, resulting in prolonged inhibition of the anaphase-promoting complex/cyclosome (APC/C) and delayed degradation of its substrates SECURIN and CYCLIN B1. Strikingly, RNA interference-mediated depletion of WT1 leads to enhanced turnover of SECURIN, decreased lag time to anaphase and defects in chromosome segregation. Our findings identify WT1 as a regulator of the mitotic checkpoint and chromosomal stability.

Exercise-induced DNA damage: Is there a relationship with inflammatory responses?

Both a systemic inflammatory response as well as DNA damage has been observed following exhaustive endurance exercise. Hypothetically, exercise-induced DNA damage might either be a consequence of inflammatory processes or causally involved in inflammation and immunological alterations after strenuous prolonged exercise (e.g. by inducing lymphocyte apoptosis and lymphocytopenia). Nevertheless, up to now only few studies have addressed this issue and there is hardly any evidence regarding a direct relationship between DNA or chromosomal damage and inflammatory responses in the context of exercise. The most conclusive picture that emerges from available data is that reactive oxygen and nitrogen species (RONS) appear to be the key effectors which link inflammation with DNA damage. Considering the time-courses of inflammatory and oxidative stress responses on the one hand and DNA effects on the other the lack of correlations between these responses might also be explained by too short observation periods. This review summarizes and discusses the recent findings on this topic. Furthermore, data from our own study are presented that aimed to verify potential associations between several endpoints of genome stability and inflammatory, immune-endocrine and muscle damage parameters in competitors of an Ironman triathlon until 19 days into recovery. The current results indicate that DNA effects in lymphocytes are not responsible for exercise-induced inflammatory responses. Furthermore, this investigation shows that inflammatory processes, vice versa, do not promote DNA damage, neither directly nor via an increased formation of RONS derived from inflammatory cells. Oxidative DNA damage might have been counteracted by training- and exercise-induced antioxidant responses. However, further studies are needed that combine advanced -omics based techniques (transcriptomics, proteomics) with state-of-the-art biochemical biomarkers to gain more insights into the underlying mechanisms.

Phylogeography of two freshwater prawn species from far-northern Queensland

The distribution of a freshwater species is often dependent on its ability to disperse within the riverine system. Species with high dispersal abilities tend to be widespread, whereas those with restricted dispersal tend to be geographically restricted and are usually given higher conservation priority. Population structure was compared between a widespread freshwater prawn species, Macrobrachium australiense, and a narrow-range endemic freshwater prawn, Macrobrachium koombooloomba. The distribution of M. australiense and M. koombooloomba did not overlap, although suggested historical river-boundary rearrangements indicate that there has been the potential for dispersal into neighbouring catchments. A fragment of the mtDNA CO1 gene was analysed and a Mantel test revealed a significant isolation by distance effect for both species. Significant overall FST values confirmed that both species exhibited low levels of dispersal, a prediction for populations inhabiting a fragmented upland environment. The level of structure in M. australiense is surprising for a widely distributed species. Not all M. australiense populations conformed to the stream-hierarchy model, with results being best explained by historical river realignment or cross-catchment dispersal. The fact that both species show limited dispersal highlights the importance of conservation in highland areas for both endemic and widely spread species.

Super DNAging - New insights into DNA integrity, genome stability and telomeres in the oldest old

Reductions in DNA integrity, genome stability, and telomere length are strongly associated with the aging process, age-related diseases as well as the age-related loss of muscle mass. However, in people reaching an age far beyond their statistical life expectancy the prevalence of diseases, such as cancer, cardiovascular disease, diabetes or dementia, is much lower compared to “averagely” aged humans. These inverse observations in nonagenarians (90–99 years), centenarians (100–109 years) and super-centenarians (110 years and older) require a closer look into dynamics underlying DNA damage within the oldest old of our society. Available data indicate improved DNA repair and antioxidant defense mechanisms in “super old” humans, which are comparable with much younger cohorts. Partly as a result of these enhanced endogenous repair and protective mechanisms, the oldest old humans appear to cope better with risk factors for DNA damage over their lifetime compared to subjects whose lifespan coincides with the statistical life expectancy. This model is supported by study results demonstrating superior chromosomal stability, telomere dynamics and DNA integrity in “successful agers”. There is also compelling evidence suggesting that life-style related factors including regular physical activity, a well-balanced diet and minimized psycho-social stress can reduce DNA damage and improve chromosomal stability. The most conclusive picture that emerges from reviewing the literature is that reaching “super old” age appears to be primarily determined by hereditary/genetic factors, while a healthy lifestyle additionally contributes to achieving the individual maximum lifespan in humans. More research is required in this rapidly growing population of super old people. In particular, there is need for more comprehensive investigations including short- and long-term lifestyle interventions as well as investigations focusing on the mechanisms causing DNA damage, mutations, and telomere shortening.

Suggestive linkage of 2p22-25 and 11q12-13 with low bone mineral density at the lumbar spine in the Irish population

Osteoporosis is a disease characterized by low bone mineral density (BMD) and poor bone quality. Peak bone density is achieved by the third decade of life, after which bone is maintained by a balanced cycle of bone resorption and synthesis. Age-related bone loss occurs as the bone resorption phase outweighs the bone synthesis phase of bone metabolism. Heritability accounts for up to 90% of the variability in BMD. Chromosomal loci including 1p36, 2p22-25, 11q12-13, parathyroid hormone receptor type 1 (PTHR1), interleukin-6 (IL-6), interleukin 1 alpha (IL-1α) and type II collagen A1/vitamin D receptor (COL11A1/VDR) have been linked or shown suggestive linkage with BMD in other populations. To determine whether these loci predispose to low BMD in the Irish population, we investigated 24 microsatellite markers at 7 chromosomal loci by linkage studies in 175 Irish families of probands with primary low BMD (T-score ≤ -1.5). Nonparametric analysis was performed using the maximum likelihood variance estimation and traditional Haseman-Elston tests on the Mapmaker/Sibs program. Suggestive evidence of linkage was observed with lumbar spine BMD at 2p22-25 (maximum LOD score 2.76) and 11q12-13 (MLS 2.55). One region, 1p36, approached suggestive linkage with femoral neck BMD (MLS 2.17). In addition, seven markers achieved LOD scores > 1.0, D2S149, D11S1313, D11S987, D11S1314 including those encompassing the PTHR1 (D3S3559, D3S1289) for lumbar spine BMD and D2S149 for femoral neck BMD. Our data suggest that genes within a these chromosomal regions are contributing to a predisposition to low BMD in the Irish population.

Combined analysis of three whole genome linkage scans for Ankylosing Spondylitis

Objective. Ankylosing spondylitis (AS) is a debilitating chronic inflammatory condition with a high degree of familiality (λs=82) and heritability (>90%) that primarily affects spinal and sacroiliac joints. Whole genome scans for linkage to AS phenotypes have been conducted, although results have been inconsistent between studies and all have had modest sample sizes. One potential solution to these issues is to combine data from multiple studies in a retrospective meta-analysis. Methods: The International Genetics of Ankylosing Spondylitis Consortium combined data from three whole genome linkage scans for AS (n=3744 subjects) to determine chromosomal markers that show evidence of linkage with disease. Linkage markers typed in different centres were integrated into a consensus map to facilitate effective data pooling. We performed a weighted meta-analysis to combine the linkage results, and compared them with the three individual scans and a combined pooled scan. Results: In addition to the expected region surrounding the HLA-B27 gene on chromosome 6, we determined that several marker regions showed significant evidence of linkage with disease status. Regions on chromosome 10q and 16q achieved 'suggestive' evidence of linkage, and regions on chromosomes 1q, 3q, 5q, 6q, 9q, 17q and 19q showed at least nominal linkage in two or more scans and in the weighted meta-analysis. Regions previously associated with AS on chromosome 2q (the IL-1 gene cluster) and 22q (CYP2D6) exhibited nominal linkage in the meta-analysis, providing further statistical support for their involvement in susceptibility to AS. Conclusion: These findings provide a useful guide for future studies aiming to identify the genes involved in this highly heritable condition. . Published by on behalf of the British Society for Rheumatology.

Target acquired: Progress and promise of targeted therapeutics in the treatment of prostate cancer

Cancer is fundamentally a genomic disease caused by mutations or rearrangements in the DNA or epigenetic machinery of a patient. An emerging field in cancer treatment targets key aberrations arising from the mutational landscape of an individual patient’s disease rather than employing a cancer-wide cytotoxic therapy approach. In prostate cancer in particular, where there is an observed variation in response to standard treatments between patients with disease of a similar pathological stage and grade, mutationdirected treatment may grow to be a viable tool for clinicians to tailor more effective treatments. This review will describe a number of mutations across multiple forms of cancer that have been successfully antagonised by targeted therapeutics including their identification, the development of targeted compounds to combat them and the development of resistance to these therapies. This review will continue to examine these same mutations in the treatment and management of prostate cancer; the prevalence of targetable mutations in prostate cancer, recent clinical trials of targeted-agents and the potential or limitations for their use.

Clinical Escherichia coli isolates utilize alpha-hemolysin to inhibit in vitro epithelial cytokine production

Uropathogenic Escherichia coli is the primary cause of urinary tract infections, which affects over 60% of women during their lifetime. UPEC exhibits a number of virulence traits that facilitate colonization of the bladder, including inhibition of cytokine production by bladder epithelial cells. The goal of this study was to identify the mechanism of this inhibition. We observed that cytokine suppression was associated with rapid cytotoxicity toward epithelial cells. We found that cytotoxicity, cytokine suppression and alpha-hemolysin production were all tightly linked in clinical isolates. We screened a UPEC fosmid library and identified clones that gained the cytotoxicity and cytokine-suppression phenotypes. Both clones contained fosmids encoding a PAI II(J96)-like domain and expressed the alpha-hemolysin (hlyA) encoded therein. Mutation of the fosmid-encoded hly operon abolished cytotoxicity and cytokine suppression. Similarly, mutation of the chromosomal hlyCABD operon of UPEC isolate F11 also abolished these phenotypes, and they could be restored by introducing the PAI II(J96)-like domain-encoding fosmid. We also examined the role of alpha-hemolysin in cytokine production both in the murine UTI model as well as patient specimens. We conclude that E. coli utilizes alpha-hemolysin to inhibit epithelial cytokine production in vitro. Its contribution to inflammation during infection requires further study.

Genome-wide linkage analysis of multiple measures of neuroticism of 2 large cohorts from Australia and the Netherlands

CONTEXT People meeting diagnostic criteria for anxiety or depressive disorders tend to score high on the personality scale of neuroticism. Studying this personality dimension can give insights into the etiology of these important psychiatric disorders. OBJECTIVES To undertake a comprehensive genome-wide linkage study of neuroticism using large study samples that have been measured multiple times and to compare the results between countries for replication and across time within countries for consistency. DESIGN Genome-wide linkage scan. SETTING Twin individuals and their family members from Australia and the Netherlands. PARTICIPANTS Nineteen thousand six hundred thirty-five sibling pairs completed self-report questionnaires for neuroticism up to 5 times over a period of up to 22 years. Five thousand sixty-nine sibling pairs were genotyped with microsatellite markers. METHODS Nonparametric linkage analyses were conducted in MERLIN-REGRESS for the mean neuroticism scores averaged across time. Additional analyses were conducted for the time-specific measures of neuroticism from each country to investigate consistency of linkage results. RESULTS Three chromosomal regions exceeded empirically derived thresholds for suggestive linkage using mean neuroticism scores: 10p 5 Kosambi cM (cM) (Dutch study sample), 14q 103 cM (Dutch study sample), and 18q 117 cM (combined Australian and Dutch study sample), but only 14q retained significance after correction for multiple testing. These regions all showed evidence for linkage in individual time-specific measures of neuroticism and 1 (18q) showed some evidence for replication between countries. Linkage intervals for these regions all overlap with regions identified in other studies of neuroticism or related traits and/or in studies of anxiety in mice. CONCLUSIONS Our results demonstrate the value of the availability of multiple measures over time and add to the optimism reported in recent reviews for replication of linkage regions for neuroticism. These regions are likely to harbor causal variants for neuroticism and its related psychiatric disorders and can inform prioritization of results from genome-wide association studies.