ALBERT

Description: The Cancer Genome Atlas (TCGA) research network has made public a large collection of clinical and molecular phenotypes of more than 10 000 tumor patients across 33 different tumor types. Using this cohort, TCGA has published over 20 marker papers detailing the genomic and epigenomic alterations associated with these tumor types. Although many important discoveries have been made by TCGA's research network, opportunities still exist to implement novel methods, thereby elucidating new biological pathways and diagnostic markers. However, mining the TCGA data presents several bioinformatics challenges, such as data retrieval and integration with clinical data and other molecular data types (e.g. RNA and DNA methylation). We developed an R/Bioconductor package called TCGAbiolinks to address these challenges and offer bioinformatics solutions by using a guided workflow to allow users to query, download and perform integrative analyses of TCGA data. We combined methods from computer science and statistics into the pipeline and incorporated methodologies developed in previous TCGA marker studies and in our own group. Using four different TCGA tumor types (Kidney, Brain, Breast and Colon) as examples, we provide case studies to illustrate examples of reproducibility, integrative analysis and utilization of different Bioconductor packages to advance and accelerate novel discoveries.

Description: Coenzyme Q10 (CoQ10) is essential for the energy production of the cells and as an electron transporter in the mitochondrial respiratory chain. CoQ10 links the mitochondrial fatty acid β-oxidation to the respiratory chain by accepting electrons from electron transfer flavoprotein-ubiquinone oxidoreductase (ETF-QO). Recently, it was shown that a group of patients with the riboflavin responsive form of multiple acyl-CoA dehydrogenation deficiency (RR-MADD) carrying inherited amino acid variations in ETF-QO also had secondary CoQ10 deficiency with beneficial effects of CoQ10 treatment, thus adding RR-MADD to an increasing number of diseases involving secondary CoQ10 deficiency. In this study, we show that moderately decreased CoQ10 levels in fibroblasts from six unrelated RR-MADD patients were associated with increased levels of mitochondrial reactive oxygen species (ROS). Treatment with CoQ10, but not with riboflavin, could normalize the CoQ10 level and decrease the level of ROS in the patient cells. Additionally, riboflavin-depleted control fibroblasts showed moderate CoQ10 deficiency, but not increased mitochondrial ROS, indicating that variant ETF-QO proteins and not CoQ10 deficiency are the causes of mitochondrial ROS production in the patient cells. Accordingly, the corresponding variant Rhodobacter sphaeroides ETF-QO proteins, when overexpressed in vitro , bind a CoQ10 pseudosubstrate, Q10Br, less tightly than the wild-type ETF-QO protein, suggesting that molecular oxygen can get access to the electrons in the misfolded ETF-QO protein, thereby generating superoxide and oxidative stress, which can be reversed by CoQ10 treatment.

Description: : Pharmacogenomics holds great promise for the development of biomarkers of drug response and the design of new therapeutic options, which are key challenges in precision medicine. However, such data are scattered and lack standards for efficient access and analysis, consequently preventing the realization of the full potential of pharmacogenomics. To address these issues, we implemented PharmacoGx , an easy-to-use, open source package for integrative analysis of multiple pharmacogenomic datasets. We demonstrate the utility of our package in comparing large drug sensitivity datasets, such as the Genomics of Drug Sensitivity in Cancer and the Cancer Cell Line Encyclopedia. Moreover, we show how to use our package to easily perform Connectivity Map analysis. With increasing availability of drug-related data, our package will open new avenues of research for meta-analysis of pharmacogenomic data. Availability and implementation : PharmacoGx is implemented in R and can be easily installed on any system. The package is available from CRAN and its source code is available from GitHub. Contact : bhaibeka@uhnresearch.ca or benjamin.haibe.kains@utoronto.ca Supplementary information: Supplementary data are available at Bioinformatics online.

Description: The large repABC plasmids of the order Rhizobiales with Class I quorum-regulated conjugative transfer systems often define the nature of the bacterium that harbors them. These otherwise diverse plasmids contain a core of highly conserved genes for replication and conjugation raising the question of their evolutionary relationships. In an analysis of 18 such plasmids these elements fall into two organizational classes, Group I and Group II, based on the sites at which cargo DNA is located. Cladograms constructed from proteins of the transfer and quorum-sensing components indicated that those of the Group I plasmids, while coevolving, have diverged from those coevolving proteins of the Group II plasmids. Moreover, within these groups the phylogenies of the proteins usually occupy similar, if not identical, tree topologies. Remarkably, such relationships were not seen among proteins of the replication system; although RepA and RepB coevolve, RepC does not. Nor do the replication proteins coevolve with the proteins of the transfer and quorum-sensing systems. Functional analysis was mostly consistent with phylogenies. TraR activated promoters from plasmids within its group, but not between groups and dimerized with TraR proteins from within but not between groups. However, oriT sequences, which are highly conserved, were processed by the transfer system of plasmids regardless of group. We conclude that these plasmids diverged into two classes based on the locations at which cargo DNA is inserted, that the quorum-sensing and transfer functions are coevolving within but not between the two groups, and that this divergent evolution extends to function.

Description: Monthly means of the magnetic field measurements at ground observatories are a key data source for studying temporal changes of the core magnetic field. However, when they are calculated in the usual way, contributions of external (magnetospheric and ionospheric) origin may remain, which make them less favourable for studying the field generated by dynamo action in the core. We remove external field predictions, including a new way of characterizing the magnetospheric ring current, from the data and then calculate revised monthly means using robust methods. The geomagnetic secular variation (SV) is calculated as the first annual differences of these monthly means, which also removes the static crustal field. SV time-series based on revised monthly means are much less scattered than those calculated from ordinary monthly means, and their variances and correlations between components are smaller. On the annual to decadal timescale, the SV is generated primarily by advection in the fluid outer core. We demonstrate the utility of the revised monthly means by calculating models of the core surface advective flow between 1997 and 2013 directly from the SV data. One set of models assumes flow that is constant over three months; such models exhibit large and rapid temporal variations. For models of this type, less complex flows achieve the same fit to the SV derived from revised monthly means than those from ordinary monthly means. However, those obtained from ordinary monthly means are able to follow excursions in SV that are likely to be external field contamination rather than core signals. Having established that we can find models that fit the data adequately, we then assess how much temporal variability is required. Previous studies have suggested that the flow is consistent with torsional oscillations (TO), solid body-like oscillations of fluid on concentric cylinders with axes aligned along the Earth's rotation axis. TO have been proposed to explain decadal timescale changes in the length-of-day. We invert for flow models where the only temporal changes are consistent with TO, but such models have an unacceptably large data misfit. However, if we relax the TO constraint to allow a little more temporal variability, we can fit the data as well as with flows assumed constant over three months, demonstrating that rapid SV changes can be reproduced by rather small flow changes. Although the flow itself changes slowly, its time derivative can be locally (temporally and spatially) large, in particular when and where core surface secular acceleration peaks. Spherical harmonic expansion coefficients of the flows are not well resolved, and many of them are strongly correlated. Averaging functions, a measure of our ability to determine the flow at a given location from the data distribution available, are poor approximations to the ideal, even when centred on points of the core surface below areas of high observatory density. Both resolution and averaging functions are noticeably worse for the toroidal flow component, which dominates the flow, than the poloidal flow component, except around the magnetic equator where averaging functions for both components are poor.

Description: The bacteriophage P1 Ref (recombination enhancement function) protein is a RecA-dependent, HNH endonuclease. It can be directed to create targeted double-strand breaks within a displacement loop formed by RecA. The 76 amino acid N-terminal region of Ref is positively charged (25/76 amino acid residues) and inherently unstructured in solution. Our investigation of N-terminal truncation variants shows this region is required for DNA binding, contains a Cys involved in incidental dimerization and is necessary for efficient Ref-mediated DNA cleavage. Specifically, Ref N-terminal truncation variants lacking between 21 and 47 amino acids are more effective RecA-mediated targeting nucleases. We propose a more refined set of options for the Ref-mediated cleavage mechanism, featuring the N-terminal region as an anchor for at least one of the DNA strand cleavage events.

Description: Many studies have investigated genotoxic effects of high Se diets but very few have addressed the genotoxicity of Se deprivation and its consequences in germ cells and none in somatic cells. To address these data gaps, C57BL/6 male mice were subjected to Se deprivation starting in the parental generation, i.e. before conception. Mice were given a diet of either low (0.01mg Se/kg diet) or normal (0.23mg Se/kg diet) Se content. Ogg1-deficient ( Ogg1 –/– ) mice were used as a sensitive model towards oxidative stress due to their reduced capacity to repair oxidised purines. Ogg1 –/– mice also mimic the repair characteristics of human post-meiotic male germ cells which have a reduced ability to repair such lesions. The genotoxicity of Se deficiency was addressed by measuring DNA lesions with the alkaline single cell gel electrophoresis (+ Fpg to detect oxidised DNA lesions) in somatic cells (nucleated blood cells and lung cells) and male germ cells (testicular cells). Total Se concentration in liver and GPx activity in plasma and testicular cells were measured. Gene mutation was evaluated by an erythrocyte-based Pig-a assay. We found that Se deprivation of F1 from their conception and until early adulthood led to the induction of DNA lesions in testicular and lung cells expressed as significantly increased levels of DNA lesions, irrespective of the mouse genotype. In blood cells, Se levels did not appear to affect DNA lesions or mutant cell frequencies. The results suggest that the testis was the most sensitive tissue. Thus, genotoxicity induced by the low Se diet in the spermatozoal genome has potential implications for the offspring.