Assessment of LD matrix measures for the analysis of biological pathway association.

Complex diseases will have multiple functional sites, and it will be invaluable to understand the cross-locus interaction in terms of linkage disequilibrium (LD) between those sites (epistasis) in addition to the haplotype-LD effects. We investigated the statistical properties of a class of matrix-based statistics to assess this epistasis. These statistical methods include two LD contrast tests (Zaykin et al., 2006) and partial least squares regression (Wang et al., 2008). To estimate Type 1 error rates and power, we simulated multiple two-variant disease models using the SIMLA software package. SIMLA allows for the joint action of up to two disease genes in the simulated data with all possible multiplicative interaction effects between them. Our goal was to detect an interaction between multiple disease-causing variants by means of their linkage disequilibrium (LD) patterns with other markers. We measured the effects of marginal disease effect size, haplotype LD, disease prevalence and minor allele frequency have on cross-locus interaction (epistasis). In the setting of strong allele effects and strong interaction, the correlation between the two disease genes was weak (r=0.2). In a complex system with multiple correlations (both marginal and interaction), it was difficult to determine the source of a significant result. Despite these complications, the partial least squares and modified LD contrast methods maintained adequate power to detect the epistatic effects; however, for many of the analyses we often could not separate interaction from a strong marginal effect. While we did not exhaust the entire parameter space of possible models, we do provide guidance on the effects that population parameters have on cross-locus interaction.

Barriers to and facilitators of interventions to counter publication bias: thematic analysis of scholarly articles and stakeholder interviews.

BACKGROUND: When the nature and direction of research results affect their chances of publication, a distortion of the evidence base - termed publication bias - results. Despite considerable recent efforts to implement measures to reduce the non-publication of trials, publication bias is still a major problem in medical research. The objective of our study was to identify barriers to and facilitators of interventions to prevent or reduce publication bias. METHODS: We systematically reviewed the scholarly literature and extracted data from articles. Further, we performed semi-structured interviews with stakeholders. We performed an inductive thematic analysis to identify barriers to and facilitators of interventions to counter publication bias. RESULTS: The systematic review identified 39 articles. Thirty-four of 89 invited interview partners agreed to be interviewed. We clustered interventions into four categories: prospective trial registration, incentives for reporting in peer-reviewed journals or research reports, public availability of individual patient-level data, and peer-review/editorial processes. Barriers we identified included economic and personal interests, lack of financial resources for a global comprehensive trial registry, and different legal systems. Facilitators identified included: raising awareness of the effects of publication bias, providing incentives to make data publically available, and implementing laws to enforce prospective registration and reporting of clinical trial results. CONCLUSIONS: Publication bias is a complex problem that reflects the complex system in which it occurs. The cooperation amongst stakeholders to increase public awareness of the problem, better tailoring of incentives to publish, and ultimately legislative regulations have the greatest potential for reducing publication bias.

Morphological and genomic characterization of Filobasidiella depauperata: a homothallic sibling species of the pathogenic cryptococcus species complex.

The fungal species Cryptococcus neoformans and Cryptococcus gattii cause respiratory and neurological disease in animals and humans following inhalation of basidiospores or desiccated yeast cells from the environment. Sexual reproduction in C. neoformans and C. gattii is controlled by a bipolar system in which a single mating type locus (MAT) specifies compatibility. These two species are dimorphic, growing as yeast in the asexual stage, and producing hyphae, basidia, and basidiospores during the sexual stage. In contrast, Filobasidiella depauperata, one of the closest related species, grows exclusively as hyphae and it is found in association with decaying insects. Examination of two available strains of F. depauperata showed that the life cycle of this fungal species shares features associated with the unisexual or same-sex mating cycle in C. neoformans. Therefore, F. depauperata may represent a homothallic and possibly an obligately sexual fungal species. RAPD genotyping of 39 randomly isolated progeny from isolate CBS7855 revealed a new genotype pattern in one of the isolated basidiospores progeny, therefore suggesting that the homothallic cycle in F. depauperata could lead to the emergence of new genotypes. Phylogenetic analyses of genes linked to MAT in C. neoformans indicated that two of these genes in F. depauperata, MYO2 and STE20, appear to form a monophyletic clade with the MATa alleles of C. neoformans and C. gattii, and thus these genes may have been recruited to the MAT locus before F. depauperata diverged. Furthermore, the ancestral MATa locus may have undergone accelerated evolution prior to the divergence of the pathogenic Cryptococcus species since several of the genes linked to the MATa locus appear to have a higher number of changes and substitutions than their MATalpha counterparts. Synteny analyses between C. neoformans and F. depauperata showed that genomic regions on other chromosomes displayed conserved gene order. In contrast, the genes linked to the MAT locus of C. neoformans showed a higher number of chromosomal translocations in the genome of F. depauperata. We therefore propose that chromosomal rearrangements appear to be a major force driving speciation and sexual divergence in these closely related pathogenic and saprobic species.

Structure, function, and phylogeny of the mating locus in the Rhizopus oryzae complex.

The Rhizopus oryzae species complex is a group of zygomycete fungi that are common, cosmopolitan saprotrophs. Some strains are used beneficially for production of Asian fermented foods but they can also act as opportunistic human pathogens. Although R. oryzae reportedly has a heterothallic (+/-) mating system, most strains have not been observed to undergo sexual reproduction and the genetic structure of its mating locus has not been characterized. Here we report on the mating behavior and genetic structure of the mating locus for 54 isolates of the R. oryzae complex. All 54 strains have a mating locus similar in overall organization to Phycomyces blakesleeanus and Mucor circinelloides (Mucoromycotina, Zygomycota). In all of these fungi, the minus (-) allele features the SexM high mobility group (HMG) gene flanked by an RNA helicase gene and a TP transporter gene (TPT). Within the R. oryzae complex, the plus (+) mating allele includes an inserted region that codes for a BTB/POZ domain gene and the SexP HMG gene. Phylogenetic analyses of multiple genes, including the mating loci (HMG, TPT, RNA helicase), ITS1-5.8S-ITS2 rDNA, RPB2, and LDH genes, identified two distinct groups of strains. These correspond to previously described sibling species R. oryzae sensu stricto and R. delemar. Within each species, discordant gene phylogenies among multiple loci suggest an outcrossing population structure. The hypothesis of random-mating is also supported by a 50:50 ratio of plus and minus mating types in both cryptic species. When crossed with tester strains of the opposite mating type, most isolates of R. delemar failed to produce zygospores, while isolates of R. oryzae produced sterile zygospores. In spite of the reluctance of most strains to mate in vitro, the conserved sex locus structure and evidence for outcrossing suggest that a normal sexual cycle occurs in both species.

Real-Time and Data-Driven Operation Optimization and Knowledge Discovery for an Enterprise Information System

An enterprise information system (EIS) is an integrated data-applications platform characterized by diverse, heterogeneous, and distributed data sources. For many enterprises, a number of business processes still depend heavily on static rule-based methods and extensive human expertise. Enterprises are faced with the need for optimizing operation scheduling, improving resource utilization, discovering useful knowledge, and making data-driven decisions.

This thesis research is focused on real-time optimization and knowledge discovery that addresses workflow optimization, resource allocation, as well as data-driven predictions of process-execution times, order fulfillment, and enterprise service-level performance. In contrast to prior work on data analytics techniques for enterprise performance optimization, the emphasis here is on realizing scalable and real-time enterprise intelligence based on a combination of heterogeneous system simulation, combinatorial optimization, machine-learning algorithms, and statistical methods.

On-demand digital-print service is a representative enterprise requiring a powerful EIS.We use real-life data from Reischling Press, Inc. (RPI), a digit-print-service provider (PSP), to evaluate our optimization algorithms.

In order to handle the increase in volume and diversity of demands, we first present a high-performance, scalable, and real-time production scheduling algorithm for production automation based on an incremental genetic algorithm (IGA). The objective of this algorithm is to optimize the order dispatching sequence and balance resource utilization. Compared to prior work, this solution is scalable for a high volume of orders and it provides fast scheduling solutions for orders that require complex fulfillment procedures. Experimental results highlight its potential benefit in reducing production inefficiencies and enhancing the productivity of an enterprise.

We next discuss analysis and prediction of different attributes involved in hierarchical components of an enterprise. We start from a study of the fundamental processes related to real-time prediction. Our process-execution time and process status prediction models integrate statistical methods with machine-learning algorithms. In addition to improved prediction accuracy compared to stand-alone machine-learning algorithms, it also performs a probabilistic estimation of the predicted status. An order generally consists of multiple series and parallel processes. We next introduce an order-fulfillment prediction model that combines advantages of multiple classification models by incorporating flexible decision-integration mechanisms. Experimental results show that adopting due dates recommended by the model can significantly reduce enterprise late-delivery ratio. Finally, we investigate service-level attributes that reflect the overall performance of an enterprise. We analyze and decompose time-series data into different components according to their hierarchical periodic nature, perform correlation analysis,

and develop univariate prediction models for each component as well as multivariate models for correlated components. Predictions for the original time series are aggregated from the predictions of its components. In addition to a significant increase in mid-term prediction accuracy, this distributed modeling strategy also improves short-term time-series prediction accuracy.

In summary, this thesis research has led to a set of characterization, optimization, and prediction tools for an EIS to derive insightful knowledge from data and use them as guidance for production management. It is expected to provide solutions for enterprises to increase reconfigurability, accomplish more automated procedures, and obtain data-driven recommendations or effective decisions.

The Innate Immune System in Acute and Chronic Wounds.

Significance: This review article provides an overview of the critical roles of the innate immune system to wound healing. It explores aspects of dysregulation of individual innate immune elements known to compromise wound repair and promote nonhealing wounds. Understanding the key mechanisms whereby wound healing fails will provide seed concepts for the development of new therapeutic approaches. Recent Advances: Our understanding of the complex interactions of the innate immune system in wound healing has significantly improved, particularly in our understanding of the role of antimicrobials and peptides and the nature of the switch from inflammatory to reparative processes. This takes place against an emerging understanding of the relationship between human cells and commensal bacteria in the skin. Critical Issues: It is well established and accepted that early local inflammatory mediators in the wound bed function as an immunological vehicle to facilitate immune cell infiltration and microbial clearance upon injury to the skin barrier. Both impaired and excessive innate immune responses can promote nonhealing wounds. It appears that the switch from the inflammatory to the proliferative phase is tightly regulated and mediated, at least in part, by a change in macrophages. Defining the factors that initiate the switch in such macrophage phenotypes and functions is the subject of multiple investigations. Future Directions: The review highlights processes that may be useful targets for further investigation, particularly the switch from M1 to M2 macrophages that appears to be critical as dysregulation of this switch occurs during defective wound healing.

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.

Evaluating functional network inference using simulations of complex biological systems.

MOTIVATION: Although many network inference algorithms have been presented in the bioinformatics literature, no suitable approach has been formulated for evaluating their effectiveness at recovering models of complex biological systems from limited data. To overcome this limitation, we propose an approach to evaluate network inference algorithms according to their ability to recover a complex functional network from biologically reasonable simulated data. RESULTS: We designed a simulator to generate data representing a complex biological system at multiple levels of organization: behaviour, neural anatomy, brain electrophysiology, and gene expression of songbirds. About 90% of the simulated variables are unregulated by other variables in the system and are included simply as distracters. We sampled the simulated data at intervals as one would sample from a biological system in practice, and then used the sampled data to evaluate the effectiveness of an algorithm we developed for functional network inference. We found that our algorithm is highly effective at recovering the functional network structure of the simulated system-including the irrelevance of unregulated variables-from sampled data alone. To assess the reproducibility of these results, we tested our inference algorithm on 50 separately simulated sets of data and it consistently recovered almost perfectly the complex functional network structure underlying the simulated data. To our knowledge, this is the first approach for evaluating the effectiveness of functional network inference algorithms at recovering models from limited data. Our simulation approach also enables researchers a priori to design experiments and data-collection protocols that are amenable to functional network inference.

Investigation and Development of a Fully 3D Tilt Capable Hybrid SPECT - CT System for Dedicated Breast Imaging

X-ray mammography has been the gold standard for breast imaging for decades, despite the significant limitations posed by the two dimensional (2D) image acquisitions. Difficulty in diagnosing lesions close to the chest wall and axilla, high amount of structural overlap and patient discomfort due to compression are only some of these limitations. To overcome these drawbacks, three dimensional (3D) breast imaging modalities have been developed including dual modality single photon emission computed tomography (SPECT) and computed tomography (CT) systems. This thesis focuses on the development and integration of the next generation of such a device for dedicated breast imaging. The goals of this dissertation work are to: [1] understand and characterize any effects of fully 3-D trajectories on reconstructed image scatter correction, absorbed dose and Hounsifeld Unit accuracy, and [2] design, develop and implement the fully flexible, third generation hybrid SPECT-CT system capable of traversing complex 3D orbits about a pendant breast volume, without interference from the other. Such a system would overcome artifacts resulting from incompletely sampled divergent cone beam imaging schemes and allow imaging closer to the chest wall, which other systems currently under research and development elsewhere cannot achieve.

The dependence of x-ray scatter radiation on object shape, size, material composition and the CT acquisition trajectory, was investigated with a well-established beam stop array (BSA) scatter correction method. While the 2D scatter to primary ratio (SPR) was the main metric used to characterize total system scatter, a new metric called ‘normalized scatter contribution’ was developed to compare the results of scatter correction on 3D reconstructed volumes. Scatter estimation studies were undertaken with a sinusoidal saddle (±15° polar tilt) orbit and a traditional circular (AZOR) orbit. Clinical studies to acquire data for scatter correction were used to evaluate the 2D SPR on a small set of patients scanned with the AZOR orbit. Clinical SPR results showed clear dependence of scatter on breast composition and glandular tissue distribution, otherwise consistent with the overall phantom-based size and density measurements. Additionally, SPR dependence was also observed on the acquisition trajectory where 2D scatter increased with an increase in the polar tilt angle of the system.

The dose delivered by any imaging system is of primary importance from the patient’s point of view, and therefore trajectory related differences in the dose distribution in a target volume were evaluated. Monte Carlo simulations as well as physical measurements using radiochromic film were undertaken using saddle and AZOR orbits. Results illustrated that both orbits deliver comparable dose to the target volume, and only slightly differ in distribution within the volume. Simulations and measurements showed similar results, and all measured dose values were within the standard screening mammography-specific, 6 mGy dose limit, which is used as a benchmark for dose comparisons.

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is very important, especially when using non-traditional trajectories. Uniform phantoms filled with various uniform density fluids were used to investigate differences in HU accuracy between saddle and AZOR orbits. Results illustrate the considerably better performance of the saddle orbit, especially close to the chest and nipple region of what would clinically be a pedant breast volume. The AZOR orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values in reconstructed volumes.

Finally, the third generation, fully-suspended SPECT-CT system was designed in and developed in our lab. A novel mechanical method using a linear motor was developed for tilting the CT system. A new x-ray source and a custom made 40 x 30 cm2 detector were integrated on to this system. The SPECT system was nested, in the center of the gantry, orthogonal to the CT source-detector pair. The SPECT system tilts on a goniometer, and the newly developed CT tilting mechanism allows ±15° maximum polar tilting of the CT system. The entire gantry is mounted on a rotation stage, allowing complex arbitrary trajectories for each system, without interference from the other, while having a common field of view. This hybrid system shows potential to be used clinically as a diagnostic tool for dedicated breast imaging.