Eco-evolutionary feedbacks in community and ecosystem ecology: interactions between the ecological theatre and the evolutionary play.

Interactions between natural selection and environmental change are well recognized and sit at the core of ecology and evolutionary biology. Reciprocal interactions between ecology and evolution, eco-evolutionary feedbacks, are less well studied, even though they may be critical for understanding the evolution of biological diversity, the structure of communities and the function of ecosystems. Eco-evolutionary feedbacks require that populations alter their environment (niche construction) and that those changes in the environment feed back to influence the subsequent evolution of the population. There is strong evidence that organisms influence their environment through predation, nutrient excretion and habitat modification, and that populations evolve in response to changes in their environment at time-scales congruent with ecological change (contemporary evolution). Here, we outline how the niche construction and contemporary evolution interact to alter the direction of evolution and the structure and function of communities and ecosystems. We then present five empirical systems that highlight important characteristics of eco-evolutionary feedbacks: rotifer-algae chemostats; alewife-zooplankton interactions in lakes; guppy life-history evolution and nutrient cycling in streams; avian seed predators and plants; and tree leaf chemistry and soil processes. The alewife-zooplankton system provides the most complete evidence for eco-evolutionary feedbacks, but other systems highlight the potential for eco-evolutionary feedbacks in a wide variety of natural systems.

Influence of Increased Human Presence in the Mills River Basin on Water Availability and Drought

Periods of drought and low streamflow can have profound impacts on both human and natural systems. People depend on a reliable source of water for numerous reasons including potable water supply and to produce economic value through agriculture or energy production. Aquatic ecosystems depend on water in addition to the economic benefits they provide to society through ecosystem services. Given that periods of low streamflow may become more extreme and frequent in the future, it is important to study the factors that control water availability during these times. In the absence of precipitation the slower hydrological response of groundwater systems will play an amplified role in water supply. Understanding the variability of the fraction of streamflow contribution from baseflow or groundwater during periods of drought provides insight into what future water availability may look like and how it can best be managed. The Mills River Basin in North Carolina is chosen as a case-study to test this understanding. First, obtaining a physically meaningful estimation of baseflow from USGS streamflow data via computerized hydrograph analysis techniques is carried out. Then applying a method of time series analysis including wavelet analysis can highlight signals of non-stationarity and evaluate the changes in variance required to better understand the natural variability of baseflow and low flows. In addition to natural variability, human influence must be taken into account in order to accurately assess how the combined system reacts to periods of low flow. Defining a combined demand that consists of both natural and human demand allows us to be more rigorous in assessing the level of sustainable use of a shared resource, in this case water. The analysis of baseflow variability can differ based on regional location and local hydrogeology, but it was found that baseflow varies from multiyear scales such as those associated with ENSO (3.5, 7 years) up to multi decadal time scales, but with most of the contributing variance coming from decadal or multiyear scales. It was also found that the behavior of baseflow and subsequently water availability depends a great deal on overall precipitation, the tracks of hurricanes or tropical storms and associated climate indices, as well as physiography and hydrogeology. Evaluating and utilizing the Duke Combined Hydrology Model (DCHM), reasonably accurate estimates of streamflow during periods of low flow were obtained in part due to the model’s ability to capture subsurface processes. Being able to accurately simulate streamflow levels and subsurface interactions during periods of drought can be very valuable to water suppliers, decision makers, and ultimately impact citizens. Knowledge of future droughts and periods of low flow in addition to tracking customer demand will allow for better management practices on the part of water suppliers such as knowing when they should withdraw more water during a surplus so that the level of stress on the system is minimized when there is not ample water supply.

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

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

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

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


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

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

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

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

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

Operationalizing the social-ecological systems framework to assess sustainability.

Environmental governance is more effective when the scales of ecological processes are well matched with the human institutions charged with managing human-environment interactions. The social-ecological systems (SESs) framework provides guidance on how to assess the social and ecological dimensions that contribute to sustainable resource use and management, but rarely if ever has been operationalized for multiple localities in a spatially explicit, quantitative manner. Here, we use the case of small-scale fisheries in Baja California Sur, Mexico, to identify distinct SES regions and test key aspects of coupled SESs theory. Regions that exhibit greater potential for social-ecological sustainability in one dimension do not necessarily exhibit it in others, highlighting the importance of integrative, coupled system analyses when implementing spatial planning and other ecosystem-based strategies.

Constructal Design of Energy Systems

This dissertation shows the use of Constructal law to find the relation between the morphing of the system configuration and the improvements in the global performance of the complex flow system. It shows that the better features of both flow and heat transfer architecture can be found and predicted by using the constructal law in energy systems. Chapter 2 shows the effect of flow configuration on the heat transfer performance of a spiral shaped pipe embedded in a cylindrical conducting volume. Several configurations were considered. The optimal spacings between the spiral turns and spire planes exist, such that the volumetric heat transfer rate is maximal. The optimized features of the heat transfer architecture are robust. Chapter 3 shows the heat transfer performance of a helically shaped pipe embedded in a cylindrical conducting volume. It shows that the optimized features of the heat transfer architecture are robust with respect to changes in several physical parameters. Chapter 4 reports analytically the formulas for effective permeability in several configurations of fissured systems, using the closed-form description of tree networks designed to provide flow access. The permeability formulas do not vary much from one tree design to the next, suggesting that similar formulas may apply to naturally fissured porous media with unknown precise details, which occur in natural reservoirs. Chapter 5 illustrates a counterflow heat exchanger consists of two plenums with a core. The results show that the overall flow and thermal resistance are lowest when the core is absent. Overall, the constructal design governs the evolution of flow configuration in nature and energy systems.

Effects of Natural Organic Matter on the Dissolution Kinetics and Bioavailability of Metal Oxide Nanoparticles

The rapid development of nanotechnology and wider applications of engineered nanomaterials (ENMs) in the last few decades have generated concerns regarding their environmental and health risks. After release into the environment, ENMs undergo aggregation, transformation, and, for metal-based nanomaterials, dissolution processes, which together determine their fate, bioavailability and toxicity to living organisms in the ecosystems. The rates of these processes are dependent on nanomaterial characteristics as well as complex environmental factors, including natural organic matter (NOM). As a ubiquitous component of aquatic systems, NOM plays a key role in the aggregation, dissolution and transformation of metal-based nanomaterials and colloids in aquatic environments.

The goal of this dissertation work is to investigate how NOM fractions with different chemical and molecular properties affect the dissolution kinetics of metal oxide ENMs, such as zinc oxide (ZnO) and copper oxide (CuO) nanoparticles (NPs), and consequently their bioavailability to aquatic vertebrate, with Gulf killifish (Fundulus grandis) embryos as model organisms.

ZnO NPs are known to dissolve at relatively fast rates, and the rate of dissolution is influenced by water chemistry, including the presence of Zn-chelating ligands. A challenge, however, remains in quantifying the dissolution of ZnO NPs, particularly for time scales that are short enough to determine rates. This dissertation assessed the application of anodic stripping voltammetry (ASV) with a hanging mercury drop electrode to directly measure the concentration of dissolved Zn in ZnO NP suspensions, without separation of the ZnO NPs from the aqueous phase. Dissolved zinc concentration measured by ASV ([Zn]ASV) was compared with that measured by inductively coupled plasma mass spectrometry (ICP-MS) after ultracentrifugation ([Zn]ICP-MS), for four types of ZnO NPs with different coatings and primary particle diameters. For small ZnO NPs (4-5 nm), [Zn]ASV was 20% higher than [Zn]ICP-MS, suggesting that these small NPs contributed to the voltammetric measurement. For larger ZnO NPs (approximately 20 nm), [Zn]ASV was (79±19)% of [Zn]ICP-MS, despite the high concentrations of ZnO NPs in suspension, suggesting that ASV can be used to accurately measure the dissolution kinetics of ZnO NPs of this primary particle size.

Using the ASV technique to directly measure dissolved zinc concentration, we examined the effects of 16 different NOM isolates on the dissolution kinetics of ZnO NPs in buffered potassium chloride solution. The observed dissolution rate constants (kobs) and dissolved zinc concentrations at equilibrium increased linearly with NOM concentration (from 0 to 40 mg-C L-1) for Suwannee River humic acid (SRHA), Suwannee River fulvic acid and Pony Lake fulvic acid. When dissolution rates were compared for the 16 NOM isolates, kobs was positively correlated with certain properties of NOM, including specific ultraviolet absorbance (SUVA), aromatic and carbonyl carbon contents, and molecular weight. Dissolution rate constants were negatively correlated to hydrogen/carbon ratio and aliphatic carbon content. The observed correlations indicate that aromatic carbon content is a key factor in determining the rate of NOM-promoted dissolution of ZnO NPs. NOM isolates with higher SUVA were also more effective at enhancing the colloidal stability of the NPs; however, the NOM-promoted dissolution was likely due to enhanced interactions between surface metal ions and NOM rather than smaller aggregate size.

Based on the above results, we designed experiments to quantitatively link the dissolution kinetics and bioavailability of CuO NPs to Gulf killifish embryos under the influence of NOM. The CuO NPs dissolved to varying degrees and at different rates in diluted 5‰ artificial seawater buffered to different pH (6.3-7.5), with or without selected NOM isolates at various concentrations (0.1-10 mg-C L-1). NOM isolates with higher SUVA and aromatic carbon content (such as SRHA) were more effective at promoting the dissolution of CuO NPs, as with ZnO NPs, especially at higher NOM concentrations. On the other hand, the presence of NOM decreased the bioavailability of dissolved Cu ions, with the uptake rate constant negatively correlated to dissolved organic carbon concentration ([DOC]) multiplied by SUVA, a combined parameter indicative of aromatic carbon concentration in the media. When the embryos were exposed to CuO NP suspension, changes in their Cu content were due to the uptake of both dissolved Cu ions and nanoparticulate CuO. The uptake rate constant of nanoparticulate CuO was also negatively correlated to [DOC]×SUVA, in a fashion roughly proportional to changes in dissolved Cu uptake rate constant. Thus, the ratio of uptake rate constants from dissolved Cu and nanoparticulate CuO (ranging from 12 to 22, on average 17±4) were insensitive to NOM type or concentration. Instead, the relative contributions of these two Cu forms were largely determined by the percentage of CuO NP that was dissolved.

Overall, this dissertation elucidated the important role that dissolved NOM plays in affecting the environmental fate and bioavailability of soluble metal-based nanomaterials. This dissertation work identified aromatic carbon content and its indicator SUVA as key NOM properties that influence the dissolution, aggregation and biouptake kinetics of metal oxide NPs and highlighted dissolution rate as a useful functional assay for assessing the relative contributions of dissolved and nanoparticulate forms to metal bioavailability. Findings of this dissertation work will be helpful for predicting the environmental risks of engineered nanomaterials.

Biological and Physical Factors Affecting the Natural History and Evolution of Encapsulated Development

The evolution of reproductive strategies involves a complex calculus of costs and benefits to both parents and offspring. Many marine animals produce embryos packaged in tough egg capsules or gelatinous egg masses attached to benthic surfaces. While these egg structures can protect against environmental stresses, the packaging is energetically costly for parents to produce. In this series of studies, I examined a variety of ecological factors affecting the evolution of benthic development as a life history strategy. I used marine gastropods as my model system because they are incredibly diverse and abundant worldwide, and they exhibit a variety of reproductive and developmental strategies.

The first study examines predation on benthic egg masses. I investigated: 1) behavioral mechanisms of predation when embryos are targeted (rather than the whole egg mass); 2) the specific role of gelatinous matrix in predation. I hypothesized that gelatinous matrix does not facilitate predation. One study system was the sea slug Olea hansineensis, an obligate egg mass predator, feeding on the sea slug Haminoea vesicula. Olea fed intensely and efficiently on individual Haminoea embryos inside egg masses but showed no response to live embryos removed from gel, suggesting that gelatinous matrix enables predation. This may be due to mechanical support of the feeding predator by the matrix. However, Haminoea egg masses outnumber Olea by two orders of magnitude in the field, and each egg mass can contain many tens of thousands of embryos, so predation pressure on individuals is likely not strong. The second system involved the snail Nassarius vibex, a non-obligate egg mass predator, feeding on the polychaete worm Clymenella mucosa. Gel neither inhibits nor promotes embryo predation for Nassarius, but because it cannot target individual embryos inside an egg mass, its feeding is slow and inefficient, and feeding rates in the field are quite low. However, snails that compete with Nassarius for scavenged food have not been seen to eat egg masses in the field, leaving Nassarius free to exploit the resource. Overall, egg mass predation in these two systems likely benefits the predators much more than it negatively affects the prey. Thus, selection for environmentally protective aspects of egg mass production may be much stronger than selection for defense against predation.

In the second study, I examined desiccation resistance in intertidal egg masses made by Haminoea vesicula, which preferentially attaches its flat, ribbon-shaped egg masses to submerged substrata. Egg masses occasionally detach and become stranded on exposed sand at low tide. Unlike adults, the encased embryos cannot avoid desiccation by selectively moving about the habitat, and the egg mass shape has high surface-area-to-volume ratio that should make it prone to drying out. Thus, I hypothesized that the embryos would not survive stranding. I tested this by deploying individual egg masses of two age classes on exposed sand bars for the duration of low tide. After rehydration, embryos midway through development showed higher rates of survival than newly-laid embryos, though for both stages survival rates over 25% were frequently observed. Laboratory desiccation trials showed that >75% survival is possible in an egg mass that has lost 65% of its water weight, and some survival (<25%) was observed even after 83% water weight lost. Although many surviving embryos in both experiments showed damage, these data demonstrate that egg mass stranding is not necessarily fatal to embryos. They may be able to survive a far greater range of conditions than they normally encounter, compensating for their lack of ability to move. Also, desiccation tolerance of embryos may reduce pressure on parents to find optimal laying substrata.

The third study takes a big-picture approach to investigating the evolution of different developmental strategies in cone snails, the largest genus of marine invertebrates. Cone snail species hatch out of their capsules as either swimming larvae or non-dispersing forms, and their developmental mode has direct consequences for biogeographic patterns. Variability in life history strategies among taxa may be influenced by biological, environmental, or phylogenetic factors, or a combination of these. While most prior research has examined these factors singularly, my aim was to investigate the effects of a host of intrinsic, extrinsic, and historical factors on two fundamental aspects of life history: egg size and egg number. I used phylogenetic generalized least-squares regression models to examine relationships between these two egg traits and a variety of hypothesized intrinsic and extrinsic variables. Adult shell morphology and spatial variability in productivity and salinity across a species geographic range had the strongest effects on egg diameter and number of eggs per capsule. Phylogeny had no significant influence. Developmental mode in Conus appears to be influenced mostly by species-level adaptations and niche specificity rather than phylogenetic conservatism. Patterns of egg size and egg number appear to reflect energetic tradeoffs with body size and specific morphologies as well as adaptations to variable environments. Overall, this series of studies highlights the importance of organism-scale biotic and abiotic interactions in evolutionary patterns.

Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.