Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within- family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

Molecular analysis of the pathway for the synthesis of thiol tripeptides in the model legume Lotus japonicus

The thiol tripeptides, glutathione (GSH) and homoglutathione (hGSH), perform multiple roles in legumes, including protection against toxicity of free radicals and heavy metals. The three genes involved in the synthesis of GSH and hGSH in the model legume, Lotus japonicus, have been fully characterized and appear to be present as single copies in the genome. The gamma-glutamylcysteine synthetase (gammaecs) gene was mapped on the long arm of chromosome 4 (70.0 centimorgans [cM]) and consists of 15 exons, whereas the glutathione synthetase (gshs) and homoglutathione synthetase (hgshs) genes were mapped on the long arm of chromosome 1 (81.3 cM) and found to be arranged in tandem, with a separation of approximately 8 kb. Both genes consist of 12 exons of exactly the same size (except exon 1, which is similar). Two types of transcripts were detected for the gshs gene, which putatively encode proteins localized in the plastids and cytosol. Promoter regions contain cis-acting regulatory elements that may be involved in the plant's response to light, hormones, and stress. Determination of transcript levels, enzyme activities, and thiol contents in nodules, roots, and leaves revealed that gammaecs and hgshs are expressed in all three plant organs, whereas gshs is significantly functional only in nodules. This strongly suggests an important role of GSH in the rhizobia-legume symbiosis.

Ratiometric and nonratiometric Ca2+ indicators for the assessment of intracellular free Ca2+ in a breast cancer cell line using a fluorescence microplate reader

Transporters of Ca2+ are potential drug targets and Ca2+ is a useful signal in the assessment of G-protein-coupled receptor activation. Assays involving the assessment of intracellular Ca2+ using microplate readers most often use Ca2+ indicators which do not exhibit a spectra shift on Ca2+ binding (e.g. fluo-3). Indicators that do exhibit a spectral shift upon Ca2+ binding (e.g. fura-2) offer potential advantages for the calibration of intracellular Ca2+ levels. However, experimental limitations may limit the use of ratiometric dyes in microplate readers capable of screening. In this study, we compared the assessment of intracellular Ca2+ in adherent breast cancer cells using ratiometric and nonratiometric Ca2+ indicators. Our results demonstrate that both fluo-3 and fura-2 detect ATP dose-dependent increases in intracellular Ca2+ in the MCF-7 breast cancer cell line and that some of the limitations in the use of fura-2 appear to be overcome by the use of glass bottom microplates. The calibrated intracellular Ca2+ levels derived using fura-2 are consistent with those from microscopy and cuvette-based studies. Fura-2 may be useful in microplate studies, where cell lines with different properties are compared or where screening treatments lead to differences in the number of cells or dye loading. (C) 2003 Elsevier B.V. All rights reserved.

The preparation and characterization of tin(IV) complexes of 2-quinolinecarboxaldehyde Schiff bases of S-methyl- and S-benzyldithiocarbazates and the X-ray crystal and molecular structures of the 2-quinolinecarboxaldehyde Schiff base of S-benzyldithi

New tin(IV) complexes of empirical formula, Sn(NNS)I-3 (NNS = anionic forms of the 2-quinolinecarboxaldehyde Schiff bases of S-methyl- and S-benzyldithiocarbazate) have been prepared and characterized by a variety of physico-chemical techniques. In the solid state, the Schiff bases exist as the thione tautomer but in solution and in the presence of tin(IV) iodide they convert to the thiol tautomer and coordinate to the tin atom in their deprotonated thiolate forms. The structures of the free ligand, Hqaldsbz and its triiodotin(IV) complex, [Sn(qaldsbz)I-3] have been determined by X-ray diffraction. The complex, [Sn(qaldsbz)I-3] has a distorted octahedral structure with the Schiff base coordinated to the tin atom as a uninegatively charged tridentate chelating agent via the quinoline nitrogen atom, the azomethine nitrogen atom and the thiolate sulfur atom. The three iodo ligands are coordinated meridionally to the tin atom. The distortion from an ideal octahedral geometry of [Sn(qaldsbz)I-3] is attributed to the restricted bite size of the tridentate Schiff base ligand. (C) 2004 Elsevier Ltd. All rights reserved.

Photodissociation of benzene under collision-free conditions: An ab initio/Rice-Ramsperger-Kassel-Marcus study

The ab initio/Rice-Ramsperger-Kassel-Marcus (RRKM) approach has been applied to investigate the photodissociation mechanism of benzene at various wavelengths upon absorption of one or two UV photons followed by internal conversion into the ground electronic state. Reaction pathways leading to various decomposition products have been mapped out at the G2M level and then the RRKM and microcanonical variational transition state theories have been applied to compute rate constants for individual reaction steps. Relative product yields (branching ratios) for C6H5+H, C6H4+H-2, C4H4+C2H2, C4H2+C2H4, C3H3+C3H3, C5H3+CH3, and C4H3+C2H3 have been calculated subsequently using both numerical integration of kinetic master equations and the steady-state approach. The results show that upon absorption of a 248 nm photon dissociation is too slow to be observable in molecular beam experiments. In photodissociation at 193 nm, the dominant dissociation channel is H atom elimination (99.6%) and the minor reaction channel is H-2 elimination, with the branching ratio of only 0.4%. The calculated lifetime of benzene at 193 nm is about 11 mus, in excellent agreement with the experimental value of 10 mus. At 157 nm, the H loss remains the dominant channel but its branching ratio decreases to 97.5%, while that for H-2 elimination increases to 2.1%. The other channels leading to C3H3+C3H3, C5H3+CH3, C4H4+C2H2, and C4H3+C2H3 play insignificant role but might be observed. For photodissociation upon absorption of two UV photons occurring through the neutral hot benzene mechanism excluding dissociative ionization, we predict that the C6H5+H channel should be less dominant, while the contribution of C6H4+H-2 and the C3H3+C3H3, CH3+C5H3, and C4H3+C2H3 radical channels should significantly increase. (C) 2004 American Institute of Physics.

Postsynthesis stabilization of free-standing mesoporous silica films

Mixed ammonia-water vapor postsynthesis treatment provides a simple and convenient method for stabilizing mesostructured silica films. X-ray diffraction, transmission electron microscopy, nitrogen adsorption/desorption, and solid-state NMR (C-13, Si-29) were applied to study the effects of mixed ammonia-water vapor at 90 degreesC on the mesostructure of the films. An increased cross-linking of the silica network was observed. Subsequent calcination of the silica films was seen to cause a bimodal pore-size distribution, with an accompanying increase in the volume and surface area ratios of the primary (d = 3 nm) to secondary (d = 5-30 nm) pores. Additionally, mixed ammonia-water treatment was observed to cause a narrowing of the primary pore-size distribution. These findings have implications for thin film based applications and devices, such as sensors, membranes, or surfaces for heterogeneous catalysis.

Modeling the molecular weight distribution of block copolymer formation in a reversible addition-fragmentation chain transfer mediated living radical polymerization

Block copolymers have become an integral part of the preparation of complex architectures through self-assembly. The use of reversible addition-fragmentation chain transfer (RAFT) allows blocks ranging from functional to nonfunctional polymers to be made with predictable molecular weight distributions. This article models block formation by varying many of the kinetic parameters. The simulations provide insight into the overall polydispersities (PDIs) that will be obtained when the chain-transfer constants in the main equilibrium steps are varied from 100 to 0.5. When the first dormant block [polymer-S-C(Z)=S] has a PDI of 1 and the second propagating radical has a low reactivity to the RAFT moiety, the overall PDI will be greater than 1 and dependent on the weight fraction of each block. When the first block has a PDI of 2 and the second propagating radical has a low reactivity to the RAFT moiety, the PDI will decrease to around 1.5 because of random coupling of two broad distributions. It is also shown how we can in principle use only one RAFT agent to obtain block copolymers with any desired molecular weight distribution. We can accomplish this by maintaining the monomer concentration at a constant level in the reactor over the course of the reaction. (c) 2005 Wiley Periodicals, Inc.

Adsorption of N-2, Ar, acetone, chloroform and acetone-chloroform mixture on carbon black in the framework of molecular layer structure theory (MLST)

Adsorption of pure nitrogen, argon, acetone, chloroform and acetone-chloroform mixture on graphitized thermal carbon black is considered at sub-critical conditions by means of molecular layer structure theory (MLST). In the present version of the MLST an adsorbed fluid is considered as a sequence of 2D molecular layers, whose Helmholtz free energies are obtained directly from the analysis of experimental adsorption isotherm of pure components. The interaction of the nearest layers is accounted for in the framework of mean field approximation. This approach allows quantitative correlating of experimental nitrogen and argon adsorption isotherm both in the monolayer region and in the range of multi-layer coverage up to 10 molecular layers. In the case of acetone and chloroform the approach also leads to excellent quantitative correlation of adsorption isotherms, while molecular approaches such as the non-local density functional theory (NLDFT) fail to describe those isotherms. We extend our new method to calculate the Helmholtz free energy of an adsorbed mixture using a simple mixing rule, and this allows us to predict mixture adsorption isotherms from pure component adsorption isotherms. The approach, which accounts for the difference in composition in different molecular layers, is tested against the experimental data of acetone-chloroform mixture (non-ideal mixture) adsorption on graphitized thermal carbon black at 50 degrees C. (C) 2005 Elsevier Ltd. All rights reserved.

Design strategies for controlling the molecular weight and rate using reversible addition-fragmentation chain transfer mediated living radical polymerization

Living radical polymerization has allowed complex polymer architectures to be synthesized in bulk, solution, and water. The most versatile of these techniques is reversible addition-fragmentation chain transfer (RAFT), which allows a wide range of functional and nonfunctional polymers to be made with predictable molecular weight distributions (MWDs), ranging from very narrow to quite broad. The great complexity of the RAFT mechanism and how the kinetic parameters affect the rate of polymerization and MWD are not obvious. Therefore, the aim of this article is to provide useful insights into the important kinetic parameters that control the rate of polymerization and the evolution of the MWD with conversion. We discuss how a change in the chain-transfer constant can affect the evolution of the MWD. It is shown how we can, in principle, use only one RAFT agent to obtain a poly-mer with any MWD. Retardation and inhibition are discussed in terms of (1) the leaving R group reactivity and (2) the intermediate radical termination model versus the slow fragmentation model. (c) 2005 Wiley Periodicals, Inc.

Steepness expansion for free surface flows in coastal aquifers

Free surface flow of groundwater in aquifers has been studied since the early 1960s. Previous investigations have been based on the Boussinesq equation, derived from the non-linear kinematic boundary condition. In fact, the Boussinesq equation is the zeroth-order equation in the shallow-water expansion. A key assumption in this expansion is that the mean thickness of the aquifer is small compared with a reference length, normally taken to be the linear decay length. In this study, we re-examine the expansion scheme for free surface groundwater flows, and propose a new expansion wherein the shallow-water assumption is replaced by a steepness assumption. A comparison with experimental data shows that the new model provides a better prediction of water table levels than the conventional shallow-water expansion. The applicable ranges of the two expansions are exhibited. (c) 2004 Elsevier B.V. All rights reserved.

Macrophages overexpressing tartrate-resistant acid phosphatase show altered profile of free radical production and enhanced capacity of bacterial killing

Activated macrophages and osteoclasts express high amounts of tartrate-resistant acid phosphatase (TRACP, acp5). TRACP has a binuclear iron center with a redox-active iron that has been shown to catalyze the formation of reactive oxygen species (ROS) by Fenton's reaction. Previous Studies Suggest that ROS generated by TRACP may participate in degradation of endocytosed bone matrix products in resorbing osteoclasts and degradation of foreign Compounds during. antigen presentation in activated macrophages. Here we have compared free radical production in macrophages of TRACP overexpressing (TRACP +) and wild-type (WT) mice. TRACP overexpression increased both ROS levels and Superoxide production. Nitric oxide production was increased in activated macrophages or WT mice, but not in TRACP+ mice, Macrophages from TRACP+ mice showed increased capacity or bacterial killing. Recombinant TRACP enzyme was capable of bacterial killing in the presence of hydrogen peroxide. These results suggest that TRACP has an important biological function in immune defense systern.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Mimicking the action of GroEL in molecular dynamics simulations: Application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.