Qualitative Analysis of the Interdisciplinary Interaction between Data Analysis Specialists and Novice Clinical Researchers

Background: The inherent complexity of statistical methods and clinical phenomena compel researchers with diverse domains of expertise to work in interdisciplinary teams, where none of them have a complete knowledge in their counterpart's field. As a result, knowledge exchange may often be characterized by miscommunication leading to misinterpretation, ultimately resulting in errors in research and even clinical practice. Though communication has a central role in interdisciplinary collaboration and since miscommunication can have a negative impact on research processes, to the best of our knowledge, no study has yet explored how data analysis specialists and clinical researchers communicate over time. Methods/Principal Findings: We conducted qualitative analysis of encounters between clinical researchers and data analysis specialists (epidemiologist, clinical epidemiologist, and data mining specialist). These encounters were recorded and systematically analyzed using a grounded theory methodology for extraction of emerging themes, followed by data triangulation and analysis of negative cases for validation. A policy analysis was then performed using a system dynamics methodology looking for potential interventions to improve this process. Four major emerging themes were found. Definitions using lay language were frequently employed as a way to bridge the language gap between the specialties. Thought experiments presented a series of ""what if'' situations that helped clarify how the method or information from the other field would behave, if exposed to alternative situations, ultimately aiding in explaining their main objective. Metaphors and analogies were used to translate concepts across fields, from the unfamiliar to the familiar. Prolepsis was used to anticipate study outcomes, thus helping specialists understand the current context based on an understanding of their final goal. Conclusion/Significance: The communication between clinical researchers and data analysis specialists presents multiple challenges that can lead to errors.

Integration and visualization of systems biology data in context of the genome

Background: High-density tiling arrays and new sequencing technologies are generating rapidly increasing volumes of transcriptome and protein-DNA interaction data. Visualization and exploration of this data is critical to understanding the regulatory logic encoded in the genome by which the cell dynamically affects its physiology and interacts with its environment. Results: The Gaggle Genome Browser is a cross-platform desktop program for interactively visualizing high-throughput data in the context of the genome. Important features include dynamic panning and zooming, keyword search and open interoperability through the Gaggle framework. Users may bookmark locations on the genome with descriptive annotations and share these bookmarks with other users. The program handles large sets of user-generated data using an in-process database and leverages the facilities of SQL and the R environment for importing and manipulating data. A key aspect of the Gaggle Genome Browser is interoperability. By connecting to the Gaggle framework, the genome browser joins a suite of interconnected bioinformatics tools for analysis and visualization with connectivity to major public repositories of sequences, interactions and pathways. To this flexible environment for exploring and combining data, the Gaggle Genome Browser adds the ability to visualize diverse types of data in relation to its coordinates on the genome. Conclusions: Genomic coordinates function as a common key by which disparate biological data types can be related to one another. In the Gaggle Genome Browser, heterogeneous data are joined by their location on the genome to create information-rich visualizations yielding insight into genome organization, transcription and its regulation and, ultimately, a better understanding of the mechanisms that enable the cell to dynamically respond to its environment.

Suppression subtractive hybridization analysis reveals expression of conserved and novel genes in male accessory glands of the ant Leptothorax gredleri

Background: During mating, insect males eject accessory gland proteins (Acps) into the female genital tract. These substances are known to affect female post-mating behavior and physiology. In addition, they may harm the female, e. g., in reducing its lifespan. This is interpreted as a consequence of sexual antagonistic co-evolution. Whereas sexual conflict abounds in non-social species, the peculiar life history of social insects (ants, bees, wasps) with lifelong pair-bonding and no re-mating aligns the reproductive interests of the sexes. Harming the female during mating would negatively affect male fitness and sexual antagonism is therefore not expected. Indeed, mating appears to increase female longevity in at least one ant species. Acps are presumed to play a role in this phenomenon, but the underlying mechanisms are unknown. In this study, we investigated genes, which are preferentially expressed in male accessory glands of the ant Leptothorax gredleri, to determine which proteins might be transferred in the seminal fluid. Results: By a suppression subtractive hybridization protocol we obtained 20 unique sequences (USs). Twelve had mutual best matches with genes predicted for Apis mellifera and Nasonia vitripennis. Functional information (Gene Ontology) was available only for seven of these, including intracellular signaling, energy-dependent transport and metabolic enzyme activities. The remaining eight USs did not match sequences from other species. Six genes were further analyzed by quantitative RT-PCR in three life cycle stages of male ants. A gene with carboxy-lyase activity and one of unpredicted function were significantly overexpressed in accessory glands of sexually mature males. Conclusions: Our study is the first one to investigate differential gene expression in ants in a context related to mating. Our findings indicate that male accessory glands of L. gredleri express a series of genes that are unique to this species, possibly representing novel genes, in addition to conserved ones for which functions can be predicted. Identifying differentially expressed genes might help to better understand molecular mechanisms involved in reproductive processes in eusocial Hymenoptera. While the novel genes could account for rapidly evolving ones driven by intra-sexual conflict between males, conserved genes imply that rather beneficial traits might get fixed by a process described as inter-sexual cooperation between males and females.

VLT-FLAMES analysis of 8 giants in the bulge metal-poor globular cluster NGC 6522: oldest cluster in the Galaxy? Analysis of 8 giants in NGC 6522

Context. NGC 6522 has been the first metal-poor globular cluster identified in the bulge by Baade. Despite its importance, very few high-resolution abundance analyses of stars in this cluster are available. The bulge metal-poor clusters may be important tracers of the early chemical enrichment of the Galaxy. Aims. The main purpose of this study is to determine metallicity and elemental ratios in individual stars of NGC 6522. Methods. High-resolution spectra of 8 giants of the bulge's globular cluster NGC 6522 were obtained at the 8m VLT UT2-Kueyen telescope with the FLAMES+GIRAFFE spectrograph. Multiband V, I, J, K(s) photometry was used to derive effective temperatures as reference values. Spectroscopic parameters were derived from Fe I and Fe II lines, and adopted for the derivation of abundance ratios. Results. The present analysis provides a metallicity [Fe/H] = -1.0 +/- 0.2. The alpha-elements oxygen, magnesium and silicon show [O/Fe] = +0.4 +/- 0.3, [Mg/Fe] = [Si/Fe] = +0.25 +/- 0.15, whereas calcium and titanium show shallower ratios of [Ca/Fe] = [Ti/Fe] = +0.15 +/- 0.15. The neutron-capture r-process element europium appears to be overabundant by [Eu/Fe] = +0.4 +/- 0.4. The neutron-capture s-elements lanthanum and barium are enhanced by [La/Fe] = +0.35 +/- 0.2 and [Ba/Fe] = +0.5 +/- 0.5. The large internal errors, indicating the large star-to-star variation in the barium and europium abundances, are also discussed. Conclusions. The moderate metallicity combined to a blue horizontal branch (BHB), are characteristics similar to those of HP 1 and NGC 6558, pointing to a population of very old globular clusters in the Galactic bulge. Also, the abundance ratios in NGC 6522 resemble those in HP 1 and NGC 6558. The ultimate conclusion is that the bulge is old, and went through an early prompt chemical enrichment.

The B0.5IVe CoRoT target HD 49330 I. Photometric analysis from CoRoT data

Context. Be stars undergo outbursts producing a circumstellar disk from the ejected material. The beating of non-radial pulsations has been put forward as a possible mechanism of ejection. Aims. We analyze the pulsational behavior of the early B0.5IVe star HD 49330 observed during the first CoRoT long run towards the Galactical anticenter (LRA1). This Be star is located close to the lower edge of the beta Cephei instability strip in the HR diagram and showed a 0.03 mag outburst during the CoRoT observations. It is thus an ideal case for testing the aforementioned hypothesis. Methods. We analyze the CoRoT light curve of HD 49330 using Fourier methods and non-linear least square fitting. Results. In this star, we find pulsation modes typical of beta Cep stars (p modes) and SPB stars (g modes) with amplitude variations along the run directly correlated with the outburst. These results provide new clues about the origin of the Be phenomenon as well as strong constraints on the seismic modelling of Be stars.

Kinematic analysis and membership status of TWA22 AB

Context. TWA22 was initially regarded as a member of the TW Hydrae association (TWA). In addition to being one of the youngest (approximate to 8 Myr) and nearest (approximate to 20 pc) stars to Earth, TWA22 has proven to be very interesting after being resolved as a tight, very low-mass binary. This binary can serve as a very useful dynamical calibrator for pre-main sequence evolutionary models. However, its membership in the TWA has been recently questioned despite due to the lack of accurate kinematic measurements. Aims. Based on proper motion, radial velocity, and trigonometric parallax measurements, we aim here to re-analyze the membership of TWA22 to young, nearby associations. Methods. Using the ESO NTT/SUSI2 telescope, we observed TWA22 AB during 5 different observing runs over 1.2 years to measure its trigonometric parallax and proper motion. This is a part of a larger project measuring trigonometric parallaxes and proper motions of most known TWA members at a sub-milliarcsec level. HARPS at the ESO 3.6 m telescope was also used to measure the system's radial velocity over 2 years. Results. We report an absolute trigonometric parallax of TWA22 AB, pi = 57.0 +/- 0.7 mas, corresponding to a distance 17.5 +/- 0.2 pc from Earth. Measured proper motions of TWA 22AB are mu(alpha) cos(delta) = -175.8 +/- 0.8 mas/yr and mu delta = -21.3 +/- 0.8 mas/yr. Finally, from HARPS measurements, we obtain a radial velocity V(rad) = 14.8 +/- 2.1 km s(-1). Conclusions. A kinematic analysis of TWA22 AB space motion and position implies that a membership of TWA22 AB to known young, nearby associations can be excluded except for the beta Pictoris and TW Hydrae associations. Membership probabilities based on the system's Galactic space motion and/or the trace-back technique support a higher chance of being a member to the beta Pictoris association. Membership of TWA22 in the TWA cannot be fully excluded because of large uncertainties in parallax measurements and radial velocities and to the uncertain internal velocity dispersion of its members.

PARTITIONED ANALYSIS FOR DIMENSIONALLY-HETEROGENEOUS HYDRAULIC NETWORKS

In this work an iterative strategy is developed to tackle the problem of coupling dimensionally-heterogeneous models in the context of fluid mechanics. The procedure proposed here makes use of a reinterpretation of the original problem as a nonlinear interface problem for which classical nonlinear solvers can be applied. Strong coupling of the partitions is achieved while dealing with different codes for each partition, each code in black-box mode. The main application for which this procedure is envisaged arises when modeling hydraulic networks in which complex and simple subsystems are treated using detailed and simplified models, correspondingly. The potentialities and the performance of the strategy are assessed through several examples involving transient flows and complex network configurations.

Topology Studies of Hydrodynamics Using Two-Particle Correlation Analysis

The effects of fluctuating initial conditions are studied in the context of relativistic heavy ion collisions where a rapidly evolving system is formed. Two-particle correlation analysis is applied to events generated with the NEXSPHERIO hydrodynamic code, starting with fluctuating nonsmooth initial conditions (IC). The results show that the nonsmoothness in the IC survives the hydroevolution and can be seen as topological features of the angular correlation function of the particles emerging from the evolving system. A long range correlation is observed in the longitudinal direction and in the azimuthal direction a double peak structure is observed in the opposite direction to the trigger particle. This analysis provides clear evidence that these are signatures of the combined effect of tubular structures present in the IC and the proceeding collective dynamics of the hot and dense medium.

alpha plus alpha scattering reexamined in the context of the Sao Paulo potential

We have analyzed a large set of alpha + alpha elastic scattering data for bombarding energies ranging from 0.6 to 29.5 MeV. Because of the complete lack of open reaction channels, the optical interaction at these energies must have a vanishing imaginary part. Thus, this system is particularly important because the corresponding elastic scattering cross sections are very sensitive to the real part of the interaction. The data were analyzed in the context of the velocity-dependent Sao Paulo potential, which is a successful theoretical model for the description of heavy-ion reactions from sub-barrier to intermediate energies. We have verified that, even in this low-energy region, the velocity dependence of the model is quite important for describing the data of the alpha + alpha system.

Structure of an Odorant-Binding Protein from the Mosquito Aedes aegypti Suggests a Binding Pocket Covered by a pH-Sensitive ""Lid""

Background: The yellow fever mosquito, Aedes aegypti, is the primary vector for the viruses that cause yellow fever, mostly in tropical regions of Africa and in parts of South America, and human dengue, which infects 100 million people yearly in the tropics and subtropics. A better understanding of the structural biology of olfactory proteins may pave the way for the development of environmentally-friendly mosquito attractants and repellents, which may ultimately contribute to reduction of mosquito biting and disease transmission. Methodology: Previously, we isolated and cloned a major, female-enriched odorant-binding protein (OBP) from the yellow fever mosquito, AaegOBP1, which was later inadvertently renamed AaegOBP39. We prepared recombinant samples of AaegOBP1 by using an expression system that allows proper formation of disulfide bridges and generates functional OBPs, which are indistinguishable from native OBPs. We crystallized AaegOBP1 and determined its three-dimensional structure at 1.85 angstrom resolution by molecular replacement based on the structure of the malaria mosquito OBP, AgamOBP1, the only mosquito OBP structure known to date. Conclusion: The structure of AaegOBP1 (= AaegOBP39) shares the common fold of insect OBPs with six alpha-helices knitted by three disulfide bonds. A long molecule of polyethylene glycol (PEG) was built into the electron-density maps identified in a long tunnel formed by a crystallographic dimer of AaegOBP1. Circular dichroism analysis indicated that delipidated AaegOBP1 undergoes a pH-dependent conformational change, which may lead to release of odorant at low pH (as in the environment in the vicinity of odorant receptors). A C-terminal loop covers the binding cavity and this ""lid"" may be opened by disruption of an array of acid-labile hydrogen bonds thus explaining reduced or no binding affinity at low pH.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

GPIomics: global analysis of glycosylphosphatidylinositol-anchored molecules of Trypanosoma cruzi

Glycosylphosphatidylinositol (GPI) anchoring is a common, relevant posttranslational modification of eukaryotic surface proteins. Here, we developed a fast, simple, and highly sensitive (high attomole-low femtomole range) method that uses liquid chromatography-tandem mass spectrometry (LC-MS(n)) for the first large-scale analysis of GPI-anchored molecules (i.e., the GPIome) of a eukaryote, Trypanosoma cruzi, the etiologic agent of Chagas disease. Our genome-wise prediction analysis revealed that approximately 12% of T. cruzi genes possibly encode GPI-anchored proteins. By analyzing the GPIome of T. cruzi insect-dwelling epimastigote stage using LC-MS(n), we identified 90 GPI species, of which 79 were novel. Moreover, we determined that mucins coded by the T. cruzi small mucin-like gene (TcSMUG S) family are the major GPI-anchored proteins expressed on the epimastigote cell surface. TcSMUG S mucin mature sequences are short (56-85 amino acids) and highly O-glycosylated, and contain few proteolytic sites, therefore, less likely susceptible to proteases of the midgut of the insect vector. We propose that our approach could be used for the high throughput GPIomic analysis of other lower and higher eukaryotes. Molecular Systems Biology 7 April 2009; doi:10.1038/msb.2009.13

A greener and highly sensitive flow-based procedure for carbaryl determination exploiting long pathlength spectrophotometry and photochemical waste degradation

An environmentally friendly analytical procedure with high sensitivity for determination of carbaryl pesticide in natural waters was developed. The flow system was designed with solenoid micro-pumps in order to improve mixing conditions and minimize reagent consumption as well as waste generation. A long pathlength (100 cm) flow cell based on a liquid core waveguide (LCW) was employed to increase the sensitivity in detection of the indophenol formed from the reaction between carbaryl and p-aminophenol (PAP). A clean-up step based on cloud-point extraction was explored to remove the interfering organic matter, avoiding the use of toxic organic solvents. A linear response was observed within the range 5-200 mu g L(-1) and the detection limit, coefficient of variation and sampling rate were estimated as 1.7 mu g L(-1) (99.7% confidence level), 0.7% (n=20) and 55 determinations per hour, respectively. The reagents consumption was 1.9 mu g of PAP and 5.7 mu g of potassium metaperiodate, with volume of 2.6 mL of effluent per determination. The proposed procedure was selective for the determination of carbaryl, without interference from other carbamate pesticides. Recoveries within 84% and 104% were estimated for carbaryl spiked to water samples and the results obtained were also in agreement with those found by a batch spectrophotometric procedure at the 95% confidence level. The waste of the analytical procedure was treated with potassium persulphate and ultraviolet irradiation, yielding a colorless residue and a decrease of 94% of total organic carbon. In addition, the residue after treatment was not toxic for Vibrio fischeri bacteria. (c) 2010 Elsevier B.V. All rights reserved.

Single interface flow analysis with accuracy assessment

Single interface flow systems (SIFA) present some noteworthy advantages when compared to other flow systems, such as a simpler configuration, a more straightforward operation and control and an undemanding optimisation routine. Moreover, the plain reaction zone establishment, which relies strictly on the mutual inter-dispersion of the adjoining solutions, could be exploited to set up multiple sequential reaction schemes providing supplementary information regarding the species under determination. In this context, strategies for accuracy assessment could be favourably implemented. To this end, the sample could be processed by two quasi-independent analytical methods and the final result would be calculated after considering the two different methods. Intrinsically more precise and accurate results would be then gathered. In order to demonstrate the feasibility of the approach, a SIFA system with spectrophotometric detection was designed for the determination of lansoprazole in pharmaceutical formulations. Two reaction interfaces with two distinct pi-acceptors, chloranilic acid (CIA) and 2,3-dichloro-5,6-dicyano-p-benzoquinone (DDQ) were implemented. Linear working concentration ranges between 2.71 x 10(-4) to 8.12 x 10(-4) mol L(-1) and 2.17 x 10(-4) to 8.12 x 10(-4) mol L(-1) were obtained for DDQ and CIA methods, respectively. When compared with the results furnished by the reference procedure, the results showed relative deviations lower than 2.7%. Furthermore. the repeatability was good, with r.s.d. lower than 3.8% and 4.7% for DDQ and CIA methods, respectively. Determination rate was about 30 h(-1). (C) 2009 Elsevier B.V. All rights reserved.

Green Chemistry-Sensitive Analytical Procedure for Photometric Determination of Orthophosphate in River and Tap Water by Use of a Simple LED-Based Photometer

In the current work a Green Analytical Chemistry (GAC) procedure for photometric determination of orthophosphate in river water at mu g L-1 concentration level is described. The flow system module and the LED-based photometer were assembled together to constitute a compact unit in order to allow that a flow cell with optical path-length of 100mm was coupled to them. The photometric procedure based on the molybdenum blue method was implemented employing the multicommuted flow injection analysis approach, which provided facilities to allow reduction of reagent consumption and as well as waste generation. Aiming to prove the usefulness of the system, orthophosphate in river and tap waters was determined. Accuracy was ascertained by spiking samples with orthophosphate solution yielding recoveries ranging from 96% up to 107%. Other profitable features such as a wide linear response range between 10 to 800 mu g L-1 [image omitted]; a detection limit (3 sigma criterion) of 2.4 mu g L-1 [image omitted]; a relative standard deviation (n=7) of 2% using a typical water sample with concentration of 120 mu g L-1 [image omitted]; reagent consumption of 3.0mg ammonium molybdate, 0.3mg hydrazine sulfate, and 0.03mg stannous chloride per determination; a waste generation of 2.4mL per determination; and a sampling throughput of 20 determination per hours were also achieved.