cudaMap: a GPU accelerated program for gene expression connectivity mapping

Background: Modern cancer research often involves large datasets and the use of sophisticated statistical techniques. Together these add a heavy computational load to the analysis, which is often coupled with issues surrounding data accessibility. Connectivity mapping is an advanced bioinformatic and computational technique dedicated to therapeutics discovery and drug re-purposing around differential gene expression analysis. On a normal desktop PC, it is common for the connectivity mapping task with a single gene signature to take >2h to complete using sscMap, a popular Java application that runs on standard CPUs (Central Processing Units). Here, we describe new software, cudaMap, which has been implemented using CUDA C/C++ to harness the computational power of NVIDIA GPUs (Graphics Processing Units) to greatly reduce processing times for connectivity mapping.

Results: cudaMap can identify candidate therapeutics from the same signature in just over thirty seconds when using an NVIDIA Tesla C2050 GPU. Results from the analysis of multiple gene signatures, which would previously have taken several days, can now be obtained in as little as 10 minutes, greatly facilitating candidate therapeutics discovery with high throughput. We are able to demonstrate dramatic speed differentials between GPU assisted performance and CPU executions as the computational load increases for high accuracy evaluation of statistical significance.

Conclusion: Emerging 'omics' technologies are constantly increasing the volume of data and information to be processed in all areas of biomedical research. Embracing the multicore functionality of GPUs represents a major avenue of local accelerated computing. cudaMap will make a strong contribution in the discovery of candidate therapeutics by enabling speedy execution of heavy duty connectivity mapping tasks, which are increasingly required in modern cancer research. cudaMap is open source and can be freely downloaded from http://purl.oclc.org/NET/cudaMap.

Retroviral-mediated gene transfer of the human erythropoietin gene into a factor-dependent cell line, TF-1

The factor-dependent cell line, TF-1, established from a patient with erythroleukaemia, shows characteristics of immature erythroblasts. Addition of granulocyte-macrophage colony stimulating factor (GM-CSF) to the culture medium is required for long-term growth of the cells. Erythropoietin (Epo) can also be used to sustain TF-1 cells but for only limited periods (approximately a week). Low levels of both growth factors can act synergistically to maintain proliferation for a longer period of time than Epo alone. To eliminate the requirement of exogenous Epo for growth, TF-1 cells were co-cultured with a retroviral secreting cell line containing the human erythropoietin (hEpo) gene and a neomycin (neo) selectable marker. TF-1 cells which exhibited neo resistance (indicating infection by the retrovirus) were then grown in low concentrations of GM-CSF without the addition of Epo. Under these conditions growth of normal TF-1 cells was not sustained. The neo-resistant cells survived for more than 14 days indicating synergy between GM-CSF and the Epo synthesised by the co-cultured TF-1 cells. Radioimmunoassays performed on growth media detected concentrations up to 1 mU/ml of Epo, implying that stable integration of the retroviral vector and expression of the hEpo gene have been achieved.

SORBS1 gene, a new candidate for diabetic nephropathy: results from a multi-stage genome-wide association study in patients with type 1 diabetes

Aims/hypothesis

An arsenate tolerance gene on chromosome 6 of rice

The genetics of arsenic tolerance in plants has not been extensively studied and no arsenic tolerance gene has been genetically mapped. Screening 20 diverse genotypes of rice for reduced root growth in 13.3 μM arsenate identified marked differences in tolerance. The most sensitive variety, Dawn, is known to be highly susceptible to straighthead, a condition linked to arsenic contamination of soil. Screening 108 recombinant inbred lines of the Bala x Azucena mapping population revealed the presence of a major gene, AsTol, which mapped between markers RZ516 and RG213 on chromosome 6. This gene is a good target for further characterisation. It should prove valuable for investigations into the physiological and molecular mechanism behind arsenic tolerance in plants and may lead to strategies aimed at breeding for arsenic contaminated regions. © New Phytologist (2004).

Wilms' tumor gene 1 (WT1) expression in childhood acute lymphoblastic leukemia:a wide range of WT1 expression levels, its impact on prognosis and minimal residual disease monitoring

Wilms' tumor gene 1 (WT1) is overexpressed in the majority (70-90%) of acute leukemias and has been identified as an independent adverse prognostic factor, a convenient minimal residual disease (MRD) marker and potential therapeutic target in acute leukemia. We examined WT1 expression patterns in childhood acute lymphoblastic leukemia (ALL), where its clinical implication remains unclear. Using a real-time quantitative PCR designed according to Europe Against Cancer Program recommendations, we evaluated WT1 expression in 125 consecutively enrolled patients with childhood ALL (106 BCP-ALL, 19 T-ALL) and compared it with physiologic WT1 expression in normal and regenerating bone marrow (BM). In childhood B-cell precursor (BCP)-ALL, we detected a wide range of WT1 levels (5 logs) with a median WT1 expression close to that of normal BM. WT1 expression in childhood T-ALL was significantly higher than in BCP-ALL (P<0.001). Patients with MLL-AF4 translocation showed high WT1 overexpression (P<0.01) compared to patients with other or no chromosomal aberrations. Older children (> or =10 years) expressed higher WT1 levels than children under 10 years of age (P<0.001), while there was no difference in WT1 expression in patients with peripheral blood leukocyte count (WBC) > or =50 x 10(9)/l and lower. Analysis of relapsed cases (14/125) indicated that an abnormal increase or decrease in WT1 expression was associated with a significantly increased risk of relapse (P=0.0006), and this prognostic impact of WT1 was independent of other main risk factors (P=0.0012). In summary, our study suggests that WT1 expression in childhood ALL is very variable and much lower than in AML or adult ALL. WT1, thus, will not be a useful marker for MRD detection in childhood ALL, however, it does represent a potential independent risk factor in childhood ALL. Interestingly, a proportion of childhood ALL patients express WT1 at levels below the normal physiological BM WT1 expression, and this reduced WT1 expression appears to be associated with a higher risk of relapse.

Connectivity mapping using a combined gene signature from multiple colorectal cancer datasets identified candidate drugs including existing chemotherapies

BACKGROUND: While the discovery of new drugs is a complex, lengthy and costly process, identifying new uses for existing drugs is a cost-effective approach to therapeutic discovery. Connectivity mapping integrates gene expression profiling with advanced algorithms to connect genes, diseases and small molecule compounds and has been applied in a large number of studies to identify potential drugs, particularly to facilitate drug repurposing. Colorectal cancer (CRC) is a commonly diagnosed cancer with high mortality rates, presenting a worldwide health problem. With the advancement of high throughput omics technologies, a number of large scale gene expression profiling studies have been conducted on CRCs, providing multiple datasets in gene expression data repositories. In this work, we systematically apply gene expression connectivity mapping to multiple CRC datasets to identify candidate therapeutics to this disease.

RESULTS: We developed a robust method to compile a combined gene signature for colorectal cancer across multiple datasets. Connectivity mapping analysis with this signature of 148 genes identified 10 candidate compounds, including irinotecan and etoposide, which are chemotherapy drugs currently used to treat CRCs. These results indicate that we have discovered high quality connections between the CRC disease state and the candidate compounds, and that the gene signature we created may be used as a potential therapeutic target in treating the disease. The method we proposed is highly effective in generating quality gene signature through multiple datasets; the publication of the combined CRC gene signature and the list of candidate compounds from this work will benefit both cancer and systems biology research communities for further development and investigations.

Development of TMTP-1 targeted designer biopolymers for gene delivery to prostate cancer

Designer biopolymers (DBPs) represent state of the art genetically engineered biomacromolecules designed to condense plasmid DNA, and overcome intra- and extra- cellular barriers to gene delivery. Three DBPs were synthesized, each with the tumor molecular targeting peptide-1 (TMTP-1) motif to specifically target metastases. Each DBP was complexed with a pEGFP-N1 reporter plasmid to permit physiochemical and biological assay analysis. Results indicated that two of the biopolymers (RMHT and RM3GT) effectively condensed pEGFP-N1 into cationic nanoparticles< 100nm and were capable of transfecting PC-3 metastatic prostate cancer cells. Conversely the anionic RMGT DBP nanoparticles could not transfect PC-3 cells. RMHT and RM3GT nanoparticles were stable in the presence of serum and protected the cargo from degradation. Additionally it was concluded that cell viability could recover post-transfection with these DBPs, which were less toxic than the commercially available transfection reagent Lipofectamine® 2000. With both DBPs, a higher transfection efficacy was observed in PC-3 cells than in the moderately metastatic, DU145, and normal, PNT2-C2, cell lines. Blocking of the TMTP-1 receptors inhibited gene transfer indicating internalization via this receptor. In conclusion RMHT and RM3GT are fully functional DBPs that address major obstacles to gene delivery and target metastatic cells expressing the TMTP-1 receptor.

A gene-signature progression approach to identifying candidate small-molecule cancer therapeutics with connectivity mapping

Gene expression connectivity mapping has gained much popularity recently with a number of successful applications in biomedical research testifying its utility and promise. Previously methodological research in connectivity mapping mainly focused on two of the key components in the framework, namely, the reference gene expression profiles and the connectivity mapping algorithms. The other key component in this framework, the query gene signature, has been left to users to construct without much consensus on how this should be done, albeit it has been an issue most relevant to end users. As a key input to the connectivity mapping process, gene signature is crucially important in returning biologically meaningful and relevant results. This paper intends to formulate a standardized procedure for constructing high quality gene signatures from a user’s perspective.

Genomic structure and expression of parathyroid hormone-related protein gene (PTHrP) in a teleost, Fugu rubripes

In this study we describe the isolation and characterisation of the parathyroid hormone-related protein (PTHrP) gene from the teleost Fugu rubripes. The gene has a relatively simple structure, compared with tetrapod PTHrP genes, composed of three exons and two introns, encompassing 2.25 kb of genomic DNA. The gene encodes a protein of 163 amino acids, with a putative signal peptide of 37 amino acids and a mature peptide of 126 amino acids. The overall homology with known tetrapod PTHrP proteins is low (36%), with a novel sequence inserted between positions 38 and 65, the absence of the conserved pentapeptide (TRSAW) and shortened C-terminal domain. The N-terminus shows greater conservation (62%), suggesting that it may have a hypercalcaemic function similar to that of tetrapod PTHrP. In situ localisation and RT–PCR have demonstrated the presence of PTHrP in a wide range of tissues with varying levels of expression. Sequence scanning of overlapping cosmids has identified three additional genes, TMPO, LDHB and KCNA1, which map to human chromosome 12, with the latter two mapping to 12p12-11.2. PTHrP in human also maps to this chromosome 12 sub-region, thus demonstrating conservation of synteny between human and Fugu.

Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links.

Metabolic traits are molecular phenotypes that can drive clinical phenotypes and may predict disease progression. Here, we report results from a metabolome- and genome-wide association study on (1)H-NMR urine metabolic profiles. The study was conducted within an untargeted approach, employing a novel method for compound identification. From our discovery cohort of 835 Caucasian individuals who participated in the CoLaus study, we identified 139 suggestively significant (P<5×10(-8)) and independent associations between single nucleotide polymorphisms (SNP) and metabolome features. Fifty-six of these associations replicated in the TasteSensomics cohort, comprising 601 individuals from São Paulo of vastly diverse ethnic background. They correspond to eleven gene-metabolite associations, six of which had been previously identified in the urine metabolome and three in the serum metabolome. Our key novel findings are the associations of two SNPs with NMR spectral signatures pointing to fucose (rs492602, P = 6.9×10(-44)) and lysine (rs8101881, P = 1.2×10(-33)), respectively. Fine-mapping of the first locus pinpointed the FUT2 gene, which encodes a fucosyltransferase enzyme and has previously been associated with Crohn's disease. This implicates fucose as a potential prognostic disease marker, for which there is already published evidence from a mouse model. The second SNP lies within the SLC7A9 gene, rare mutations of which have been linked to severe kidney damage. The replication of previous associations and our new discoveries demonstrate the potential of untargeted metabolomics GWAS to robustly identify molecular disease markers.

HIV-1 drug resistance transmission networks in southwest Switzerland.

To determine viral subtypes and resistance mutations to antiretroviral treatment (ART) in untreated HIV-1 acutely infected subjects from Southwest Switzerland. Clinical samples were obtained from the HIV primary infection cohort from Lausanne. Briefly, pol gene was amplified by nested PCR and sequenced to generate a 1?kb sequence spanning protease and reverse transcriptase key protein regions. Nucleotide sequences were used to assess viral genotype and ART resistance mutations. Blood specimens and medical information were obtained from 30 patients. Main viral subtypes corresponded to clade B, CRF02_AG, and F1. Resistant mutations to PIs consisted of L10V and accessory mutations 16E and 60E present in all F1 clades. The NNRTI major resistant mutation 103N was detected in all F1 viruses and in other 2 clades. Additionally, we identified F1 sequences from other 6 HIV infected and untreated individuals from Southwest Switzerland, harboring nucleotide motifs and resistance mutations to ART as observed in the F1 strains from the cohort. These data reveal a high transmission rate (16.6%) for NNRTI resistant mutation 103N in a cohort of HIV acute infection. Three of the 5 resistant strains were F1 clades closely related to other F1 isolates from HIV-1 infection untreated patients also coming from Southwest Switzerland. Overall, we provide strong evidence towards an HIV-1 resistant transmission network in Southwest Switzerland. These findings have relevant implications for the local molecular mapping of HIV-1 and future ART surveillance studies in the region.