Advanced Techniques for Image Quality Assessment of Modern X-ray Computed Tomography Systems

X-ray computed tomography (CT) imaging constitutes one of the most widely used diagnostic tools in radiology today with nearly 85 million CT examinations performed in the U.S in 2011. CT imparts a relatively high amount of radiation dose to the patient compared to other x-ray imaging modalities and as a result of this fact, coupled with its popularity, CT is currently the single largest source of medical radiation exposure to the U.S. population. For this reason, there is a critical need to optimize CT examinations such that the dose is minimized while the quality of the CT images is not degraded. This optimization can be difficult to achieve due to the relationship between dose and image quality. All things being held equal, reducing the dose degrades image quality and can impact the diagnostic value of the CT examination.

A recent push from the medical and scientific community towards using lower doses has spawned new dose reduction technologies such as automatic exposure control (i.e., tube current modulation) and iterative reconstruction algorithms. In theory, these technologies could allow for scanning at reduced doses while maintaining the image quality of the exam at an acceptable level. Therefore, there is a scientific need to establish the dose reduction potential of these new technologies in an objective and rigorous manner. Establishing these dose reduction potentials requires precise and clinically relevant metrics of CT image quality, as well as practical and efficient methodologies to measure such metrics on real CT systems. The currently established methodologies for assessing CT image quality are not appropriate to assess modern CT scanners that have implemented those aforementioned dose reduction technologies.

Thus the purpose of this doctoral project was to develop, assess, and implement new phantoms, image quality metrics, analysis techniques, and modeling tools that are appropriate for image quality assessment of modern clinical CT systems. The project developed image quality assessment methods in the context of three distinct paradigms, (a) uniform phantoms, (b) textured phantoms, and (c) clinical images.

The work in this dissertation used the “task-based” definition of image quality. That is, image quality was broadly defined as the effectiveness by which an image can be used for its intended task. Under this definition, any assessment of image quality requires three components: (1) A well defined imaging task (e.g., detection of subtle lesions), (2) an “observer” to perform the task (e.g., a radiologists or a detection algorithm), and (3) a way to measure the observer’s performance in completing the task at hand (e.g., detection sensitivity/specificity).

First, this task-based image quality paradigm was implemented using a novel multi-sized phantom platform (with uniform background) developed specifically to assess modern CT systems (Mercury Phantom, v3.0, Duke University). A comprehensive evaluation was performed on a state-of-the-art CT system (SOMATOM Definition Force, Siemens Healthcare) in terms of noise, resolution, and detectability as a function of patient size, dose, tube energy (i.e., kVp), automatic exposure control, and reconstruction algorithm (i.e., Filtered Back-Projection– FPB vs Advanced Modeled Iterative Reconstruction– ADMIRE). A mathematical observer model (i.e., computer detection algorithm) was implemented and used as the basis of image quality comparisons. It was found that image quality increased with increasing dose and decreasing phantom size. The CT system exhibited nonlinear noise and resolution properties, especially at very low-doses, large phantom sizes, and for low-contrast objects. Objective image quality metrics generally increased with increasing dose and ADMIRE strength, and with decreasing phantom size. The ADMIRE algorithm could offer comparable image quality at reduced doses or improved image quality at the same dose (increase in detectability index by up to 163% depending on iterative strength). The use of automatic exposure control resulted in more consistent image quality with changing phantom size.

Based on those results, the dose reduction potential of ADMIRE was further assessed specifically for the task of detecting small (<=6 mm) low-contrast (<=20 HU) lesions. A new low-contrast detectability phantom (with uniform background) was designed and fabricated using a multi-material 3D printer. The phantom was imaged at multiple dose levels and images were reconstructed with FBP and ADMIRE. Human perception experiments were performed to measure the detection accuracy from FBP and ADMIRE images. It was found that ADMIRE had equivalent performance to FBP at 56% less dose.

Using the same image data as the previous study, a number of different mathematical observer models were implemented to assess which models would result in image quality metrics that best correlated with human detection performance. The models included naïve simple metrics of image quality such as contrast-to-noise ratio (CNR) and more sophisticated observer models such as the non-prewhitening matched filter observer model family and the channelized Hotelling observer model family. It was found that non-prewhitening matched filter observers and the channelized Hotelling observers both correlated strongly with human performance. Conversely, CNR was found to not correlate strongly with human performance, especially when comparing different reconstruction algorithms.

The uniform background phantoms used in the previous studies provided a good first-order approximation of image quality. However, due to their simplicity and due to the complexity of iterative reconstruction algorithms, it is possible that such phantoms are not fully adequate to assess the clinical impact of iterative algorithms because patient images obviously do not have smooth uniform backgrounds. To test this hypothesis, two textured phantoms (classified as gross texture and fine texture) and a uniform phantom of similar size were built and imaged on a SOMATOM Flash scanner (Siemens Healthcare). Images were reconstructed using FBP and a Sinogram Affirmed Iterative Reconstruction (SAFIRE). Using an image subtraction technique, quantum noise was measured in all images of each phantom. It was found that in FBP, the noise was independent of the background (textured vs uniform). However, for SAFIRE, noise increased by up to 44% in the textured phantoms compared to the uniform phantom. As a result, the noise reduction from SAFIRE was found to be up to 66% in the uniform phantom but as low as 29% in the textured phantoms. Based on this result, it clear that further investigation was needed into to understand the impact that background texture has on image quality when iterative reconstruction algorithms are used.

To further investigate this phenomenon with more realistic textures, two anthropomorphic textured phantoms were designed to mimic lung vasculature and fatty soft tissue texture. The phantoms (along with a corresponding uniform phantom) were fabricated with a multi-material 3D printer and imaged on the SOMATOM Flash scanner. Scans were repeated a total of 50 times in order to get ensemble statistics of the noise. A novel method of estimating the noise power spectrum (NPS) from irregularly shaped ROIs was developed. It was found that SAFIRE images had highly locally non-stationary noise patterns with pixels near edges having higher noise than pixels in more uniform regions. Compared to FBP, SAFIRE images had 60% less noise on average in uniform regions for edge pixels, noise was between 20% higher and 40% lower. The noise texture (i.e., NPS) was also highly dependent on the background texture for SAFIRE. Therefore, it was concluded that quantum noise properties in the uniform phantoms are not representative of those in patients for iterative reconstruction algorithms and texture should be considered when assessing image quality of iterative algorithms.

The move beyond just assessing noise properties in textured phantoms towards assessing detectability, a series of new phantoms were designed specifically to measure low-contrast detectability in the presence of background texture. The textures used were optimized to match the texture in the liver regions actual patient CT images using a genetic algorithm. The so called “Clustured Lumpy Background” texture synthesis framework was used to generate the modeled texture. Three textured phantoms and a corresponding uniform phantom were fabricated with a multi-material 3D printer and imaged on the SOMATOM Flash scanner. Images were reconstructed with FBP and SAFIRE and analyzed using a multi-slice channelized Hotelling observer to measure detectability and the dose reduction potential of SAFIRE based on the uniform and textured phantoms. It was found that at the same dose, the improvement in detectability from SAFIRE (compared to FBP) was higher when measured in a uniform phantom compared to textured phantoms.

The final trajectory of this project aimed at developing methods to mathematically model lesions, as a means to help assess image quality directly from patient images. The mathematical modeling framework is first presented. The models describe a lesion’s morphology in terms of size, shape, contrast, and edge profile as an analytical equation. The models can be voxelized and inserted into patient images to create so-called “hybrid” images. These hybrid images can then be used to assess detectability or estimability with the advantage that the ground truth of the lesion morphology and location is known exactly. Based on this framework, a series of liver lesions, lung nodules, and kidney stones were modeled based on images of real lesions. The lesion models were virtually inserted into patient images to create a database of hybrid images to go along with the original database of real lesion images. ROI images from each database were assessed by radiologists in a blinded fashion to determine the realism of the hybrid images. It was found that the radiologists could not readily distinguish between real and virtual lesion images (area under the ROC curve was 0.55). This study provided evidence that the proposed mathematical lesion modeling framework could produce reasonably realistic lesion images.

Based on that result, two studies were conducted which demonstrated the utility of the lesion models. The first study used the modeling framework as a measurement tool to determine how dose and reconstruction algorithm affected the quantitative analysis of liver lesions, lung nodules, and renal stones in terms of their size, shape, attenuation, edge profile, and texture features. The same database of real lesion images used in the previous study was used for this study. That database contained images of the same patient at 2 dose levels (50% and 100%) along with 3 reconstruction algorithms from a GE 750HD CT system (GE Healthcare). The algorithms in question were FBP, Adaptive Statistical Iterative Reconstruction (ASiR), and Model-Based Iterative Reconstruction (MBIR). A total of 23 quantitative features were extracted from the lesions under each condition. It was found that both dose and reconstruction algorithm had a statistically significant effect on the feature measurements. In particular, radiation dose affected five, three, and four of the 23 features (related to lesion size, conspicuity, and pixel-value distribution) for liver lesions, lung nodules, and renal stones, respectively. MBIR significantly affected 9, 11, and 15 of the 23 features (including size, attenuation, and texture features) for liver lesions, lung nodules, and renal stones, respectively. Lesion texture was not significantly affected by radiation dose.

The second study demonstrating the utility of the lesion modeling framework focused on assessing detectability of very low-contrast liver lesions in abdominal imaging. Specifically, detectability was assessed as a function of dose and reconstruction algorithm. As part of a parallel clinical trial, images from 21 patients were collected at 6 dose levels per patient on a SOMATOM Flash scanner. Subtle liver lesion models (contrast = -15 HU) were inserted into the raw projection data from the patient scans. The projections were then reconstructed with FBP and SAFIRE (strength 5). Also, lesion-less images were reconstructed. Noise, contrast, CNR, and detectability index of an observer model (non-prewhitening matched filter) were assessed. It was found that SAFIRE reduced noise by 52%, reduced contrast by 12%, increased CNR by 87%. and increased detectability index by 65% compared to FBP. Further, a 2AFC human perception experiment was performed to assess the dose reduction potential of SAFIRE, which was found to be 22% compared to the standard of care dose.

In conclusion, this dissertation provides to the scientific community a series of new methodologies, phantoms, analysis techniques, and modeling tools that can be used to rigorously assess image quality from modern CT systems. Specifically, methods to properly evaluate iterative reconstruction have been developed and are expected to aid in the safe clinical implementation of dose reduction technologies.

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Coupling of beta2-adrenoceptor to Gi proteins and its physiological relevance in murine cardiac myocytes.

-Transgenic mouse models have been developed to manipulate beta-adrenergic receptor (betaAR) signal transduction. Although several of these models have altered betaAR subtypes, the specific functional sequelae of betaAR stimulation in murine heart, particularly those of beta2-adrenergic receptor (beta2AR) stimulation, have not been characterized. In the present study, we investigated effects of beta2AR stimulation on contraction, [Ca2+]i transient, and L-type Ca2+ currents (ICa) in single ventricular myocytes isolated from transgenic mice overexpressing human beta2AR (TG4 mice) and wild-type (WT) littermates. Baseline contractility of TG4 heart cells was increased by 3-fold relative to WT controls as a result of the presence of spontaneous beta2AR activation. In contrast, beta2AR stimulation by zinterol or isoproterenol plus a selective beta1-adrenergic receptor (beta1AR) antagonist CGP 20712A failed to enhance the contractility in TG4 myocytes, and more surprisingly, beta2AR stimulation was also ineffective in increasing contractility in WT myocytes. Pertussis toxin (PTX) treatment fully rescued the ICa, [Ca2+]i, and contractile responses to beta2AR agonists in both WT and TG4 cells. The PTX-rescued murine cardiac beta2AR response is mediated by cAMP-dependent mechanisms, because it was totally blocked by the inhibitory cAMP analog Rp-cAMPS. These results suggest that PTX-sensitive G proteins are responsible for the unresponsiveness of mouse heart to agonist-induced beta2AR stimulation. This was further corroborated by an increased incorporation of the photoreactive GTP analog [gamma-32P]GTP azidoanilide into alpha subunits of Gi2 and Gi3 after beta2AR stimulation by zinterol or isoproterenol plus the beta1AR blocker CGP 20712A. This effect to activate Gi proteins was abolished by a selective beta2AR blocker ICI 118,551 or by PTX treatment. Thus, we conclude that (1) beta2ARs in murine cardiac myocytes couple to concurrent Gs and Gi signaling, resulting in null inotropic response, unless the Gi signaling is inhibited; (2) as a special case, the lack of cardiac contractile response to beta2AR agonists in TG4 mice is not due to a saturation of cell contractility or of the cAMP signaling cascade but rather to an activation of beta2AR-coupled Gi proteins; and (3) spontaneous beta2AR activation may differ from agonist-stimulated beta2AR signaling.

Enhancement of cardiac function after adenoviral-mediated in vivo intracoronary beta2-adrenergic receptor gene delivery.

Exogenous gene delivery to alter the function of the heart is a potential novel therapeutic strategy for treatment of cardiovascular diseases such as heart failure (HF). Before gene therapy approaches to alter cardiac function can be realized, efficient and reproducible in vivo gene techniques must be established to efficiently transfer transgenes globally to the myocardium. We have been testing the hypothesis that genetic manipulation of the myocardial beta-adrenergic receptor (beta-AR) system, which is impaired in HF, can enhance cardiac function. We have delivered adenoviral transgenes, including the human beta2-AR (Adeno-beta2AR), to the myocardium of rabbits using an intracoronary approach. Catheter-mediated Adeno-beta2AR delivery produced diffuse multichamber myocardial expression, peaking 1 week after gene transfer. A total of 5 x 10(11) viral particles of Adeno-beta2AR reproducibly produced 5- to 10-fold beta-AR overexpression in the heart, which, at 7 and 21 days after delivery, resulted in increased in vivo hemodynamic function compared with control rabbits that received an empty adenovirus. Several physiological parameters, including dP/dtmax as a measure of contractility, were significantly enhanced basally and showed increased responsiveness to the beta-agonist isoproterenol. Our results demonstrate that global myocardial in vivo gene delivery is possible and that genetic manipulation of beta-AR density can result in enhanced cardiac performance. Thus, replacement of lost receptors seen in HF may represent novel inotropic therapy.

Human dimensions of bycatch reduction technology: Current assumptions and directions for future research

Bycatch reduction technology (BRT) modifies fishing gear to increase selectivity and avoid capture of non-target species, or to facilitate their non-lethal release. As a solution to fisheries-related mortality of non-target species, BRT is an attractive option; effectively implemented, BRT presents a technical 'fix' that can reduce pressure for politically contentious and economically detrimental interventions, such as fisheries closures. While a number of factors might contribute to effective implementation, our review of BRT literature finds that research has focused on technical design and experimental performance of individual technologies. In contrast, and with a few notable exceptions, research on the human and institutional context of BRT, and more specifically on how fishers respond to BRT, is limited. This is not to say that fisher attitudes are ignored or overlooked, but that incentives for fisher uptake of BRT are usually assumed rather than assessed or demonstrated. Three assumptions about fisher incentives dominate: (1) economic incentives will generate acceptance of BRT; (2) enforcement will generate compliance with BRT; and (3) 'participation' by fishers will increase acceptance and compliance, and overall support for BRT. In this paper, we explore evidence for and against these assumptions and situate our analysis in the wider social science literature on fisheries. Our goal is to highlight the need and suggest focal areas for further research. © Inter-Research 2008.

Bayesian inference of the number of factors in gene-expression analysis: application to human virus challenge studies.

BACKGROUND: Nonparametric Bayesian techniques have been developed recently to extend the sophistication of factor models, allowing one to infer the number of appropriate factors from the observed data. We consider such techniques for sparse factor analysis, with application to gene-expression data from three virus challenge studies. Particular attention is placed on employing the Beta Process (BP), the Indian Buffet Process (IBP), and related sparseness-promoting techniques to infer a proper number of factors. The posterior density function on the model parameters is computed using Gibbs sampling and variational Bayesian (VB) analysis. RESULTS: Time-evolving gene-expression data are considered for respiratory syncytial virus (RSV), Rhino virus, and influenza, using blood samples from healthy human subjects. These data were acquired in three challenge studies, each executed after receiving institutional review board (IRB) approval from Duke University. Comparisons are made between several alternative means of per-forming nonparametric factor analysis on these data, with comparisons as well to sparse-PCA and Penalized Matrix Decomposition (PMD), closely related non-Bayesian approaches. CONCLUSIONS: Applying the Beta Process to the factor scores, or to the singular values of a pseudo-SVD construction, the proposed algorithms infer the number of factors in gene-expression data. For real data the "true" number of factors is unknown; in our simulations we consider a range of noise variances, and the proposed Bayesian models inferred the number of factors accurately relative to other methods in the literature, such as sparse-PCA and PMD. We have also identified a "pan-viral" factor of importance for each of the three viruses considered in this study. We have identified a set of genes associated with this pan-viral factor, of interest for early detection of such viruses based upon the host response, as quantified via gene-expression data.

Bioengineered human myobundles mimic clinical responses of skeletal muscle to drugs.

Existing in vitro models of human skeletal muscle cannot recapitulate the organization and function of native muscle, limiting their use in physiological and pharmacological studies. Here, we demonstrate engineering of electrically and chemically responsive, contractile human muscle tissues ('myobundles') using primary myogenic cells. These biomimetic constructs exhibit aligned architecture, multinucleated and striated myofibers, and a Pax7(+) cell pool. They contract spontaneously and respond to electrical stimuli with twitch and tetanic contractions. Positive correlation between contractile force and GCaMP6-reported calcium responses enables non-invasive tracking of myobundle function and drug response. During culture, myobundles maintain functional acetylcholine receptors and structurally and functionally mature, evidenced by increased myofiber diameter and improved calcium handling and contractile strength. In response to diversely acting drugs, myobundles undergo dose-dependent hypertrophy or toxic myopathy similar to clinical outcomes. Human myobundles provide an enabling platform for predictive drug and toxicology screening and development of novel therapeutics for muscle-related disorders.

Improving Indwelling Glucose Sensor Performance: Porous, Dexamethasone-Releasing Coatings that Modulate the Foreign Body Response

Inflammation and the formation of an avascular fibrous capsule have been identified as the key factors controlling the wound healing associated failure of implantable glucose sensors. Our aim is to guide advantageous tissue remodeling around implanted sensor leads by the temporal release of dexamethasone (Dex), a potent anti-inflammatory agent, in combination with the presentation of a stable textured surface.

First, Dex-releasing polyurethane porous coatings of controlled pore size and thickness were fabricated using salt-leaching/gas-foaming technique. Porosity, pore size, thickness, drug release kinetics, drug loading amount, and drug bioactivity were evaluated. In vitro sensor functionality test were performed to determine if Dex-releasing porous coatings interfered with sensor performance (increased signal attenuation and/or response times) compared to bare sensors. Drug release from coatings monitored over two weeks presented an initial fast release followed by a slower release. Total release from coatings was highly dependent on initial drug loading amount. Functional in vitro testing of glucose sensors deployed with porous coatings against glucose standards demonstrated that highly porous coatings minimally affected signal strength and response rate. Bioactivity of the released drug was determined by monitoring Dex-mediated, dose-dependent apoptosis of human peripheral blood derived monocytes in culture.

The tissue modifying effects of Dex-releasing porous coatings were accessed by fully implanting Tygon® tubing in the subcutaneous space of healthy and diabetic rats. Based on encouraging results from these studies, we deployed Dex-releasing porous coatings from the tips of functional sensors in both diabetic and healthy rats. We evaluated if the tissue modifying effects translated into accurate, maintainable and reliable sensor signals in the long-term. Sensor functionality was accessed by continuously monitoring glucose levels and performing acute glucose challenges at specified time points.

Sensors treated with porous Dex-releasing coatings showed diminished inflammation and enhanced vascularization of the tissue surrounding the implants in healthy rats. Functional sensors with Dex-releasing porous coatings showed enhanced sensor sensitivity over a 21-day period when compared to controls. Enhanced sensor sensitivity was accompanied with an increase in sensor signal lag and MARD score. These results indicated that Dex-loaded porous coatings were able to elicit a favorable tissue response, and that such tissue microenvironment could be conducive towards extending the performance window of glucose sensors in vivo.

The diabetic pilot animal study showed differences in wound healing patters between healthy and diabetic subjects. Diabetic rats showed lower levels of inflammation and vascularization of the tissue surrounding implants when compared to their healthy counterparts. Also, functional sensors treated with Dex-releasing porous coatings did not show enhanced sensor sensitivity over a 21-day period. Moreover, increased in sensor signal lag and MARD scores were present in porous coated sensors regardless of Dex-loading when compared to bare implants. These results suggest that the altered wound healing patterns presented in diabetic tissues may lead to premature sensor failure when compared to sensors implanted in healthy rats.

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.

Human Vascular Microphysiological System for in vitro Drug Screening.

In vitro human tissue engineered human blood vessels (TEBV) that exhibit vasoactivity can be used to test human toxicity of pharmaceutical drug candidates prior to pre-clinical animal studies. TEBVs with 400-800 μM diameters were made by embedding human neonatal dermal fibroblasts or human bone marrow-derived mesenchymal stem cells in dense collagen gel. TEBVs were mechanically strong enough to allow endothelialization and perfusion at physiological shear stresses within 3 hours after fabrication. After 1 week of perfusion, TEBVs exhibited endothelial release of nitric oxide, phenylephrine-induced vasoconstriction, and acetylcholine-induced vasodilation, all of which were maintained up to 5 weeks in culture. Vasodilation was blocked with the addition of the nitric oxide synthase inhibitor L-N(G)-Nitroarginine methyl ester (L-NAME). TEBVs elicited reversible activation to acute inflammatory stimulation by TNF-α which had a transient effect upon acetylcholine-induced relaxation, and exhibited dose-dependent vasodilation in response to caffeine and theophylline. Treatment of TEBVs with 1 μM lovastatin for three days prior to addition of Tumor necrosis factor - α (TNF-α) blocked the injury response and maintained vasodilation. These results indicate the potential to develop a rapidly-producible, endothelialized TEBV for microphysiological systems capable of producing physiological responses to both pharmaceutical and immunological stimuli.