Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.

Identification of stimulatory and inhibitory inputs to the hypothalamic-pituitary-adrenal axis during hypoglycaemia or transport in ewes

This study used the novel approach of statistical modelling to investigate the control of hypothalamic-pituitary-adrenal (HPA) axis and quantify temporal relationships between hormones. Two experimental paradigms were chosen, insulin-induced hypoglycaemia and 2 h transport, to assess differences in control between noncognitive and cognitive stimuli. Vasopressin and corticotropin-releasing hormone (CRH) were measured in hypophysial portal plasma, and adrenocorticotropin hormone (ACTH) and cortisol in jugular plasma of conscious sheep, and deconvolution analysis was used to calculate secretory rates, before modelling. During hypoglycaemia, the relationship between plasma glucose and vasopressin or CRH was best described by log(10) transforming variables (i.e. a positive power-curve relationship). A negative-feedback relationship with log(10) cortisol concentration 2 h previously was detected. Analysis of the 'transport' stimulus suggested that the strength of the perceived stimulus decreased over time after accounting for cortisol facilitation and negative-feedback. The time course of vasopressin and CRH responses to each stimulus were different However, at the pituitary level, the data suggested that log(10) ACTH secretion rate was related to log(10) vasopressin and CRH concentrations with very similar regression coefficients and an identical ratio of actions (2.3 : 1) for both stimuli. Similar magnitude negative-feedback effects of log(10) cortisol at -110 min (hypoglycaemia) or -40 min (transport) were detected, and both models contained a stimulatory relationship with cortisol at 0 min (facilitation). At adrenal gland level, cortisol secretory rates were related to simultaneously measured untransformed ACTH concentration but the regression coefficient for the hypoglycaemia model was 2.5-fold greater than for transport. No individual sustained maximum cortisol secretion for longer than 20 min during hypoglycaemia and 40 min during transport. These unique models demonstrate that corticosteroid negative-feedback is a significant control mechanism at both the pituitary and hypothalamus. The amplitude of HPA response may be related to stimulus intensity and corticosteroid negative-feedback, while duration depended on feedback alone.

Extraction behavior of the synergistic system 2,6-bis-(benzoxazolyl)-4-dodecyloxylpyridine and 2-bromodecanoic acid using Am and Eu as radioactive tracers

In this study, the extraction properties of a synergistic system consisting of 2,6-bis-(benzoxazolyl)-4-dodecyloxylpyridine (BODO) and 2-bromodecanoic acid (HA) in tert-butyl benzene (TBB) have been investigated as a function of ionic strength by varying the nitrate ion and perchlorate ion concentrations. The influence of the hydrogen ion concentration has also been investigated. Distribution ratios between 0.03-12 and 0.003-0.8 have been found for Am(III) and Eu(HI), respectively, but there were no attempts to maximize these values. It has been shown that the distribution ratios decrease with increasing amounts of ClO4-, NO3-, and H+. The mechanisms, however, by which the decrease occurs, are different. In the case of increasing perchlorate ion concentration, the decrease in extraction is linear in a log-log plot of the distribution ratio vs. the ionic strength, while in the nitrate case the complexation between nitrate and Am or Eu increases at high nitrate ion concentrations and thereby decreases the distribution ratio in a non-linearway. The decrease in extraction could be caused by changes in activity coefficients that can be explained with specific ion interaction theory (SIT); shielding of the metal ions, and by nitrate complexation with Am and Eu as competing mechanism at high ionic strengths. The separation factor between Am and Eu reaches a maximum at similar to1 M nitrate ion concentration. Thereafter the values decrease with increasing nitrate ion concentrations.

A comparison of the reactivity of germylene and dimethylgermylene with some methylgermanes. Direct kinetic and quantum chemical studies

Time-resolved studies of germylene, GeH2, and dimethygermylene, GeMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to try to obtain rate coefficients for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas-phase. GeH2 + Me2GeH2 was studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 296-553 K. Only slight pressure dependences were found (at 386, 447 and 553 K). RRKM modelling was carried out to fit these pressure dependences. The high pressure rate coefficients gave the Arrhenius parameters: log(A/cm(3) molecule(-1)s(-1)) = -10.99 +/- 0.07 and E-a = -(7.35 +/- 0.48) kJ mol(-1). No reaction could be found between GeMe2 + Me2GeH2 at any temperature up to 549 K, and upper limits of ca. 10(-14) cm(3) molecule(-1)s(-1) were set for the rate coefficients. A rate coefficient of (1.33 +/- 0.04) x 10(-11)cm(3) molecule(-1)s(-1) was also obtained for GeH2 + MeGeH3 at 296 K. No reaction was found between GeMe2 and MeGeH3. Rate coefficient comparisons showed, inter alia, that in the substrate germane Me-for-H substitution increased the magnitudes of rate coefficients significantly, while in the germylene Me-for-H substitution decreased the magnitudes of rate coefficients by at least four orders of magnitude. Quantum chemical calculations (G2(MP2,SVP)// B3LYP level) supported these findings and showed that the lack of reactivity of GeMe2 is caused by a positive energy barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper.

Surprisingly slow reaction of dimethylsilylene with dimethylgermane: time-resolved kinetic studies and related quantum chemical calculations

Time-resolved studies of silylene, SiH2, and dimethylsilylene, SiMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to obtain rate constants for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas phase. SiMe2 + Me2GeH2 was studied at five temperatures in the range 299-555 K. Problems of substrate UV absorption at 193 nm at temperatures above 400 K meant that only three temperatures could be used reliably for rate constant measurement. These rate constants gave the Arrhenius parameters log(A/cm(3) molecule(-1) s(-1)) = -13.25 +/- 0.16 and E-a = -(5.01 +/- 1.01) kJ mol(-1). Only room temperature studies of SiH2 were carried out. These gave values of (4.05 +/- 0.06) x 10(-10) cm(3) molecule(-1) s(-1) (SiH2 + Me2GeH2 at 295 K) and also (4.41 +/- 0.07) x 10(-10) cm(3) molecule(-1) s(-1) (SiH2 + MeGeH3 at 296 K). Rate constant comparisons show the surprising result that SiMe2 reacts 12.5 times slower with Me2GeH2 than with Me2SiH2. Quantum chemical calculations (G2(MP2,SVP)//B3LYP level) of the model Si-H and Ge-H insertion processes of SiMe2 with SiH4/MeSiH3 and GeH4/MeGeH3 support these findings and show that the lower reactivity of SiMe2 with Ge-H bonds is caused by a higher secondary barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper. Other, related, comparisons of silylene reactivity are also presented.

Direct time-resolved study of the gas-phase reactions of germylene with ethyl- and diethylgermane: absolute rate constants, temperature dependences, and mechanism

Time-resolved studies of germylene, GeH2, generated by the 193 nm laser flash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reactions with ethyl- and diethylgermanes in the gas phase. The reactions were studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 297-564 K. Only slight pressure dependences were found for GeH2 + EtGeH3 (399, 486, and 564 K). The high pressure rate constants gave the following Arrhenius parameters: for GeH2 + EtGeH3, log A = -10.75 +/- 0.08 and E-a = -6.7 +/- 0.6 kJ mol(-1); for GeH2 + Et2GeH2, log A = -10.68 +/- 0.11 and E-a = -6.95 +/- 0.80 kJ mol(-1). These are consistent with fast, near collision-controlled, association processes at 298 K. RRKM modeling calculations are, for the most part, consistent with the observed pressure dependence of GeH2 + EtGeH3. The ethyl substituent effects have been extracted from these results and are much larger than the analogous methyl substituent effects in the SiH2 + methylsilane reaction series. This is consistent with a mechanistic model for Ge-H insertion in which the intermediate complex has a sizable secondary barrier to rearrangement.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with oxygen

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with O-2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second order rate constants at 10 Torr were fitted to the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.08 +/- 0.04) + (1.57 +/- 0.32 kJ mol(-1))/RT ln10 The decrease in rate constant values with increasing temperature, although systematic is very small. The rate constants showed slight increases in value with pressure at each temperature, but this was scarcely beyond experimental uncertainty. From estimates of Lennard-Jones collision rates, this reaction is occurring at ca. 1 in 20 collisions, almost independent of pressure and temperature. Ab initio calculations at the G3 level backed further by multi-configurational (MC) SCF calculations, augmented by second order perturbation theory (MRMP2), support a mechanism in which the initial adduct, H2SiOO, formed in the triplet state (T), undergoes intersystem crossing to the more stable singlet state (S) prior to further low energy isomerisation processes leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are H2O + SiO. The decomposition of the intermediate cyclo-siladioxirane, via O-O bond fission, plays an important role in the overall process. The bottleneck for the overall process appears to be the T -> S process in H2SiOO. This process has a small spin orbit coupling matrix element, consistent with an estimate of its rate constant of 1 x 10(9) s(-1) obtained with the aid of RRKM theory. This interpretation preserves the idea that, as in its reactions in general, SiH2 initially reacts at the encounter rate with O-2. The low values for the secondary reaction barriers on the potential energy surface account for the lack of an observed pressure dependence. Some comparisons are drawn with the reactions of CH2 + O-2 and SiCl2 + O-2.

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with nitric oxide

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with NO. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 299-592 K. The second-order rate constants at 10 Torr fitted the Arrhenius equation log(k/cm(3) molecule(-1) s(-1)) = (- 11.66 +/- 0.01) + (6.20 +/- 0.10 kJ mol(-1))IRT In 10 The rate constants showed a variation with pressure of a factor of ca. 2 over the available range, almost independent of temperature. The data could not be fitted by RRKM calculations to a simple third body assisted association reaction alone. However, a mechanistic model with an additional (pressure independent) side channel gave a reasonable fit to the data. Ab initio calculations at the G3 level supported a mechanism in which the initial adduct, bent H2SiNO, can ring close to form cyclo-H2SiNO, which is partially collisionally stabilized. In addition, bent H2SiNO can undergo a low barrier isomerization reaction leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are NH2 + SiO. The rate controlling barrier for this latter pathway is only 16 kJ mol(-1) below the energy of SiH2 + NO. This is consistent with the kinetic findings. A particular outcome of this work is that, despite the pressure dependence and the effects of the secondary barrier (in the side reaction), the initial encounter of SiH2 with NO occurs at the collision rate. Thus, silylene can be as reactive with odd electron molecules as with many even electron species. Some comparisons are drawn with the reactions of CH2 + NO and SiCl2 + NO.

The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD2, generated by laser flash photolysis of phenylsilane-d(3), have been carried out to obtain rate constants for its bimolecular reaction with C2H2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol(-1))/RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH2 with C2H2 and C2D2. Additionally, pressure-dependent rate constants for the reaction of SiH2 with C2H2 in the presence of He (1-100 Tort) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC2H4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH2 + C2H2 - SWPO + C2H4. The results are compared and contrasted with previous studies of this reaction system.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 1. HCl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with HCL The reaction was studied in the gas phase at 10 Torr total pressure in SF6 bath gas, at five temperatures in the range of 296-611 K. The second-order rate constants fitted the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.51 +/- 0.06) + (1.92 +/- 0.47 kJ mol(-1))/RTIn10 Experiments at other pressures showed that these rate constants were unaffected by pressure in the range of 10-100 Torr, but showed small decreases in value of no more than 20% ( +/- 10%) at I Toff, at both the highest and lowest temperatures. The data are consistent with formation of an initial weakly bound donor-acceptor complex, which reacts by two parallel pathways. The first is by chlorine-to-silicon H-shift to make vibrationally excited chlorosilane, SiH3Cl*, which yields HSiCl by H-2 elimination from silicon. In the second pathway, the complex proceeds via H-2 elimination (4-center process) to make chlorosilylene, HSiCl, directly. This interpretation is supported by ab initio quantum calculations carried out at the G3 level which reveal the direct H-2 elimination route for the first time. RRKM modeling predicts the approximate magnitude of the pressure effect but is unable to determine the proportions of each pathway. The experimental data agree with the only previous measurements at room temperature. Comparisons with other reactions of SiH2 are also drawn.

Mapping antimalarial pharmacophores as a useful tool for the rapid discovery of drugs effective in vivo: design, construction, characterization, and pharmacology of metaquine

Resistant strains of Plasmodium falciparum and the unavailability of useful antimalarial vaccines reinforce the need to develop new efficacious antimalarials. This study details a pharmacophore model that has been used to identify a potent, soluble, orally bioavailable antimalarial bisquinoline, metaquine (N,N'-bis(7-chloroquinolin-4-yl)benzene-1,3-diamine) (dihydrochloride), which is active against Plasmodium berghei in vivo (oral ID50 of 25 mu mol/kg) and multidrug-resistant Plasmodium falciparum K1 in vitro (0.17 mu M). Metaquine shows strong affinity for the putative antimalarial receptor, heme at pH 7.4 in aqueous DMSO. Both crystallographic analyses and quantum mechanical calculations (HF/6-31+G*) reveal important regions of protonation and bonding thought to persist at parasitic vacuolar pH concordant with our receptor model. Formation of drug-heme adduct in solution was confirmed using high-resolution positive ion electrospray mass spectrometry. Metaquine showed strong binding with the receptor in a 1: 1 ratio (log K = 5.7 +/- 0.1) that was predicted by molecular mechanics calculations. This study illustrates a rational multidisciplinary approach for the development of new 4-aminoquinoline antimalarials, with efficacy superior to chloroquine, based on the use of a pharmacophore model.

Survivability of a probiotic Lactobacillus casei in the gastrointestinal tract of healthy human volunteers and its impact on the faecal microflora

Aim: The aim of this study was to measure the gastrointestinal survival of Lactobacillus casei and its impact on the gut microflora in healthy human volunteers. Methods and Results: Twenty healthy volunteers took part in a double-blind placebo-controlled probiotic feeding study (10 fed probiotic, 10 fed placebo). The probiotic was delivered in two 65 ml aliquots of fermented milk drink (FMD) daily for 21 days at a dose of 8.6 +/- 0.1 Log(10)Lact. casei CFU ml(-1) FMD. Faecal samples were collected before, during and after FMD or placebo consumption, and important groups of faecal bacteria enumerated by fluorescent in situ hybridization (FISH) using oligonucleotide probes targeting the 16S rRNA. The fed Lact. casei was enumerated using selective nutrient agar and colony identity confirmed by pulsed field gel electrophoresis. Seven days after ingestion of FMD, the Lact. casei was recovered from faecal samples taken from the active treatment group at 7.1 +/- 0.4 Log(10) CFU g(-1) faeces (mean +/- SD, n = 9) and numbers were maintained at this level until day 21. Lact. casei persisted in six volunteers until day 28 at 5.0 +/- 0.9 Log(10) CFU g(-1) faeces (mean +/- SD, n = 6). Numbers of faecal lactobacilli increased significantly upon FMD ingestion. In addition, the numbers of bifidobacteria were higher on days 7 and 21 than on days 0 and 28 in both FMD fed and placebo fed groups. Consumption of Lact. casei had little discernible effect on other bacterial groups enumerated. Conclusions: Daily consumption of FMD enabled a probiotic Lact. casei strain to be maintained in the gastrointestinal tract of volunteers at a stable relatively high population level during the probiotic feeding period. Significance and Impact of the Study: The study has confirmed that this probiotic version of Lact. casei survives well within the human gastrointestinal tract.

Relationship between sublethal injury and microbial inactivation by the combination of high hydrostatic pressure and citral or tert-butyl hydroquinone

The aim was to investigate (i) the occurrence of sublethal injury in Listeria monocytogenes, Escherichia coli, and Saccharomyces cerevisiae after high hydrostatic pressure (HHP) treatment as a function of the treatment medium pH and composition and (ii) the relationship between the occurrence of sublethal injury and the inactivating effect of a combination of HHP and two antimicrobial compounds, tert-butyl hydroquinone (TBHQ) and citral. The three microorganisms showed a high proportion of sublethally injured cells (up to 99.99% of the surviving population) after HHP. In E. coli and L. monocytogenes, the extent of inactivation and sublethal injury depended on the pH and the composition of the treatment medium, whereas in S. cerevisiae, inactivation and sublethal injury were independent of medium pH or composition under the conditions tested. TBHQ alone was not lethal to E. coli or L. monocytogenes but acted synergistically with HHP and 24-h refrigeration, resulting in a viability decrease of >5 log(10) cycles of both organisms. The antimicrobial effect of citral depended on the microorganism and the treatment medium pH. Acting alone for 24 h under refrigeration, 1,000 ppm of citral caused a reduction of 5 log(10) cycles of E. coli at pH 7.0 and almost 3 log(10) cycles of L. monocytogenes at pH 4.0. The combination of citral and HHP also showed a synergistic effect. Our results have confirmed that the detection of sublethal injury after HHP may contribute to the identification of those treatment conditions under which HHP may act synergistically with other preserving processes.

An in vitro study of the probiotic potential of a bile-salt-hydrolyzing Lactobacillus fermentum strain, and determination of its cholesterol-lowering properties

This study evaluated the use of a bile-salt-hydrolyzing Lactobacillus fermentum strain as a probiotic with potential hypocholesterolemic properties. The effect of L. fermentum on representative microbial populations and overall metabolic activity of the human intestinal microbiota was investigated using a three-stage continuous culture system. Also, the use of galactooligosaccharides as a prebiotic to enhance growth and/or activity of the Lactobacillus strain was evaluated. Administration of L. fermentum resulted in a decrease in the overall bifidobacterial population (ca. 1 log unit). In the in vitro system, no significant changes were observed in the total bacterial, Lactobacillus, Bacteroides, and clostridial populations through L. fermentum supplementation. Acetate production decreased by 9 to 27%, while the propionate and butyrate concentrations increased considerably (50 to 90% and 52 to 157%, respectively). A general, although lesser, increase in the production of lactate was observed with the administration of the L. fermentum strain. Supplementation of the prebiotic to the culture medium did not cause statistically significant changes in either the numbers or the activity of the microbiota, although an increase in the butyrate production was seen (29 to 39%). Results from this in vitro study suggest that L.Fermentum KC5b is a candidate probiotic which may affect cholesterol metabolism. The short-chain fatty acid concentrations, specifically the molar proportion of propionate and/or bile salt deconjugation, are probably the major mechanism involved in the purported cholesterol-lowering properties of this strain.