Integration Analysis of Three Omics Data Using Penalized Regression Methods: An Application to Bladder Cancer

Omics data integration is becoming necessary to investigate the genomic mechanisms involved in complex diseases. During the integration process, many challenges arise such as data heterogeneity, the smaller number of individuals in comparison to the number of parameters, multicollinearity, and interpretation and validation of results due to their complexity and lack of knowledge about biological processes. To overcome some of these issues, innovative statistical approaches are being developed. In this work, we propose a permutation-based method to concomitantly assess significance and correct by multiple testing with the MaxT algorithm. This was applied with penalized regression methods (LASSO and ENET) when exploring relationships between common genetic variants, DNA methylation and gene expression measured in bladder tumor samples. The overall analysis flow consisted of three steps: (1) SNPs/CpGs were selected per each gene probe within 1Mb window upstream and downstream the gene; (2) LASSO and ENET were applied to assess the association between each expression probe and the selected SNPs/CpGs in three multivariable models (SNP, CPG, and Global models, the latter integrating SNPs and CPGs); and (3) the significance of each model was assessed using the permutation-based MaxT method. We identified 48 genes whose expression levels were significantly associated with both SNPs and CPGs. Importantly, 36 (75%) of them were replicated in an independent data set (TCGA) and the performance of the proposed method was checked with a simulation study. We further support our results with a biological interpretation based on an enrichment analysis. The approach we propose allows reducing computational time and is flexible and easy to implement when analyzing several types of omics data. Our results highlight the importance of integrating omics data by applying appropriate statistical strategies to discover new insights into the complex genetic mechanisms involved in disease conditions.     

Evolution of the prevalence and incidence of consumption of antidepressants in a Spanish region (2002-2007)

Background The treatment of depressive disorders involves the administration of drugs of proven efficacy at the correct doses and for specific periods of time, in conjunction with psychotherapeutic support.Aim To assess the evolution of the consumption of antidepressants in the Health Region of Lleida (Spain).Method A retrospective cohort study of the antidepressant medication prescribed via the Spanish National Health System in the Health Region of Lleida between 2002 and 2007. The variables recorded in the study were age, sex, number of patients in antidepressant treatment in the Health Region of Lleida, length of treatment and type of drug. The prevalence of the population of the health region who were receiving antidepressant drugs and the incidence for each particular year was calculated.Results The mean prevalence of patients in treatment with antidepressant drugs was 8.5% (5% in men and 12.1% in women). The highest prevalence was observed in the higher age groups. By therapeutic groups, selective serotonin reuptake inhibitors (SSRIs) were the most frequently prescribed, five times more than the next group, tricyclics/heterocyclics. The follow-up assessment of the medication prescribed showed that one out of every four patients did not continue treatment after the first month, and 38.4% did not continue after three months. Very few were treated for more than six months.Conclusion This study stresses the high rate of antidepressant treatment in the older women's group. One of every four treatments initiated did not last more than one month. Over the six-year period, 16 506 patients dropped out of treatment.     

Crystal structure, magnetic properties and Mössbauer studies of [Fe(qsal)2][Ni(dmit)2]

A new compound of formula [Fe(qsal)₂][Ni(dmit)₂] (1) has been synthesised, structurally and magnetically characterised (qsalH = N-(8-quinolyl)salicylaldimine, dmit²⁻ = 1,3-dithiol-2-thione-4,5-dithiolato). Its structural features and its magnetic behaviour were compared with those of [Fe(qsal)₂]-based complexes, and more particularly [Fe(qsal)₂][Ni(dmit)₂] · 2CH₃CN.     

Age-related increase in xanthine oxidase activity in human plasma and rat tissues

This study assessed the role of xanthine oxidase in vascular ageing. A positive correlation between xanthine oxidase activity and age was found in human plasma. Similar results were found in rat plasma. Xanthine oxidase expression and activity in homogenates from the aortic wall were significantly higher in samples from old rats than in their young counterparts (p < 0.01). In rat skeletal muscle homogenates both xanthine oxidase expression and activity showed a similar age-related profile. Superoxide production by xanthine oxidase in aortic rings was higher in aged rats. Uric acid, the final product of xanthine oxidase has been proposed as a risk factor for coronary heart disease and an independent marker of worse prognosis in patients with moderate-to-severe chronic heart failure. These results give a possible explanation for this correlation and underscore the role of xanthine oxidase in ageing.     

Optimisation of a quantitative polymerase chain reaction-based strategy for the detection and quantification of human herpesvirus 6 DNA in patients undergoing allogeneic haematopoietic stem cell transplantation

Human herpesvirus 6 (HHV-6) may cause severe complications after haematopoietic stem cell transplantation (HSCT). Monitoring this virus and providing precise, rapid and early diagnosis of related clinical diseases, constitute essential measures to improve outcomes. A prospective survey on the incidence and clinical features of HHV-6 infections after HSCT has not yet been conducted in Brazilian patients and the impact of this infection on HSCT outcome remains unclear. A rapid test based on real-time quantitative polymerase chain reaction (qPCR) has been optimised to screen and quantify clinical samples for HHV-6. The detection step was based on reaction with TaqMan^® hydrolysis probes. A set of previously described primers and probes have been tested to evaluate efficiency, sensitivity and reproducibility. The target efficiency range was 91.4% with linearity ranging from 10-10⁶ copies/reaction and a limit of detection of five copies/reaction or 250 copies/mL of plasma. The qPCR assay developed in the present study was simple, rapid and sensitive, allowing the detection of a wide range of HHV-6 loads. In conclusion, this test may be useful as a practical tool to help elucidate the clinical relevance of HHV-6 infection and reactivation in different scenarios and to determine the need for surveillance.     

The UBC-40 Urothelial Bladder Cancer cell line index: a genomic resource for functional studies

Background

Urothelial bladder cancer is a highly heterogeneous disease. Cancer cell lines are useful tools for its study. This is a comprehensive genomic characterization of 40 urothelial bladder carcinoma (UBC) cell lines including information on origin, mutation status of genes implicated in bladder cancer (FGFR3, PIK3CA, TP53, and RAS), copy number alterations assessed using high density SNP arrays, uniparental disomy (UPD) events, and gene expression.

Results

Based on gene mutation patterns and genomic changes we identify lines representative of the FGFR3-driven tumor pathway and of the TP53/RB tumor suppressor-driven pathway. High-density array copy number analysis identified significant focal gains (1q32, 5p13.1-12, 7q11, and 7q33) and losses (i.e. 6p22.1) in regions altered in tumors but not previously described as affected in bladder cell lines. We also identify new evidence for frequent regions of UPD, often coinciding with regions reported to be lost in tumors. Previously undescribed chromosome X losses found in UBC lines also point to potential tumor suppressor genes. Cell lines representative of the FGFR3-driven pathway showed a lower number of UPD events.

Conclusions

Overall, there is a predominance of more aggressive tumor subtypes among the cell lines. We provide a cell line classification that establishes their relatedness to the major molecularly-defined bladder tumor subtypes. The compiled information should serve as a useful reference to the bladder cancer research community and should help to select cell lines appropriate for the functional analysis of bladder cancer genes, for example those being identified through massive parallel sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1450-3) contains supplementary material, which is available to authorized users.     

Finessing atlas data for species distribution models

Aim The spatial resolution of species atlases and therefore resulting model predictions are often too coarse for local applications. Collecting distribution data at a finer resolution for large numbers of species requires a comprehensive sampling effort, making it impractical and expensive. This study outlines the incorporation of existing knowledge into a conventional approach to predict the distribution of Bonelli’s eagle (Aquila fasciata) at a resolution 100 times finer than available atlas data. Location Malaga province, Andalusia, southern Spain. Methods A Bayesian expert system was proposed to utilize the knowledge from distribution models to yield the probability of a species being recorded at a finer resolution (1 × 1 km) than the original atlas data (10 × 10 km). The recorded probability was then used as a weight vector to generate a sampling scheme from the species atlas to enhance the accuracy of the modelling procedure. The maximum entropy for species distribution modelling (MaxEnt) was used as the species distribution model. A comparison was made between the results of the MaxEnt using the enhanced and, the random sampling scheme, based on four groups of environmental variables: topographic, climatic, biological and anthropogenic. Results The models with the sampling scheme enhanced by an expert system had a higher discriminative capacity than the baseline models. The downscaled (i.e. finer scale) species distribution maps using a hybrid MaxEnt/expert system approach were more specific to the nest locations and were more contrasted than those of the baseline model. Main conclusions The proposed method is a feasible substitute for comprehensive field work. The approach developed in this study is applicable for predicting the distribution of Bonelli’s eagle at a local scale from a national‐level occurrence data set; however, the usefulness of this approach may be limited to well‐known species.