Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping. © 2013 McArt et al.

Metamaterials Nanosensor for Label-free Analysis of Conformation and Molecular Interaction of Biomolecules

Analysis of molecular interaction and conformational dynamics of biomolecules is of paramount importance in understanding of their vital functions in complex biological systems, disease detection, and new drug development. Plasmonic biosensors based upon surface plasmon resonance and localized surface plasmon resonance have become the predominant workhorse for detecting accumulated biomass caused by molecular binding events. However, unlike surface-enhanced Raman spectroscopy (SERS), the plasmonic biosensors indeed are not suitable tools to interrogate vibrational signatures of conformational transitions required for biomolecules to interact. Here, we show that plasmonic metamaterials can offer two transducing channels for parallel acquisition of optical transmission and sensitive SERS spectra at the biointerface, simultaneously probing the conformational states and binding affinity of biomolecules, e.g. G-quadruplexes, in different environments (Fig. 1). We further demonstrate the use of the metamaterials for fingerprinting and detection of arginine-glycine-glycine domain of nucleolin, a cancer biomarker which specifically binds to a G-quadruplex, with the picomolar sensitivity. The dual-mode nanosensor will significantly contribute to unraveling the complexes of the conformational dynamics of biomolecules as well as to improving specificity of biodetection assays.

Finding exoplanets orbiting young active stars - I. Technique

Stellar activity, such as starspots, can induce radial velocity (RV) variations that can mask or even mimic the RV signature of orbiting exoplanets. For this reason RV exoplanet surveys have been unsuccessful when searching for planets around young, active stars and are therefore failing to explore an important regime which can help to reveal how planets form and migrate. This paper describes a new technique to remove spot signatures from the stellar line-profiles of moderately rotating, active stars (v sin i ranging from 10 to 50 km s(-1)). By doing so it allows planetary RV signals to be uncovered. We used simulated models of a G5V type star with differing dark spots on its surface along with archive data of the known active star HD 49933 to validate our method. The results showed that starspots could be effectively cleaned from the line-profiles so that the stellar RV jitter was reduced by more than 80 per cent. Applying this procedure to the same models and HD 49933 data, but with fake planets injected, enabled the effective removal of starspots so that Jupiter mass planets on short orbital periods were successfully recovered. These results show that this approach can be useful in the search for hot-Jupiter planets that orbit around young, active stars with a v sin i of similar to 10-50 km/s.

Metamaterials-Based Label-Free Nanosensor for Conformation and Affinity Biosensing

Analysis of molecular interaction and conformational dynamics of biomolecules is of paramount importance in understanding of their vital functions in complex biological systems, disease detection, and new drug development. Plasmonic biosensors based upon surface plasmon resonance and localized surface plasmon resonance have become the predominant workhorse for detecting accumulated biomass caused by molecular binding events. However, unlike surface-enhanced Raman spectroscopy (SERS), the plasmonic biosensors indeed are not suitable tools to interrogate vibrational signatures of conformational transitions required for biomolecules to interact. Here, we show that highly tunable plasmonic metamaterials can offer two transducing channels for parallel acquisition of optical transmission and sensitive SERS spectra at the biointerface, simultaneously probing the conformational states and binding affinity of biomolecules, e.g. G-quadruplexes, in different environments. We further demonstrate the use of the metamaterials for fingerprinting and detection of arginine-glycine-glycine domain of nucleolin, a cancer biomarker which specifically binds to a G-quadruplex, with the picomolar sensitivity.

Killer immunoglobulin-like Receptors (KIR) haplogroups A and B track with Natural Killer Cells and Cytokine Profile in Aged Subjects: Observations from Octo/Nonagenarians in the Belfast Elderly Longitudinal Free-living Aging STudy (BELFAST)

Background: Natural Killer Cells (NK) play an important role in detection and elimination of virus-infected, damaged or cancer cells. NK cell function is guided by expression of Killer Immunoglobulin-like Receptors (KIRs) and contributed to by the cytokine milieu. KIR molecules are grouped on NK cells into stimulatory and inhibitory KIR haplotypes A and B, through which NKs sense and tolerate HLA self-antigens or up-regulate the NK-cytotoxic response to cells with altered HLA self-antigens, damaged by viruses or tumours. We have previously described increased numbers of NK and NK-related subsets in association with sIL-2R cytokine serum levels in BELFAST octo/nonagenarians. We hypothesised that changes in KIR A and B haplotype gene frequencies could explain the increased cytokine profiles and NK compartments previously described in Belfast Elderly Longitudinal Free-living Aging STudy (BELFAST) octo/nonagenarians, who show evidence of ageing well.

Results: In the BELFAST study, 24% of octo/nonagenarians carried the KIR A haplotype and 76% KIR B haplotype with no differences for KIR A haplogroup frequency between male or female subjects (23% v 24%; p=0.88) or for KIR B haplogroup (77% v 76%; p=0.99). Octo/nonagenarian KIR A haplotype carriers showed increased NK numbers and percentage compared to Group B KIR subjects (p=0.003; p=0.016 respectively). There were no KIR A/ B haplogroup-associated changes for related CD57+CD8 (high or low) subsets. Using logistic regression, KIR B carriers were predicted to have higher IL-12 cytokine levels compared to KIR A carriers by about 3% (OR 1.03, confidence limits CI 0.99–1.09; p=0.027) and 14% higher levels for TGF-ß (active), a cytokine with an anti-inflammatory role, (OR 1.14, confidence limits CI 0.99–1.09; p=0.002).

Conclusion: In this observational study, BELFAST octo/nonagenarians carrying KIR A haplotype showed higher NK cell numbers and percentage compared to KIR B carriers. Conversely, KIR B haplotype carriers, with genes encoding for activating KIRs, showed a tendency for higher serum pro-inflammatory cytokines compared to KIR A carriers. While the findings in this study should be considered exploratory they may serve to stimulate debate about the immune signatures of those who appear to age slowly and who represent a model for good quality survivor-hood.© 2013 Rea et al.; licensee BioMed Central Ltd.

Τesting models for the beginnings of the Aurignacian and the advent of figurative art and music: The radiocarbon chronology of Geißenklösterle

The German site of Geißenklösterle is crucial to debates concerning the European Middle to Upper Palaeolithic transition and the origins of the Aurignacian in Europe. Previous dates from the site are
central to an important hypothesis, the Kulturpumpe model, which posits that the Swabian Jura was an area where crucial behavioural developments took place and then spread to other parts of Europe. The previous chronology (critical to the model), is based mainly on radiocarbon dating, but remains poorly constrained due to the dating resolution and the variability of dates. The cause of these problems is disputed, but two principal explanations have been proposed: a) larger than expected variations in the production of atmospheric radiocarbon, and b) taphonomic in?uences in the site mixing the bones that were dated into different parts of the site. We reinvestigate the chronology using a new series of radiocarbon determinations obtained from the Mousterian, Aurignacian and Gravettian levels. The results strongly imply that the previous dates were affected by insuf?cient decontamination of the bone collagen prior to dating. Using an ultra?ltration protocol the chronometric picture becomes much clearer. Comparison of the results against other recently dated sites in other parts of Europe suggests the Early Aurignacian levels are earlier than other sites in the south of France and Italy, but not as early as recently dated sites which suggest a pre-Aurignacian dispersal of modern humans to Italy byw45000 cal BP. They are consistent with the importance of the Danube Corridor as a key route for the movement of people and ideas. The new dates fail to refute the Kulturpumpe model and suggest that Swabian Jura is a region that contributed signi?cantly to the evolution of symbolic behaviour as indicated by early evidence for ?gurative art, music and mythical imagery. © 2012 Elsevier Ltd. All rights reserved.

The impact of personalised therapies on respiratory medicine

Stratified approaches to treating disease are very attractive, as efficacy is maximised by identifying responders using a companion diagnostic or by careful phenotyping. This approach will spare non-responders form potential side-effects. This has been pioneered in oncology where single genes or gene signatures indicate tumours that will respond to specific chemotherapies. Stratified approaches to the treatment of asthma with biological therapies are currently being extensively studied. In cystic fibrosis (CF), therapies have been developed that are targeted at specific functional classes of mutations. Ivacaftor, the first of such therapies, potentiates dysfunctional cystic fibrosis transmembrane conductance regulator (CFTR) protein Class III mutations and is now available in the USA and some European countries. Pivotal studies in patients with a G551D mutation, the most common Class III mutation, have demonstrated significant improvements in clinically important outcomes such as spirometry and exacerbations. Sweat chloride was significantly reduced demonstrating a functional effect on the dysfunctional CFTR protein produced by the G551D mutation. Symptom scores are also greatly improved to a level that indicates that this is a transformational treatment for many patients. This stratified approach to the development of therapies based on the functional class of the mutations in CF is likely to lead to new drugs or combinations that will correct the basic defect in many patients with CF. © ERS 2013.

Fast recognition of musical sounds based on timbre

Human listeners seem to have an impressive ability to recognize a wide variety of natural sounds. However, there is surprisingly little quantitative evidence to characterize this fundamental ability. Here the speed and accuracy of musical-sound recognition were measured psychophysically with a rich but acoustically balanced stimulus set. The set comprised recordings of notes from musical instruments and sung vowels. In a first experiment, reaction times were collected for three target categories: voice, percussion, and strings. In a go/no-go task, listeners reacted as quickly as possible to members of a target category while withholding responses to distractors (a diverse set of musical instruments). Results showed near-perfect accuracy and fast reaction times, particularly for voices. In a second experiment, voices were recognized among strings and vice-versa. Again, reaction times to voices were faster. In a third experiment, auditory chimeras were created to retain only spectral or temporal features of the voice. Chimeras were recognized accurately, but not as quickly as natural voices. Altogether, the data suggest rapid and accurate neural mechanisms for musical-sound recognition based on selectivity to complex spectro-temporal signatures of sound sources.

Infrared non-detection of fomalhaut b: Implications for the planet interpretation

The nearby A4-type star Fomalhaut hosts a debris belt in the form of an eccentric ring, which is thought to be caused by dynamical influence from a giant planet companion. In 2008, a detection of a point source inside the inner edge of the ring was reported and was interpreted as a direct image of the planet, named Fomalhaut b. The detection was made at 600-800nm, but no corresponding signatures were found in the near-infrared range, where the bulk emission of such a planet should be expected. Here, we present deep observations of Fomalhaut with Spitzer/IRAC at 4.5 µm, using a novel point-spread function subtraction technique based on angular differential imaging and Locally Optimized Combination of Images, in order to substantially improve the Spitzer contrast at small separations. The results provide more than an order ofmagnitude improvement in the upper flux limit of Fomalhaut b and exclude the possibility that any flux from a giant planet surface contributes to the observed flux at visible wavelengths. This renders any direct connection between the observed light source and the dynamically inferred giant planet highly unlikely. We discuss several possible interpretations of the total body of observations of the Fomalhaut system and find that the interpretation that best matches the available data for the observed source is scattered light from a transient or semi-transient dust cloud. © 2012 The American Astronomical Society. All rights reserved.

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.

Model selection for prognostic time-to-event gene signature discovery with applications in early breast cancer data

Model selection between competing models is a key consideration in the discovery of prognostic multigene signatures. The use of appropriate statistical performance measures as well as verification of biological significance of the signatures is imperative to maximise the chance of external validation of the generated signatures. Current approaches in time-to-event studies often use only a single measure of performance in model selection, such as logrank test p-values, or dichotomise the follow-up times at some phase of the study to facilitate signature discovery. In this study we improve the prognostic signature discovery process through the application of the multivariate partial Cox model combined with the concordance index, hazard ratio of predictions, independence from available clinical covariates and biological enrichment as measures of signature performance. The proposed framework was applied to discover prognostic multigene signatures from early breast cancer data. The partial Cox model combined with the multiple performance measures were used in both guiding the selection of the optimal panel of prognostic genes and prediction of risk within cross validation without dichotomising the follow-up times at any stage. The signatures were successfully externally cross validated in independent breast cancer datasets, yielding a hazard ratio of 2.55 [1.44, 4.51] for the top ranking signature.

Geochemical records of palaeoenvironmental controls on peat forming processes in the Mfabeni peatland, Kwazulu Natal, South Africa since the Late Pleistocene

The Mfabeni peatland is the only known sub-tropical coastal fen that transcends the Last Glacial Maximum (LGM). This ca. 10m thick peat sequence provides a continuous sedimentation record spanning from the late Pleistocene to present (basal age c. 47kcalyr BP). We investigated the paleaeoenvironmental controls on peat formation and organic matter source input at the Mfabeni fen by: 1) exploring geochemical records (mass accumulation rate, total organic carbon, carbon accumulation rate, δC, δN and C/N ratio) to delineate primary production, organic matter source input, preservation and diagenetic processes, and 2) employ these geochemical signatures to reconstruct the palaeoenvironmental conditions and prevailing climate that drove carbon accumulation in the peatland. We established that the Mfabeni peat sediments have undergone minimal diagenetic alteration. The peat sequence was divided into 5 linear sedimentation rate (LSR) stages indicating distinct changes in climate and hydrological conditions: LSR stage 1 (c. 47 to c. 32.2kcalyr BP): predominantly cool and wet climate with C4 plant assemblages, interrupted by two short warming events. LSR stage 2 (c. 32.2 to c. 27.6kcalyr BP): dry and windy climate followed by a brief warm and wet period with increased C4 sedge swamp vegetation. LSR stage 3 (c. 27.6 to c. 20.3kcalyr BP): initial cool and wet period with prevailing C4 sedge plant assemblage until c. 23kcalyr BP; then an abrupt change to dry and cool glacial conditions and steady increases in C3 grasses. LSR stage 4 (c. 20.3 to c. 10.4kcalyr BP): continuation of cool and dry conditions and strong C3 grassland signature until c. 15kcalyr BP, after which precipitation increases. LSR stage 5 (c. 10.4kcalyr BP to present): characterised by extreme fluctuations between pervasive wet and warm to cool interglacial conditions with intermittent abrupt millennial-scale cooling/drying events and oscillations between C3 and C4 plant assemblages. In this study we reconstructed a high-resolution record of local hydrology, bulk plant assemblage and inferred climate since the Late Pleistocene, which suggest an anti-phase link between Southern African and the Northern Hemisphere, most notably during Heinrich (5 to 2) and Younger Dryas events. © 2013 Elsevier B.V.

Theory of phonon dissipation in the conduction of stressed Au nanowires

Recent experiments on Au break junctions [Phys. Rev. Lett. 88 (2002) 216803] have characterized the nonlinear conductance of stretched short Au nanowires. They reveal in the voltage range 10-20 meV the signatures of dissipation effects, likely due to phonons in the nanowire, reducing the conductance below the quantized value of 2e(2)/h. We present here a theory, based on a model tight-binding Hamiltonian and on non-equilibrium Green's function techniques, which accounts for the main features of the experiment. The theory helps in revealing details of the experiment which need to be addressed with a more realistic, less idealized, theoretical framework. (C) 2004 Elsevier B.V. All rights reserved.

Exploring Vertex Interactions in Ferroelectric Flux-Closure Domains

Using piezoresponse force microscopy, we have observed the progressive development of ferroelectric flux-closure domain structures and Landau−Kittel-type domain patterns, in 300 nm thick single-crystal BaTiO3 platelets. As the microstructural development proceeds, the rate of change of the domain configuration is seen to decrease exponentially. Nevertheless, domain wall velocities throughout are commensurate with creep processes in oxide ferroelectrics. Progressive screening of macroscopic destabilizing fields, primarily the surface-related depolarizing field, successfully describes the main features of the observed kinetics. Changes in the separation of domain-wall vertex junctions prompt a consideration that vertex−vertex interactions could be influencing the measured kinetics. However, the expected dynamic signatures associated with direct vertex−vertex interactions are not resolved. If present, our measurements confine the length scale for interaction between vertices to the order of a few hundred nanometers.

Interacting supernovae from photoionization-confined shells around red supergiant stars

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.