Transcriptional analysis of persistent Chlamydia pneumoniae infection in vitro

Chlamydia pneumoniae can cause acute respiratory infections including pneumonia. Repeated and persistent Chlamydia infections occur and persistent C. pneumoniae infection may have a role in the pathogenesis of atherosclerosis and coronary heart disease and may also contribute to the development of chronic inflammatory lung diseases like chronic obstructive pulmonary disease (COPD) and asthma. In this thesis in vitro models for persistent C. pneumonia infection were established in epithelial and monocyte/macrophage cell lines. Expression of host cell genes in the persistent C. pneumoniae infection model of epithelial cells was studied by microarray and RT-PCR. In the monocyte/macrophage infection model expression of selected C. pneumoniae genes were studied by RT-PCR and immunofluorescence microscopy. Chlamydia is able to modulate host cell gene expression and apoptosis of host cells, which may assist Chlamydia to evade the host cells' immune responses. This, in turn, may lead to extended survival of the organism inside epithelial cells and promote the development of persistent infection. To simulate persistent C. pneumoniae infection in vivo, we set up a persistent infection model exposing the HL cell cultures to IFN-gamma. When HL cell cultures were treated with moderate concentration of IFN-gamma, the replication of C. pneumoniae DNA was unaffected while differentiation into infectious elementary bodies (EB) was strongly inhibited. By transmission electron microscopy small atypical inclusions were identified in IFN-gamma treated cultures. No second cycle of infection was observed in cells exposed to IFN-gamma , whereas C. pneumoniae was able to undergo a second cycle of infection in unexposed HL cells. Although monocytic cells can naturally restrict chlamydial growth, IFN-gamma further reduced production of infectious C. pneumoniae in Mono Mac 6 cells. Under both studied conditions no second cycle of infection could be detected in monocytic cell line suggesting persistent infection in these cells. As a step toward understanding the role of host genes in the development and pathogenesis of persistent C. pneumoniae infection, modulation of host cell gene expression during IFN-gamma induced persistent infection was examined and compared to that seen during active C. pneumoniae infection or IFN-gamma treatment. Total RNA was collected at 6 to 150 h after infection of an epithelial cell line (HL) and analyzed by a cDNA array (available at that time) representing approximately 4000 human transcripts. In initial analysis 250 of the 4000 genes were identified as differentially expressed upon active and persistent chlamydial infection and IFN-gamma treatment. In persistent infection more potent up-regulation of many genes was observed in IFN-gamma induced persistent infection than in active infection or in IFN-gamma treated cell cultures. Also sustained up-regulation was observed for some genes. In addition, we could identify nine host cell genes whose transcription was specifically altered during the IFN-gamma induced persistent C. pneumoniae infection. Strongest up-regulation in persistent infection in relation to controls was identified for insulin like growth factor binding protein 6, interferon-stimulated protein 15 kDa, cyclin D1 and interleukin 7 receptor. These results suggest that during persistent infection, C. pneumoniae reprograms the host transcriptional machinery regulating a variety of cellular processes including adhesion, cell cycle regulation, growth and inflammatory response, all of which may play important roles in the pathogenesis of persistent C. pneumoniae infection. C. pneumoniae DNA can be detected in peripheral blood mononuclear cells indicating that the bacterium can also infect monocytic cells in vivo and thereby monocytes can assist the spread of infection from the lungs to other anatomical sites. Persistent infection established at these sites could promote inflammation and enhance pathology. Thus, the mononuclear cells are in a strategic position in the development of persistent infection. To investigate the intracellular replication and fate of C. pneumoniae in mononuclear cells we analyzed the transcription of 11 C. pneumoniae genes in Mono Mac 6 cells during infection by real time RT-PCR. Our results suggest that the transcriptional profile of the studied genes in monocytes is different from that seen in epithelial cells and that IFN-gamma has a less significant effect on C. pneumoniae transcription in monocytes. Furthermore, our study shows that type III secretion system (T3SS) related genes are transcribed and that Chlamydia possesses a functional T3SS during infection in monocytes. Since C. pneumoniae infection in monocytes has been implicated to have reduced antibiotic susceptibility, this creates opportunities for novel therapeutics targeting T3SS in the management of chlamydial infection in monocytes.

Studies of the Human Transcriptome

Gene expression is one of the most critical factors influencing the phenotype of a cell. As a result of several technological advances, measuring gene expression levels has become one of the most common molecular biological measurements to study the behaviour of cells. The scientific community has produced enormous and constantly increasing collection of gene expression data from various human cells both from healthy and pathological conditions. However, while each of these studies is informative and enlighting in its own context and research setup, diverging methods and terminologies make it very challenging to integrate existing gene expression data to a more comprehensive view of human transcriptome function. On the other hand, bioinformatic science advances only through data integration and synthesis. The aim of this study was to develop biological and mathematical methods to overcome these challenges and to construct an integrated database of human transcriptome as well as to demonstrate its usage. Methods developed in this study can be divided in two distinct parts. First, the biological and medical annotation of the existing gene expression measurements needed to be encoded by systematic vocabularies. There was no single existing biomedical ontology or vocabulary suitable for this purpose. Thus, new annotation terminology was developed as a part of this work. Second part was to develop mathematical methods correcting the noise and systematic differences/errors in the data caused by various array generations. Additionally, there was a need to develop suitable computational methods for sample collection and archiving, unique sample identification, database structures, data retrieval and visualization. Bioinformatic methods were developed to analyze gene expression levels and putative functional associations of human genes by using the integrated gene expression data. Also a method to interpret individual gene expression profiles across all the healthy and pathological tissues of the reference database was developed. As a result of this work 9783 human gene expression samples measured by Affymetrix microarrays were integrated to form a unique human transcriptome resource GeneSapiens. This makes it possible to analyse expression levels of 17330 genes across 175 types of healthy and pathological human tissues. Application of this resource to interpret individual gene expression measurements allowed identification of tissue of origin with 92.0% accuracy among 44 healthy tissue types. Systematic analysis of transcriptional activity levels of 459 kinase genes was performed across 44 healthy and 55 pathological tissue types and a genome wide analysis of kinase gene co-expression networks was done. This analysis revealed biologically and medically interesting data on putative kinase gene functions in health and disease. Finally, we developed a method for alignment of gene expression profiles (AGEP) to perform analysis for individual patient samples to pinpoint gene- and pathway-specific changes in the test sample in relation to the reference transcriptome database. We also showed how large-scale gene expression data resources can be used to quantitatively characterize changes in the transcriptomic program of differentiating stem cells. Taken together, these studies indicate the power of systematic bioinformatic analyses to infer biological and medical insights from existing published datasets as well as to facilitate the interpretation of new molecular profiling data from individual patients.

Integration of miRNA and mRNA expression with DNA copy number of Ewing sarcoma cell lines

Ewing sarcoma is an aggressive and poorly differentiated malignancy of bone and soft tissue. It primarily affects children, adolescents, and young adults, with a slight male predominance. It is characterized by a translocation between chromosomes 11 and 22 resulting in the EWSR1-FLI1fusion transcription factor. The aim of this study is to identify putative Ewing sarcoma target genes through an integrative analysis of three microarray data sets. Array comparative genomic hybridization is used to measure changes in DNA copy number, and analyzed to detect common chromosomal aberrations. mRNA and miRNA microarrays are used to measure expression of protein-coding and miRNA genes, and these results integrated with the copy number data. Chromosomal aberrations typically contain also bystanders in addition to the driving tumor suppressor and oncogenes, and integration with expression helps to identify the true targets. Correlation between expression of miRNAs and their predicted target mRNAs is also evaluated to assess the results of post-transcriptional miRNA regulation on mRNA levels. The highest frequencies of copy number gains were identified in chromosome 8, 1q, and X. Losses were most frequent in 9p21.3, which also showed an enrichment of copy number breakpoints relative to the rest of the genome. Copy number losses in 9p21.3 were found have a statistically significant effect on the expression of MTAP, but not on CDKN2A, which is a known tumor-suppressor in the same locus. MTAP was also down-regulated in the Ewing sarcoma cell lines compared to mesenchymal stem cells. Genes exhibiting elevated expression in association with copy number gains and up-regulation compared to the reference samples included DCAF7, ENO2, MTCP1, andSTK40. Differentially expressed miRNAs were detected by comparing Ewing sarcoma cell lines against mesenchymal stem cells. 21 up-regulated and 32 down-regulated miRNAs were identified, includingmiR-145, which has been previously linked to Ewing sarcoma. The EWSR1-FLI1 fusion gene represses miR-145, which in turn targets FLI1 forming a mutually repressive feedback loop. In addition higher expression linked to copy number gains and compared to mesenchymal stem cells, STK40 was also found to be a target of four different miRNAs that were all down-regulated in Ewing sarcoma cell lines compared to the reference samples. SLCO5A1 was identified as the only up-regulated gene within a frequently gained region in chromosome 8. This region was gained in over 90 % of the cell lines, and also with a higher frequency than the neighboring regions. In addition, SLCO5A1 was found to be a target of three miRNAs that were down-regulated compared to the mesenchymal stem cells.

Regulation and Function of GATA Transcription Factors in Adrenocortical Tumors and Granulosa Cells

Transcription factors play a key role in tumor development, in which dysfunction of genes regulating tissue growth and differentiation is a central phenomenon. The GATA family of transcription factors consists of six members that bind to a consensus DNA sequence (A/T)GATA(A/G) in gene promoters and enhancers. The two GATA factors expressed in the adrenal cortex are GATA-4 and GATA-6. In both mice and humans, GATA-4 can be detected only during the fetal period, whereas GATA-6 expression is abundant both throughout development and in the adult. It is already established that GATA factors are important in both normal development and tumorigenesis of several endocrine organs, and expression of GATA-4 and GATA-6 is detected in adrenocortical tumors. The aim of this study was to elucidate the function of these factors in adrenocortical tumor growth. In embryonal development, the adrenocortical cells arise and differentiate from a common pool with gonadal steroidogenic cells, the urogenital ridge. As the adult adrenal cortex undergoes constant renewal, it is hypothesized that undifferentiated adrenocortical progenitor cells reside adjacent to the adrenal capsule and give rise to daughter cells that differentiate and migrate centripetally. A diverse array of hormones controls the differentiation, growth and survival of steroidogenic cells in the adrenal gland and the gonads. Factors such as luteinizing hormone and inhibins, traditionally associated with gonadal steroidogenic cells, can also influence the function of adrenocortical cells in physiological and pathophysiological states. Certain inbred strains of mice develop subcapsular adrenocortical tumors in response to gonadectomy. In this study, we found that these tumors express GATA-4, normally absent from the adult adrenal cortex, while GATA-6 expression is downregulated. Gonadal markers such as luteinizing hormone receptor, anti-Müllerian hormone and P450c17 are also expressed in the neoplastic cells, and the tumors produce gonadal hormones. The tumor cells have lost the expression of melanocortin-2 receptor and the CYP enzymes necessary for the synthesis of corticosterone and aldosterone. By way of xenograft studies utilizing NU/J nude mice, we confirmed that chronic gonadotropin elevation is sufficient to induce adrenocortical tumorigenesis in susceptible inbred strains. Collectively, these studies suggest that subcapsular adrenocortical progenitor cells can, under certain conditions, adopt a gonadal fate. We studied the molecular mechanisms involved in gene regulation in endocrine cells in order to elucidate the role of GATA factors in endocrine tissues. Ovarian granulosa cells express both GATA-4 and GATA-6, and the TGF-β signaling pathway is active in these cells. Inhibin-α is both a target gene for, and an atypical or antagonistic member of the TGF-β growth factor superfamily. In this study, we show that GATA-4 is required for TGF-β-mediated inhibin-α promoter activation in granulosa cells, and that GATA-4 physically interacts with Smad3, a TGF-β downstream protein. Apart from the regulation of steroidogenesis and other events in normal tissues, TGF-β signaling is implicated in tumors of multiple organs, including the adrenal cortex. Another signaling pathway found often to be aberrantly active in adrenocortical tumors is the Wnt pathway. As both of these pathways regulate the expression of inhibin-α, a transcriptional target for GATA-4 and GATA-6, we wanted to investigate whether GATA factors are associated with the components of these signaling cascades in human adrenocortical tumors. We found that the expression of Wnt co-receptors LRP5 and LRP6, Smad3, GATA-6 and SF-1 was diminished in adrenocortical carcinomas with poor outcome. All of these factors drive inhibin-α expression, and their expression in adrenocortical tumors correlated with that of inhibin-α. The results support a tumor suppressor role previously suggested for inhibin-α in the mouse adrenal cortex, and offer putative pathways associated with adrenocortical tumor aggressiveness. Unraveling the role of GATA factors and associated molecules in human and mouse adrenocortical tumors could ultimately contribute to the development of diagnostic tools and future therapies for these diseases.

Implementation of a phase array diffuse optical tomographic imager

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180° along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.

Studying the resistivity imaging of chicken tissue phantoms with different current patterns in Electrical Impedance Tomography (EIT)

A current injection pattern in Electrical Impedance Tomography (EIT) has its own current distribution profile within the domain under test. Hence, different current patterns have different sensitivity, spatial resolution and distinguishability. Image reconstruction studies with practical phantoms are essential to assess the performance of EIT systems for their validation, calibration and comparison purposes. Impedance imaging of real tissue phantoms with different current injection methods is also essential for better assessment of the biomedical EIT systems. Chicken tissue paste phantoms and chicken tissue block phantoms are developed and the resistivity image reconstruction is studied with different current injection methods. A 16-electrode array is placed inside the phantom tank and the tank is filled with chicken muscle tissue paste or chicken tissue blocks as the background mediums. Chicken fat tissue, chicken bone, air hole and nylon cylinders are used as the inhomogeneity to obtained different phantom configurations. A low magnitude low frequency constant sinusoidal current is injected at the phantom boundary with opposite and neighboring current patterns and the boundary potentials are measured. Resistivity images are reconstructed from the boundary data using EIDORS and the reconstructed images are analyzed with the contrast parameters calculated from their elemental resistivity profiles. Results show that the resistivity profiles of all the phantom domains are successfully reconstructed with a proper background resistivity and high inhomogeneity resistivity for both the current injection methods. Reconstructed images show that, for all the chicken tissue phantoms, the inhomogeneities are suitably reconstructed with both the current injection protocols though the chicken tissue block phantom and opposite method are found more suitable. It is observed that the boundary potentials of the chicken tissue block phantoms are higher than the chicken tissue paste phantom. SNR of the chicken tissue block phantoms are found comparatively more and hence the chicken tissue block phantom is found more suitable for its lower noise performance. The background noise is found less in opposite method for all the phantom configurations which yields the better resistivity images with high PCR and COC and proper IRMean and IRMax neighboring method showed higher noise level for both the chicken tissue paste phantoms and chicken tissue block phantoms with all the inhomogeneities. Opposite method is found more suitable for both the chicken tissue phantoms, and also, chicken tissue block phantoms are found more suitable compared to the chicken tissue paste phantom. (C) 2012 Elsevier Ltd. All rights reserved.

Genomic mapping of cAMP receptor protein (CRPMt) in Mycobacterium tuberculosis: relation to transcriptional start sites and the role of CRPMt as a transcription factor

Chromatin immunoprecipitation identified 191 binding sites of Mycobacterium tuberculosis cAMP receptor protein (CRPMt) at endogenous expression levels using a specific alpha-CRPMt antibody. Under these native conditions an equal distribution between intragenic and intergenic locations was observed. CRPMt binding overlapped a palindromic consensus sequence. Analysis by RNA sequencing revealed widespread changes in transcriptional profile in a mutant strain lacking CRPMt during exponential growth, and in response to nutrient starvation. Differential expression of genes with a CRPMt-binding site represented only a minor portion of this transcriptional reprogramming with similar to 19% of those representing transcriptional regulators potentially controlled by CRPMt. The subset of genes that are differentially expressed in the deletion mutant under both culture conditions conformed to a pattern resembling canonical CRP regulation in Escherichia coli, with binding close to the transcriptional start site associated with repression and upstream binding with activation. CRPMt can function as a classical transcription factor in M. tuberculosis, though this occurs at only a subset of CRPMt-binding sites.

The mycobacterial iron-dependent regulator IdeR induces ferritin (bfrB) by alleviating Lsr2 repression

Emerging evidence indicates that precise regulation of iron (Fe) metabolism and maintenance of Fe homeostasis in Mycobacterium tuberculosis (Mtb) are essential for its survival and proliferation in the host. IdeR is a central transcriptional regulator of Mtb genes involved in Fe metabolism. While it is well understood how IdeR functions as a repressor, how it induces transcription of a subset of its targets is still unclear. We investigated the molecular mechanism of IdeR-mediated positive regulation of bfrB, the gene encoding the major Fe-storage protein of Mtb. We found that bfrB induction by Fe required direct interaction of IdeR with a DNA sequence containing four tandem IdeR-binding boxes located upstream of the bfrB promoter. Results of in vivo and in vitro transcription assays identified a direct repressor of bfrB, the histone-like protein Lsr2. IdeR counteracted Lsr2-mediated repression in vitro, suggesting that IdeR induces bfrB transcription by antagonizing the repressor activity of Lsr2. Together, these results elucidate the main mechanism of bfrB positive regulation by IdeR and identify Lsr2 as a new factor contributing to Fe homeostasis in mycobacteria.

DNA Damage and Transcriptional Changes in the Gills of Mytilus galloprovincialis Exposed to Nanomolar Doses of Combined Metal Salts (Cd, Cu, Hg)

[ENG]Aiming at an integrated and mechanistic view of the early biological effects of selected metals in the marine sentinel organism Mytilus galloprovincialis, we exposed mussels for 48 hours to 50, 100 and 200 nM solutions of equimolar Cd, Cu and Hg salts and measured cytological and molecular biomarkers in parallel. Focusing on the mussel gills, first target of toxic water contaminants and actively proliferating tissue, we detected significant dose-related increases of cells with micronuclei and other nuclear abnormalities in the treated mussels, with differences in the bioconcentration of the three metals determined in the mussel flesh by atomic absorption spectrometry. Gene expression profiles, determined in the same individual gills in parallel, revealed some transcriptional changes at the 50 nM dose, and substantial increases of differentially expressed genes at the 100 and 200 nM doses, with roughly similar amounts of up- and down-regulated genes. The functional annotation of gill transcripts with consistent expression trends and significantly altered at least in one dose point disclosed the complexity of the induced cell response. The most evident transcriptional changes concerned protein synthesis and turnover, ion homeostasis, cell cycle regulation and apoptosis, and intracellular trafficking (transcript sequences denoting heat shock proteins, metal binding thioneins, sequestosome 1 and proteasome subunits, and GADD45 exemplify up-regulated genes while transcript sequences denoting actin, tubulins and the apoptosis inhibitor 1 exemplify down-regulated genes). Overall, nanomolar doses of co-occurring free metal ions have induced significant structural and functional changes in the mussel gills: the intensity of response to the stimulus measured in laboratory supports the additional validation of molecular markers of metal exposure to be used in Mussel Watch programs

Sequence-function relationships in E. coli transcriptional regulation

Understanding how transcriptional regulatory sequence maps to regulatory function remains a difficult problem in regulatory biology. Given a particular DNA sequence for a bacterial promoter region, we would like to be able to say which transcription factors bind there, how strongly they bind, and whether they interact with each other and/or RNA polymerase, with the ultimate objective of integrating knowledge of these parameters into a prediction of gene expression levels. The theoretical framework of statistical thermodynamics provides a useful framework for doing so, enabling us to predict how gene expression levels depend on transcription factor binding energies and concentrations. We used thermodynamic models, coupled with models of the sequence-dependent binding energies of transcription factors and RNAP, to construct a genotype to phenotype map for the level of repression exhibited by the lac promoter, and tested it experimentally using a set of promoter variants from E. coli strains isolated from different natural environments. For this work, we sought to ``reverse engineer'' naturally occurring promoter sequences to understand how variations in promoter sequence affects gene expression. The natural inverse of this approach is to ``forward engineer'' promoter sequences to obtain targeted levels of gene expression. We used a high precision model of RNAP-DNA sequence dependent binding energy, coupled with a thermodynamic model relating binding energy to gene expression, to predictively design and verify a suite of synthetic E. coli promoters whose expression varied over nearly three orders of magnitude.

However, although thermodynamic models enable predictions of mean levels of gene expression, it has become evident that cell-to-cell variability or ``noise'' in gene expression can also play a biologically important role. In order to address this aspect of gene regulation, we developed models based on the chemical master equation framework and used them to explore the noise properties of a number of common E. coli regulatory motifs; these properties included the dependence of the noise on parameters such as transcription factor binding strength and copy number. We then performed experiments in which these parameters were systematically varied and measured the level of variability using mRNA FISH. The results showed a clear dependence of the noise on these parameters, in accord with model predictions.

Finally, one shortcoming of the preceding modeling frameworks is that their applicability is largely limited to systems that are already well-characterized, such as the lac promoter. Motivated by this fact, we used a high throughput promoter mutagenesis assay called Sort-Seq to explore the completely uncharacterized transcriptional regulatory DNA of the E. coli mechanosensitive channel of large conductance (MscL). We identified several candidate transcription factor binding sites, and work is continuing to identify the associated proteins.

A Detailed Analysis of Transcriptional Regulators Affecting the Saccharomyces cerevisiae Heat-shock Gene SSA1

The yeast Saccharomyces cerevisiae contains a family of hsp70 related genes. One member of this family, SSA1, encodes a 70kD heat-shock protein which in addition to its heat inducible expression has a significant basal level of expression. The first 500 bp upstream of the SSA1 start point of transcription was examined by DNAse I protection analysis. The results reveal the presence of at least 14 factor binding sites throughout the upstream promoter region. The function of these binding sites has been examined using a series of 5' promoter deletions fused to the recorder gene lacZ in a centromere-containing yeast shuttle vector. The following sites have been identified in the promoter and their activity in yeast determined individually with a centromere-based recorder plasmid containing a truncated CYC1 /lacZ fusion: a heat-shock element or HSE which is sufficient to convey heat-shock response on the recorder plasmid; a homology to the SV40 'core' sequence which can repress the GCN4 recognition element (GCRE) and the yAP1 recognition element (ARE), and has been designated a upstream repression element or URE; a 'G'-rich region named G-box which can also convey heatshock response on the recorder plasmid; and a purine-pyrimidine alternating sequence name GT-box which is an activator of transcription. A series of fusion constructs were made to identify a putative silencer-like element upstream of SSA1. This element is position dependent and has been localized to a region containing both an ABF1 binding site and a RAP1 binding site. Five site-specific DNA-binding factors are identified and their purification is presented: the heat-shock transcription factor or HSTF, which recognizes the HSE; the G-box binding factor or GBF; the URE recognition factor or URF; the GT-box binding factor; and the GC-box binding factor or yeast Sp1.

Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis.

Cytosine methylation is important for transposon silencing and epigenetic regulation of endogenous genes, although the extent to which this DNA modification functions to regulate the genome is still unknown. Here we report the first comprehensive DNA methylation map of an entire genome, at 35 base pair resolution, using the flowering plant Arabidopsis thaliana as a model. We find that pericentromeric heterochromatin, repetitive sequences, and regions producing small interfering RNAs are heavily methylated. Unexpectedly, over one-third of expressed genes contain methylation within transcribed regions, whereas only approximately 5% of genes show methylation within promoter regions. Interestingly, genes methylated in transcribed regions are highly expressed and constitutively active, whereas promoter-methylated genes show a greater degree of tissue-specific expression. Whole-genome tiling-array transcriptional profiling of DNA methyltransferase null mutants identified hundreds of genes and intergenic noncoding RNAs with altered expression levels, many of which may be epigenetically controlled by DNA methylation.

Nonlinear Acoustics in Underwater and Biomedical Applications: Array Performance Degradation and Time Reversal Invariance

This dissertation describes a model for acoustic propagation in inhomogeneous flu- ids, and explores the focusing by arrays onto targets under various conditions. The work explores the use of arrays, in particular the time reversal array, for underwater and biomedical applications. Aspects of propagation and phasing which can lead to reduced focusing effectiveness are described. An acoustic wave equation was derived for the propagation of finite-amplitude waves in lossy time-varying inhomogeneous fluid media. The equation was solved numerically in both Cartesian and cylindrical geometries using the finite-difference time-domain (FDTD) method. It was found that time reversal arrays are sensitive to several debilitating factors. Focusing ability was determined to be adequate in the presence of temporal jitter in the time reversed signal only up to about one-sixth of a period. Thermoviscous absorption also had a debilitating effect on focal pressure for both linear and nonlinear propagation. It was also found that nonlinearity leads to degradation of focal pressure through amplification of the received signal at the array, and enhanced absorption in the shocked waveforms. This dissertation also examined the heating effects of focused ultrasound in a tissue-like medium. The application considered is therapeutic heating for hyperther- mia. The acoustic model and a thermal model for tissue were coupled to solve for transient and steady temperature profiles in tissue-like media. The Pennes bioheat equation was solved using the FDTD method to calculate the temperature fields in tissue-like media from focused acoustic sources. It was found that the temperature-dependence of the medium's background prop- erties can play an important role in the temperature predictions. Finite-amplitude effects contributed excess heat when source conditions were provided for nonlinear ef- fects to manifest themselves. The effect of medium heterogeneity was also found to be important in redistributing the acoustic and temperature fields, creating regions with hotter and colder temperatures than the mean by local scattering and lensing action. These temperature excursions from the mean were found to increase monotonically with increasing contrast in the medium's properties.

A host transcriptional signature for presymptomatic detection of infection in humans exposed to influenza H1N1 or H3N2.

There is great potential for host-based gene expression analysis to impact the early diagnosis of infectious diseases. In particular, the influenza pandemic of 2009 highlighted the challenges and limitations of traditional pathogen-based testing for suspected upper respiratory viral infection. We inoculated human volunteers with either influenza A (A/Brisbane/59/2007 (H1N1) or A/Wisconsin/67/2005 (H3N2)), and assayed the peripheral blood transcriptome every 8 hours for 7 days. Of 41 inoculated volunteers, 18 (44%) developed symptomatic infection. Using unbiased sparse latent factor regression analysis, we generated a gene signature (or factor) for symptomatic influenza capable of detecting 94% of infected cases. This gene signature is detectable as early as 29 hours post-exposure and achieves maximal accuracy on average 43 hours (p = 0.003, H1N1) and 38 hours (p-value = 0.005, H3N2) before peak clinical symptoms. In order to test the relevance of these findings in naturally acquired disease, a composite influenza A signature built from these challenge studies was applied to Emergency Department patients where it discriminates between swine-origin influenza A/H1N1 (2009) infected and non-infected individuals with 92% accuracy. The host genomic response to Influenza infection is robust and may provide the means for detection before typical clinical symptoms are apparent.

Acoustic radiation force impulse imaging (ARFI) on an IVUS circular array.

Our long-term goal is the detection and characterization of vulnerable plaque in the coronary arteries of the heart using intravascular ultrasound (IVUS) catheters. Vulnerable plaque, characterized by a thin fibrous cap and a soft, lipid-rich necrotic core is a precursor to heart attack and stroke. Early detection of such plaques may potentially alter the course of treatment of the patient to prevent ischemic events. We have previously described the characterization of carotid plaques using external linear arrays operating at 9 MHz. In addition, we previously modified circular array IVUS catheters by short-circuiting several neighboring elements to produce fixed beamwidths for intravascular hyperthermia applications. In this paper, we modified Volcano Visions 8.2 French, 9 MHz catheters and Volcano Platinum 3.5 French, 20 MHz catheters by short-circuiting portions of the array for acoustic radiation force impulse imaging (ARFI) applications. The catheters had an effective transmit aperture size of 2 mm and 1.5 mm, respectively. The catheters were connected to a Verasonics scanner and driven with pushing pulses of 180 V p-p to acquire ARFI data from a soft gel phantom with a Young's modulus of 2.9 kPa. The dynamic response of the tissue-mimicking material demonstrates a typical ARFI motion of 1 to 2 microns as the gel phantom displaces away and recovers back to its normal position. The hardware modifications applied to our IVUS catheters mimic potential beamforming modifications that could be implemented on IVUS scanners. Our results demonstrate that the generation of radiation force from IVUS catheters and the development of intravascular ARFI may be feasible.