Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

The induced innovation hypothesis and energy-saving technological change

We develop a methodology for testing Hicks's induced innovation hypothesis by estimating a product-characteristics model of energy-using consumer durables, augmenting the hypothesis to allow for the influence of government regulations. For the products we explored, the evidence suggests that (i) the rate of overall innovation was independent of energy prices and regulations; (ii) the direction of innovation was responsive to energy price changes for some products but not for others; (iii) energy price changes induced changes in the subset of technically feasible models that were offered for sale; (iv) this responsiveness increased substantially during the period after energy-efficiency product labeling was required; and (v) nonetheless, a sizable portion of efficiency improvements were autonomous.

Influence of Outlet Boundary Conditions on Cerebrovascular Aneurysm Hemodynamics

Computational fluid dynamic (CFD) studies of blood flow in cerebrovascular aneurysms have potential to improve patient treatment planning by enabling clinicians and engineers to model patient-specific geometries and compute predictors and risks prior to neurovascular intervention. However, the use of patient-specific computational models in clinical settings is unfeasible due to their complexity, computationally intensive and time-consuming nature. An important factor contributing to this challenge is the choice of outlet boundary conditions, which often involves a trade-off between physiological accuracy, patient-specificity, simplicity and speed. In this study, we analyze how resistance and impedance outlet boundary conditions affect blood flow velocities, wall shear stresses and pressure distributions in a patient-specific model of a cerebrovascular aneurysm. We also use geometrical manipulation techniques to obtain a model of the patient’s vasculature prior to aneurysm development, and study how forces and stresses may have been involved in the initiation of aneurysm growth. Our CFD results show that the nature of the prescribed outlet boundary conditions is not as important as the relative distributions of blood flow through each outlet branch. As long as the appropriate parameters are chosen to keep these flow distributions consistent with physiology, resistance boundary conditions, which are simpler, easier to use and more practical than their impedance counterparts, are sufficient to study aneurysm pathophysiology, since they predict very similar wall shear stresses, time-averaged wall shear stresses, time-averaged pressures, and blood flow patterns and velocities. The only situations where the use of impedance boundary conditions should be prioritized is if pressure waveforms are being analyzed, or if local pressure distributions are being evaluated at specific time points, especially at peak systole, where the use of resistance boundary conditions leads to unnaturally large pressure pulses. In addition, we show that in this specific patient, the region of the blood vessel where the neck of the aneurysm developed was subject to abnormally high wall shear stresses, and that regions surrounding blebs on the aneurysmal surface were subject to low, oscillatory wall shear stresses. Computational models using resistance outlet boundary conditions may be suitable to study patient-specific aneurysm progression in a clinical setting, although several other challenges must be addressed before these tools can be applied clinically.

Reference groups and product line decisions: An experimental investigation of limited editions and product proliferation

Some luxury goods manufacturers offer limited editions of their products, whereas some others market multiple product lines. Researchers have found that reference groups shape consumer evaluations of these product categories. Yet little empirical research has examined how reference groups affect the product line decisions of firms. Indeed, in a field setting it is quite a challenge to isolate reference group effects from contextual effects and correlated effects. In this paper, we propose a parsimonious model that allows us to study how reference groups influence firm behavior and that lends itself to experimental analysis. With the aid of the model we investigate the behavior of consumers in a laboratory setting where we can focus on the reference group effects after controlling for the contextual and correlated effects. The experimental results show that in the presence of strong reference group effects, limited editions and multiple products can help improve firms' profits. Furthermore, the trends in the purchase decisions of our participants point to the possibility that they are capable of introspecting close to two steps of thinking at the outset of the game and then learning through reinforcement mechanisms. © 2010 INFORMS.

Altered gene expression and DNA damage in peripheral blood cells from Friedreich's ataxia patients: cellular model of pathology.

The neurodegenerative disease Friedreich's ataxia (FRDA) is the most common autosomal-recessively inherited ataxia and is caused by a GAA triplet repeat expansion in the first intron of the frataxin gene. In this disease, transcription of frataxin, a mitochondrial protein involved in iron homeostasis, is impaired, resulting in a significant reduction in mRNA and protein levels. Global gene expression analysis was performed in peripheral blood samples from FRDA patients as compared to controls, which suggested altered expression patterns pertaining to genotoxic stress. We then confirmed the presence of genotoxic DNA damage by using a gene-specific quantitative PCR assay and discovered an increase in both mitochondrial and nuclear DNA damage in the blood of these patients (p<0.0001, respectively). Additionally, frataxin mRNA levels correlated with age of onset of disease and displayed unique sets of gene alterations involved in immune response, oxidative phosphorylation, and protein synthesis. Many of the key pathways observed by transcription profiling were downregulated, and we believe these data suggest that patients with prolonged frataxin deficiency undergo a systemic survival response to chronic genotoxic stress and consequent DNA damage detectable in blood. In conclusion, our results yield insight into the nature and progression of FRDA, as well as possible therapeutic approaches. Furthermore, the identification of potential biomarkers, including the DNA damage found in peripheral blood, may have predictive value in future clinical trials.

Selective enhancement of donor hematopoietic cell engraftment by the CXCR4 antagonist AMD3100 in a mouse transplantation model.

The interaction between stromal cell-derived factor-1 (SDF-1) with CXCR4 chemokine receptors plays an important role in hematopoiesis following hematopoietic stem cell transplantation. We examined the efficacy of post transplant administration of a specific CXCR4 antagonist (AMD3100) in improving animal survival and in enhancing donor hematopoietic cell engraftment using a congeneic mouse transplantation model. AMD3100 was administered subcutaneously at 5 mg/kg body weight 3 times a week beginning at day +2 post-transplant. Post-transplant administration of AMD3100 significantly improves animal survival. AMD3100 reduces pro-inflammatory cytokine/chemokine production. Furthermore, post transplant administration of AMD3100 selectively enhances donor cell engraftment and promotes recovery of all donor cell lineages (myeloid cells, T and B lymphocytes, erythrocytes and platelets). This enhancement results from a combined effect of increased marrow niche availability and greater cell division induced by AMD3100. Our studies shed new lights into the biological roles of SDF-1/CXCR4 interaction in hematopoietic stem cell engraftment following transplantation and in transplant-related mortality. Our results indicate that AMD3100 provides a novel approach for enhancing hematological recovery following transplantation, and will likely benefit patients undergoing transplantation.

Discharge competence and pattern formation in peatlands: a meta-ecosystem model of the Everglades ridge-slough landscape.

Regular landscape patterning arises from spatially-dependent feedbacks, and can undergo catastrophic loss in response to changing landscape drivers. The central Everglades (Florida, USA) historically exhibited regular, linear, flow-parallel orientation of high-elevation sawgrass ridges and low-elevation sloughs that has degraded due to hydrologic modification. In this study, we use a meta-ecosystem approach to model a mechanism for the establishment, persistence, and loss of this landscape. The discharge competence (or self-organizing canal) hypothesis assumes non-linear relationships between peat accretion and water depth, and describes flow-dependent feedbacks of microtopography on water depth. Closed-form model solutions demonstrate that 1) this mechanism can produce spontaneous divergence of local elevation; 2) divergent and homogenous states can exhibit global bi-stability; and 3) feedbacks that produce divergence act anisotropically. Thus, discharge competence and non-linear peat accretion dynamics may explain the establishment, persistence, and loss of landscape pattern, even in the absence of other spatial feedbacks. Our model provides specific, testable predictions that may allow discrimination between the self-organizing canal hypotheses and competing explanations. The potential for global bi-stability suggested by our model suggests that hydrologic restoration may not re-initiate spontaneous pattern establishment, particularly where distinct soil elevation modes have been lost. As a result, we recommend that management efforts should prioritize maintenance of historic hydroperiods in areas of conserved pattern over restoration of hydrologic regimes in degraded regions. This study illustrates the value of simple meta-ecosystem models for investigation of spatial processes.

Antisense Gene Silencing and Bacteriophages as Novel Disinfection Processes for Engineered Systems

The growth and proliferation of invasive bacteria in engineered systems is an ongoing problem. While there are a variety of physical and chemical processes to remove and inactivate bacterial pathogens, there are many situations in which these tools are no longer effective or appropriate for the treatment of a microbial target. For example, certain strains of bacteria are becoming resistant to commonly used disinfectants, such as chlorine and UV. Additionally, the overuse of antibiotics has contributed to the spread of antibiotic resistance, and there is concern that wastewater treatment processes are contributing to the spread of antibiotic resistant bacteria.

Due to the continually evolving nature of bacteria, it is difficult to develop methods for universal bacterial control in a wide range of engineered systems, as many of our treatment processes are static in nature. Still, invasive bacteria are present in many natural and engineered systems, where the application of broad acting disinfectants is impractical, because their use may inhibit the original desired bioprocesses. Therefore, to better control the growth of treatment resistant bacteria and to address limitations with the current disinfection processes, novel tools that are both specific and adaptable need to be developed and characterized.

In this dissertation, two possible biological disinfection processes were investigated for use in controlling invasive bacteria in engineered systems. First, antisense gene silencing, which is the specific use of oligonucleotides to silence gene expression, was investigated. This work was followed by the investigation of bacteriophages (phages), which are viruses that are specific to bacteria, in engineered systems.


For the antisense gene silencing work, a computational approach was used to quantify the number of off-targets and to determine the effects of off-targets in prokaryotic organisms. For the organisms of Escherichia coli K-12 MG1655 and Mycobacterium tuberculosis H37Rv the mean number of off-targets was found to be 15.0 + 13.2 and 38.2 + 61.4, respectively, which results in a reduction of greater than 90% of the effective oligonucleotide concentration. It was also demonstrated that there was a high variability in the number of off-targets over the length of a gene, but that on average, there was no general gene location that could be targeted to reduce off-targets. Therefore, this analysis needs to be performed for each gene in question. It was also demonstrated that the thermodynamic binding energy between the oligonucleotide and the mRNA accounted for 83% of the variation in the silencing efficiency, compared to the number of off-targets, which explained 43% of the variance of the silencing efficiency. This suggests that optimizing thermodynamic parameters must be prioritized over minimizing the number of off-targets. In conclusion for the antisense work, these results suggest that off-target hybrids can account for a greater than 90% reduction in the concentration of the silencing oligonucleotides, and that the effective concentration can be increased through the rational design of silencing targets by minimizing off-target hybrids.

Regarding the work with phages, the disinfection rates of bacteria in the presence of phages was determined. The disinfection rates of E. coli K12 MG1655 in the presence of coliphage Ec2 ranged up to 2 h-1, and were dependent on both the initial phage and bacterial concentrations. Increasing initial phage concentrations resulted in increasing disinfection rates, and generally, increasing initial bacterial concentrations resulted in increasing disinfection rates. However, disinfection rates were found to plateau at higher bacterial and phage concentrations. A multiple linear regression model was used to predict the disinfection rates as a function of the initial phage and bacterial concentrations, and this model was able to explain 93% of the variance in the disinfection rates. The disinfection rates were also modeled with a particle aggregation model. The results from these model simulations suggested that at lower phage and bacterial concentrations there are not enough collisions to support active disinfection rates, which therefore, limits the conditions and systems where phage based bacterial disinfection is possible. Additionally, the particle aggregation model over predicted the disinfection rates at higher phage and bacterial concentrations of 108 PFU/mL and 108 CFU/mL, suggesting other interactions were occurring at these higher concentrations. Overall, this work highlights the need for the development of alternative models to more accurately describe the dynamics of this system at a variety of phage and bacterial concentrations. Finally, the minimum required hydraulic residence time was calculated for a continuous stirred-tank reactor and a plug flow reactor (PFR) as a function of both the initial phage and bacterial concentrations, which suggested that phage treatment in a PFR is theoretically possible.

In addition to determining disinfection rates, the long-term bacterial growth inhibition potential was determined for a variety of phages with both Gram-negative and Gram-positive bacteria. It was determined, that on average, phages can be used to inhibit bacterial growth for up to 24 h, and that this effect was concentration dependent for various phages at specific time points. Additionally, it was found that a phage cocktail was no more effective at inhibiting bacterial growth over the long-term than the best performing phage in isolation.

Finally, for an industrial application, the use of phages to inhibit invasive Lactobacilli in ethanol fermentations was investigated. It was demonstrated that phage 8014-B2 can achieve a greater than 3-log inactivation of Lactobacillus plantarum during a 48 h fermentation. Additionally, it was shown that phages can be used to protect final product yields and maintain yeast viability. Through modeling the fermentation system with differential equations it was determined that there was a 10 h window in the beginning of the fermentation run, where the addition of phages can be used to protect final product yields, and after 20 h no additional benefit of the phage addition was observed.

In conclusion, this dissertation improved the current methods for designing antisense gene silencing targets for prokaryotic organisms, and characterized phages from an engineering perspective. First, the current design strategy for antisense targets in prokaryotic organisms was improved through the development of an algorithm that minimized the number of off-targets. For the phage work, a framework was developed to predict the disinfection rates in terms of the initial phage and bacterial concentrations. In addition, the long-term bacterial growth inhibition potential of multiple phages was determined for several bacteria. In regard to the phage application, phages were shown to protect both final product yields and yeast concentrations during fermentation. Taken together, this work suggests that the rational design of phage treatment is possible and further work is needed to expand on this foundation.

The Function and Regulation of Photobodies in Phytochrome Signaling

Light is a critical environmental signal that regulates every phase of the plant life cycle, from germination to floral initiation. Of the many light receptors in the model plant Arabidopsis thaliana, the red- and far-red light-sensing phytochromes (phys) are arguably the best studied, but the earliest events in the phy signaling pathway remain poorly understood. One of the earliest phy signaling events is the translocation of photoactivated phys from the cytoplasm to the nucleus, where they localize to subnuclear foci termed photobodies; in continuous light, photobody localization correlates closely with the light-dependent inhibition of embryonic stem growth. Despite a growing body of evidence supporting the biological significance of photobodies in light signaling, photobodies have also been shown to be dispensable for seedling growth inhibition in continuous light, so their physiological importance remains controversial; additionally, the molecular components that are required for phy localization to photobodies are largely unknown. The overall goal of my dissertation research was to gain insight into the early steps of phy signaling by further defining the role of photobodies in this process and identifying additional intragenic and extragenic requirements for phy localization to photobodies.

Even though the domain structure of phys has been extensively studied, not all of the intramolecular requirements for phy localization to photobodies are known. Previous studies have shown that the entire C-terminus of phys is both necessary and sufficient for their localization to photobodies. However, the importance of the individual subdomains of the C-terminus is still unclear. For example a truncation lacking part of the most C-terminal domain, the histidine kinase-related domain (HKRD), can still localize to small photobodies in the light and behaves like a weak allele. However, a point mutation within the HKRD renders the entire molecule completely inactive. To resolve this discrepancy, I explored the hypothesis that this point mutation might impair the dimerization of the HKRD; dimerization has been shown to occur via the C-terminus of phy and is required for more efficient signaling. I show that this point mutation impairs nuclear localization of phy as well as its subnuclear localization to photobodies. Additionally, yeast-two-hybrid analysis shows that the wild-type HKRD can homodimerize but that the HKRD containing the point mutation fails to dimerize with both itself and with wild-type HKRD. These results demonstrate that dimerization of the HKRD is required for both nuclear and photobody localization of phy.

Studies of seedlings grown in diurnal conditions show that photoactivated phy can persist into darkness to repress seedling growth; a seedling's growth rate is therefore fastest at the end of the night. To test the idea that photobodies could be involved in regulating seedling growth in the dark, I compared the growth of two transgenic Arabidopsis lines, one in which phy can localize to photobodies (PBG), and one in which it cannot (NGB). Despite these differences in photobody morphology, both lines are capable of transducing light signals and inhibiting seedling growth in continuous light. After the transition from red light to darkness, the PBG line was able to repress seedling growth, as well as the accumulation of the growth-promoting, light-labile transcription factor PHYTOCHROME INTERACTING FACTOR 3 (PIF3), for eighteen hours, and this correlated perfectly with the presence of photobodies. Reducing the amount of active phy by either reducing the light intensity or adding a phy-inactivating far-red pulse prior to darkness led to faster accumulation of PIF3 and earlier seedling growth. In contrast, the NGB line accumulated PIF3 even in the light, and seedling growth was only repressed for six hours; this behavior was similar in NGB regardless of the light treatment. These results suggest that photobodies are required for the degradation of PIF3 and for the prolonged stabilization of active phy in darkness. They also support the hypothesis that photobody localization of phys could serve as an instructive cue during the light-to-dark transition, thereby fine-tuning light-dependent responses in darkness.

In addition to determining an intragenic requirement for photobody localization and further exploring the significance of photobodies in phy signaling, I wanted to identify extragenic regulators of photobody localization. A recent study identified one such factor, HEMERA (HMR); hmr mutants do not form large photobodies, and they are tall and albino in the light. To identify other components in the HMR-mediated branch of the phy signaling pathway, I performed a forward genetic screen for suppressors of a weak hmr allele. Surprisingly, the first three mutants isolated from the screen were alleles of the same novel gene, SON OF HEMERA (SOH). The soh mutations rescue all of the phenotypes associated with the weak hmr allele, and they do so in an allele-specific manner, suggesting a direct interaction between SOH and HMR. Null soh alleles, which were isolated in an independent, tall, albino screen, are defective in photobody localization, demonstrating that SOH is an extragenic regulator of phy localization to photobodies that works in the same genetic pathway as HMR.

In this work, I show that dimerization of the HKRD is required for both the nuclear and photobody localization of phy. I also demonstrate a tight correlation between photobody localization and PIF3 degradation, further establishing the significance of photobodies in phy signaling. Finally, I identify a novel gene, SON OF HEMERA, whose product is necessary for phy localization to photobodies in the light, thereby isolating a new extragenic determinant of photobody localization. These results are among the first to focus exclusively on one of the earliest cellular responses to light - photobody localization of phys - and they promise to open up new avenues into the study of a poorly understood facet of the phy signaling pathway.

HDAC inhibitors cause site-specific chromatin remodeling at PU.1-bound enhancers in K562 cells

BACKGROUND: Small molecule inhibitors of histone deacetylases (HDACi) hold promise as anticancer agents for particular malignancies. However, clinical use is often confounded by toxicity, perhaps due to indiscriminate hyperacetylation of cellular proteins. Therefore, elucidating the mechanisms by which HDACi trigger differentiation, cell cycle arrest, or apoptosis of cancer cells could inform development of more targeted therapies. We used the myelogenous leukemia line K562 as a model of HDACi-induced differentiation to investigate chromatin accessibility (DNase-seq) and expression (RNA-seq) changes associated with this process. RESULTS: We identified several thousand specific regulatory elements [~10 % of total DNase I-hypersensitive (DHS) sites] that become significantly more or less accessible with sodium butyrate or suberanilohydroxamic acid treatment. Most of the differential DHS sites display hallmarks of enhancers, including being enriched for non-promoter regions, associating with nearby gene expression changes, and increasing luciferase reporter expression in K562 cells. Differential DHS sites were enriched for key hematopoietic lineage transcription factor motifs, including SPI1 (PU.1), a known pioneer factor. We found PU.1 increases binding at opened DHS sites with HDACi treatment by ChIP-seq, but PU.1 knockdown by shRNA fails to block the chromatin accessibility and expression changes. A machine-learning approach indicates H3K27me3 initially marks PU.1-bound sites that open with HDACi treatment, suggesting these sites are epigenetically poised. CONCLUSIONS: We find HDACi treatment of K562 cells results in site-specific chromatin remodeling at epigenetically poised regulatory elements. PU.1 shows evidence of a pioneer role in this process by marking poised enhancers but is not required for transcriptional activation.