Temporal dynamics of host molecular responses differentiate symptomatic and asymptomatic influenza a infection.

Exposure to influenza viruses is necessary, but not sufficient, for healthy human hosts to develop symptomatic illness. The host response is an important determinant of disease progression. In order to delineate host molecular responses that differentiate symptomatic and asymptomatic Influenza A infection, we inoculated 17 healthy adults with live influenza (H3N2/Wisconsin) and examined changes in host peripheral blood gene expression at 16 timepoints over 132 hours. Here we present distinct transcriptional dynamics of host responses unique to asymptomatic and symptomatic infections. We show that symptomatic hosts invoke, simultaneously, multiple pattern recognition receptors-mediated antiviral and inflammatory responses that may relate to virus-induced oxidative stress. In contrast, asymptomatic subjects tightly regulate these responses and exhibit elevated expression of genes that function in antioxidant responses and cell-mediated responses. We reveal an ab initio molecular signature that strongly correlates to symptomatic clinical disease and biomarkers whose expression patterns best discriminate early from late phases of infection. Our results establish a temporal pattern of host molecular responses that differentiates symptomatic from asymptomatic infections and reveals an asymptomatic host-unique non-passive response signature, suggesting novel putative molecular targets for both prognostic assessment and ameliorative therapeutic intervention in seasonal and pandemic influenza.

Stochastic E2F activation and reconciliation of phenomenological cell-cycle models.

The transition of the mammalian cell from quiescence to proliferation is a highly variable process. Over the last four decades, two lines of apparently contradictory, phenomenological models have been proposed to account for such temporal variability. These include various forms of the transition probability (TP) model and the growth control (GC) model, which lack mechanistic details. The GC model was further proposed as an alternative explanation for the concept of the restriction point, which we recently demonstrated as being controlled by a bistable Rb-E2F switch. Here, through a combination of modeling and experiments, we show that these different lines of models in essence reflect different aspects of stochastic dynamics in cell cycle entry. In particular, we show that the variable activation of E2F can be described by stochastic activation of the bistable Rb-E2F switch, which in turn may account for the temporal variability in cell cycle entry. Moreover, we show that temporal dynamics of E2F activation can be recast into the frameworks of both the TP model and the GC model via parameter mapping. This mapping suggests that the two lines of phenomenological models can be reconciled through the stochastic dynamics of the Rb-E2F switch. It also suggests a potential utility of the TP or GC models in defining concise, quantitative phenotypes of cell physiology. This may have implications in classifying cell types or states.

Exploring Mechanisms of Bacterial Adaptation to Seasonal Temperature Change

This research examines three potential mechanisms by which bacteria can adapt to different temperatures: changes in strain-level population structure, gene regulation and particle colonization. For the first two mechanisms, I utilize bacterial strains from the Vibrionaceae family due to their ease of culturability, ubiquity in coastal environments and status as a model system for marine bacteria. I first examine vibrio seasonal dynamics in temperate, coastal water and compare the thermal performance of strains that occupy different thermal environments. Our results suggest that there are tradeoffs in adaptation to specific temperatures and that thermal specialization can occur at a very fine phylogenetic scale. The observed thermal specialization over relatively short evolutionary time-scales indicates that few genes or cellular processes may limit expansion to a different thermal niche. I then compare the genomic and transcriptional changes associated with thermal adaptation in closely-related vibrio strains under heat and cold stress. The two vibrio strains have very similar genomes and overall exhibit similar transcriptional profiles in response to temperature stress but their temperature preferences are determined by differential transcriptional responses in shared genes as well as temperature-dependent regulation of unique genes. Finally, I investigate the temporal dynamics of particle-attached and free-living bacterial community in coastal seawater and find that microhabitats exert a stronger forcing on microbial communities than environmental variability, suggesting that particle-attachment could buffer the impacts of environmental changes and particle-associated communities likely respond to the presence of distinct eukaryotes rather than commonly-measured environmental parameters. Integrating these results will offer new perspectives on the mechanisms by which bacteria respond to seasonal temperature changes as well as potential adaptations to climate change-driven warming of the surface oceans.

Inheritance of the CENP-A chromatin domain is spatially and temporally constrained at human centromeres.

BACKGROUND: Chromatin containing the histone variant CENP-A (CEN chromatin) exists as an essential domain at every centromere and heritably marks the location of kinetochore assembly. The size of the CEN chromatin domain on alpha satellite DNA in humans has been shown to vary according to underlying array size. However, the average amount of CENP-A reported at human centromeres is largely consistent, implying the genomic extent of CENP-A chromatin domains more likely reflects variations in the number of CENP-A subdomains and/or the density of CENP-A nucleosomes within individual subdomains. Defining the organizational and spatial properties of CEN chromatin would provide insight into centromere inheritance via CENP-A loading in G1 and the dynamics of its distribution between mother and daughter strands during replication. RESULTS: Using a multi-color protein strategy to detect distinct pools of CENP-A over several cell cycles, we show that nascent CENP-A is equally distributed to sister centromeres. CENP-A distribution is independent of previous or subsequent cell cycles in that centromeres showing disproportionately distributed CENP-A in one cycle can equally divide CENP-A nucleosomes in the next cycle. Furthermore, we show using extended chromatin fibers that maintenance of the CENP-A chromatin domain is achieved by a cycle-specific oscillating pattern of new CENP-A nucleosomes next to existing CENP-A nucleosomes over multiple cell cycles. Finally, we demonstrate that the size of the CENP-A domain does not change throughout the cell cycle and is spatially fixed to a similar location within a given alpha satellite DNA array. CONCLUSIONS: We demonstrate that most human chromosomes share similar patterns of CENP-A loading and distribution and that centromere inheritance is achieved through specific placement of new CENP-A near existing CENP-A as assembly occurs each cell cycle. The loading pattern fixes the location and size of the CENP-A domain on individual chromosomes. These results suggest that spatial and temporal dynamics of CENP-A are important for maintaining centromere identity and genome stability.

Epigenetic Basis of Centromere Maintenance and Inheritance

Centromeres are essential chromosomal loci at which kinetochore formation occurs for spindle fiber attachment during mitosis and meiosis, guiding proper segregation of chromosomes. In humans, centromeres are located at large arrays of alpha satellite DNA, contributing to but not defining centromere function. The histone variant CENP-A assembles at alpha satellite DNA, epigenetically defining the centromere. CENP-A containing chromatin exists as an essential domain composed of blocks of CENP-A nucleosomes interspersed with blocks of H3 nucleosomes, and is surrounded by pericentromeric heterochromatin. In order to maintain genomic stability, the CENP-A domain is propagated epigenetically over each cell division; disruption of propagation is associated with chromosome instabilities such as aneuploidy, found in birth defects and in cancer.

The CENP-A chromatin domain occupies 30-45% of the alpha satellite array, varying in genomic distance according to the underlying array size. However, the molecular mechanisms that control assembly and organization of CENP-A chromatin within its genomic context remain unclear. The domain may shift, expand, or contract, as CENP-A is loaded and dispersed each cell cycle. We hypothesized that in order to maintain genome stability, the centromere is inherited as static chromatin domains, maintaining size and position within the pericentric heterochromatin. Utilizing stretched chromatin fibers, I found that CENP-A chromatin is limited to a sub-region of the alpha satellite array that is fixed in size and location through the cell cycle and across populations.

The average amount of CENP-A at human centromeres is largely consistent, implying that the variation in size of CENP-A domains reflects variations in the number of CENP-A subdomains and/or the density of CENP-A nucleosomes. Multi-color nascent protein labeling experiments were utilized to examine the distribution and incorporation of distinct pools of CENP-A over several cell cycles. I found that in each cell cycle there is independent CENP-A distribution, occurring equally between sister centromeres across all chromosomes, in similar quantities. Furthermore, centromere inheritance is achieved through specific placement of CENP-A, following an oscillating pattern that fixes the location and size of the CENP-A domain. These results suggest that spatial and temporal dynamics of CENP-A are important for maintaining centromere and genome stability.

Spike avalanches exhibit universal dynamics across the sleep-wake cycle.

BACKGROUND: Scale-invariant neuronal avalanches have been observed in cell cultures and slices as well as anesthetized and awake brains, suggesting that the brain operates near criticality, i.e. within a narrow margin between avalanche propagation and extinction. In theory, criticality provides many desirable features for the behaving brain, optimizing computational capabilities, information transmission, sensitivity to sensory stimuli and size of memory repertoires. However, a thorough characterization of neuronal avalanches in freely-behaving (FB) animals is still missing, thus raising doubts about their relevance for brain function. METHODOLOGY/PRINCIPAL FINDINGS: To address this issue, we employed chronically implanted multielectrode arrays (MEA) to record avalanches of action potentials (spikes) from the cerebral cortex and hippocampus of 14 rats, as they spontaneously traversed the wake-sleep cycle, explored novel objects or were subjected to anesthesia (AN). We then modeled spike avalanches to evaluate the impact of sparse MEA sampling on their statistics. We found that the size distribution of spike avalanches are well fit by lognormal distributions in FB animals, and by truncated power laws in the AN group. FB data surrogation markedly decreases the tail of the distribution, i.e. spike shuffling destroys the largest avalanches. The FB data are also characterized by multiple key features compatible with criticality in the temporal domain, such as 1/f spectra and long-term correlations as measured by detrended fluctuation analysis. These signatures are very stable across waking, slow-wave sleep and rapid-eye-movement sleep, but collapse during anesthesia. Likewise, waiting time distributions obey a single scaling function during all natural behavioral states, but not during anesthesia. Results are equivalent for neuronal ensembles recorded from visual and tactile areas of the cerebral cortex, as well as the hippocampus. CONCLUSIONS/SIGNIFICANCE: Altogether, the data provide a comprehensive link between behavior and brain criticality, revealing a unique scale-invariant regime of spike avalanches across all major behaviors.

Dynamics of visual receptive fields in the macaque frontal eye field.

Neuronal receptive fields (RFs) provide the foundation for understanding systems-level sensory processing. In early visual areas, investigators have mapped RFs in detail using stochastic stimuli and sophisticated analytical approaches. Much less is known about RFs in prefrontal cortex. Visual stimuli used for mapping RFs in prefrontal cortex tend to cover a small range of spatial and temporal parameters, making it difficult to understand their role in visual processing. To address these shortcomings, we implemented a generalized linear model to measure the RFs of neurons in the macaque frontal eye field (FEF) in response to sparse, full-field stimuli. Our high-resolution, probabilistic approach tracked the evolution of RFs during passive fixation, and we validated our results against conventional measures. We found that FEF neurons exhibited a surprising level of sensitivity to stimuli presented as briefly as 10 ms or to multiple dots presented simultaneously, suggesting that FEF visual responses are more precise than previously appreciated. FEF RF spatial structures were largely maintained over time and between stimulus conditions. Our results demonstrate that the application of probabilistic RF mapping to FEF and similar association areas is an important tool for clarifying the neuronal mechanisms of cognition.

Assessment of the Spatial and Temporal Distribution of Functional Connectivity in Resting-State BOLD fMRI

Recent research into resting-state functional magnetic resonance imaging (fMRI) has shown that the brain is very active during rest. This thesis work utilizes blood oxygenation level dependent (BOLD) signals to investigate the spatial and temporal functional network information found within resting-state data, and aims to investigate the feasibility of extracting functional connectivity networks using different methods as well as the dynamic variability within some of the methods. Furthermore, this work looks into producing valid networks using a sparsely-sampled sub-set of the original data.

In this work we utilize four main methods: independent component analysis (ICA), principal component analysis (PCA), correlation, and a point-processing technique. Each method comes with unique assumptions, as well as strengths and limitations into exploring how the resting state components interact in space and time.

Correlation is perhaps the simplest technique. Using this technique, resting-state patterns can be identified based on how similar the time profile is to a seed region’s time profile. However, this method requires a seed region and can only identify one resting state network at a time. This simple correlation technique is able to reproduce the resting state network using subject data from one subject’s scan session as well as with 16 subjects.

Independent component analysis, the second technique, has established software programs that can be used to implement this technique. ICA can extract multiple components from a data set in a single analysis. The disadvantage is that the resting state networks it produces are all independent of each other, making the assumption that the spatial pattern of functional connectivity is the same across all the time points. ICA is successfully able to reproduce resting state connectivity patterns for both one subject and a 16 subject concatenated data set.

Using principal component analysis, the dimensionality of the data is compressed to find the directions in which the variance of the data is most significant. This method utilizes the same basic matrix math as ICA with a few important differences that will be outlined later in this text. Using this method, sometimes different functional connectivity patterns are identifiable but with a large amount of noise and variability.

To begin to investigate the dynamics of the functional connectivity, the correlation technique is used to compare the first and second halves of a scan session. Minor differences are discernable between the correlation results of the scan session halves. Further, a sliding window technique is implemented to study the correlation coefficients through different sizes of correlation windows throughout time. From this technique it is apparent that the correlation level with the seed region is not static throughout the scan length.

The last method introduced, a point processing method, is one of the more novel techniques because it does not require analysis of the continuous time points. Here, network information is extracted based on brief occurrences of high or low amplitude signals within a seed region. Because point processing utilizes less time points from the data, the statistical power of the results is lower. There are also larger variations in DMN patterns between subjects. In addition to boosted computational efficiency, the benefit of using a point-process method is that the patterns produced for different seed regions do not have to be independent of one another.

This work compares four unique methods of identifying functional connectivity patterns. ICA is a technique that is currently used by many scientists studying functional connectivity patterns. The PCA technique is not optimal for the level of noise and the distribution of the data sets. The correlation technique is simple and obtains good results, however a seed region is needed and the method assumes that the DMN regions is correlated throughout the entire scan. Looking at the more dynamic aspects of correlation changing patterns of correlation were evident. The last point-processing method produces a promising results of identifying functional connectivity networks using only low and high amplitude BOLD signals.