Comparison of pattern detection methods in microarray time series of the segmentation clock.

While genome-wide gene expression data are generated at an increasing rate, the repertoire of approaches for pattern discovery in these data is still limited. Identifying subtle patterns of interest in large amounts of data (tens of thousands of profiles) associated with a certain level of noise remains a challenge. A microarray time series was recently generated to study the transcriptional program of the mouse segmentation clock, a biological oscillator associated with the periodic formation of the segments of the body axis. A method related to Fourier analysis, the Lomb-Scargle periodogram, was used to detect periodic profiles in the dataset, leading to the identification of a novel set of cyclic genes associated with the segmentation clock. Here, we applied to the same microarray time series dataset four distinct mathematical methods to identify significant patterns in gene expression profiles. These methods are called: Phase consistency, Address reduction, Cyclohedron test and Stable persistence, and are based on different conceptual frameworks that are either hypothesis- or data-driven. Some of the methods, unlike Fourier transforms, are not dependent on the assumption of periodicity of the pattern of interest. Remarkably, these methods identified blindly the expression profiles of known cyclic genes as the most significant patterns in the dataset. Many candidate genes predicted by more than one approach appeared to be true positive cyclic genes and will be of particular interest for future research. In addition, these methods predicted novel candidate cyclic genes that were consistent with previous biological knowledge and experimental validation in mouse embryos. Our results demonstrate the utility of these novel pattern detection strategies, notably for detection of periodic profiles, and suggest that combining several distinct mathematical approaches to analyze microarray datasets is a valuable strategy for identifying genes that exhibit novel, interesting transcriptional patterns.

Features of programmed cell death in intact Xenopus oocytes and early embryos revealed by near-infrared fluorescence and real-time monitoring.

Factors influencing apoptosis of vertebrate eggs and early embryos have been studied in cell-free systems and in intact embryos by analyzing individual apoptotic regulators or caspase activation in static samples. A novel method for monitoring caspase activity in living Xenopus oocytes and early embryos is described here. The approach, using microinjection of a near-infrared caspase substrate that emits fluorescence only after its proteolytic cleavage by active effector caspases, has enabled the elucidation of otherwise cryptic aspects of apoptotic regulation. In particular, we show that brief caspase activity (10 min) is sufficient to cause apoptotic death in this system. We illustrate a cytochrome c dose threshold in the oocyte, which is lowered by Smac, a protein that binds thereby neutralizing the inhibitor of apoptosis proteins. We show that meiotic oocytes develop resistance to cytochrome c, and that the eventual death of oocytes arrested in meiosis is caspase-independent. Finally, data acquired through imaging caspase activity in the Xenopus embryo suggest that apoptosis in very early development is not cell-autonomous. These studies both validate this assay as a useful tool for apoptosis research and reveal subtleties in the cell death program during early development. Moreover, this method offers a potentially valuable screening modality for identifying novel apoptotic regulators.

Transcription factors MYOCD, SRF, Mesp1 and SMARCD3 enhance the cardio-inducing effect of GATA4, TBX5, and MEF2C during direct cellular reprogramming.

Transient overexpression of defined combinations of master regulator genes can effectively induce cellular reprogramming: the acquisition of an alternative predicted phenotype from a differentiated cell lineage. This can be of particular importance in cardiac regenerative medicine wherein the heart lacks the capacity to heal itself, but simultaneously contains a large pool of fibroblasts. In this study we determined the cardio-inducing capacity of ten transcription factors to actuate cellular reprogramming of mouse embryonic fibroblasts into cardiomyocyte-like cells. Overexpression of transcription factors MYOCD and SRF alone or in conjunction with Mesp1 and SMARCD3 enhanced the basal but necessary cardio-inducing effect of the previously reported GATA4, TBX5, and MEF2C. In particular, combinations of five or seven transcription factors enhanced the activation of cardiac reporter vectors, and induced an upregulation of cardiac-specific genes. Global gene expression analysis also demonstrated a significantly greater cardio-inducing effect when the transcription factors MYOCD and SRF were used. Detection of cross-striated cells was highly dependent on the cell culture conditions and was enhanced by the addition of valproic acid and JAK inhibitor. Although we detected Ca(2+) transient oscillations in the reprogrammed cells, we did not detect significant changes in resting membrane potential or spontaneously contracting cells. This study further elucidates the cardio-inducing effect of the transcriptional networks involved in cardiac cellular reprogramming, contributing to the ongoing rational design of a robust protocol required for cardiac regenerative therapies.

Radioactive in situ hybridization for detecting diverse gene expression patterns in tissue.

Knowing the timing, level, cellular localization, and cell type that a gene is expressed in contributes to our understanding of the function of the gene. Each of these features can be accomplished with in situ hybridization to mRNAs within cells. Here we present a radioactive in situ hybridization method modified from Clayton et al. (1988)(1) that has been working successfully in our lab for many years, especially for adult vertebrate brains(2-5). The long complementary RNA (cRNA) probes to the target sequence allows for detection of low abundance transcripts(6,7). Incorporation of radioactive nucleotides into the cRNA probes allows for further detection sensitivity of low abundance transcripts and quantitative analyses, either by light sensitive x-ray film or emulsion coated over the tissue. These detection methods provide a long-term record of target gene expression. Compared with non-radioactive probe methods, such as DIG-labeling, the radioactive probe hybridization method does not require multiple amplification steps using HRP-antibodies and/or TSA kit to detect low abundance transcripts. Therefore, this method provides a linear relation between signal intensity and targeted mRNA amounts for quantitative analysis. It allows processing 100-200 slides simultaneously. It works well for different developmental stages of embryos. Most developmental studies of gene expression use whole embryos and non-radioactive approaches(8,9), in part because embryonic tissue is more fragile than adult tissue, with less cohesion between cells, making it difficult to see boundaries between cell populations with tissue sections. In contrast, our radioactive approach, due to the larger range of sensitivity, is able to obtain higher contrast in resolution of gene expression between tissue regions, making it easier to see boundaries between populations. Using this method, researchers could reveal the possible significance of a newly identified gene, and further predict the function of the gene of interest.

Public goods and procreation.

Procreation is the ultimate public goods problem. Each new child affects the welfare of many other people, and some (but not all) children produce uncompensated value that future people will enjoy. This essay addresses challenges that arise if we think of procreation and parenting as public goods. These include whether individual choices are likely to lead to a socially desirable outcome, and whether changes in laws, social norms, or access to genetic engineering and embryo selection might improve the aggregate outcome of our reproductive choices.

Evolutionary characters, phenotypes and ontologies: curating data from the systematic biology literature.

BACKGROUND: The wealth of phenotypic descriptions documented in the published articles, monographs, and dissertations of phylogenetic systematics is traditionally reported in a free-text format, and it is therefore largely inaccessible for linkage to biological databases for genetics, development, and phenotypes, and difficult to manage for large-scale integrative work. The Phenoscape project aims to represent these complex and detailed descriptions with rich and formal semantics that are amenable to computation and integration with phenotype data from other fields of biology. This entails reconceptualizing the traditional free-text characters into the computable Entity-Quality (EQ) formalism using ontologies. METHODOLOGY/PRINCIPAL FINDINGS: We used ontologies and the EQ formalism to curate a collection of 47 phylogenetic studies on ostariophysan fishes (including catfishes, characins, minnows, knifefishes) and their relatives with the goal of integrating these complex phenotype descriptions with information from an existing model organism database (zebrafish, http://zfin.org). We developed a curation workflow for the collection of character, taxonomic and specimen data from these publications. A total of 4,617 phenotypic characters (10,512 states) for 3,449 taxa, primarily species, were curated into EQ formalism (for a total of 12,861 EQ statements) using anatomical and taxonomic terms from teleost-specific ontologies (Teleost Anatomy Ontology and Teleost Taxonomy Ontology) in combination with terms from a quality ontology (Phenotype and Trait Ontology). Standards and guidelines for consistently and accurately representing phenotypes were developed in response to the challenges that were evident from two annotation experiments and from feedback from curators. CONCLUSIONS/SIGNIFICANCE: The challenges we encountered and many of the curation standards and methods for improving consistency that we developed are generally applicable to any effort to represent phenotypes using ontologies. This is because an ontological representation of the detailed variations in phenotype, whether between mutant or wildtype, among individual humans, or across the diversity of species, requires a process by which a precise combination of terms from domain ontologies are selected and organized according to logical relations. The efficiencies that we have developed in this process will be useful for any attempt to annotate complex phenotypic descriptions using ontologies. We also discuss some ramifications of EQ representation for the domain of systematics.

A unified anatomy ontology of the vertebrate skeletal system.

The skeleton is of fundamental importance in research in comparative vertebrate morphology, paleontology, biomechanics, developmental biology, and systematics. Motivated by research questions that require computational access to and comparative reasoning across the diverse skeletal phenotypes of vertebrates, we developed a module of anatomical concepts for the skeletal system, the Vertebrate Skeletal Anatomy Ontology (VSAO), to accommodate and unify the existing skeletal terminologies for the species-specific (mouse, the frog Xenopus, zebrafish) and multispecies (teleost, amphibian) vertebrate anatomy ontologies. Previous differences between these terminologies prevented even simple queries across databases pertaining to vertebrate morphology. This module of upper-level and specific skeletal terms currently includes 223 defined terms and 179 synonyms that integrate skeletal cells, tissues, biological processes, organs (skeletal elements such as bones and cartilages), and subdivisions of the skeletal system. The VSAO is designed to integrate with other ontologies, including the Common Anatomy Reference Ontology (CARO), Gene Ontology (GO), Uberon, and Cell Ontology (CL), and it is freely available to the community to be updated with additional terms required for research. Its structure accommodates anatomical variation among vertebrate species in development, structure, and composition. Annotation of diverse vertebrate phenotypes with this ontology will enable novel inquiries across the full spectrum of phenotypic diversity.

Transgenerational Effects of Early Life Starvation on Growth, Reproduction, and Stress Resistance in Caenorhabditis elegans.

Starvation during early development can have lasting effects that influence organismal fitness and disease risk. We characterized the long-term phenotypic consequences of starvation during early larval development in Caenorhabditis elegans to determine potential fitness effects and develop it as a model for mechanistic studies. We varied the amount of time that larvae were developmentally arrested by starvation after hatching ("L1 arrest"). Worms recovering from extended starvation grew slowly, taking longer to become reproductive, and were smaller as adults. Fecundity was also reduced, with the smallest individuals most severely affected. Feeding behavior was impaired, possibly contributing to deficits in growth and reproduction. Previously starved larvae were more sensitive to subsequent starvation, suggesting decreased fitness even in poor conditions. We discovered that smaller larvae are more resistant to heat, but this correlation does not require passage through L1 arrest. The progeny of starved animals were also adversely affected: Embryo quality was diminished, incidence of males was increased, progeny were smaller, and their brood size was reduced. However, the progeny and grandprogeny of starved larvae were more resistant to starvation. In addition, the progeny, grandprogeny, and great-grandprogeny were more resistant to heat, suggesting epigenetic inheritance of acquired resistance to starvation and heat. Notably, such resistance was inherited exclusively from individuals most severely affected by starvation in the first generation, suggesting an evolutionary bet-hedging strategy. In summary, our results demonstrate that starvation affects a variety of life-history traits in the exposed animals and their descendants, some presumably reflecting fitness costs but others potentially adaptive.

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.

Interception of host angiogenic signalling limits mycobacterial growth.

Pathogenic mycobacteria induce the formation of complex cellular aggregates called granulomas that are the hallmark of tuberculosis. Here we examine the development and consequences of vascularization of the tuberculous granuloma in the zebrafish-Mycobacterium marinum infection model, which is characterized by organized granulomas with necrotic cores that bear striking resemblance to those of human tuberculosis. Using intravital microscopy in the transparent larval zebrafish, we show that granuloma formation is intimately associated with angiogenesis. The initiation of angiogenesis in turn coincides with the generation of local hypoxia and transcriptional induction of the canonical pro-angiogenic molecule Vegfaa. Pharmacological inhibition of the Vegf pathway suppresses granuloma-associated angiogenesis, reduces infection burden and limits dissemination. Moreover, anti-angiogenic therapies synergize with the first-line anti-tubercular antibiotic rifampicin, as well as with the antibiotic metronidazole, which targets hypoxic bacterial populations. Our data indicate that mycobacteria induce granuloma-associated angiogenesis, which promotes mycobacterial growth and increases spread of infection to new tissue sites. We propose the use of anti-angiogenic agents, now being used in cancer regimens, as a host-targeting tuberculosis therapy, particularly in extensively drug-resistant disease for which current antibiotic regimens are largely ineffective.

CPAG: software for leveraging pleiotropy in GWAS to reveal similarity between human traits links plasma fatty acids and intestinal inflammation.

Meta-analyses of genome-wide association studies (GWAS) have demonstrated that the same genetic variants can be associated with multiple diseases and other complex traits. We present software called CPAG (Cross-Phenotype Analysis of GWAS) to look for similarities between 700 traits, build trees with informative clusters, and highlight underlying pathways. Clusters are consistent with pre-defined groups and literature-based validation but also reveal novel connections. We report similarity between plasma palmitoleic acid and Crohn's disease and find that specific fatty acids exacerbate enterocolitis in zebrafish. CPAG will become increasingly powerful as more genetic variants are uncovered, leading to a deeper understanding of complex traits. CPAG is freely available at www.sourceforge.net/projects/CPAG/.

Structure of the polyisoprenyl-phosphate glycosyltransferase GtrB and insights into the mechanism of catalysis.

The attachment of a sugar to a hydrophobic polyisoprenyl carrier is the first step for all extracellular glycosylation processes. The enzymes that perform these reactions, polyisoprenyl-glycosyltransferases (PI-GTs) include dolichol phosphate mannose synthase (DPMS), which generates the mannose donor for glycosylation in the endoplasmic reticulum. Here we report the 3.0 Å resolution crystal structure of GtrB, a glucose-specific PI-GT from Synechocystis, showing a tetramer in which each protomer contributes two helices to a membrane-spanning bundle. The active site is 15 Å from the membrane, raising the question of how water-soluble and membrane-embedded substrates are brought into apposition for catalysis. A conserved juxtamembrane domain harbours disease mutations, which compromised activity in GtrB in vitro and in human DPM1 tested in zebrafish. We hypothesize a role of this domain in shielding the polyisoprenyl-phosphate for transport to the active site. Our results reveal the basis of PI-GT function, and provide a potential molecular explanation for DPM1-related disease.

NgBR is essential for endothelial cell glycosylation and vascular development.

NgBR is a transmembrane protein identified as a Nogo-B-interacting protein and recently has been shown to be a subunit required for cis-prenyltransferase (cisPTase) activity. To investigate the integrated role of NgBR in vascular development, we have characterized endothelial-specific NgBR knockout embryos. Here, we show that endothelial-specific NgBR knockout results in embryonic lethality due to vascular development defects in yolk sac and embryo proper. Loss of NgBR in endothelial cells reduces proliferation and promotes apoptosis of the cells largely through defects in the glycosylation of key endothelial proteins including VEGFR2, VE-cadherin, and CD31, and defective glycosylation can be rescued by treatment with the end product of cisPTase activity, dolichol phosphate. Moreover, NgBR functions in endothelial cells during embryogenesis are Nogo-B independent. These data uniquely show the importance of NgBR and protein glycosylation during vascular development.