Transcriptional Signatures of Sleep Duration Discordance in Monozygotic Twins.

Introduction: Habitual short sleep duration is associated with adverse metabolic, cardiovascular, and inflammatory effects. Co-twin study methodologies account for familial (eg, genetics and shared environmental) confounding, allowing assessment of subtle environmental effects, such as the effect of habitual short sleep duration on gene expression. Therefore, we investigated gene expression in monozygotic twins discordant for actigraphically phenotyped habitual sleep duration. Methods: Eleven healthy monozygotic twin pairs (82% female; mean age 42.7 years; SD = 18.1), selected based on subjective sleep duration discordance, were objectively phenotyped for habitual sleep duration with 2 weeks of wrist actigraphy. Peripheral blood leukocyte (PBL) RNA from fasting blood samples was obtained on the final day of actigraphic measurement and hybridized to Illumina humanHT-12 microarrays. Differential gene expression was determined between paired samples and mapped to functional categories using Gene Ontology. Finally, a more comprehensive gene set enrichment analysis was performed based on the entire PBL transcriptome. Results: The mean 24-hour sleep duration of the total sample was 439.2 minutes (SD = 46.8 minutes; range 325.4-521.6 minutes). Mean within-pair sleep duration difference per 24 hours was 64.4 minutes (SD = 21.2; range 45.9-114.6 minutes). The twin cohort displayed distinctive pathway enrichment based on sleep duration differences. Habitual short sleep was associated with up-regulation of genes involved in transcription, ribosome, translation, and oxidative phosphorylation. Unexpectedly, genes down-regulated in short sleep twins were highly enriched in immuno-inflammatory pathways such as interleukin signaling and leukocyte activation, as well as developmental programs, coagulation cascade, and cell adhesion. Conclusions: Objectively assessed habitual sleep duration in monozygotic twin pairs appears to be associated with distinct patterns of differential gene expression and pathway enrichment. By accounting for familial confounding and measuring real life sleep duration, our study shows the transcriptomic effects of habitual short sleep on dysregulated immune response and provides a potential link between sleep deprivation and adverse metabolic, cardiovascular, and inflammatory outcomes.

Effect of temperature on DNA secondary structure in the absence and presence of 0.5 M tetramethylammonium chloride.

Changes in the average secondary structures of three different linear DNAs over the premelting region from 5 to 60 degrees C were investigated by measuring their CD spectra and also their torsion elastic constants () by time-resolved fluorescence polarization anisotropy. For one of these DNAs, the Haell fragment of pBR322, the apparent diffusion coefficients [Dapp(k)] at small and large scattering vectors (k) were also measured by dynamic light scattering. With increasing temperature, all three DNAs exhibited typical premelting changes in their CD spectra, and these were accompanied by 1.4- to 1.7-fold decreases in . Also for the 1876 base pair fragment, Dapp(k) at large scattering vectors, which is sensitive to the dynamic bending rigidity, decreased by 17%, even though there was no change at small scattering vectors, where Dapp(k) = D0 is the translational diffusion coefficient of the center-of-mass. These observations demonstrate conclusively that the premelting CD changes of these DNAs are associated with a significant change in average secondary structure and mechanical properties, though not in persistence length. In the presence of 0.5 M tetramethylammonium chloride (TMA-Cl) the premelting change in CD is largely suppressed, and the corresponding changes in and Dapp(k) at large scattering vectors are substantially diminished. These observations suggest that TMA-Cl, which binds preferentially to A.T-rich regions and stabilizes those regions (relative to G.C-rich regions) against melting, effectively stabilizes the prevailing low-temperature secondary structure sufficiently that the DNA is effectively trapped in that state over the temperature range observed.

On the origin of the temperature dependence of the supercoiling free energy.

Monte Carlo simulations using temperature-invariant torsional and bending rigidities fail to predict the rather steep decline of the experimental supercoiling free energy with increasing temperature, and consequently fail to predict the correct sign and magnitude of the supercoiling entropy. To illustrate this problem, values of the twist energy parameter (E(T)), which governs the supercoiling free energy, were simulated using temperature-invariant torsion and bending potentials and compared to experimental data on pBR322 over a range of temperatures. The slope, -dE(T)/dT, of the simulated values is also compared to the slope derived from previous calorimetric data. The possibility that the discrepancies arise from some hitherto undetected temperature dependence of the torsional rigidity was investigated. The torsion elastic constant of an 1876-bp restriction fragment of pBR322 was measured by time-resolved fluorescence polarization anisotropy of intercalated ethidium over the range 278-323 K, and found to decline substantially over that interval. Simulations of a 4349-bp model DNA were performed using these measured temperature-dependent torsional rigidities. The slope, -dE(T)/dT, of the simulated data agrees satisfactorily with the slope derived from previous calorimetric measurements, but still lies substantially below that of Duguet's data. Models that involve an equilibrium between different secondary structure states with different intrinsic twists and torsion constants provide the most likely explanation for the variation of the torsion constant with T and other pertinent observations.

Comparison of hard-cylinder and screened Coulomb interactions in the modeling of supercoiled DNAs.

A 1000 base pair (bp) model supercoiled DNA is simulated using spherical screened Coulomb interactions between subunits on one hand and equivalent hard-cylinder interactions on the other. The amplitudes, or effective charges, of the spherical screened Coulomb electrostatic potentials are chosen so that the electrostatic potential surrounding the middle of a linear array of 2001 subunits (31.8 A diameter) closely matches the solution of the nonlinear Poisson-Boltzmann equation for a cylinder with 12 A radius and the full linear charge density of DNA at all distances beyond the 24 A hard-core diameter. This superposition of spherical screened Coulomb potentials is practically identical to the particular solution of the cylindrical linearized Poisson-Boltzmann equation that matches the solution of the nonlinear Poisson-Boltzmann equation at large distances. The interaction energy between subunits is reckoned from the effective charges according to the standard DLVO expression. The equivalent hard-cylinder diameter is chosen following Stigter's protocol for matching second virial coefficients, but for the full linear charge density of DNA. The electrostatic persistence length of the model with screened Coulomb interactions is extremely sensitive to the (arbitrarily) chosen subunit length at the higher salt concentrations. The persistence length of the hard-cylinder model is adjusted to match that of the screened Coulomb model for each ionic condition. Simulations for a superhelix density sigma = -0.05 using a spherical screened Coulomb interaction plus a 24 A hard-cylinder core (SCPHC) potential indicate that the radius of gyration of this 1000 bp DNA actually undergoes a slight increase as the NaCl concentration is raised from 0.01 to 1.0M. Thus, merely softening the potential from hard-cylinder to screened Coulomb form does not produce a large decrease in radius of gyration with increasing NaCl concentration for DNAs of this size. Radii of gyration, static structure factors, and diffusion coefficients obtained using the equivalent hard-cylinder (EHC) potential agree well with those obtained using the SCPHC potential in 1.0M NaCl, but in 0.1M NaCl the agreement is not as good, and in 0.01M NaCl the agreement is definitely unsatisfactory. These conclusions differ in significant respects from those obtained in previous studies.

The question of long-range allosteric transitions in DNA.

The question of long-range allosteric transitions of DNA secondary structure and their possible involvement in transcriptional activation is discussed in the light of new results. A variety of recent evidence strongly supports a fluctuating long-range description of DNA secondary structure. Balanced equilibria between two or more different secondary structures, and the occurrence of very large domain sizes, have been documented in several instances. Long-range allosteric effects stemming from changes in sequence or secondary structure over a small region of the DNA have been observed to extend over distances up to hundreds of base pairs in some cases. The discovery that coherent bending strain beyond a threshold level in small (N < or = 250 base pairs (bp)] circular DNAs significantly alters the DNA secondary structure has important implications, especially for transcriptional activators that either bend the DNA directly or are involved in the formation of DNA loops of sufficiently small size (N < or = 250 bp). Whether the RNA polymerase is activated primarily via protein: protein contacts, as is widely believed, or instead via a bend-induced allosteric transition of the DNA in such a small loop, is now an open question. Binding of the transcriptional activator Sp1 to linear DNA induces a remarkably long-range change in its secondary structure, and catabolite activator protein binding to a supercoiled DNA behaves similarly, though possibly for different reasons. Compelling evidence for a bend-induced long-range structural transmission effect of the transcriptional activator integration host factor on RNA polymerase activity was recently reported. These results may augur a new paradigm in which allosteric transitions of duplex DNA, as well as of the proteins, are involved in the regulation of transcription.

Effects of Na+ and Mg2+ on the structures of supercoiled DNAs: comparison of simulations with experiments.

Recent cryo-electron microscopy (cryo-EM) results suggest that sufficient NaCl concentration (> or approximately 0.1 M) and superhelix density (> or approximately-0.05) cause circular DNAs to adopt highly extended, tightly interwound configurations, in which the strands are laterally contiguous along almost their entire length. Millimolar levels of MgCl2 reportedly act synergistically with NaCl to produce similar conformations. However, Monte Carlo simulations with purely repulsive interduplex forces failed to reproduce such structures. In the present work, solution measurements of particular physical properties were performed both to characterize the effects of Na+ and Mg2+ on DNA structure and to provide quantitative tests of Monte Carlo simulations of circular DNAs. Supercoiled p30 delta DNAs in 10 mM Tris plus 0, 0.122, and 0.1 M NaCl, and 0.1 M NaCl plus 4 mM Mg2+ were examined by static and dynamic light scattering (LS and DLS), time-resolved fluorescence polarization anisotropy (FPA) of intercalated ethidium, and circular dichroism (CD) spectroscopy. Upon addition of 0.122 M NaCl, the radius of gyration (Rg) decreased substantially, which indicates that p30 delta adopts a more compact structure. This contradicts the cryo-EM studies, where molecular extension and Rg both increase upon adding 0.1 M NaCl. In 0.1 M NaCl, the torsion constant measured by FPA is practically invariant to superhelix density, and the plateau diffusion coefficient at large scattering vector (Dplat) is likewise nearly the same at both relaxed and native superhelix densities. Such invariance is difficult to reconcile with any transition from relaxed circles to tightly interwound structures with laterally contiguous strands. Metropolis Monte Carlo simulations were performed to generate canonically distributed sets of structures, from which average Do values and scattered intensity ratios, [symbol: see text]I (zero) [symbol: see text]/[symbol: see text] l(k) [symbol: see text], were calculated. Agreement between simulations and experiments in regard to [symbol: see text] I(O) [symbol: see text] /[symbol: see text] I(k) [symbol: see text], D(zero) and the supercoiling free energy, delta Gsc (delta l), is remarkably good for the most extensively studied p30 delta samples. The simulated structures exhibit no sign of very tight interwinding with extensive lateral contacts, but instead exhibit most probable superhelix diameters of 85 to 90 A. When 4 mM Mg2+ was added to native supercoiled p30 delta in 0.1 M NaCl, Rg decreased, D(zero) increased, and the longest internal relaxation rate (1/tau 2(zero)) increased, all of which indicate a further overall contraction of the molecular envelope. The torsion constant exhibited a slight increase that is hardly statistically significant. In this case, agreement between the simulations and experiments was only semi-quantitative for most samples investigated, although the predicted contraction was exhibited by all five samples of p30 delta and one of pBR322 DNA. The simulated structures in 0.1 M NaCl plus 4 mM Mg2+ again showed no sign of extensive lateral contacts. A plausible explanation is proposed for the highly extended, tightly interwound structures seen in cryo-EM, and explicitly tested by Monte Carlo simulations of a 1000 bp circular DNA at +25 and -50 degrees C. Structures identical to those seen in cryo-EM are in fact the equilibrium structures in the simulations at -50 degrees C, and the estimated time for equilibration (2.3 x 10(-6) second) is much smaller than the estimated time for vitrification (1 x 10(-4) second).

Rapid analysis of mitotic histone H1 phosphorylation by cationic disc electrophoresis at neutral pH in minigels.

A new method is described for analysis of histone H1 and other basic proteins by cationic disc electrophoresis in polyacrylamide gels at neutral pH. The multiphasic buffer (disc) system uses Na+ as leading ion, L-histidine as trailing ion, and Hepes as buffering counterion. These "Hepes/histidine gels" have three advantages over conventional acid-urea gels for studies of H1 phosphorylation and dephosphorylation: speed, convenience, and the need for only small amounts of cells or chromatin. Core histones and their acetylated forms can also be separated in gels containing 0.4% Triton X-100. The difference in electrophoretic mobility between mitotic (superphosphorylated) and interphase H1 from HeLa cells is approximately twice as great at neutral pH as at pH 4.5, making it possible to separate these two H1 forms rapidly and easily in Hepes/histidine "minigels" only 5-cm long. Total histones can be rapidly prepared by simply neutralizing 0.2 N HCl extracts, and the entire analysis, from harvesting cells to destaining gels, can be carried out in 1 day. The stacking effect of the disc system produces sharp bands and high resolution even with relatively dilute samples.