A dimensionless number for understanding the evolutionary dynamics of antigenically variable RNA viruses.

Antigenically variable RNA viruses are significant contributors to the burden of infectious disease worldwide. One reason for their ubiquity is their ability to escape herd immunity through rapid antigenic evolution and thereby to reinfect previously infected hosts. However, the ways in which these viruses evolve antigenically are highly diverse. Some have only limited diversity in the long-run, with every emergence of a new antigenic variant coupled with a replacement of the older variant. Other viruses rapidly accumulate antigenic diversity over time. Others still exhibit dynamics that can be considered evolutionary intermediates between these two extremes. Here, we present a theoretical framework that aims to understand these differences in evolutionary patterns by considering a virus's epidemiological dynamics in a given host population. Our framework, based on a dimensionless number, probabilistically anticipates patterns of viral antigenic diversification and thereby quantifies a virus's evolutionary potential. It is therefore similar in spirit to the basic reproduction number, the well-known dimensionless number which quantifies a pathogen's reproductive potential. We further outline how our theoretical framework can be applied to empirical viral systems, using influenza A/H3N2 as a case study. We end with predictions of our framework and work that remains to be done to further integrate viral evolutionary dynamics with disease ecology.

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.

Dynamics of the Disk-Pendulum Coupled System With Vertical Excitation

This paper investigates the static and dynamic characteristics of the semi-elliptical rocking disk on which a pendulum pinned. This coupled system’s response is also analyzed analytically and numerically when a vertical harmonic excitation is applied to the bottom of the rocking disk. Lagrange’s Equation is used to derive the motion equations of the disk-pendulum coupled system. The second derivative test for the system’s potential energy shows how the location of the pendulum’s pivotal point affects the number and stability of equilibria, and the change of location presents different bifurcation diagrams for different geometries of the rocking disk. For both vertically excited and unforced cases, the coupled system shows chaos easily, but the proper chosen parameters can still help the system reach and keep the steady state. For the steady state of the vertically excited rocking disk without a pendulum, the variation of the excitation’s amplitude and frequency result in the hysteresis for the amplitude of the response. When a pendulum is pinned on the rocking disk, three major categories of steady states are presently in the numerical way.

Role of Vinculin in Regulating Force-Sensitive Dynamics of Adherens Junctions During Collective Cell Migration

Dynamic processes such as morphogenesis and tissue patterning require the precise control of many cellular processes, especially cell migration. Historically, these processes are thought to be mediated by genetic and biochemical signaling pathways. However, recent advances have unraveled a previously unappreciated role of mechanical forces in regulating these homeostatic processes in of multicellular systems. In multicellular systems cells adhere to both deformable extracellular matrix (ECM) and other cells, which are sources of applied forces and means of mechanical support. Cells detect and respond to these mechanical signals through a poorly understood process called mechanotransduction, which can have profound effects on processes such as cell migration. These effects are largely mediated by the sub cellular structures that link cells to the ECM, called focal adhesions (FAs), or cells to other cells, termed adherens junctions (AJs).

Overall this thesis is comprised of my work on identifying a novel force dependent function of vinculin, a protein which resides in both FAs and AJs - in dynamic process of collective migration. Using a collective migration assay as a model for collective cell behavior and a fluorescence resonance energy transfer (FRET) based molecular tension sensor for vinculin I demonstrated a spatial gradient of tension across vinculin in the direction of migration. To define this novel force-dependent role of vinculin in collective migration I took advantage of previously established shRNA based vinculin knock down Marin-Darby Canine Kidney (MDCK) epithelial cells.

The first part of my thesis comprises of my work demonstrating the mechanosensitive role of vinculin at AJ’s in collectively migrating cells. Using vinculin knockdown cells and vinculin mutants, which specifically disrupt vinculin’s ability to bind actin (VinI997A) or disrupt its ability to localize to AJs without affecting its localization at FAs (VinY822F), I establish a role of force across vinculin in E-cadherin internalization and clipping. Furthermore by measuring E-cadherin dynamics using fluorescence recovery after bleaching (FRAP) analysis I show that vinculin inhibition affects the turnover of E-cadherin at AJs. Together these data reveal a novel mechanosensitive role of vinculin in E-cadherin internalization and turnover in a migrating cell layer, which is contrary to the previously identified role of vinculin in potentiating E-cadherin junctions in a static monolayer.

For the last part of my thesis I designed a novel tension sensor to probe tension across N-cadherin (NTS). N-cadherin plays a critical role in cardiomyocytes, vascular smooth muscle cells, neurons and neural crest cells. Similar to E-cadherin, N-cadherin is also believed to bear tension and play a role in mechanotransduction pathways. To identify the role of tension across N-cadherin I designed a novel FRET-based molecular tension sensor for N-cadherin. I tested the ability of NTS to sense molecular tension in vascular smooth muscle cells, cardiomyocytes and cancer cells. Finally in collaboration with the Horwitz lab we have been able to show a role of tension across N-cadherin in synaptogenesis of neurons.

Dynamics of Delta/Notch signaling on endomesoderm segregation in the sea urchin embryo.

Endomesoderm is the common progenitor of endoderm and mesoderm early in the development of many animals. In the sea urchin embryo, the Delta/Notch pathway is necessary for the diversification of this tissue, as are two early transcription factors, Gcm and FoxA, which are expressed in mesoderm and endoderm, respectively. Here, we provide a detailed lineage analysis of the cleavages leading to endomesoderm segregation, and examine the expression patterns and the regulatory relationships of three known regulators of this cell fate dichotomy in the context of the lineages. We observed that endomesoderm segregation first occurs at hatched blastula stage. Prior to this stage, Gcm and FoxA are co-expressed in the same cells, whereas at hatching these genes are detected in two distinct cell populations. Gcm remains expressed in the most vegetal endomesoderm descendant cells, while FoxA is downregulated in those cells and activated in the above neighboring cells. Initially, Delta is expressed exclusively in the micromeres, where it is necessary for the most vegetal endomesoderm cell descendants to express Gcm and become mesoderm. Our experiments show a requirement for a continuous Delta input for more than two cleavages (or about 2.5 hours) before Gcm expression continues in those cells independently of further Delta input. Thus, this study provides new insights into the timing mechanisms and the molecular dynamics of endomesoderm segregation during sea urchin embryogenesis and into the mode of action of the Delta/Notch pathway in mediating mesoderm fate.

Effect of different jump distributions on the dynamics of jump processes

The paper investigates stochastic processes forced by independent and identically distributed jumps occurring according to a Poisson process. The impact of different distributions of the jump amplitudes are analyzed for processes with linear drift. Exact expressions of the probability density functions are derived when jump amplitudes are distributed as exponential, gamma, and mixture of exponential distributions for both natural and reflecting boundary conditions. The mean level-crossing properties are studied in relation to the different jump amplitudes. As an example of application of the previous theoretical derivations, the role of different rainfall-depth distributions on an existing stochastic soil water balance model is analyzed. It is shown how the shape of distribution of daily rainfall depths plays a more relevant role on the soil moisture probability distribution as the rainfall frequency decreases, as predicted by future climatic scenarios. © 2010 The American Physical Society.

Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses--using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

Long-term dynamics of death rates of emphysema, asthma, and pneumonia and improving air quality.

BACKGROUND: The respiratory tract is a major target of exposure to air pollutants, and respiratory diseases are associated with both short- and long-term exposures. We hypothesized that improved air quality in North Carolina was associated with reduced rates of death from respiratory diseases in local populations. MATERIALS AND METHODS: We analyzed the trends of emphysema, asthma, and pneumonia mortality and changes of the levels of ozone, sulfur dioxide (SO2), nitrogen dioxide (NO2), carbon monoxide (CO), and particulate matters (PM2.5 and PM10) using monthly data measurements from air-monitoring stations in North Carolina in 1993-2010. The log-linear model was used to evaluate associations between air-pollutant levels and age-adjusted death rates (per 100,000 of population) calculated for 5-year age-groups and for standard 2000 North Carolina population. The studied associations were adjusted by age group-specific smoking prevalence and seasonal fluctuations of disease-specific respiratory deaths. RESULTS: Decline in emphysema deaths was associated with decreasing levels of SO2 and CO in the air, decline in asthma deaths-with lower SO2, CO, and PM10 levels, and decline in pneumonia deaths-with lower levels of SO2. Sensitivity analyses were performed to study potential effects of the change from International Classification of Diseases (ICD)-9 to ICD-10 codes, the effects of air pollutants on mortality during summer and winter, the impact of approach when only the underlying causes of deaths were used, and when mortality and air-quality data were analyzed on the county level. In each case, the results of sensitivity analyses demonstrated stability. The importance of analysis of pneumonia as an underlying cause of death was also highlighted. CONCLUSION: Significant associations were observed between decreasing death rates of emphysema, asthma, and pneumonia and decreases in levels of ambient air pollutants in North Carolina.

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.

The spatiotemporal dynamics of autobiographical memory: neural correlates of recall, emotional intensity, and reliving.

We sought to map the time course of autobiographical memory retrieval, including brain regions that mediate phenomenological experiences of reliving and emotional intensity. Participants recalled personal memories to auditory word cues during event-related functional magnetic resonance imaging (fMRI). Participants pressed a button when a memory was accessed, maintained and elaborated the memory, and then gave subjective ratings of emotion and reliving. A novel fMRI approach based on timing differences capitalized on the protracted reconstructive process of autobiographical memory to segregate brain areas contributing to initial access and later elaboration and maintenance of episodic memories. The initial period engaged hippocampal, retrosplenial, and medial and right prefrontal activity, whereas the later period recruited visual, precuneus, and left prefrontal activity. Emotional intensity ratings were correlated with activity in several regions, including the amygdala and the hippocampus during the initial period. Reliving ratings were correlated with activity in visual cortex and ventromedial and inferior prefrontal regions during the later period. Frontopolar cortex was the only brain region sensitive to emotional intensity across both periods. Results were confirmed by time-locked averages of the fMRI signal. The findings indicate dynamic recruitment of emotion-, memory-, and sensory-related brain regions during remembering and their dissociable contributions to phenomenological features of the memories.

The impact of host immune status on the within-host and population dynamics of antigenic immune escape.

Antigenically evolving pathogens such as influenza viruses are difficult to control owing to their ability to evade host immunity by producing immune escape variants. Experimental studies have repeatedly demonstrated that viral immune escape variants emerge more often from immunized hosts than from naive hosts. This empirical relationship between host immune status and within-host immune escape is not fully understood theoretically, nor has its impact on antigenic evolution at the population level been evaluated. Here, we show that this relationship can be understood as a trade-off between the probability that a new antigenic variant is produced and the level of viraemia it reaches within a host. Scaling up this intra-host level trade-off to a simple population level model, we obtain a distribution for variant persistence times that is consistent with influenza A/H3N2 antigenic variant data. At the within-host level, our results show that target cell limitation, or a functional equivalent, provides a parsimonious explanation for how host immune status drives the generation of immune escape mutants. At the population level, our analysis also offers an alternative explanation for the observed tempo of antigenic evolution, namely that the production rate of immune escape variants is driven by the accumulation of herd immunity. Overall, our results suggest that disease control strategies should be further assessed by considering the impact that increased immunity--through vaccination--has on the production of new antigenic variants.

In vitro fluid dynamics of the Ahmed glaucoma valve modified with expanded polytetrafluoroethylene.

PURPOSE: Long-term intraocular pressure reduction by glaucoma drainage devices (GDDs) is often limited by the fibrotic capsule that forms around them. Prior work demonstrates that modifying a GDD with a porous membrane promotes a vascularized and more permeable capsule. This work examines the in vitro fluid dynamics of the Ahmed valve after enclosing the outflow tract with a porous membrane of expanded polytetrafluoroethylene (ePTFE). MATERIALS AND METHODS: The control and modified Ahmed implants (termed porous retrofitted implant with modified enclosure or PRIME-Ahmed) were submerged in saline and gelatin and perfused in a system that monitored flow (Q) and pressure (P). Flow rates of 1-50 μl/min were applied and steady state pressure recorded. Resistance was calculated by dividing pressure by flow. RESULTS: Modifying the Ahmed valve implant outflow with expanded ePTFE increased pressure and resistance. Pressure at a flow of 2 μl/min was increased in the PRIME-Ahmed (11.6 ± 1.5 mm Hg) relative to the control implant (6.5 ± 1.2 mm Hg). Resistance at a flow of 2 μl/min was increased in the PRIME-Ahmed (5.8 ± 0.8 mm Hg/μl/min) when compared to the control implant (3.2 ± 0.6 mm Hg/μl/min). CONCLUSIONS: Modifying the outflow tract of the Ahmed valve with a porous membrane adds resistance that decreases with increasing flow. The Ahmed valve implant behaves as a variable resistor. It is partially open at low pressures and provides reduced resistance at physiologic flow rates.

The East Asian Collection at Duke University: Dynamics of Change.

Traces the history of Duke's East Asian Studies program and associated library collections from the beginning of the twentieth century to the present. Describes the strengths of the Japanese, Chinese and Korean collections, materials in special collections and cooperation with the University of North Carolina.

The neural dynamics of stimulus and response conflict processing as a function of response complexity and task demands.

Both stimulus and response conflict can disrupt behavior by slowing response times and decreasing accuracy. Although several neural activations have been associated with conflict processing, it is unclear how specific any of these are to the type of stimulus conflict or the amount of response conflict. Here, we recorded electrical brain activity, while manipulating the type of stimulus conflict in the task (spatial [Flanker] versus semantic [Stroop]) and the amount of response conflict (two versus four response choices). Behaviorally, responses were slower to incongruent versus congruent stimuli across all task and response types, along with overall slowing for higher response-mapping complexity. The earliest incongruency-related neural effect was a short-duration frontally-distributed negativity at ~200 ms that was only present in the Flanker spatial-conflict task. At longer latencies, the classic fronto-central incongruency-related negativity 'N(inc)' was observed for all conditions, but was larger and ~100 ms longer in duration with more response options. Further, the onset of the motor-related lateralized readiness potential (LRP) was earlier for the two vs. four response sets, indicating that smaller response sets enabled faster motor-response preparation. The late positive complex (LPC) was present in all conditions except the two-response Stroop task, suggesting this late conflict-related activity is not specifically related to task type or response-mapping complexity. Importantly, across tasks and conditions, the LRP onset at or before the conflict-related N(inc), indicating that motor preparation is a rapid, automatic process that interacts with the conflict-detection processes after it has begun. Together, these data highlight how different conflict-related processes operate in parallel and depend on both the cognitive demands of the task and the number of response options.