Gene expression in large pedigrees: analytic approaches.

BackgroundWe currently have the ability to quantify transcript abundance of messenger RNA (mRNA), genome-wide, using microarray technologies. Analyzing genotype, phenotype and expression data from 20 pedigrees, the members of our Genetic Analysis Workshop (GAW) 19 gene expression group published 9 papers, tackling some timely and important problems and questions. To study the complexity and interrelationships of genetics and gene expression, we used established statistical tools, developed newer statistical tools, and developed and applied extensions to these tools.MethodsTo study gene expression correlations in the pedigree members (without incorporating genotype or trait data into the analysis), 2 papers used principal components analysis, weighted gene coexpression network analysis, meta-analyses, gene enrichment analyses, and linear mixed models. To explore the relationship between genetics and gene expression, 2 papers studied expression quantitative trait locus allelic heterogeneity through conditional association analyses, and epistasis through interaction analyses. A third paper assessed the feasibility of applying allele-specific binding to filter potential regulatory single-nucleotide polymorphisms (SNPs). Analytic approaches included linear mixed models based on measured genotypes in pedigrees, permutation tests, and covariance kernels. To incorporate both genotype and phenotype data with gene expression, 4 groups employed linear mixed models, nonparametric weighted U statistics, structural equation modeling, Bayesian unified frameworks, and multiple regression.Results and discussionRegarding the analysis of pedigree data, we found that gene expression is familial, indicating that at least 1 factor for pedigree membership or multiple factors for the degree of relationship should be included in analyses, and we developed a method to adjust for familiality prior to conducting weighted co-expression gene network analysis. For SNP association and conditional analyses, we found FaST-LMM (Factored Spectrally Transformed Linear Mixed Model) and SOLAR-MGA (Sequential Oligogenic Linkage Analysis Routines -Major Gene Analysis) have similar type 1 and type 2 errors and can be used almost interchangeably. To improve the power and precision of association tests, prior knowledge of DNase-I hypersensitivity sites or other relevant biological annotations can be incorporated into the analyses. On a biological level, eQTL (expression quantitative trait loci) are genetically complex, exhibiting both allelic heterogeneity and epistasis. Including both genotype and phenotype data together with measurements of gene expression was found to be generally advantageous in terms of generating improved levels of significance and in providing more interpretable biological models.ConclusionsPedigrees can be used to conduct analyses of and enhance gene expression studies

A new geometric invariant on initial data for Einstein equations

For a given asymptotically flat initial data set for Einstein equations a new geometric invariant is constructed. This invariant measure the departure of the data set from the stationary regime, it vanishes if and only if the data is stationary. In vacuum, it can be interpreted as a measure of the total amount of radiation contained in the data.Comment: 5 pages. Important corrections regarding the generalization to the non-time symmetric cas

Enhanced Integrin α4β1-Mediated Adhesion Contributes to a Mobilization Defect of Endothelial Progenitor Cells in Diabetes.

Diabetes is associated with a deficit of circulating endothelial progenitor cells (EPCs), which has been attributed to their defective mobilization from the bone marrow. The basis for this mobilization defect is not completely understood, and we sought to determine if hyperglycemic conditions enhanced EPC adhesion. We found that culturing EPCs in high glucose media increased adhesion to bone marrow stromal cells. This enhanced adhesion was associated with decreased expression of protein kinase A regulatory subunit 1β (PRKAR1β), activation of protein kinase A (PKA), and phosphorylation of α4-integrin on serine 988. This potentiated adhesion was reversed by treatment with a PKA inhibitor, overexpression of PRKAR1β, or expression of a phosphorylation-defective α4-integrin variant (α4[S988A]). Using a model of type 1 diabetes, we showed that α4(S988A)-expressing mice have more circulating EPCs than their wild-type counterparts. Moreover, diabetic α4(S988A) mice demonstrate enhanced revascularization after hind limb ischemia. Thus, we have identified a novel signaling mechanism activating PKA in diabetes (downregulation of an inhibitory regulatory subunit) that leads to deficits of circulating EPCs and impaired vascular repair, which could be reversed by α4-integrin mutation

Stretching Semiflexible Polymer Chains: Evidence for the Importance of Excluded Volume Effects from Monte Carlo Simulation

Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled by self-avoiding walks on the simple cubic lattice ($d=3$ dimensions) and square lattice ($d=2$ dimensions), varying chain stiffness by an energy penalty $\epsilon_b$ for chain bending. In the absence of excluded volume interactions, the persistence length $\ell_p$ of the polymers would then simply be $\ell_p=\ell_b(2d-2)^{-1}q_b^{-1}$ with $q_b= \exp(-\epsilon_b/k_BT)$, the bond length $\ell_b$ being the lattice spacing, and $k_BT$ is the thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method (PERM), both $q_b$ and the chain length $N$ are varied over a wide range $(0.005 \leq q_b \leq 1, \; N \leq 50000$), and also a stretching force $f$ is applied to one chain end (fixing the other end at the origin). In the absence of this force, in $d=2$ a single crossover from rod-like behavior (for contour lengths less than $\ell_p$) to swollen coils occurs, invalidating the Kratky-Porod model, while in $d=3$ a double crossover occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that are swollen due to the excluded volume interaction. If the stretching force is applied, excluded volume interactions matter for the force versus extension relation irrespective of chain stiffness in $d=2$, while theories based on the Kratky-Porod model are found to work in $d=3$ for stiff chains in an intermediate regime of chain extensions. While for $q_b \ll 1$ in this model a persistence length can be estimated from the initial decay of bond-orientational correlations, it is argued that this is not possible for more complex wormlike chains (e.g. bottle-brush polymers). Consequences for the proper interpretation of experiments are briefly discussed.Comment: 23 pages, 17 figures, 2 tables, to be published in J. Chem. Phys. (2011

The Effects of Stacking on the Configurations and Elasticity of Single Stranded Nucleic Acids

Stacking interactions in single stranded nucleic acids give rise to configurations of an annealed rod-coil multiblock copolymer. Theoretical analysis identifies the resulting signatures for long homopolynucleotides: A non monotonous dependence of size on temperature, corresponding effects on cyclization and a plateau in the extension force law. Explicit numerical results for poly(dA) and poly(rU) are presented.Comment: 4 pages and 2 figures. Accepted in Phys. Rev. E Rapid Com

On the existence of initial data containing isolated black holes

We present a general construction of initial data for Einstein's equations containing an arbitrary number of black holes, each of which is instantaneously in equilibrium. Each black hole is taken to be a marginally trapped surface and plays the role of the inner boundary of the Cauchy surface. The black hole is taken to be instantaneously isolated if its outgoing null rays are shear-free. Starting from the choice of a conformal metric and the freely specifiable part of the extrinsic curvature in the bulk, we give a prescription for choosing the shape of the inner boundaries and the boundary conditions that must be imposed there. We show rigorously that with these choices, the resulting non-linear elliptic system always admits solutions.Comment: 11 pages, 2 figures, RevTeX

A variational principle for stationary, axisymmetric solutions of Einstein's equations

Stationary, axisymmetric, vacuum, solutions of Einstein's equations are obtained as critical points of the total mass among all axisymmetric and $(t,\phi)$ symmetric initial data with fixed angular momentum. In this variational principle the mass is written as a positive definite integral over a spacelike hypersurface. It is also proved that if absolute minimum exists then it is equal to the absolute minimum of the mass among all maximal, axisymmetric, vacuum, initial data with fixed angular momentum. Arguments are given to support the conjecture that this minimum exists and is the extreme Kerr initial data.Comment: 21 page

Uniqueness and Non-uniqueness in the Einstein Constraints

The conformal thin sandwich (CTS) equations are a set of four of the Einstein equations, which generalize the Laplace-Poisson equation of Newton's theory. We examine numerically solutions of the CTS equations describing perturbed Minkowski space, and find only one solution. However, we find {\em two} distinct solutions, one even containing a black hole, when the lapse is determined by a fifth elliptic equation through specification of the mean curvature. While the relationship of the two systems and their solutions is a fundamental property of general relativity, this fairly simple example of an elliptic system with non-unique solutions is also of broader interest.Comment: 4 pages, 4 figures; abstract and introduction rewritte

3+1 Approach to the Long Wavelength Iteration Scheme

Large-scale inhomogeneities and anisotropies are modeled using the Long Wavelength Iteration Scheme. In this scheme solutions are obtained as expansions in spatial gradients, which are taken to be small. It is shown that the choice of foliation for spacetime can make the iteration scheme more effective in two respects: (i) the shift vector can be chosen so as to dilute the effect of anisotropy on the late-time value of the extrinsic curvature of the spacelike hypersurfaces of the foliation; and (ii) pure gauge solutions present in a similar calculation using the synchronous gauge vanish when the spacelike hypersurfaces have extrinsic curvature with constant trace. We furthermore verify the main conclusion of the synchronous gauge calculation which is large-scale inhomogeneity decays if the matter--considered to be that of a perfect-fluid with a barotropic equation of state--violates the strong-energy condition. Finally, we obtain the solution for the lapse function and discuss its late-time behaviour. It is found that the lapse function is well-behaved when the matter violates the strong energy condition.Comment: 21 pages, TeX file, already publishe