New data on programmed risks of death in normal mice and mutants with growth delay

Study of the lifespans of normal (non mutant) mice and growth delay mutants has shown that mortality rates for both kinds of animals exhibit reproducible fluctuations. In the case of the mutant mice, the positions of peaks on the differential mortality curves (mortality rate plotted against lifespan) coincided in different-sex groups of animals and in same-sex subgroups of animals. Differential mortality curves of the mutant mice also had a peak at 1 month of age that was absent from the differential mortality curves of the normal mice. In the case of normal animals, positions of most peaks were the same in the studied independent subgroups of males, and to a lesser extent – independent subgroups of females, which might be explained by a shift in mortality peak positions due to the reproductive activity of females. Similar positions of mortality rate peaks in the differential mortality curves for animals from independent groups and subgroups indicate the existence of increased risks of death at specific ages. The observed pattern could be due to the programming in the genome of both the periods of increased risk of death and the intermitting intervals of stable development.     

Salt-concentration dependence of thermal denaturation of restriction fragment DNAs from phi X174

High-resolution differential melting curves of phi X174 Y1 and Y2 restriction fragment DNAs, for which the base sequences were known, were measured at various sodium ion concentrations ranging from 195 to 2.3 mM. The curves were resolved into component peaks, and the change in the melting temperature, the change in the area, and the change in the breadth of each peak with change in salt concentration were examined. The locations of the melting regions corresponding to the peaks in the melting curves were assigned based on theoretical calculations of melting curves and stability maps. It was found that as the salt concentration was decreased from the high to the intermediate range, the breadths of the peaks on the low-temperature side decreased whereas those on the high-temperature side remained almost constant, and also the separation between the peaks along the temperature axis increased. Changes in the positions of peaks relative to one another were interpreted in terms of the difference in the free energy increase between a loop state and an end-coil state as the salt concentration decreased.     

Peak quantification in surface-enhanced laser desorption/ionization by using mixture models

Surface-enhanced laser desorption/ionization (SELDI) time of flight (TOF) is a mass spectrometry technology for measuring the composition of a sampled protein mixture. A mass spectrum contains peaks corresponding to proteins in the sample. The peak areas are proportional to the measured concentrations of the corresponding proteins. Quantifying peak areas is difficult for existing methods because peak shapes are not constant across a spectrum and because peaks often overlap. We present a new method for quantifying peak areas. Our method decomposes a spectrum into peaks and a baseline using so-called statistical finite mixture models. We illustrate our method in detail on 8 samples from culture media of adipose tissue and globally on 64 samples from serum to compare our method to the standard Ciphergen method. Both methods give similar estimates for singleton peaks, but not for overlapping peaks. The Ciphergen method overestimates the heights of such peaks while our method still gives appropriate estimates. Peak quantification is an important step in pre-processing SELDI-TOF data and improvements therein will pay off in the later biomarker discovery phase.     

Robust protein quantitation in chromatographic fractions using MALDI-MS of tryptic peptides

A method is introduced to evaluate protein concentrations using the height sum of all MALDI-MS peaks that unambiguously match theoretic tryptic peptide masses of the protein sought after. The method uses native chromatographic protein fractionation prior to digestion but does not require any depletion, labeling, derivatization, or preparation of a compound similar to the analyte. All peak heights of tryptic peptides are normalized with the peak height of a unique standard peptide added to the MALDI-MS samples. The sum of normalized peak heights, S(n), or the normalized mean peak height, M(n), reflects the concentration of the respective protein. For fractions containing various proteins, S(n) and M(n) can be used to compare concentrations of a protein between different fractions. For fractions with one predominating protein, they can be used to estimate concentration ratios between fractions, or to quantify the fractional protein concentration after calibration with pure protein solutions. Initial native fractionation retains the possibility to apply all conventional analytic procedures. Moreover, it renders the method relatively robust to MS mass accuracy. The method was validated with albumin, transferrin, alpha1-antitrypsin, and immunoglobulin G within highly complex chromatographic fractions of pathological and normal sera, which contained the respective intact native protein in dominating as well as minor concentrations. The correlation found between S(n) and the protein concentration as determined with ELISA showed that the method can be applied to select markers for distinguishing between normal and pathological serum samples.     

Evolutionary genetics of lifespan and mortality rates in two populations of the seed beetle, Callosobruchus maculatus

The age at which individuals die varies substantially within and between species, but we still have little understanding of why there is such variation in life expectancy. We examined sex-specific and genetic variation in adult lifespan and the shape of mortality curves both within and between two populations of the seed beetle, Callosobruchus maculatus, that differ in a suite of life history characters associated with adaptation to different host species. Mean adult lifespan and the shape of the logistic mortality curves differed substantially between males and females (males had lower initial mortality rates, but a faster increase in the rate of mortality with increasing age) and between populations (they differed in the rate of increase in mortality with age). Larger individuals lived longer than smaller individuals, both because they had lower initial mortality rates and a slower increase in the rate of mortality with increasing age. However, differences in body size were not adequate to explain the differences in mortality between the sexes or populations. Both lifespan and mortality rates were genetically variable within populations and genetic variance/covariance matrices for lifespan differed between the populations and sexes. This study thus demonstrated substantial genetic variation in lifespan and mortality rates within and between populations of C. maculatus.     

pParse: a method for accurate determination of monoisotopic peaks in high-resolution mass spectra

Determining the monoisotopic peak of a precursor is a first step in interpreting mass spectra, which is basic but non-trivial. The reason is that in the isolation window of a precursor, other peaks interfere with the determination of the monoisotopic peak, leading to wrong mass-to-charge ratio or charge state. Here we propose a method, named pParse, to export the most probable monoisotopic peaks for precursors, including co-eluted precursors. We use the relationship between the position of the highest peak and the mass of the first peak to detect candidate clusters. Then, we extract three features to sort the candidate clusters: (i) the sum of the intensity, (ii) the similarity of the experimental and the theoretical isotopic distribution, and (iii) the similarity of elution profiles. We showed that the recall of pParse, MaxQuant, and BioWorks was 98-98.8%, 0.5-17%, and 1.8-36.5% at the same precision, respectively. About 50% of tandem mass spectra are triggered by multiple precursors which are difficult to identify. Then we design a new scoring function to identify the co-eluted precursors. About 26% of all identified peptides were exclusively from co-eluted peptides. Therefore, accurately determining monoisotopic peaks, including co-eluted precursors, can greatly increase peptide identification rate.     

Complex genetic architecture of population differences in adult lifespan of a beetle: nonadditive inheritance, gender differences, body size and a large maternal effect

Evolutionary responses to selection can be complicated when there is substantial nonadditivity, which limits our ability to extrapolate from simple models of selection to population differentiation and speciation. Studies of Drosophila melanogaster indicate that lifespan and the rate of senescence are influenced by many genes that have environment- and sex-specific effects. These studies also demonstrate that interactions among alleles (dominance) and loci (epistasis) are common, with the degree of interaction differing between the sexes and among environments. However, little is known about the genetic architecture of lifespan or mortality rates for organisms other than D. melanogaster. We studied genetic architecture of differences in lifespan and shapes of mortality curves between two populations of the seed beetle, Callosobruchus maculatus (South India and Burkina Faso populations). These two populations differ in various traits (such as body size and adult lifespan) that have likely evolved via host-specific selection. We found that the genetic architecture of lifespan differences between populations differs substantially between males and females; there was a large maternal effect on male lifespan (but not on female lifespan), and substantial dominance of long-life alleles in females (but not males). The large maternal effect in males was genetically based (there was no significant cytoplasmic effect) likely due to population differences in maternal effects genes that influence lifespan of progeny. Rearing host did not affect the genetic architecture of lifespan, and there was no evidence that genes on the Y-chromosome influence the population differences in lifespan. Epistatic interactions among loci were detectable for the mortality rate of both males and females, but were detectable for lifespan only after controlling for body size variation among lines. The detection of epistasis, dominance, and sex-specific genetic effects on C. maculatus lifespan is consistent with results from line cross and quantitative trait locus studies of D. melanogaster.     

A computer analysis of electrophoresis and ultracentrifugation patterns.

A program DIANA is designed to fit the scanning curves of electrophoresis and ultracentrifugation patterns with 10-15 overlapping peaks by a sum of Gaussian-like distribitions. The parameters of the distributions are adjusted by minimizing x, using a problem oriented minimizing procedure. The number of molecular components as well as starting values of the positions and the widths of the peaks must be choosen as input for each type of pattern. Programs are written in FORTRAN IV.     

Studies on the microstructure of bile salt aggregates in aqueous n-alkanol solutions by Small Angle Neutron Scattering

Small Angle Neutron Scattering (SANS) measurements on two bile salt micelle systems sodium cholate (NaC) and sodium deoxycholate (NaDC) each 0.1 M in the absence and presence of 0.05 M 1-butanol, 1-pentanol, 1-hexanol and 1-octanol were carried out in D20 at ambient conditions (25 degrees C). The scattering cross-section as a function of wave vector transfer (Q) showed the presence of correlation peaks characteristic of strong inter micellar interactions. The correlation peak positions were shifted to low Q values in the presence of n-alkanols for both NaC and NaDC micelles. Best fit curves of SANS were produced by applying Hayter-Penfold analysis considering monodisperse ellipsoid model for the micelles. The best fits were found for constant semi-minor axis = 8A with an increasing semi-major axis for increasing n-alkanol chain lengths. The micellar growths in presence of n-alkanols were studied using the ESR correlation time measurements on suitable spin probe incorporated micelles. The growth parameter and the hydrodynamic radii values were found to be agreeable with SANS data. Intermicellar interaction potentials seem to increase with the Cn of n-alkanols.     

A novel peak detection approach with chemical noise removal using short‐time FFT for prOTOF MS data

Peak detection is a pivotal first step in biomarker discovery from MS data and can significantly influence the results of downstream data analysis steps. We developed a novel automatic peak detection method for prOTOF MS data, which does not require a priori knowledge of protein masses. Random noise is removed by an undecimated wavelet transform and chemical noise is attenuated by an adaptive short‐time discrete Fourier transform. Isotopic peaks corresponding to a single protein are combined by extracting an envelope over them. Depending on the S/N, the desired peaks in each individual spectrum are detected and those with the highest intensity among their peak clusters are recorded. The common peaks among all the spectra are identified by choosing an appropriate cut‐off threshold in the complete linkage hierarchical clustering. To remove the 1 Da shifting of the peaks, the peak corresponding to the same protein is determined as the detected peak with the largest number among its neighborhood. We validated this method using a data set of serial peptide and protein calibration standards. Compared with MoverZ program, our new method detects more peaks and significantly enhances S/N of the peak after the chemical noise removal. We then successfully applied this method to a data set from prOTOF MS spectra of albumin and albumin‐bound proteins from serum samples of 59 patients with carotid artery disease compared to vascular disease‐free patients to detect peaks with S/N≥2. Our method is easily implemented and is highly effective to define peaks that will be used for disease classification or to highlight potential biomarkers.     

Presence of immunoreactive dynorphin in human cerebrospinal fluid

Immunoreactive dynorphin (I-dynorphin) was measured by radioimmunoassay in human cerebrospinal fluid (CSF). I-dynorphin concentration in CSF was 30 +/- 2 pg/ml. Sephadex G-25 column chromatography showed the main peak eluted at the position of dynorphin-(1-17). HPLC elution profile of this major peak from gel filtration showed a large peak corresponding to the position of dynorphin-(1-17) and small peaks corresponding to the positions of dynorphin-(1-13), dynorphin-(1-12) and other unknown peptides.     

Separation and partial characterization of two alkaline peptidases from cotyledons of resting kidney beans, Phaseolus vulagaris

Extracts prepared from cotyledons of resting kidney beans ( Phaseolus vulgaris L., cv. Processor) rapidly hydrolyzed two dipeptides, Leu-Tyr and Ala-Gly, with pH optima at 9.2 and 8.5, respectively. On ion exchange chromatography on DEAE-Sephacel the two activities eluted as separate peaks, showing that they were due to two different peptidases. The extracts also hydrolyzed Leu-β-naphthylamide optimally at pH 6.4; this activity eluted as a third peak berween the other peaks. The activity peak acting on Leu-Tyr and Ala-Gly rapidly hydrolized two tripeptides, showing that it was an aminopeptidase, whereas the Ala-Gly hydrolyzing peak acted only on dipeptides. The activities against Leu-Tyr and Ala-Gly were also separated by gel chromatography on Sephacryl S-300 with elution positions corresponding to M, values of about 360 000 and 105 000. The aminopeptidase was inhibited by bestatin, and the dipeptidase was inhibited by p-hydroxymercuribenzoate. Both enzymes were inhibited by o-phenanthroline. In most of their properties the two kidneys bean enzymes resembled the alkaline aminopeptidase and the dipeptidase earlier characterized from barley grains.     

Efficient unfolding pattern recognition in single molecule force spectroscopy data

Background

Single-molecule force spectroscopy (SMFS) is a technique that measures the force necessary to unfold a protein. SMFS experiments generate Force-Distance (F-D) curves. A statistical analysis of a set of F-D curves reveals different unfolding pathways. Information on protein structure, conformation, functional states, and inter- and intra-molecular interactions can be derived.

Results

In the present work, we propose a pattern recognition algorithm and apply our algorithm to datasets from SMFS experiments on the membrane protein bacterioRhodopsin (bR). We discuss the unfolding pathways found in bR, which are characterised by main peaks and side peaks. A main peak is the result of the pairwise unfolding of the transmembrane helices. In contrast, a side peak is an unfolding event in the alpha-helix or other secondary structural element. The algorithm is capable of detecting side peaks along with main peaks. Therefore, we can detect the individual unfolding pathway as the sequence of events labeled with their occurrences and co-occurrences special to bR's unfolding pathway. We find that side peaks do not co-occur with one another in curves as frequently as main peaks do, which may imply a synergistic effect occurring between helices. While main peaks co-occur as pairs in at least 50% of curves, the side peaks co-occur with one another in less than 10% of curves. Moreover, the algorithm runtime scales well as the dataset size increases.

Conclusions

Our algorithm satisfies the requirements of an automated methodology that combines high accuracy with efficiency in analyzing SMFS datasets. The algorithm tackles the force spectroscopy analysis bottleneck leading to more consistent and reproducible results.     

High-sensitivity scanning calorimetric study of mixtures of cholesterol with dimyristoyl- and dipalmitoylphosphatidylcholines

A highly sensitive and stable scanning microcalorimeter is employed in a reinvestigation of the effect of cholesterol on multilamellar suspensions of dimyristoylphosphatidylcholine (DMPC) and dipalmitoylphosphatidylcholine (DPPC). Below 20 mol % cholesterol the DPPC mixtures give heat-capacity curves each of which can be resolved into a narrow and a broad peak, suggesting the coexistence of two immiscible solid phases; above 20 mol % only the broad peak is observed and this disappears at about 50 mol %. The DMPC mixtures show a more complicated behavior; from about 13.5 to 20 mol % cholesterol the observed curves appear to be the sum of three component peaks. As with the DPPC mixtures, only a single broad peak is observed above 20 mol % cholesterol, and this broad peak becomes undetectable above about 50 mol %. These results are discussed.     

Potential bias in NMR relaxation data introduced by peak intensity analysis and curve fitting methods

We present an evaluation of the accuracy and precision of relaxation rates calculated using a variety of methods, applied to data sets obtained for several very different protein systems. We show that common methods of data evaluation, such as the determination of peak heights and peak volumes, may be subject to bias, giving incorrect values for quantities such as R₁ and R₂. For example, one common method of peak-height determination, using a search routine to obtain the peak-height maximum in successive spectra, may be a source of significant systematic error in the relaxation rate. The alternative use of peak volumes or of a fixed coordinate position for the peak height in successive spectra gives more accurate results, particularly in cases where the signal/noise is low, but these methods have inherent problems of their own. For example, volumes are difficult to quantitate for overlapped peaks. We show that with any method of sampling the peak intensity, the choice of a 2- or 3-parameter equation to fit the exponential relaxation decay curves can dramatically affect both the accuracy and precision of the calculated relaxation rates. In general, a 2-parameter fit of relaxation decay curves is preferable. However, for very low intensity peaks a 3 parameter fit may be more appropriate.     

Daily activity profiles over the lifespan of female medflies as biomarkers of aging and longevity

The relationship between the early-age activity of Mediterranean fruit flies (medflies) or other fruit flies and their lifespan has not been much studied, in contrast to the connections between lifespan and diet, sexual signaling, and reproduction. The objective of this study is to assess intra-day and day-to-day activity profiles of female Mediterranean fruit flies and their role as biomarker of longevity as well as to explore the relationships between these activity profiles, diet, and age-at-death throughout the lifespan. We use advanced statistical methods from functional data analysis (FDA). Three distinct patterns of activity variations in early-age activity profiles can be distinguished. A low-caloric diet is associated with a delayed activity peak, while a high-caloric diet is linked with an earlier activity peak. We find that age-at-death of individual medflies is connected to their activity profiles in early life. An increased risk of mortality is associated with increased activity in early age, as well as with a higher contrast between daytime and nighttime activity. Conversely, medflies are more likely to have a longer lifespan when they are fed a medium-caloric diet and when their daily activity is more evenly distributed across the early-age span and between daytime and nighttime. The before-death activity profile of medflies displays two characteristic before-death patterns, where one pattern is characterized by slowly declining daily activity and the other by a sudden decline in activity that is followed by death.     

Changes in the longevity of rats under the influence of enzyme imprinting by phenobarbital

Experiments on the Wistar rats showed that neonatal daily injections of phenobarbital (35 mg/kg) during first 3 days of life resulted in the enzyme imprinting of liver microsomal monooxygenases. Rise in the activities of liver microsomal oxidation enzymes is constantly maintained during all the life leading to increase in average lifespan of rats. Analysis of the survival curves in Gompertz equation co-ordinates showed that enzyme imprinting by phenobarbital caused changes in mortality patterns at different stages of ontogenesis. The phenomenon of enzyme imprinting by phenobarbital and lifespan prolongation was registered only in females but not in males. An inverse correlation was found between the duration of phenobarbital sleeping time and lifespan of rats.     

Improved peak detection and quantification of mass spectrometry data acquired from surface-enhanced laser desorption and ionization by denoising spectra with the undecimated discrete wavelet transform

Mass spectrometry is being used to find disease-related patterns in mixtures of proteins derived from biological fluids. Questions have been raised about the reproducibility and reliability of peak quantifications using this technology. We collected nipple aspirate fluid from breast cancer patients and healthy women, pooled them into a quality control sample, and produced 24 replicate SELDI spectra. We developed a novel algorithm to process the spectra, denoising with the undecimated discrete wavelet transform (UDWT), and evaluated it for consistency and reproducibility. UDWT efficiently decomposes spectra into noise and signal. The noise is consistent and uncorrelated. Baseline correction produces isolated peak clusters separated by flat regions. Our method reproducibly detects more peaks than the method implemented in Ciphergen software. After normalization and log transformation, the mean coefficient of variation of peak heights is 10.6%. Our method to process spectra provides improvements over existing methods. Denoising using the UDWT appears to be an important step toward obtaining results that are more accurate. It improves the reproducibility of quantifications and supplies tools for investigation of the variations in the technology more carefully. Further study will be required, because we do not have a gold standard providing an objective assessment of which peaks are present in the samples.     

Parallel inductive kinetics of fluorescence and photoacoustic signal in dark-adapted thalli of Bryopsis maxima

The inductive kinetics of fluorescence and photoacoustic signal were measured simultaneously in dark-adapted thalli of the green coenocytic alga Bryopsis maxima. Under illumination with weak red light modulated at 60 Hz, the fluorescence yield varied, showing three maxima P, M₁ and M₂ almost immediately, 10 s and 6 min after the onset of the illumination, respectively (Yamagishi, A., Satoh, K. and Katoh, S. (1978) Plant Cell Physiol. 19, 17–25). The photoacoustic signal also showed inductive transients which parallel well those of the fluorescence up to the M₂ stage. After M₂, the photoacoustic signal remained at a constant level, while the emission yield gradually decreased. The first peak of the fluorescence induction and a corresponding peak of the photoacoustic transients were selectively eliminated by prior illumination or methyl viologen treatment of the dark-adapted thalli. The second peaks of the two induction curves were abolished by carbonylcyanide-m-chlorophenylhydrazone, whereas dicyclohexylcarbodiimide enhanced their peak heights and suppressed the subsequent decreases. The results indicate that the fluorescence yield is mainly determined by the redox state of the Photosystem II reaction center throughout the induction period except the last phase. Mechanisms underlying inductive transients of fluorescence are discussed in the light of the present findings.     

Multi-temporal ecological analysis of Jeffrey pine beetle outbreak dynamics within the Lake Tahoe Basin

From 1991 to 1996, Jeffrey pine beetles (Dendroctonus jeffreyi Hopkins) (JPB) caused tree mortality throughout the Lake Tahoe Basin during a severe drought. Census data were collected annually on 10,721 trees to assess patterns of JPB-caused mortality. This represents the most extensive tree-level, spatiotemporal dataset collected to-date documenting bark beetle activity. Our study was an exploratory assessment of characteristics associated with the probability of successful JPB mass-attack (P_JPB) and group aggregation behavior that occurred throughout various outbreak phases. Numerous characteristics associated with P_JPB varied by outbreak phase although population pressure and forest density had positive associations during all phases. During the incipient phase, JPBs caused mortality in individual trees and small groups within toeslope topographic positions and P_JPB had a negative relationship with stem diameter. In the epidemic phase, JPB activity occurred in all topographic positions and caused mortality in spatially expanding clusters. P_JPB had a curvilinear relationship with tree diameter and a negative relationship with proximity to nearest brood tree. Majority (92–96 %) of mass-attacked trees were within 30 m of a brood tree during the peak epidemic years. During the post-epidemic phase, mortality clusters progressively decreased while dispersal distances between mass-attacked and brood trees increased. Post-epidemic P_JPB had a negative relationship with stem diameter and mortality was concentrated in the mid and upper-slope topographic positions. Results indicate mortality predictions are reasonable for the epidemic phase but not for incipient and post-epidemic phases. Ecological factors influencing JPB-caused tree mortality, clustered mortality patterns, and transitions from environmental to dynamic determinism are discussed.