Detecting separate time scales in genetic expression data.

BACKGROUND: Biological processes occur on a vast range of time scales, and many of them occur concurrently. As a result, system-wide measurements of gene expression have the potential to capture many of these processes simultaneously. The challenge however, is to separate these processes and time scales in the data. In many cases the number of processes and their time scales is unknown. This issue is particularly relevant to developmental biologists, who are interested in processes such as growth, segmentation and differentiation, which can all take place simultaneously, but on different time scales. RESULTS: We introduce a flexible and statistically rigorous method for detecting different time scales in time-series gene expression data, by identifying expression patterns that are temporally shifted between replicate datasets. We apply our approach to a Saccharomyces cerevisiae cell-cycle dataset and an Arabidopsis thaliana root developmental dataset. In both datasets our method successfully detects processes operating on several different time scales. Furthermore we show that many of these time scales can be associated with particular biological functions. CONCLUSIONS: The spatiotemporal modules identified by our method suggest the presence of multiple biological processes, acting at distinct time scales in both the Arabidopsis root and yeast. Using similar large-scale expression datasets, the identification of biological processes acting at multiple time scales in many organisms is now possible.

Dramatic variability of the carbonate system at a temperate coastal ocean site (Beaufort, North Carolina, USA) is regulated by physical and biogeochemical processes on multiple timescales.

Increasing atmospheric carbon dioxide (CO2) from anthropogenic sources is acidifying marine environments resulting in potentially dramatic consequences for the physical, chemical and biological functioning of these ecosystems. If current trends continue, mean ocean pH is expected to decrease by ~0.2 units over the next ~50 years. Yet, there is also substantial temporal variability in pH and other carbon system parameters in the ocean resulting in regions that already experience change that exceeds long-term projected trends in pH. This points to short-term dynamics as an important layer of complexity on top of long-term trends. Thus, in order to predict future climate change impacts, there is a critical need to characterize the natural range and dynamics of the marine carbonate system and the mechanisms responsible for observed variability. Here, we present pH and dissolved inorganic carbon (DIC) at time intervals spanning 1 hour to >1 year from a dynamic, coastal, temperate marine system (Beaufort Inlet, Beaufort NC USA) to characterize the carbonate system at multiple time scales. Daily and seasonal variation of the carbonate system is largely driven by temperature, alkalinity and the balance between primary production and respiration, but high frequency change (hours to days) is further influenced by water mass movement (e.g. tides) and stochastic events (e.g. storms). Both annual (~0.3 units) and diurnal (~0.1 units) variability in coastal ocean acidity are similar in magnitude to 50 year projections of ocean acidity associated with increasing atmospheric CO2. The environmental variables driving these changes highlight the importance of characterizing the complete carbonate system rather than just pH. Short-term dynamics of ocean carbon parameters may already exert significant pressure on some coastal marine ecosystems with implications for ecology, biogeochemistry and evolution and this shorter term variability layers additive effects and complexity, including extreme values, on top of long-term trends in ocean acidification.

Mechanism of Shear Thickening in Reversibly Cross-linked Supramolecular Polymer Networks.

We report here the nonlinear rheological properties of metallo-supramolecular networks formed by the reversible cross-linking of semi-dilute unentangled solutions of poly(4-vinylpyridine) (PVP) in dimethyl sulfoxide (DMSO). The reversible cross-linkers are bis-Pd(II) or bis-Pt(II) complexes that coordinate to the pyridine functional groups on the PVP. Under steady shear, shear thickening is observed above a critical shear rate, and that critical shear rate is experimentally correlated with the lifetime of the metal-ligand bond. The onset and magnitude of the shear thickening depend on the amount of cross-linkers added. In contrast to the behavior observed in most transient networks, the time scale of network relaxation is found to increase during shear thickening. The primary mechanism of shear thickening is ascribed to the shear-induced transformation of intrachain cross-linking to interchain cross-linking, rather than nonlinear high tension along polymer chains that are stretched beyond the Gaussian range.

Deformation of stem cell nuclei by nanotopographical cues.

Cells sense cues in their surrounding microenvironment. These cues are converted into intracellular signals and transduced to the nucleus in order for the cell to respond and adapt its function. Within the nucleus, structural changes occur that ultimately lead to changes in the gene expression. In this study, we explore the structural changes of the nucleus of human mesenchymal stem cells as an effect of topographical cues. We use a controlled nanotopography to drive shape changes to the cell nucleus, and measure the changes with both fluorescence microscopy and a novel light scattering technique. The nucleus changes shape dramatically in response to the nanotopography, and in a manner dependent on the mechanical properties of the substrate. The kinetics of the nuclear deformation follows an unexpected trajectory. As opposed to a gradual shape change in response to the topography, once the cytoskeleton attains an aligned and elongation morphology on the time scale of several hours, the nucleus changes shape rapidly and intensely.

Bayesian Analysis of Latent Threshold Dynamic Models

We discuss a general approach to dynamic sparsity modeling in multivariate time series analysis. Time-varying parameters are linked to latent processes that are thresholded to induce zero values adaptively, providing natural mechanisms for dynamic variable inclusion/selection. We discuss Bayesian model specification, analysis and prediction in dynamic regressions, time-varying vector autoregressions, and multivariate volatility models using latent thresholding. Application to a topical macroeconomic time series problem illustrates some of the benefits of the approach in terms of statistical and economic interpretations as well as improved predictions. Supplementary materials for this article are available online. © 2013 Copyright Taylor and Francis Group, LLC.

The evolution of airplanes

The prevailing view is that we cannot witness biological evolution because it occurred on a time scale immensely greater than our lifetime. Here, we show that we can witness evolution in our lifetime by watching the evolution of the flying human-and-machine species: the airplane. We document this evolution, and we also predict it based on a physics principle: the constructal law. We show that the airplanes must obey theoretical allometric rules that unite them with the birds and other animals. For example, the larger airplanes are faster, more efficient as vehicles, and have greater range. The engine mass is proportional to the body size: this scaling is analogous to animal design, where the mass of the motive organs (muscle, heart, lung) is proportional to the body size. Large or small, airplanes exhibit a proportionality between wing span and fuselage length, and between fuel load and body size. The animal-design counterparts of these features are evident. The view that emerges is that the evolution phenomenon is broader than biological evolution. The evolution of technology, river basins, and animal design is one phenomenon, and it belongs in physics. © 2014 AIP Publishing LLC.

The Ontogeny of Mucosal and Systemic Antibody Responses to HIV-1 Infection

The humoral immune system plays a critical role in the clearance of numerous pathogens. In the setting of HIV-1 infection, the virus infects, integrates its genome into the host's cells, replicates, and establishes a reservoir of virus-infected cells. The initial antibody response to HIV-1 infection is targeted to non-neutralizing epitopes on HIV-1 Env gp41, and when a neutralizing response does develop months after transmission, it is specific for the autologous founder virus and the virus escapes rapidly. After continuous waves of antibody mediated neutralization and viral escape, a small subset of infected individuals eventually develop broad and potent heterologous neutralizing antibodies years after infection. In this dissertation, I have studied the ontogeny of mucosal and systemic antibody responses to HIV-1 infection by means of three distinct aims: 1. Determine the origin of the initial antibody response to HIV-1 infection. 2. Characterize the role of restricted VH and VL gene segment usage in shaping the antibody response to HIV-1 infection. 3. Determine the role of persistence of B cell clonal lineages in shaping the mutation frequencies of HIV-1 reactive antibodies.

After the introduction (Chapter 1) and methods (Chapter 2), Chapter 3 of this dissertation describes a study of the antibody response of terminal ileum B cells to HIV-1 envelope (Env) in early and chronic HIV-1 infection and provides evidence for the role of environmental antigens in shaping the repertoire of B cells that respond to HIV-1 infection. Previous work by Liao et al. demonstrated that the initial plasma cell response in the blood to acute HIV-1 infection is to gp41 and is derived from a polyreactive memory B cell pool. Many of these antibodies cross-reacted with commensal bacteria, Therefore, in Chapter 3, the relationship of intestinal B cell reactivity with commensal bacteria to HIV-1 infection-induced antibody response was probed using single B cell sorting, reverse transcription and nested polymerase chain reaction (RT- PCR) methods, and recombinant antibody technology. The dominant B cell response in the terminal ileum was to HIV-1 envelope (Env) gp41, and 82% of gp41- reactive antibodies cross-reacted with commensal bacteria whole cell lysates. Pyrosequencing of blood B cells revealed HIV-1 antibody clonal lineages shared between ileum and blood. Mutated IgG antibodies cross-reactive with both Env gp41 and commensal bacteria could also be isolated from the terminal ileum of HIV-1 uninfected individuals. Thus, the antibody response to HIV-1 can be shaped by intestinal B cells stimulated by commensal bacteria prior to HIV-1 infection to develop a pre-infection pool of memory B cells cross-reactive with HIV-1 gp41.

Chapter 4 details the study of restricted VH and VL gene segment usage for gp41 and gp120 antibody induction following acute HIV-1 infection; mutations in gp41 lead to virus enhanced neutralization sensitivity. The B cell repertoire of antibodies induced in a HIV-1 infected African individual, CAP206, who developed broadly neutralizing antibodies (bnAbs) directed to the HIV-1 envelope gp41 membrane proximal external region (MPER), is characterized. Understanding the selection of virus mutants by neutralizing antibodies is critical to understanding the role of antibodies in control of HIV-1 replication and prevention from HIV-1 infection. Previously, an MPER neutralizing antibody, CAP206-CH12, with the binding footprint identical to that of MPER broadly neutralizing antibody 4E10, that like 4E10 utilized the VH1-69 and VK3-20 variable gene segments was isolated from this individual (Morris et al., 2011). Using single B cell sorting, RT- PCR methods, and recombinant antibody technology, Chapter 4 describes the isolation of a VH1-69, Vk3-20 glycan-dependent clonal lineage from CAP206, targeted to gp120, that has the property of neutralizing a neutralization sensitive CAP206 transmitted/founder (T/F) and heterologous viruses with mutations at amino acids 680 or 681 in the MPER 4E10/CH12 binding site. These data demonstrate sites within the MPER bnAb epitope (aa 680-681) in which mutations can be selected that lead to viruses with enhanced sensitivity to autologous and heterologous neutralizing antibodies.

In Chapter 5, I have completed a comparison of evolution of B cell clonal lineages in two HIV-1 infected individuals who have a predominant VH1-69 response to HIV-1 infection--one who produces broadly neutralizing MPER-reactive mAbs and one who does not. Autologous neutralization in the plasma takes ~12 weeks to develop (Gray et al., 2007; Tomaras et al., 2008b). Only a small subset of HIV-1 infected individuals develops high plasma levels of broad and potent heterologous neutralization, and when it does occur, it typically takes 3-4 years to develop (Euler et al., 2010; Gray et al., 2007; 2011; Tomaras et al., 2011). The HIV-1 bnAbs that have been isolated to date have a number of unusual characteristics including, autoreactivity and high levels of somatic hypermutations, which are typically tightly regulated by immune control mechanisms (Haynes et al., 2005; 2012b; Kwong and Mascola, 2012; Scheid et al., 2009a). The VH mutation frequencies of bnAbs average ~15% but have been shown to be as high as 32% (reviewed in Mascola and Haynes, 2013; Kwong and Mascola, 2012). The high frequency of somatic hypermutations suggests that the B cell clonal lineages that eventually produce bnAbs undergo high-levels of affinity maturation, implying prolonged germinal center (GC) reactions and high levels of T cell help. To study the duration of HIV-1- reactive B cell clonal persistence, HIV-1 reactive and non HIV-1- reactive B cell clonal lineages were isolated from an HIV-1 infected individual that produces bnAbs, CAP206, and an HIV-1 infected individual who does not produce bnAbs, 004-0. Single B cell sorting, RT-PCR and recombinant antibody technology was used to isolate and produce monoclonal antibodies from multiple time points from each individual. B cell sequences clonally related to mAbs isolated by single cell PCR were identified within pyrosequences of longitudinal samples of these two individuals. Both individuals produced long-lived B cell clones that persisted from 0-232 weeks in CAP206, and 0-238 weeks in 004-0. The average length of persistence of clones containing members isolated from two separate time points was 91.5 weeks both individuals. Examples of the continued evolution of clonal lineages were observed in both the bnAb and non-bnAb individual. These data indicated that the ability to generate persistent and evolving B cell clonal lineages occurs in both bnAb and non-bnAb individuals, suggesting that some alternative host or viral factor is critical for the generation of highly mutated broadly neutralizing antibodies.

Together the studies described in Chapter 3-5 show that multiple factors influence the antibody response to HIV-1 infection. The initial antibody response to HIV-1 Env gp41 can be shaped by a B cell response to intestinal commensal bacteria prior to HIV-1 infection. VH and VL gene segment restriction can impact the B cell response to multiple HIV-1 antigens, and virus escape mutations in the MPER can confer enhanced neutralization sensitivity to autologous and heterologous antibodies. Finally, the ability to generate long-lived HIV-1 clonal lineages in and of itself does not confer on the host the ability to produce bnAbs.

Later Life Consequences of Developmental Mitochondrial DNA Damage in C. elegans

Mitochondria are responsible for producing the vast majority of cellular ATP, and are therefore critical to organismal health [1]. They contain thir own genomes (mtDNA) which encode 13 proteins that are all subunits of the mitochondrial respiratory chain (MRC) and are essential for oxidative phosphorylation [2]. mtDNA is present in multiple copies per cell, usually between 103 and 104 , though this number is reduced during certain developmental stages [3, 4]. The health of the mitochondrial genome is also important to the health of the organism, as mutations in mtDNA lead to human diseases that collectively affect approximately 1 in 4000 people [5, 6]. mtDNA is more susceptible than nuclear DNA (nucDNA) to damage by many environmental pollutants, for reasons including the absence of Nucleotide Excision Repair (NER) in the mitochondria [7]. NER is a highly functionally conserved DNA repair pathway that removes bulky, helix distorting lesions such as those caused by ultraviolet C (UVC) radiation and also many environmental toxicants, including benzo[a]pyrene (BaP) [8]. While these lesions cannot be repaired, they are slowly removed through a process that involves mitochondrial dynamics and autophagy [9, 10]. However, when present during development in C. elegans, this damage reduces mtDNA copy number and ATP levels [11]. We hypothesize that this damage, when present during development, will result in mitochondrial dysfunction and increase the potential for adverse outcomes later in life.

To test this hypothesis, 1st larval stage (L1) C. elegans are exposed to 3 doses of 7.5J/m2 ultraviolet C radiation 24 hours apart, leading to the accumulation of mtDNA damage [9, 11]. After exposure, many mitochondrial endpoints are assessed at multiple time points later in life. mtDNA and nucDNA damage levels and genome copy numbers are measured via QPCR and real-time PCR , respectively, every 2 day for 10 days. Steady state ATP levels are measured via luciferase expressing reporter strains and traditional ATP extraction methods. Oxygen consumption is measured using a Seahorse XFe24 extra cellular flux analyzer. Gene expression changes are measured via real time PCR and targeted metabolomics via LC-MS are used to investigate changes in organic acid, amino acid and acyl-carnitine levels. Lastly, nematode developmental delay is assessed as growth, and measured via imaging and COPAS biosort.

I have found that despite being removed, UVC induced mtDNA damage during development leads to persistent deficits in energy production later in life. mtDNA copy number is permanently reduced, as are ATP levels, though oxygen consumption is increased, indicating inefficient or uncoupled respiration. Metabolomic data and mutant sensitivity indicate a role for NADPH and oxidative stress in these results, and exposed nematodes are more sensitive to the mitochondrial poison rotenone later in life. These results fit with the developmental origin of health and disease hypothesis, and show the potential for environmental exposures to have lasting effects on mitochondrial function.

Lastly, we are currently working to investigate the potential for irreparable mtDNA lesions to drive mutagenesis in mtDNA. Mutations in mtDNA lead to a wide range of diseases, yet we currently do not understand the environmental component of what causes them. In vitro evidence suggests that UVC induced thymine dimers can be mutagenic [12]. We are using duplex sequencing of C. elegans mtDNA to determine mutation rates in nematodes exposed to our serial UVC protocol. Furthermore, by including mutant strains deficient in mitochondrial fission and mitophagy, we hope to determine if deficiencies in these processes will further increase mtDNA mutation rates, as they are implicated in human diseases.

Later Life Consequences of Subteratogenic Exposure to a Complex PAH Mixture in the Atlantic Killifish (Fundulus heteroclitus)

Subteratogenic and other low-level chronic exposures to toxicant mixtures are an understudied threat to environmental and human health. It is especially important to understand the effects of these exposures for contaminants, such as polycyclic aromatic hydrocarbons (PAHs) a large group of more than 100 individual compounds, which are important environmental (including aquatic) contaminants. Aquatic sediments constitute a major sink for hydrophobic pollutants, and studies show PAHs can persist in sediments over time. Furthermore, estuarine systems (namely breeding grounds) are of particular concern, as they are highly impacted by a wide variety of pollutants, and estuarine fishes are often exposed to some of the highest levels of contaminants of any vertebrate taxon. Acute embryonic exposure to PAHs results in cardiac teratogenesis in fish, and early life exposure to certain individual PAHs and PAH mixtures cause heart alterations with decreased swimming capacity in adult fish. Consequently, the heart and cardiorespiratory system are thought to be targets of PAH mixture exposure. While many studies have investigated acute, teratogenic PAH exposures, few studies have longitudinally examined the impacts of subtle, subteratogenic PAH mixture exposures, which are arguably more broadly applicable to environmental contamination scenarios. The goal of this dissertation was to highlight the later-life consequences of early-life exposure to subteratogenic concentrations of a complex, environmentally relevant PAH mixture.

A unique population of Fundulus heteroclitus (the Atlantic killifish or mummichog, hereafter referred to as killifish), has adapted to creosote-based polycyclic aromatic hydrocarbons (PAHs) found at the Atlantic Wood Industries (AW) Superfund site in the southern branch of the Elizabeth River, VA, USA. This killifish population survives in a site heavily contaminated with a mixture of PAHs from former creosote operations. They have developed resistance to the acute toxicity and teratogenic effects caused by the mixture of PAHs in sediment from the site. The primary goal of this dissertation was to compare and contrast later-life outcomes of early-life, subteratogenic PAH mixture exposure in both the Atlantic Wood killifish (AW) and a naïve reference population of killifish from King’s Creek (KC; a relatively uncontaminated tributary of the Severn River, VA). Killifish from both populations were exposed to subteratogenic concentrations of a complex PAH-sediment extract, Elizabeth River Sediment Extract (ERSE), made by collecting sediment from the AW site. Fish were reared over a 5-month period in the laboratory, during which they were examined for a variety of molecular, physiological and behavioral responses.

The central aims of my dissertation were to determine alterations to embryonic gene expression, larval swimming activity, adult behavior, heart structure, enzyme activity, and swimming/cardiorespiratory performance following subteratogenic exposure to ERSE. I hypothesized that subteratogenic exposure to ERSE would impair cardiac ontogenic processes in a way that would be detectable via gene expression in embryos, and that the misregulation of cardiac genes would help to explain activity changes, behavioral deficits, and later-life swimming deficiencies. I also hypothesized that fish heart structure would be altered. In addition, I hypothesized that the AW killifish population would be resistant to developmental exposures and perform normally in later life challenges. To investigate these hypotheses, a series of experiments were carried out in PAH-adapted killifish from Elizabeth River and in reference killifish. As an ancillary project to the primary aims of the dissertation, I examined the toxicity of weaker aryl hydrocarbon receptor (AHR) agonists in combination with fluoranthene (FL), an inhibitor of cytochrome P4501A1 (CYP1A1). This side project was conducted in both Danio rerio (zebrafish) and the KC and AW killifish.

Embryonic gene expression was measured in both killifish populations over an ERSE dose response with multiple time points (12, 24, 48, and 144 hours post exposure). Genes known to play critical roles in cardiac structure/development, cardiac function, and angiogenesis were elevated, indicating cardiac damage and activation of cardiovascular repair mechanisms. These data helped to inform later-life swimming performance and cardiac histology studies. Behavior was assessed during light and dark cycles in larvae of both populations following developmental exposure to ERSE. While KC killifish showed activity differences following exposure, AW killifish showed no significant changes even at concentrations that would cause overt cardiac toxicity in KC killifish. Juvenile behavior experiments demonstrated hyperactivity following ERSE exposure in KC killifish, but no significant behavioral changes in AW killifish. Adult swimming performance via prolonged critical swimming capacity (Ucrit) demonstrated performance costs in the AW killifish. Furthermore, swimming performance decline was observed in KC killifish following exposure to increasing dilutions of ERSE. Lastly, cardiac histology suggested that early-life exposure to ERSE could result in cardiac structural alteration and extravasation of blood into the pericardial cavity.

Responses to AHR agonists resulted in a ranking of relative potency for agonists, and determined which agonists, when combined with FL, caused cardiac teratogenesis. These experiments showed interesting species differences for zebrafish and killifish. To probe mechanisms responsible for cardiotoxicity, a CYP1A-morpholino and a AHR2-morpholino were used to mimic FL effects or attempt to rescue cardiac deformities respectively. Findings suggested that the cardiac toxicity elicited by weak agonist + FL exposure was likely driven by AHR-independent mechanisms. These studies stand in contrast to previous research from our lab showing that moderate AHR agonist + FL caused cardiac toxicity that can be partially rescued by AHR-morpholino knockdown.

My findings will form better characterization of mechanisms of PAH toxicity, and advance our understanding of how subteratogenic mixtures of PAHs exert their toxic action in naïve killifish. Furthermore, these studies will provide a framework for investigating how subteratogenic exposures to PAH mixtures can impact aquatic organismal health and performance. Most importantly, these experiments have the potential to help inform risk assessment in fish, mammals, and potentially humans. Ultimately, this research will help protect populations exposed to subtle PAH-contamination.

Activity in descending dopaminergic neurons represents but is not required for leg movements in the fruit fly Drosophila.

Modulatory descending neurons (DNs) that link the brain to body motor circuits, including dopaminergic DNs (DA-DNs), are thought to contribute to the flexible control of behavior. Dopamine elicits locomotor-like outputs and influences neuronal excitability in isolated body motor circuits over tens of seconds to minutes, but it remains unknown how and over what time scale DA-DN activity relates to movement in behaving animals. To address this question, we identified DA-DNs in the Drosophila brain and developed an electrophysiological preparation to record and manipulate the activity of these cells during behavior. We find that DA-DN spike rates are rapidly modulated during a subset of leg movements and scale with the total speed of ongoing leg movements, whether occurring spontaneously or in response to stimuli. However, activating DA-DNs does not elicit leg movements in intact flies, nor do acute bidirectional manipulations of DA-DN activity affect the probability or speed of leg movements over a time scale of seconds to minutes. Our findings indicate that in the context of intact descending control, changes in DA-DN activity are not sufficient to influence ongoing leg movements and open the door to studies investigating how these cells interact with other descending and local neuromodulatory inputs to influence body motor output.

Stress Relaxation for Granular Materials near Jamming under Cyclic Compression.

We have explored isotropically jammed states of semi-2D granular materials through cyclic compression. In each compression cycle, systems of either identical ellipses or bidisperse disks transition between jammed and unjammed states. We determine the evolution of the average pressure P and structure through consecutive jammed states. We observe a transition point ϕ_{m} above which P persists over many cycles; below ϕ_{m}, P relaxes slowly. The relaxation time scale associated with P increases with packing fraction, while the relaxation time scale for collective particle motion remains constant. The collective motion of the ellipses is hindered compared to disks because of the rotational constraints on elliptical particles.