Metabolomics Pilot Study Identifies Desynchronization of 24-H Rhythms and Distinct Intra-patient Variability Patterns in Critical Illness: A Preliminary Report.

Background: Synchronized circadian rhythms play a key role in coordinating physiologic health. Desynchronized circadian rhythms may predispose individuals to disease or be indicative of underlying disease. Intensive care unit (ICU) patients likely experience desynchronized circadian rhythms due to disruptive environmental conditions in the ICU and underlying pathophysiology. This observational pilot study was undertaken to determine if 24-h rhythms are altered in ICU patients relative to healthy controls by profiling 24-h rhythms in vital signs and plasma metabolites. Methods: We monitored daily rhythms in 5 healthy controls and 5 ICU patients for 24 h. Heart rate and blood pressure were measured every 30 min, temperature was measured every hour, and blood was sampled for mass spectrometry-based plasma metabolomics every 4 h. Bedside sound levels were measured every minute. Twenty-four hours rhythms were evaluated in vitals and putatively identified plasma metabolites individually and in each group using the cosinor method. Results: ICU patient rooms were significantly louder than healthy controls' rooms and average noise levels were above EPA recommendations. Healthy controls generally had significant 24-h rhythms individually and as a group. While a few ICU patients had significant 24-h rhythms in isolated variables, no significant rhythms were identified in ICU patients as a group, except in cortisol. This indicates a lack of coherence in phases and amplitudes among ICU patients. Finally, principal component analysis of metabolic profiles showed surprising patterns in plasma sample clustering. Each ICU patient's samples were clearly discernable in individual clusters, separate from a single cluster of healthy controls. Conclusions: In this pilot study, ICU patients' 24-h rhythms show significant desynchronization compared to healthy controls. Clustering of plasma metabolic profiles suggests that metabolomics could be used to track individual patients' clinical courses longitudinally. Our results show global disordering of metabolism and the circadian system in ICU patients which should be characterized further in order to determine implications for patient care.

"High-Throughput Characterization of Region-Specific Mitochondrial Function and Morphology".

The tissue-specific etiology of aging and stress has been elusive due to limitations in data processing of current techniques. Despite that many techniques are high-throughput, they usually use singular features of the data (e.g. whole fluorescence). One technology at the nexus of fluorescence-based screens is large particle flow cytometry ("biosorter"), capable of recording positional fluorescence and object granularity information from many individual live animals. Current processing of biosorter data, however, do not integrate positional information into their analysis and data visualization. Here, we present a bioanalytical platform for the quantification of positional information ("longitudinal profiling") of C. elegans, which we posit embodies the benefits of both high-throughput screening and high-resolution microscopy. We show the use of these techniques in (1) characterizing distinct responses of a transcriptional reporter to various stresses in defined anatomical regions, (2) identifying regions of high mitochondrial membrane potential in live animals, (3) monitoring regional mitochondrial activity in aging models and during development, and (4) screening for regulators of muscle mitochondrial dynamics in a high-throughput format. This platform offers a significant improvement in the quality of high-throughput biosorter data analysis and visualization, opening new options for region-specific phenotypic screening of complex physiological phenomena and mitochondrial biology.

Probing subunit-subunit interactions in the yeast vacuolar ATPase by peptide arrays.

BACKGROUND: Vacuolar (H(+))-ATPase (V-ATPase; V(1)V(o)-ATPase) is a large multisubunit enzyme complex found in the endomembrane system of all eukaryotic cells where its proton pumping action serves to acidify subcellular organelles. In the plasma membrane of certain specialized tissues, V-ATPase functions to pump protons from the cytoplasm into the extracellular space. The activity of the V-ATPase is regulated by a reversible dissociation mechanism that involves breaking and re-forming of protein-protein interactions in the V(1)-ATPase - V(o)-proton channel interface. The mechanism responsible for regulated V-ATPase dissociation is poorly understood, largely due to a lack of detailed knowledge of the molecular interactions that are responsible for the structural and functional link between the soluble ATPase and membrane bound proton channel domains. METHODOLOGY/PRINCIPAL FINDINGS: To gain insight into where some of the stator subunits of the V-ATPase associate with each other, we have developed peptide arrays from the primary sequences of V-ATPase subunits. By probing the peptide arrays with individually expressed V-ATPase subunits, we have identified several key interactions involving stator subunits E, G, C, H and the N-terminal domain of the membrane bound a subunit. CONCLUSIONS: The subunit-peptide interactions identified from the peptide arrays complement low resolution structural models of the eukaryotic vacuolar ATPase obtained from transmission electron microscopy. The subunit-subunit interaction data are discussed in context of our current model of reversible enzyme dissociation.

Engineering an artificial zymogen by alternate frame protein folding.

Alternate frame folding (AFF) is a novel mechanism by which allostery can be introduced into a protein where none may have existed previously. We employ this technology to convert the cytotoxic ribonuclease barnase into an artificial zymogen that is activated by HIV-1 protease. The AFF modification entails partial duplication of the polypeptide chain and mutation of a key catalytic residue in one of the duplicated segments. The resulting molecule can fold in one of two "frames" to yield the wild-type structure or a circularly permuted form in which the positions of the N- and C-termini are exchanged with a surface loop. It cannot take on both structures simultaneously because each competes for a shared amino acid sequence. An HIV-1 protease recognition sequence is inserted into one of the surface loops in the nonpermuted frame, and cleavage induces a shift from the nonpermuted fold to the permuted fold. Using the AFF mechanism, we were able to suppress k(cat)/K(M) by 250-fold in the proenzyme relative to wild-type barnase. HIV-1 protease cleavage subsequently increases k(cat)/K(M) by 130-fold. AFF is significant because it is general and can in principle be used to control activity of many enzymes, including those whose functions are not regulated by any existing mechanism.

Identification of conserved core xylem gene sets: conifer cDNA microarray development, transcript profiling and computational analyses.

One approach for investigating the molecular basis of wood formation is to integrate microarray profiling data sets and sequence analyses, comparing tree species with model plants such as Arabidopsis. Conifers may be included in comparative studies thanks to large-scale expressed sequence tag (EST) analyses, which enable the development of cDNA microarrays with very significant genome coverage. A microarray of 10,400 low-redundancy sequences was designed starting from white spruce (Picea glauca (Moench.) Voss) cDNAs. Computational procedures that were developed to ensure broad transcriptome coverage and efficient PCR amplification were used to select cDNA clones, which were re-sequenced in the microarray manufacture process. White spruce transcript profiling experiments that compared secondary xylem to phloem and needles identified 360 xylem-preferential gene sequences. The functional annotations of all differentially expressed sequences were highly consistent with the results of similar analyses carried out in angiosperm trees and herbaceous plants. Computational analyses comparing the spruce microarray sequences and core xylem gene sets from Arabidopsis identified 31 transcripts that were highly conserved in angiosperms and gymnosperms, in terms of both sequence and xylem expression. Several other spruce sequences have not previously been linked to xylem differentiation (including genes encoding TUBBY-like domain proteins (TLPs) and a gibberellin insensitive (gai) gene sequence) or were shown to encode proteins of unknown function encompassing diverse conserved domains of unknown function.

Automated SNP detection from a large collection of white spruce expressed sequences: contributing factors and approaches for the categorization of SNPs.

BACKGROUND: High-throughput genotyping technologies represent a highly efficient way to accelerate genetic mapping and enable association studies. As a first step toward this goal, we aimed to develop a resource of candidate Single Nucleotide Polymorphisms (SNP) in white spruce (Picea glauca [Moench] Voss), a softwood tree of major economic importance. RESULTS: A white spruce SNP resource encompassing 12,264 SNPs was constructed from a set of 6,459 contigs derived from Expressed Sequence Tags (EST) and by using the bayesian-based statistical software PolyBayes. Several parameters influencing the SNP prediction were analysed including the a priori expected polymorphism, the probability score (PSNP), and the contig depth and length. SNP detection in 3' and 5' reads from the same clones revealed a level of inconsistency between overlapping sequences as low as 1%. A subset of 245 predicted SNPs were verified through the independent resequencing of genomic DNA of a genotype also used to prepare cDNA libraries. The validation rate reached a maximum of 85% for SNPs predicted with either PSNP > or = 0.95 or > or = 0.99. A total of 9,310 SNPs were detected by using PSNP > or = 0.95 as a criterion. The SNPs were distributed among 3,590 contigs encompassing an array of broad functional categories, with an overall frequency of 1 SNP per 700 nucleotide sites. Experimental and statistical approaches were used to evaluate the proportion of paralogous SNPs, with estimates in the range of 8 to 12%. The 3,789 coding SNPs identified through coding region annotation and ORF prediction, were distributed into 39% nonsynonymous and 61% synonymous substitutions. Overall, there were 0.9 SNP per 1,000 nonsynonymous sites and 5.2 SNPs per 1,000 synonymous sites, for a genome-wide nonsynonymous to synonymous substitution rate ratio (Ka/Ks) of 0.17. CONCLUSION: We integrated the SNP data in the ForestTreeDB database along with functional annotations to provide a tool facilitating the choice of candidate genes for mapping purposes or association studies.

Generation, annotation, analysis and database integration of 16,500 white spruce EST clusters.

BACKGROUND: The sequencing and analysis of ESTs is for now the only practical approach for large-scale gene discovery and annotation in conifers because their very large genomes are unlikely to be sequenced in the near future. Our objective was to produce extensive collections of ESTs and cDNA clones to support manufacture of cDNA microarrays and gene discovery in white spruce (Picea glauca [Moench] Voss). RESULTS: We produced 16 cDNA libraries from different tissues and a variety of treatments, and partially sequenced 50,000 cDNA clones. High quality 3' and 5' reads were assembled into 16,578 consensus sequences, 45% of which represented full length inserts. Consensus sequences derived from 5' and 3' reads of the same cDNA clone were linked to define 14,471 transcripts. A large proportion (84%) of the spruce sequences matched a pine sequence, but only 68% of the spruce transcripts had homologs in Arabidopsis or rice. Nearly all the sequences that matched the Populus trichocarpa genome (the only sequenced tree genome) also matched rice or Arabidopsis genomes. We used several sequence similarity search approaches for assignment of putative functions, including blast searches against general and specialized databases (transcription factors, cell wall related proteins), Gene Ontology term assignation and Hidden Markov Model searches against PFAM protein families and domains. In total, 70% of the spruce transcripts displayed matches to proteins of known or unknown function in the Uniref100 database (blastx e-value < 1e-10). We identified multigenic families that appeared larger in spruce than in the Arabidopsis or rice genomes. Detailed analysis of translationally controlled tumour proteins and S-adenosylmethionine synthetase families confirmed a twofold size difference. Sequences and annotations were organized in a dedicated database, SpruceDB. Several search tools were developed to mine the data either based on their occurrence in the cDNA libraries or on functional annotations. CONCLUSION: This report illustrates specific approaches for large-scale gene discovery and annotation in an organism that is very distantly related to any of the fully sequenced genomes. The ArboreaSet sequences and cDNA clones represent a valuable resource for investigations ranging from plant comparative genomics to applied conifer genetics.