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.

Gene expression array testing of FFPE archival breast tumor samples: an optimized protocol for WG-DASL sample preparation.

Archived formalin-fixed, paraffin embedded (FFPE) tissues constitute a vast, well-annotated, but underexploited resource for the molecular study of cancer progression, largely because degradation, chemical modification, and cross-linking, render FFPE RNA a suboptimal substrate for conventional analytical methods. We report here a modified protocol for RNA extraction from FFPE tissues which maximized the success rate (with 100% of samples) in the expression profiling of a set of 60 breast cancer samples on the WG-DASL platform; yielding data of sufficient quality such that in hierarchical clustering (a) 12/12 (100%) replicates correctly identified their respective counterparts, with a high self-correlation (r = 0.979), and (b) the overall sample set grouped with high specificity into ER+ (38/40; 95%) and ER- (18/20; 90%) subtypes. These results indicate that a large fraction of decade-old FFPE samples, of diverse institutional origins and processing histories, can yield RNA suitable for gene expression profiling experiments.

Quality assessment and data analysis for microRNA expression arrays.

MicroRNAs are small (approximately 22 nt) RNAs that regulate gene expression and play important roles in both normal and disease physiology. The use of microarrays for global characterization of microRNA expression is becoming increasingly popular and has the potential to be a widely used and valuable research tool. However, microarray profiling of microRNA expression raises a number of data analytic challenges that must be addressed in order to obtain reliable results. We introduce here a universal reference microRNA reagent set as well as a series of nonhuman spiked-in synthetic microRNA controls, and demonstrate their use for quality control and between-array normalization of microRNA expression data. We also introduce diagnostic plots designed to assess and compare various normalization methods. We anticipate that the reagents and analytic approach presented here will be useful for improving the reliability of microRNA microarray experiments.

Interferon regulatory factor-1 regulates reconstituted extracellular matrix (rECM)-mediated apoptosis in human mammary epithelial cells.

Interactions between extracellular matrix (ECM) and mammary epithelial cells are critical for mammary gland homeostasis and apoptotic signaling. Interferon regulatory factor-1 (IRF-1) is a transcriptional regulator that promotes apoptosis during mammary gland involution and p53-independent apoptosis. We have recently shown that rapid cell surface tamoxifen (Tam) signaling promotes apoptosis in normal human mammary epithelial cells that were acutely damaged by expression of human papillomavirus type-16 E6 protein (*HMEC-E6). Apoptosis was mediated by recruitment of CREB-binding protein (CBP) to the gamma-activating sequence (GAS) element of the IRF-1 promoter, induction of IRF-1 and caspase-1/-3 activation. Here, we show that growth factor-depleted, reconstituted ECM (rECM), similar to Tam, promotes apoptosis in *HMEC-E6 cells through induction of IRF-1. Apoptosis was temporally associated with recruitment of CBP to the GAS element of the IRF-1 promoter, induction of IRF-1 expression and caspase-1/-3 activation. Small interfering RNA-mediated suppression of IRF-1 protein expression in *HMEC-E6 cells blocked (1) induction of IRF-1, (2) caspase-1/-3 activation and (3) apoptosis. These observations demonstrate that IRF-1 promotes rECM-mediated apoptosis and provide evidence that both rECM and rapid Tam signaling transcriptionally activate IRF-1 through recruitment of CBP to the IRF-1 GAS promoter complex.

Genomic instability during Myc-induced lymphomagenesis in the bursa of Fabricius.

Retroviral vector-mediated overexpression of c-myc in embryonic bursal precursors induces multi-staged tumorigenesis beginning with preneoplastic-transformed follicles (TF) and progressing to clonal metastatic B-cell lymphomas. Using a 13K chicken cDNA microarray, specifically enriched for chicken immune system expressed sequence tagged (ESTs), we carried out array-based comparative genomic hybridization (array-CGH) and detected significant DNA copy number change at many loci on most or all chromosomes in both early TF and end-stage lymphomas. Formation of long palindromes, through breakage-fusion-bridge cycles, is thought to play an early role in gene amplification. Employing genome-wide analysis of palindrome formation (GAPF), we detected extensive palindrome formation in early TF and end-stage lymphomas. The population of loci showing amplification by array-CGH was enriched for palindromes detected by GAPF providing strong evidence for genetic instability early in Myc-induced tumorigenesis and further support for the role of palindromes in gene amplification. Comparing gene copy number change and RNA expression changes profiled on the same cDNA array, we detected very little consistent contribution of gene copy number change to RNA expression changes. Palindromic loci in TF and tumors, however, were expressed, many at high levels, suggesting an abundance of RNA species with long double-stranded segments generated during tumorigenesis.

Project normal: defining normal variance in mouse gene expression.

The mouse has become an indispensable and versatile model organism for the study of development, genetics, behavior, and disease. The application of comprehensive gene expression profiling technologies to compare normal and diseased tissues or to assess molecular alterations resulting from various experimental interventions has the potential to provide highly detailed qualitative and quantitative descriptions of these processes. Ideally, to interpret experimental data, the magnitude and diversity of gene expression for the system under study should be well characterized, yet little is known about the normal variation of mouse gene expression in vivo. To assess natural differences in murine gene expression, we used a 5406-clone spotted cDNA microarray to quantitate transcript levels in the kidney, liver, and testis from each of 6 normal male C57BL6 mice. We used ANOVA to compare the variance across the six mice to the variance among four replicate experiments performed for each mouse tissue. For the 6 kidney samples, 102 of 3,088 genes (3.3%) exhibited a statistically significant mouse variance at a level of 0.05. In the testis, 62 of 3,252 genes (1.9%) showed statistically significant variance, and in the liver, there were 21 of 2,514 (0.8%) genes with significantly variable expression. Immune-modulated, stress-induced, and hormonally regulated genes were highly represented among the transcripts that were most variable. The expression levels of several genes varied significantly in more than one tissue. These studies help to define the baseline level of variability in mouse gene expression and emphasize the importance of replicate microarray experiments.

Widespread collaboration of Isw2 and Sin3-Rpd3 chromatin remodeling complexes in transcriptional repression.

The yeast Isw2 chromatin remodeling complex functions in parallel with the Sin3-Rpd3 histone deacetylase complex to repress early meiotic genes upon recruitment by Ume6p. For many of these genes, the effect of an isw2 mutation is partially masked by a functional Sin3-Rpd3 complex. To identify the full range of genes repressed or activated by these factors and uncover hidden targets of Isw2-dependent regulation, we performed full genome expression analyses using cDNA microarrays. We find that the Isw2 complex functions mainly in repression of transcription in a parallel pathway with the Sin3-Rpd3 complex. In addition to Ume6 target genes, we find that many Ume6-independent genes are derepressed in mutants lacking functional Isw2 and Sin3-Rpd3 complexes. Conversely, we find that ume6 mutants, but not isw2 sin3 or isw2 rpd3 double mutants, have reduced fidelity of mitotic chromosome segregation, suggesting that one or more functions of Ume6p are independent of Sin3-Rpd3 and Isw2 complexes. Chromatin structure analyses of two nonmeiotic genes reveals increased DNase I sensitivity within their regulatory regions in an isw2 mutant, as seen previously for one meiotic locus. These data suggest that the Isw2 complex functions at Ume6-dependent and -independent loci to create DNase I-inaccessible chromatin structure by regulating the positioning or placement of nucleosomes.

Analysis of gene expression during myc oncogene-induced lymphomagenesis in the bursa of Fabricius.

The transcriptional effects of deregulated myc gene overexpression are implicated in tumorigenesis in a spectrum of experimental and naturally occurring neoplasms. In follicles of the chicken bursa of Fabricius, myc induction of B-cell neoplasia requires a target cell population present during early bursal development and progresses through preneoplastic transformed follicles to metastatic lymphomas. We developed a chicken immune system cDNA microarray to analyze broad changes in gene expression that occur during normal embryonic B-cell development and during myc-induced neoplastic transformation in the bursa. The number of mRNAs showing at least 3-fold change was greater during myc-induced lymphomagenesis than during normal development, and hierarchical cluster analysis of expression patterns revealed that levels of several hundred mRNAs varied in concert with levels of myc overexpression. A set of 41 mRNAs were most consistently elevated in myc-overexpressing preneoplastic and neoplastic cells, most involved in processes thought to be subject to regulation by Myc. The mRNAs for another cluster of genes were overexpressed in neoplasia independent of myc expression level, including a small subset with the expression signature of embryonic bursal lymphocytes. Overexpression of myc, and some of the genes overexpressed with myc, may be important for generation of preneoplastic transformed follicles. However, expression profiles of late metastatic tumors showed a large variation in concert with myc expression levels, and some showed minimal myc overexpression. Therefore, high-level myc overexpression may be more important in the early induction of these lymphomas than in maintenance of late-stage metastases.

Chromatin profiling using targeted DNA adenine methyltransferase.

Chromatin is the highly complex structure consisting of DNA and hundreds of associated proteins. Most chromatin proteins exert their regulatory and structural functions by binding to specific chromosomal loci. Knowledge of the identity of these in vivo target loci is essential for the understanding of the functions and mechanisms of action of chromatin proteins. We report here large-scale mapping of in vivo binding sites of chromatin proteins, using a novel approach based on a combination of targeted DNA methylation and microarray technology. We show that three distinct chromatin proteins in Drosophila melanogaster cells each associate with specific sets of genes. HP1 binds predominantly to pericentric genes and transposable elements. GAGA factor associates with euchromatic genes that are enriched in (GA)n motifs. A Drosophila homolog of Saccharomyces cerevisiae Sir2p is associated with several active genes and is excluded from heterochromatin. High-resolution, genome-wide maps of target loci of chromatin proteins ('chromatin profiles') provide new insights into chromatin structure and gene regulation.

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.