Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

Color-Flu Fluorescent Reporter Influenza A Viruses Allow for In Vivo Studies of Innate Immune Function in Zebrafish

Influenza virus infection can cause severe respiratory disease and is estimated to cause millions of illnesses annually. Studies on the contribution of the innate immune response to influenza A virus (IAV) to viral pathogenesis may yield new antiviral strategies. Zebrafish larvae are useful models for studying the innate immune response to pathogens, including IAV, in vivo. Here, we demonstrate how Color-flu, four fluorescent IAV strains originally developed for mice, can be used to study the host response to infection by simultaneously monitoring infected cells, neutrophils, and macrophages in vivo. Using this model, we show how the angiotensin-converting enzyme inhibitor, ramipril, and mitophagy inhibitor, MDIVI-1, improved survival, decreased viral burden, and improved the respiratory burst response to IAV infection. The Color-flu zebrafish larvae model of IAV infection is complementary to other models where the dynamics of infection and the response of innate immune cells can be visualized in a transparent host in vivo.

Multi-spectral Fluorescent Reporter Influenza A Viruses Allow for in vivo Studies of Innate Immune Function in Zebrafish

Influenza virus infection can cause severe respiratory disease and is estimated to cause millions of illnesses annually. Studies of the contribution of the innate immune response to influenza A virus (IAV) to viral pathogenesis may yield new antiviral strategies. Zebrafish larvae are useful models to study the innate immune response to pathogens, including IAV, in vivo. Here, we demonstrate how Color-flu, four fluorescent IAV strains originally developed for mice, can be used to study host-virus interactions by simultaneously monitoring virus particles, neutrophils, and macrophages in vivo. Using this model, we show how the angiotensin-converting enzyme inhibitor, ramipril, and mitophagy inhibitor, MDIVI-1, improved survival, decreased viral burden, and improved the respiratory burst response to IAV infection. The Color-flu zebrafish model of IAV infection is complementary to other models as it is the only model where interactions between virus particles and host cells in an intact vertebrate can be visualized in vivo.

medna-metadata: an open-source data management system for tracking environmental DNA samples and metadata

MOTIVATION

Environmental DNA (eDNA), as a rapidly expanding research field, stands to benefit from shared resources including sampling protocols, study designs, discovered sequences, and taxonomic assignments to sequences. High-quality community shareable eDNA resources rely heavily on comprehensive metadata documentation that captures the complex workflows covering field sampling, molecular biology lab work, and bioinformatic analyses. There are limited sources that provide documentation of database development on comprehensive metadata for eDNA and these workflows and no open-source software.

RESULTS

We present medna-metadata, an open-source, modular system that aligns with Findable, Accessible, Interoperable, and Reusable guiding principles that support scholarly data reuse and the database and application development of a standardized metadata collection structure that encapsulates critical aspects of field data collection, wet lab processing, and bioinformatic analysis. Medna-metadata is showcased with metabarcoding data from the Gulf of Maine (Polinski et al., 2019).

AVAILABILITY AND IMPLEMENTATION

The source code of the medna-metadata web application is hosted on GitHub (https://github.com/Maine-eDNA/medna-metadata). Medna-metadata is a docker-compose installable package. Documentation can be found at https://medna-metadata.readthedocs.io/en/latest/?badge=latest. The application is implemented in Python, PostgreSQL and PostGIS, RabbitMQ, and NGINX, with all major browsers supported. A demo can be found at https://demo.metadata.maine-edna.org/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cetylpyridinium chloride (CPC) reduces zebrafish mortality from influenza infection: Super-resolution microscopy reveals CPC interference with multiple protein interactions with phosphatidylinositol 4,5-bisphosphate in immune function

The COVID-19 pandemic raises significance for a potential influenza therapeutic compound, cetylpyridinium chloride (CPC), which has been extensively used in personal care products as a positively-charged quaternary ammonium antibacterial agent. CPC is currently in clinical trials to assess its effects on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) morbidity. Two published studies have provided mouse and human data indicating that CPC may alleviate influenza infection, and here we show that CPC (0.1 μM, 1 h) reduces zebrafish mortality and viral load following influenza infection. However, CPC mechanisms of action upon viral-host cell interaction are currently unknown. We have utilized super-resolution fluorescence photoactivation localization microscopy to probe the mode of CPC action. Reduction in density of influenza viral protein hemagglutinin (HA) clusters is known to reduce influenza infectivity: here, we show that CPC (at non-cytotoxic doses, 5-10 μM) reduces HA density and number of HA molecules per cluster within the plasma membrane of NIH-3T3 mouse fibroblasts. HA is known to colocalize with the negatively-charged mammalian lipid phosphatidylinositol 4,5-bisphosphate (PIP2); here, we show that nanoscale co-localization of HA with the PIP2-binding Pleckstrin homology (PH) reporter in the plasma membrane is diminished by CPC. CPC also dramatically displaces the PIP2-binding protein myristoylated alanine-rich C-kinase substrate (MARCKS) from the plasma membrane of rat RBL-2H3 mast cells; this disruption of PIP2 is correlated with inhibition of mast cell degranulation. Together, these findings offer a PIP2-focused mechanism underlying CPC disruption of influenza and suggest potential pharmacological use of this drug as an influenza therapeutic to reduce global deaths from viral disease.

Beneficial impacts of neuromuscular electrical stimulation on muscle structure and function in the zebrafish model of Duchenne muscular dystrophy

Neuromuscular electrical stimulation (NMES) allows activation of muscle fibers in the absence of voluntary force generation. NMES could have the potential to promote muscle homeostasis in the context of muscle disease, but the impacts of NMES on diseased muscle are not well understood. We used the zebrafish Duchenne muscular dystrophy (dmd) mutant and a longitudinal design to elucidate the consequences of NMES on muscle health. We designed four neuromuscular stimulation paradigms loosely based on weightlifting regimens. Each paradigm differentially affected neuromuscular structure, function, and survival. Only endurance neuromuscular stimulation (eNMES) improved all outcome measures. We found that eNMES improves muscle and neuromuscular junction morphology, swimming, and survival. Heme oxygenase and integrin alpha7 are required for eNMES-mediated improvement. Our data indicate that neuromuscular stimulation can be beneficial, suggesting that the right type of activity may benefit patients with muscle disease.

The MAPK/ERK Pathway and the Role of DUSP1 in JCPyV Infection of Primary Astrocytes

JC polyomavirus (JCPyV) is a neuroinvasive pathogen causing a fatal, demyelinating disease of the central nervous system (CNS) known as progressive multifocal leukoencephalopathy (PML). Within the CNS, JCPyV predominately targets two cell types: oligodendrocytes and astrocytes. The underlying mechanisms of astrocytic infection are poorly understood, yet recent findings suggest critical differences in JCPyV infection of primary astrocytes compared to a widely studied immortalized cell model. RNA sequencing was performed in primary normal human astrocytes (NHAs) to analyze the transcriptomic profile that emerges during JCPyV infection. Through a comparative analysis, it was validated that JCPyV requires the mitogen-activated protein kinase, extracellular signal-regulated kinase (MAPK/ERK) pathway, and additionally requires the expression of dual-specificity phosphatases (DUSPs). Specifically, the expression of DUSP1 is needed to establish a successful infection in NHAs, yet this was not observed in an immortalized cell model of JCPyV infection. Additional analyses demonstrated immune activation uniquely observed in NHAs. These results support the hypothesis that DUSPs within the MAPK/ERK pathway impact viral infection and influence potential downstream targets and cellular pathways. Collectively, this research implicates DUSP1 in JCPyV infection of primary human astrocytes, and most importantly, further resolves the signaling events that lead to successful JCPyV infection in the CNS.

Meta-analysis: Congruence of genomic and phenotypic differentiation across diverse natural study systems

Linking genotype to phenotype is a primary goal for understanding the genomic underpinnings of evolution. However, little work has explored whether patterns of linked genomic and phenotypic differentiation are congruent across natural study systems and traits. Here, we investigate such patterns with a meta-analysis of studies examining population-level differentiation at subsets of loci and traits putatively responding to divergent selection. We show that across the 31 studies (88 natural population-level comparisons) we examined, there was a moderate (R 2 = 0.39) relationship between genomic differentiation (F ST ) and phenotypic differentiation (P ST ) for loci and traits putatively under selection. This quantitative relationship between P ST and F ST for loci under selection in diverse taxa provides broad context and cross-system predictions for genomic and phenotypic adaptation by natural selection in natural populations. This context may eventually allow for more precise ideas of what constitutes "strong" differentiation, predictions about the effect size of loci, comparisons of taxa evolving in nonparallel ways, and more. On the other hand, links between P ST and F ST within studies were very weak, suggesting that much work remains in linking genomic differentiation to phenotypic differentiation at specific phenotypes. We suggest that linking genotypes to specific phenotypes can be improved by correlating genomic and phenotypic differentiation across a spectrum of diverging populations within a taxon and including wide coverage of both genomes and phenomes.

Modeling Virus-Induced Inflammation in Zebrafish: A Balance Between Infection Control and Excessive Inflammation

The inflammatory response to viral infection in humans is a dynamic process with complex cell interactions that are governed by the immune system and influenced by both host and viral factors. Due to this complexity, the relative contributions of the virus and host factors are best studied in vivo using animal models. In this review, we describe how the zebrafish (Danio rerio) has been used as a powerful model to study host-virus interactions and inflammation by combining robust forward and reverse genetic tools with in vivo imaging of transparent embryos and larvae. The innate immune system has an essential role in the initial inflammatory response to viral infection. Focused studies of the innate immune response to viral infection are possible using the zebrafish model as there is a 4-6 week timeframe during development where they have a functional innate immune system dominated by neutrophils and macrophages. During this timeframe, zebrafish lack a functional adaptive immune system, so it is possible to study the innate immune response in isolation. Sequencing of the zebrafish genome has revealed significant genetic conservation with the human genome, and multiple studies have revealed both functional conservation of genes, including those critical to host cell infection and host cell inflammatory response. In addition to studying several fish viruses, zebrafish infection models have been developed for several human viruses, including influenza A, noroviruses, chikungunya, Zika, dengue, herpes simplex virus type 1, Sindbis, and hepatitis C virus. The development of these diverse viral infection models, coupled with the inherent strengths of the zebrafish model, particularly as it relates to our understanding of macrophage and neutrophil biology, offers opportunities for far more intensive studies aimed at understanding conserved host responses to viral infection. In this context, we review aspects relating to the evolution of innate immunity, including the evolution of viral pattern recognition receptors, interferons and interferon receptors, and non-coding RNAs.

Chronic cortisol exposure in early development leads to neuroendocrine dysregulation in adulthood

Objective

Chronic early life stress can affect development of the neuroendocrine stress system, leading to its persistent dysregulation and consequently increased disease risk in adulthood. One contributing factor is thought to be epigenetic programming in response to chronic cortisol exposure during early development. We have previously shown that zebrafish embryos treated chronically with cortisol develop into adults with constitutively elevated whole-body cortisol and aberrant immune gene expression. Here we further characterize that phenotype by assessing persistent effects of the treatment on cortisol tissue distribution and dynamics, chromatin accessibility, and activities of glucocorticoid-responsive regulatory genes klf9 and fkbp5. To that end cortisol levels in different tissues of fed and fasted adults were measured using ELISA, open chromatin in adult blood cells was mapped using ATAC-seq, and gene activity in adult blood and brain cells was measured using qRT-PCR.

Results

Adults derived from cortisol-treated embryos have elevated whole-body cortisol with aberrantly regulated tissue distribution and dynamics that correlate with differential activity of klf9 and fkbp5 in blood and brain.

Modulation of TNFα Activity by the microRNA Let-7 Coordinates Zebrafish Heart Regeneration

The adult zebrafish is capable of regenerating heart muscle, resolving collagen tissue, and fully restoring heart function throughout its life. In this study, we show that the highly upregulated, epicardium-enriched microRNA let-7i functions in wound closure and cardiomyocyte proliferation. RNA sequencing experiments identified upregulated expression of members of the tumor necrosis factor (TNF) signaling pathway in the absence of let-7. Importantly, co-suppression of TNF and let-7 activity rescued epicardium migration and cardiomyocyte proliferation defects induced by depletion of let-7 alone. Sensitizing animals to low levels of TNF activity before injury culminated in repressed cardiomyocyte proliferation and wound closure defects, suggesting that levels of inflammation at the onset of injury are critical for heart regeneration. Our studies indicate that injury-induced reduction in TNF signaling by let-7 in the epicardium creates a pro-regenerative environment for cardiomyocyte proliferation during adult heart regeneration.

Profiling microRNA expression in Atlantic killifish (Fundulus heteroclitus) gill and responses to arsenic and hyperosmotic stress

The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24 h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24 h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.

A role for the CXCR4-CXCL12 axis in the little skate, Leucoraja erinacea

The interaction between C-X-C chemokine receptor type 4 (CXCR4) and its cognate ligand C-X-C motif chemokine ligand 12 (CXCL12) plays a critical role in regulating hematopoietic stem cell activation and subsequent cellular mobilization. Extensive studies of these genes have been conducted in mammals, but much less is known about the expression and function of CXCR4 and CXCL12 in non-mammalian vertebrates. In the present study, we identify simultaneous expression of CXCR4 and CXCL12 orthologs in the epigonal organ (the primary hematopoietic tissue) of the little skate, Leucoraja erinacea. Genetic and phylogenetic analyses were functionally supported by significant mobilization of leukocytes following administration of Plerixafor, a CXCR4 antagonist and clinically important drug. Our results provide evidence that, as in humans, Plerixafor disrupts CXCR4/CXCL12 binding in the little skate, facilitating release of leukocytes into the bloodstream. Our study illustrates the value of the little skate as a model organism, particularly in studies of hematopoiesis and potentially for preclinical research on hematological and vascular disorders.

RegenDbase: a comparative database of noncoding RNA regulation of tissue regeneration circuits across multiple taxa

Regeneration is an endogenous process of tissue repair that culminates in complete restoration of tissue and organ function. While regenerative capacity in mammals is limited to select tissues, lower vertebrates like zebrafish and salamanders are endowed with the capacity to regenerate entire limbs and most adult tissues, including heart muscle. Numerous profiling studies have been conducted using these research models in an effort to identify the genetic circuits that accompany tissue regeneration. Most of these studies, however, are confined to an individual injury model and/or research organism and focused primarily on protein encoding transcripts. Here we describe RegenDbase, a new database with the functionality to compare and contrast gene regulatory pathways within and across tissues and research models. RegenDbase combines pipelines that integrate analysis of noncoding RNAs in combination with protein encoding transcripts. We created RegenDbase with a newly generated comprehensive dataset for adult zebrafish heart regeneration combined with existing microarray and RNA-sequencing studies on multiple injured tissues. In this current release, we detail microRNA-mRNA regulatory circuits and the biological processes these interactions control during the early stages of heart regeneration. Moreover, we identify known and putative novel lncRNAs and identify their potential target genes based on proximity searches. We postulate that these candidate factors underscore robust regenerative capacity in lower vertebrates. RegenDbase provides a systems-level analysis of tissue regeneration genetic circuits across injury and animal models and addresses the growing need to understand how noncoding RNAs influence these changes in gene expression.

Global gene profiling of aging lungs in Atp8b1 mutant mice

Objective

Recent studies implicate cardiolipin oxidation in several age-related diseases. Atp8b1 encoding Type 4 P-type ATPases is a cardiolipin transporter. Mutation in Atp8b1 gene or inflammation of the lungs impairs the capacity of Atp8b1 to clear cardiolipin from lung fluid. However, the link between Atp8b1 mutation and age-related gene alteration is unknown. Therefore, we investigated how Atp8b1 mutation alters age-related genes.

Methods

We performed Affymetrix gene profiling of lungs isolated from young (7-9 wks, n=6) and aged (14 months, 14 M, n=6) C57BL/6 and Atp8b1 mutant mice. In addition, Ingenuity Pathway Analysis (IPA) was performed. Differentially expressed genes were validated by quantitative real-time PCR (qRT-PCR).

Results

Global transcriptome analysis revealed 532 differentially expressed genes in Atp8b1 lungs, 157 differentially expressed genes in C57BL/6 lungs, and 37 overlapping genes. IPA of age-related genes in Atp8b1 lungs showed enrichment of Xenobiotic metabolism and Nrf2-mediated signaling pathways. The increase in Adamts2 and Mmp13 transcripts in aged Atp8b1 lungs was validated by qRT-PCR. Similarly, the decrease in Col1a1 and increase in Cxcr6 transcripts was confirmed in both Atp8b1 mutant and C57BL/6 lungs.

Conclusion

Based on transcriptome profiling, our study indicates that Atp8b1 mutant mice may be susceptible to age-related lung diseases.

A retrospective comparison of neoadjuvant chemoradiotherapy regimens for locally advanced esophageal cancer

Preoperative chemoradiotherapy (CRT) with carboplatin/paclitaxel has been shown to increase survival in patients with esophageal cancer, including gastroesophageal junction (GE) junction cancer, over surgery alone; however, there have been no studies comparing the different neoadjuvant CRT regimens. We retrospectively evaluated the long-term results of trimodality therapy for patients with locally advanced esophageal cancer treated on several chemotherapy regimens. Between 1999 and 2014, 215 patients with locally advanced esophageal cancer underwent neoadjuvant CRT followed by surgical resection. The median age was 62 years (range 21-84), 80.5% were men and 86% had adenocarcinoma. The following chemotherapy regimens were administered: cisplatin/5FU (14.9%), cisplatin/irinotecan (35.8%), carboplatin/paclitaxel (35.8%), and other (9.7%). The majority of patients (92.1%) received a radiation dose of 50.4 Gy. Predictors of toxicities and surgical complications were assessed using logistic regression. Overall survival (OS) and recurrence-free survival (RFS) were estimated using the Kaplan-Meier method and proportional hazards regression was used to model time-to-event outcomes. The median follow-up among surviving patients was 4.1 years (range 0.4,13). The median OS was 3.0 years from time of diagnosis and OS was 36.8% at 5 years. RFS was 34.9% at 5 years. After neoadjuvant CRT, 34.7% of patients achieved a pathologic complete response including 60.7% of squamous cell carcinoma patients and 18.4% of adenocarcinoma patients (P < 0.001) and 66% were downstaged. Of the variables examined, pathologic stage, preoperative baseline cardiac comorbidity, postoperative cardiac or pulmonary complications, and chemotherapy regimen were associated with OS. Using cisplatin and 5FU as the reference regimen, patients treated with carboplatin/paclitaxel had significantly improved OS (HR = 0.47, P = 0.017 after adjusting for surgery type, radiation modality, baseline cardiac comorbidity, and preoperative stage) with 5-year OS rate of 66%. The most common surgical complications were cardiac in 61 patients (28.5%) and pulmonary in 52 patients (24.3%). Cardiac complications were associated with age (OR 1.05, P = 0.007) and cardiac comorbidity (OR 2.6, P = 0.02) and pulmonary complications with female gender (OR 3.98, P < 0.001). Forty-four patients (20.5%) required readmission within 30 days of discharge, and readmission was associated with cardiac comorbidity (OR 2.7, P = 0.017). Three patients died within 30 days of surgery. We observed an association between neoadjuvant carboplatin/paclitaxel and improved overall survival that requires confirmation in a prospective randomized trial.

The landscape of extreme genomic variation in the highly adaptable Atlantic killifish

Understanding and predicting the fate of populations in changing environments require knowledge about the mechanisms that support phenotypic plasticity and the adaptive value and evolutionary fate of genetic variation within populations. Atlantic killifish (Fundulus heteroclitus) exhibit extensive phenotypic plasticity that supports large population sizes in highly fluctuating estuarine environments. Populations have also evolved diverse local adaptations. To yield insights into the genomic variation that supports their adaptability, we sequenced a reference genome and 48 additional whole genomes from a wild population. Evolution of genes associated with cell cycle regulation and apoptosis is accelerated along the killifish lineage, which is likely tied to adaptations for life in highly variable estuarine environments. Genome-wide standing genetic variation, including nucleotide diversity and copy number variation, is extremely high. The highest diversity genes are those associated with immune function and olfaction, whereas genes under greatest evolutionary constraint are those associated with neurological, developmental, and cytoskeletal functions. Reduced genetic variation is detected for tight junction proteins, which in killifish regulate paracellular permeability that supports their extreme physiological flexibility. Low-diversity genes engage in more regulatory interactions than high-diversity genes, consistent with the influence of pleiotropic constraint on molecular evolution. High genetic variation is crucial for continued persistence of species given the pace of contemporary environmental change. Killifish populations harbor among the highest levels of nucleotide diversity yet reported for a vertebrate species, and thus may serve as a useful model system for studying evolutionary potential in variable and changing environments.

Advancing Exposure Science through Chemical Data Curation and Integration in the Comparative Toxicogenomics Database

BACKGROUND

Exposure science studies the interactions and outcomes between environmental stressors and human or ecological receptors. To augment its role in understanding human health and the exposome, we aimed to centralize and integrate exposure science data into the broader biological framework of the Comparative Toxicogenomics Database (CTD), a public resource that promotes understanding of environmental chemicals and their effects on human health.

OBJECTIVES

We integrated exposure data within the CTD to provide a centralized, freely available resource that facilitates identification of connections between real-world exposures, chemicals, genes/proteins, diseases, biological processes, and molecular pathways.

METHODS

We developed a manual curation paradigm that captures exposure data from the scientific literature using controlled vocabularies and free text within the context of four primary exposure concepts: stressor, receptor, exposure event, and exposure outcome. Using data from the Agricultural Health Study, we have illustrated the benefits of both centralization and integration of exposure information with CTD core data.

RESULTS

We have described our curation process, demonstrated how exposure data can be accessed and analyzed in the CTD, and shown how this integration provides a broad biological context for exposure data to promote mechanistic understanding of environmental influences on human health.

CONCLUSIONS

Curation and integration of exposure data within the CTD provides researchers with new opportunities to correlate exposures with human health outcomes, to identify underlying potential molecular mechanisms, and to improve understanding about the exposome.

CITATION

Grondin CJ, Davis AP, Wiegers TC, King BL, Wiegers JA, Reif DM, Hoppin JA, Mattingly CJ. 2016. Advancing exposure science through chemical data curation and integration in the Comparative Toxicogenomics Database. Environ Health Perspect 124:1592-1599; http://dx.doi.org/10.1289/EHP174.

The Comparative Toxicogenomics Database: update 2017

The Comparative Toxicogenomics Database (CTD; http://ctdbase.org/) provides information about interactions between chemicals and gene products, and their relationships to diseases. Core CTD content (chemical-gene, chemical-disease and gene-disease interactions manually curated from the literature) are integrated with each other as well as with select external datasets to generate expanded networks and predict novel associations. Today, core CTD includes more than 30.5 million toxicogenomic connections relating chemicals/drugs, genes/proteins, diseases, taxa, Gene Ontology (GO) annotations, pathways, and gene interaction modules. In this update, we report a 33% increase in our core data content since 2015, describe our new exposure module (that harmonizes exposure science information with core toxicogenomic data) and introduce a novel dataset of GO-disease inferences (that identify common molecular underpinnings for seemingly unrelated pathologies). These advancements centralize and contextualize real-world chemical exposures with molecular pathways to help scientists generate testable hypotheses in an effort to understand the etiology and mechanisms underlying environmentally influenced diseases.

Cortisol-treated zebrafish embryos develop into pro-inflammatory adults with aberrant immune gene regulation

Chronic early-life stress increases adult susceptibility to numerous health problems linked to chronic inflammation. One way that this may occur is via glucocorticoid-induced developmental programming. To gain insight into such programming we treated zebrafish embryos with cortisol and examined the effects on both larvae and adults. Treated larvae had elevated whole-body cortisol and glucocorticoid signaling, and upregulated genes associated with defense response and immune system processes. In adulthood the treated fish maintained elevated basal cortisol levels in the absence of exogenous cortisol, and constitutively mis-expressed genes involved in defense response and its regulation. Adults derived from cortisol-treated embryos displayed defective tailfin regeneration, heightened basal expression of pro-inflammatory genes, and failure to appropriately regulate those genes following injury or immunological challenge. These results support the hypothesis that chronically elevated glucocorticoid signaling early in life directs development of a pro-inflammatory adult phenotype, at the expense of immunoregulation and somatic regenerative capacity.

A Conserved MicroRNA Regulatory Circuit Is Differentially Controlled during Limb/Appendage Regeneration

Background

Although regenerative capacity is evident throughout the animal kingdom, it is not equally distributed throughout evolution. For instance, complex limb/appendage regeneration is muted in mammals but enhanced in amphibians and teleosts. The defining characteristic of limb/appendage regenerative systems is the formation of a dedifferentiated tissue, termed blastema, which serves as the progenitor reservoir for regenerating tissues. In order to identify a genetic signature that accompanies blastema formation, we employ next-generation sequencing to identify shared, differentially regulated mRNAs and noncoding RNAs in three different, highly regenerative animal systems: zebrafish caudal fins, bichir pectoral fins and axolotl forelimbs.

Results

These studies identified a core group of 5 microRNAs (miRNAs) that were commonly upregulated and 5 miRNAs that were commonly downregulated, as well as 4 novel tRNAs fragments with sequences conserved with humans. To understand the potential function of these miRNAs, we built a network of 1,550 commonly differentially expressed mRNAs that had functional relationships to 11 orthologous blastema-associated genes. As miR-21 was the most highly upregulated and most highly expressed miRNA in all three models, we validated the expression of known target genes, including the tumor suppressor, pdcd4, and TGFβ receptor subunit, tgfbr2 and novel putative target genes such as the anti-apoptotic factor, bcl2l13, Choline kinase alpha, chka and the regulator of G-protein signaling, rgs5.

Conclusions

Our extensive analysis of RNA-seq transcriptome profiling studies in three regenerative animal models, that diverged in evolution ~420 million years ago, reveals a common miRNA-regulated genetic network of blastema genes. These comparative studies extend our current understanding of limb/appendage regeneration by identifying previously unassociated blastema genes and the extensive regulation by miRNAs, which could serve as a foundation for future functional studies to examine the process of natural cellular reprogramming in an injury context.

Abnormal Osmotic Avoidance Behavior in C. elegans Is Associated with Increased Hypertonic Stress Resistance and Improved Proteostasis

Protein function is controlled by the cellular proteostasis network. Proteostasis is energetically costly and those costs must be balanced with the energy needs of other physiological functions. Hypertonic stress causes widespread protein damage in C. elegans. Suppression and management of protein damage is essential for optimal survival under hypertonic conditions. ASH chemosensory neurons allow C. elegans to detect and avoid strongly hypertonic environments. We demonstrate that mutations in osm-9 and osm-12 that disrupt ASH mediated hypertonic avoidance behavior or genetic ablation of ASH neurons are associated with enhanced survival during hypertonic stress. Improved survival is not due to altered systemic volume homeostasis or organic osmolyte accumulation. Instead, we find that osm-9(ok1677) mutant and osm-9(RNAi) worms exhibit reductions in hypertonicity induced protein damage in non-neuronal cells suggesting that enhanced proteostasis capacity may account for improved hypertonic stress resistance in worms with defects in osmotic avoidance behavior. RNA-seq analysis revealed that genes that play roles in managing protein damage are upregulated in osm-9(ok1677) worms. Our findings are consistent with a growing body of work demonstrating that intercellular communication between neuronal and non-neuronal cells plays a critical role in integrating cellular stress resistance with other organismal physiological demands and associated energy costs.

The Interrelationships of Placental Mammals and the Limits of Phylogenetic Inference

Placental mammals comprise three principal clades: Afrotheria (e.g., elephants and tenrecs), Xenarthra (e.g., armadillos and sloths), and Boreoeutheria (all other placental mammals), the relationships among which are the subject of controversy and a touchstone for debate on the limits of phylogenetic inference. Previous analyses have found support for all three hypotheses, leading some to conclude that this phylogenetic problem might be impossible to resolve due to the compounded effects of incomplete lineage sorting (ILS) and a rapid radiation. Here we show, using a genome scale nucleotide data set, microRNAs, and the reanalysis of the three largest previously published amino acid data sets, that the root of Placentalia lies between Atlantogenata and Boreoeutheria. Although we found evidence for ILS in early placental evolution, we are able to reject previous conclusions that the placental root is a hard polytomy that cannot be resolved. Reanalyses of previous data sets recover Atlantogenata + Boreoeutheria and show that contradictory results are a consequence of poorly fitting evolutionary models; instead, when the evolutionary process is better-modeled, all data sets converge on Atlantogenata. Our Bayesian molecular clock analysis estimates that marsupials diverged from placentals 157-170 Ma, crown Placentalia diverged 86-100 Ma, and crown Atlantogenata diverged 84-97 Ma. Our results are compatible with placental diversification being driven by dispersal rather than vicariance mechanisms, postdating early phases in the protracted opening of the Atlantic Ocean.

Calcium activated K⁺ channels in the electroreceptor of the skate confirmed by cloning. Details of subunits and splicing

Elasmobranchs detect small potentials using excitable cells of the ampulla of Lorenzini which have calcium-activated K(+) channels, first described in 1974. A distinctive feature of the outward current in voltage clamped ampullae is its apparent insensitivity to voltage. The sequence of a BK channel α isoform expressed in the ampulla of the skate was characterized. A signal peptide is present at the beginning of the gene. When compared to human isoform 1 (the canonical sequence), the largest difference was absence of a 59 amino acid region from the S8-S9 intra-cellular linker that contains the strex regulatory domain. The ampulla isoform was also compared with the isoform predicted in late skate embryos where strex was also absent. The BK voltage sensors were conserved in both skate isoforms. Differences between the skate and human BK channel included alternative splicing. Alternative splicing occurs at seven previously defined sites that are characteristic for BK channels in general and hair cells in particular. Skate BK sequences were highly similar to the Australian ghost shark and several other vertebrate species. Based on alignment of known BK sequences with the skate genome and transcriptome, there are at least two isoforms of Kcnma1α expressed in the skate. One of the β subunits (β4), which is known to decrease voltage sensitivity, was also identified in the skate genome and transcriptome and in the ampulla. These studies advance our knowledge of BK channels and suggest further studies in the ampulla and other excitable tissues.

A Uniform System for the Annotation of Vertebrate microRNA Genes and the Evolution of the Human microRNAome

Although microRNAs (miRNAs) are among the most intensively studied molecules of the past 20 years, determining what is and what is not a miRNA has not been straightforward. Here, we present a uniform system for the annotation and nomenclature of miRNA genes. We show that less than a third of the 1,881 human miRBase entries, and only approximately 16% of the 7,095 metazoan miRBase entries, are robustly supported as miRNA genes. Furthermore, we show that the human repertoire of miRNAs has been shaped by periods of intense miRNA innovation and that mature gene products show a very different tempo and mode of sequence evolution than star products. We establish a new open access database--MirGeneDB ( http://mirgenedb.org )--to catalog this set of miRNAs, which complements the efforts of miRBase but differs from it by annotating the mature versus star products and by imposing an evolutionary hierarchy upon this curated and consistently named repertoire.

A Fluorescence-Based Assay for Proteinuria Screening in Larval Zebrafish (Danio rerio)

Analysis of genes compromising the glomerular filtration barrier in rodent models using transgenic or knockdown approaches is time- and resource-consuming and often leads to unsatisfactory results. Therefore, it would be beneficial to have a selection tool indicating that your gene of interest is in fact associated with proteinuria. Zebrafish (Danio rerio) is a rapid screening tool to study effects in glomerular filtration barrier integrity after genetic manipulation. We use either injection of high-molecular-weight dextrans or a transgenic fluorescent fish line [Tg(l-fabp:DBP:EGFP)] expressing a vitamin D-binding protein fused with eGFP for indirect detection of proteinuria. A loss of high-molecular-weight proteins from the circulation of the fish into the urine can be identified by monitoring fluorescence intensity in the zebrafish eye. Paired with an optimized analysis method, this assay provides an effective screening solution to detect filtration barrier damage with proteinuria before moving to a mammalian system.

Gene Expression Changes Associated With the Developmental Plasticity of Sea Urchin Larvae in Response to Food Availability

Planktotrophic sea urchin larvae are developmentally plastic: in response to food scarcity, development of the juvenile rudiment is suspended and larvae instead develop elongated arms, thus increasing feeding capacity and extending larval life. Here, data are presented on the effect of different feeding regimes on gene expression in larvae of the green sea urchin Strongylocentrotus droebachiensis. These data indicate that during periods of starvation, larvae down-regulate genes involved in growth and metabolic activity while up-regulating genes involved in lipid transport, environmental sensing, and defense. Additionally, we show that starvation increases FoxO activity and that in well-fed larvae rapamycin treatment impedes rudiment growth, indicating that the latter requires TOR activity. These results suggest that the developmental plasticity of echinoplutei is regulated by genes known to control aging and longevity in other animals.

Natural selection canalizes expression variation of environmentally induced plasticity-enabling genes

Many organisms survive fluctuating and extreme environmental conditions by manifesting multiple distinct phenotypes during adulthood by means of developmental processes that enable phenotypic plasticity. We report on the discovery of putative plasticity-enabling genes that are involved in transforming the gill of the euryhaline teleost fish, Fundulus heteroclitus, from its freshwater to its seawater gill-type, a process that alters both morphology and function. Gene expression that normally enables osmotic plasticity is inhibited by arsenic. Gene sets defined by antagonistic interactions between arsenic and salinity show reduced transcriptional variation among individual fish, suggesting unusually accurate and precise regulatory control of these genes, consistent with the hypothesis that they participate in a canalized developmental response. We observe that natural selection acts to preserve canalized gene expression in populations of killifish that are most tolerant to abrupt salinity change and that these populations show the least variability in their transcription of genes enabling plasticity of the gill. We found that genes participating in this highly canalized and conserved plasticity-enabling response had significantly fewer and less complex associations with transcriptional regulators than genes that respond only to arsenic or salinity. Collectively these findings, which are drawn from the relationships between environmental challenge, plasticity, and canalization among populations, suggest that the selective processes that facilitate phenotypic plasticity do so by targeting the regulatory networks that gives rise to the response. These findings also provide a generalized, conceptual framework of how genes might interact with the environment and evolve toward the development of plastic traits.

The Comparative Toxicogenomics Database's 10th year anniversary: update 2015

Ten years ago, the Comparative Toxicogenomics Database (CTD; http://ctdbase.org/) was developed out of a need to formalize, harmonize and centralize the information on numerous genes and proteins responding to environmental toxic agents across diverse species. CTD's initial approach was to facilitate comparisons of nucleotide and protein sequences of toxicologically significant genes by curating these sequences and electronically annotating them with chemical terms from their associated references. Since then, however, CTD has vastly expanded its scope to robustly represent a triad of chemical–gene, chemical–disease and gene–disease interactions that are manually curated from the scientific literature by professional biocurators using controlled vocabularies, ontologies and structured notation. Today, CTD includes 24 million toxicogenomic connections relating chemicals/drugs, genes/proteins, diseases, taxa, phenotypes, Gene Ontology annotations, pathways and interaction modules. In this 10th year anniversary update, we outline the evolution of CTD, including our increased data content, new ‘Pathway View’ visualization tool, enhanced curation practices, pilot chemical–phenotype results and impending exposure data set. The prototype database originally described in our first report has transformed into a sophisticated resource used actively today to help scientists develop and test hypotheses about the etiologies of environmentally influenced diseases.

SkateBase, an elasmobranch genome project and collection of molecular resources for chondrichthyan fishes

Chondrichthyan fishes are a diverse class of gnathostomes that provide a valuable perspective on fundamental characteristics shared by all jawed and limbed vertebrates. Studies of phylogeny, species diversity, population structure, conservation, and physiology are accelerated by genomic, transcriptomic and protein sequence data. These data are widely available for many sarcopterygii (coelacanth, lungfish and tetrapods) and actinoptergii (ray-finned fish including teleosts) taxa, but limited for chondrichthyan fishes.  In this study, we summarize available data for chondrichthyes and describe resources for one of the largest projects to characterize one of these fish, Leucoraja erinacea, the little skate.  SkateBase ( http://skatebase.org) serves as the skate genome project portal linking data, research tools, and teaching resources.

Toward consilience in reptile phylogeny: miRNAs support an archosaur, not lepidosaur, affinity for turtles

Understanding the phylogenetic position of crown turtles (Testudines) among amniotes has been a source of particular contention. Recent morphological analyses suggest that turtles are sister to all other reptiles, whereas the vast majority of gene sequence analyses support turtles as being inside Diapsida, and usually as sister to crown Archosauria (birds and crocodilians). Previously, a study using microRNAs (miRNAs) placed turtles inside diapsids, but as sister to lepidosaurs (lizards and Sphenodon) rather than archosaurs. Here, we test this hypothesis with an expanded miRNA presence/absence dataset, and employ more rigorous criteria for miRNA annotation. Significantly, we find no support for a turtle + lepidosaur sister-relationship; instead, we recover strong support for turtles sharing a more recent common ancestor with archosaurs. We further test this result by analyzing a super-alignment of precursor miRNA sequences for every miRNA inferred to have been present in the most recent common ancestor of tetrapods. This analysis yields a topology that is fully congruent with our presence/absence analysis; our results are therefore in accordance with most gene sequence studies, providing strong, consilient molecular evidence from diverse independent datasets regarding the phylogenetic position of turtles.

A CTD-Pfizer collaboration: manual curation of 88,000 scientific articles text mined for drug-disease and drug-phenotype interactions

Improving the prediction of chemical toxicity is a goal common to both environmental health research and pharmaceutical drug development. To improve safety detection assays, it is critical to have a reference set of molecules with well-defined toxicity annotations for training and validation purposes. Here, we describe a collaboration between safety researchers at Pfizer and the research team at the Comparative Toxicogenomics Database (CTD) to text mine and manually review a collection of 88 629 articles relating over 1 200 pharmaceutical drugs to their potential involvement in cardiovascular, neurological, renal and hepatic toxicity. In 1 year, CTD biocurators curated 2 54 173 toxicogenomic interactions (1 52 173 chemical–disease, 58 572 chemical–gene, 5 345 gene–disease and 38 083 phenotype interactions). All chemical–gene–disease interactions are fully integrated with public CTD, and phenotype interactions can be downloaded. We describe Pfizer’s text-mining process to collate the articles, and CTD’s curation strategy, performance metrics, enhanced data content and new module to curate phenotype information. As well, we show how data integration can connect phenotypes to diseases. This curation can be leveraged for information about toxic endpoints important to drug safety and help develop testable hypotheses for drug–disease events. The availability of these detailed, contextualized, high-quality annotations curated from seven decades’ worth of the scientific literature should help facilitate new mechanistic screening assays for pharmaceutical compound survival. This unique partnership demonstrates the importance of resource sharing and collaboration between public and private entities and underscores the complementary needs of the environmental health science and pharmaceutical communities.

Database URL: http://ctdbase.org/

miRNAs: small genes with big potential in metazoan phylogenetics

microRNAs (miRNAs) are a key component of gene regulatory networks and have been implicated in the regulation of virtually every biological process found in multicellular eukaryotes. What makes them interesting from a phylogenetic perspective is the high conservation of primary sequence between taxa, their accrual in metazoan genomes through evolutionary time, and the rarity of secondary loss in most metazoan taxa. Despite these properties, the use of miRNAs as phylogenetic markers has not yet been discussed within a clear conceptual framework. Here we highlight five properties of miRNAs that underlie their utility in phylogenetics: 1) The processes of miRNA biogenesis enable the identification of novel miRNAs without prior knowledge of sequence; 2) The continuous addition of miRNA families to metazoan genomes through evolutionary time; 3) The low level of secondary gene loss in most metazoan taxa; 4) The low substitution rate in the mature miRNA sequence; and 5) The small probability of convergent evolution of two miRNAs. Phylogenetic analyses using both Bayesian and parsimony methods on a eumetazoan miRNA data set highlight the potential of miRNAs to become an invaluable new tool, especially when used as an additional line of evidence, to resolve previously intractable nodes within the tree of life.

MicroRNAs support the monophyly of enteropneust hemichordates

Understanding the evolutionary history of deuterostomes requires elucidating the phylogenetic interrelationships amongst the constituent taxa. Although the monophyly and interrelationships among the three principal groups-the chordates, the echinoderms, and the hemichordates-are well established, as are the internal relationships among the echinoderm and chordate taxa, the interrelationships among the principal groups of hemichordates-the harrimaniid enteropneusts, the ptychoderid enteropneusts, and the pterobranchs-remain unresolved. Depending on the study some find enteropneusts paraphyletic with pterobranchs (e.g., Cephalodiscus) more closely related to the harrimaniid enteropneusts (e.g., Saccoglossus) than either are to the ptychoderid enteropneusts (e.g., Ptychodera), whereas other studies support a monophyletic Enteropneusta. To try and resolve between these two competing hypotheses, we turned to microRNAs, small ∼22 nt non-coding RNA genes that have been shown to shed insight into particularly difficult phylogenetic questions. Using deep sequencing we characterized the small RNA repertoires of two hemichordate species, Cephalodiscus hodgsoni and Ptychodera flava, and the crinoid echinoderm Antedon mediterranea, and combined our results with the described complements of the hemichordate Saccoglossus kowalevskii, the sea urchin Strongylocentrotus purpuratus, and the starfish Patiria miniata. Our data unambiguously support the monophyly of Enteropneusts as S. kowalevskii shares 12 miRNA sequences with P. flava that are not present in the C. hodgsoni or A. mediterranea libraries, and have never been reported from another metazoan taxon. Thus, these data resolve the phylogenetic position of pterobranchs, ultimately allowing for a better understanding of body plan evolution throughout the deuterostomes.

Ranking transitive chemical-disease inferences using local network topology in the comparative toxicogenomics database

Exposure to chemicals in the environment is believed to play a critical role in the etiology of many human diseases. To enhance understanding about environmental effects on human health, the Comparative Toxicogenomics Database (CTD; http://ctdbase.org) provides unique curated data that enable development of novel hypotheses about the relationships between chemicals and diseases. CTD biocurators read the literature and curate direct relationships between chemicals-genes, genes-diseases, and chemicals-diseases. These direct relationships are then computationally integrated to create additional inferred relationships; for example, a direct chemical-gene statement can be combined with a direct gene-disease statement to generate a chemical-disease inference (inferred via the shared gene). In CTD, the number of inferences has increased exponentially as the number of direct chemical, gene and disease interactions has grown. To help users navigate and prioritize these inferences for hypothesis development, we implemented a statistic to score and rank them based on the topology of the local network consisting of the chemical, disease and each of the genes used to make an inference. In this network, chemicals, diseases and genes are nodes connected by edges representing the curated interactions. Like other biological networks, node connectivity is an important consideration when evaluating the CTD network, as the connectivity of nodes follows the power-law distribution. Topological methods reduce the influence of highly connected nodes that are present in biological networks. We evaluated published methods that used local network topology to determine the reliability of protein-protein interactions derived from high-throughput assays. We developed a new metric that combines and weights two of these methods and uniquely takes into account the number of common neighbors and the connectivity of each entity involved. We present several CTD inferences as case studies to demonstrate the value of this metric and the biological relevance of the inferences.

The Comparative Toxicogenomics Database: update 2013

The Comparative Toxicogenomics Database (CTD; http://ctdbase.org/) provides information about interactions between environmental chemicals and gene products and their relationships to diseases. Chemical-gene, chemical-disease and gene-disease interactions manually curated from the literature are integrated to generate expanded networks and predict many novel associations between different data types. CTD now contains over 15 million toxicogenomic relationships. To navigate this sea of data, we added several new features, including DiseaseComps (which finds comparable diseases that share toxicogenomic profiles), statistical scoring for inferred gene-disease and pathway-chemical relationships, filtering options for several tools to refine user analysis and our new Gene Set Enricher (which provides biological annotations that are enriched for gene sets). To improve data visualization, we added a Cytoscape Web view to our ChemComps feature, included color-coded interactions and created a 'slim list' for our MEDIC disease vocabulary (allowing diseases to be grouped for meta-analysis, visualization and better data management). CTD continues to promote interoperability with external databases by providing content and cross-links to their sites. Together, this wealth of expanded chemical-gene-disease data, combined with novel ways to analyze and view content, continues to help users generate testable hypotheses about the molecular mechanisms of environmental diseases.

The Alström syndrome protein, ALMS1, interacts with α-actinin and components of the endosome recycling pathway

Alström syndrome (ALMS) is a progressive multi-systemic disorder characterized by cone-rod dystrophy, sensorineural hearing loss, childhood obesity, insulin resistance and cardiac, renal, and hepatic dysfunction. The gene responsible for Alström syndrome, ALMS1, is ubiquitously expressed and has multiple splice variants. The protein encoded by this gene has been implicated in ciliary function, cell cycle control, and intracellular transport. To gain better insight into the pathways through which ALMS1 functions, we carried out a yeast two hybrid (Y2H) screen in several mouse tissue libraries to identify ALMS1 interacting partners. The majority of proteins found to interact with the murine carboxy-terminal end (19/32) of ALMS1 were α-actinin isoforms. Interestingly, several of the identified ALMS1 interacting partners (α-actinin 1, α-actinin 4, myosin Vb, rad50 interacting 1 and huntingtin associated protein1A) have been previously associated with endosome recycling and/or centrosome function. We examined dermal fibroblasts from human subjects bearing a disruption in ALMS1 for defects in the endocytic pathway. Fibroblasts from these patients had a lower uptake of transferrin and reduced clearance of transferrin compared to controls. Antibodies directed against ALMS1 N- and C-terminal epitopes label centrosomes and endosomal structures at the cleavage furrow of dividing MDCK cells, respectively, suggesting isoform-specific cellular functions. Our results suggest a role for ALMS1 variants in the recycling endosome pathway and give us new insights into the pathogenesis of a subset of clinical phenotypes associated with ALMS.

Clusters of ancestrally related genes that show paralogy in whole or in part are a major feature of the genomes of humans and other species

Arrangements of genes along chromosomes are a product of evolutionary processes, and we can expect that preferable arrangements will prevail over the span of evolutionary time, often being reflected in the non-random clustering of structurally and/or functionally related genes. Such non-random arrangements can arise by two distinct evolutionary processes: duplications of DNA sequences that give rise to clusters of genes sharing both sequence similarity and common sequence features and the migration together of genes related by function, but not by common descent [1], [2], [3]. To provide a background for distinguishing between the two, which is important for future efforts to unravel the evolutionary processes involved, we here provide a description of the extent to which ancestrally related genes are found in proximity.

Towards this purpose, we combined information from five genomic datasets, InterPro, SCOP, PANTHER, Ensembl protein families, and Ensembl gene paralogs. The results are provided in publicly available datasets (http://cgd.jax.org/datasets/clustering/paraclustering.shtml) describing the extent to which ancestrally related genes are in proximity beyond what is expected by chance (i.e. form paraclusters) in the human and nine other vertebrate genomes, as well as the D. melanogaster, C. elegans, A. thaliana, and S. cerevisiae genomes.

With the exception of Saccharomyces, paraclusters are a common feature of the genomes we examined. In the human genome they are estimated to include at least 22% of all protein coding genes. Paraclusters are far more prevalent among some gene families than others, are highly species or clade specific and can evolve rapidly, sometimes in response to environmental cues. Altogether, they account for a large portion of the functional clustering previously reported in several genomes.

Community annotation and bioinformatics workforce development in concert--Little Skate Genome Annotation Workshops and Jamborees

Recent advances in high-throughput DNA sequencing technologies have equipped biologists with a powerful new set of tools for advancing research goals. The resulting flood of sequence data has made it critically important to train the next generation of scientists to handle the inherent bioinformatic challenges. The North East Bioinformatics Collaborative (NEBC) is undertaking the genome sequencing and annotation of the little skate (Leucoraja erinacea) to promote advancement of bioinformatics infrastructure in our region, with an emphasis on practical education to create a critical mass of informatically savvy life scientists. In support of the Little Skate Genome Project, the NEBC members have developed several annotation workshops and jamborees to provide training in genome sequencing, annotation and analysis. Acting as a nexus for both curation activities and dissemination of project data, a project web portal, SkateBase (http://skatebase.org) has been developed. As a case study to illustrate effective coupling of community annotation with workforce development, we report the results of the Mitochondrial Genome Annotation Jamborees organized to annotate the first completely assembled element of the Little Skate Genome Project, as a culminating experience for participants from our three prior annotation workshops. We are applying the physical/virtual infrastructure and lessons learned from these activities to enhance and streamline the genome annotation workflow, as we look toward our continuing efforts for larger-scale functional and structural community annotation of the L. erinacea genome.

MicroRNAs support a turtle + lizard clade

Despite much interest in amniote systematics, the origin of turtles remains elusive. Traditional morphological phylogenetic analyses place turtles outside Diapsida-amniotes whose ancestor had two fenestrae in the temporal region of the skull (among the living forms the tuatara, lizards, birds and crocodilians)-and allied with some unfenestrate-skulled (anapsid) taxa. Nonetheless, some morphological analyses place turtles within Diapsida, allied with Lepidosauria (tuatara and lizards). Most molecular studies agree that turtles are diapsids, but rather than allying them with lepidosaurs, instead place turtles near or within Archosauria (crocodilians and birds). Thus, three basic phylogenetic positions for turtles with respect to extant Diapsida are currently debated: (i) sister to Diapsida, (ii) sister to Lepidosauria, or (iii) sister to, or within, Archosauria. Interestingly, although these three alternatives are consistent with a single unrooted four-taxon tree for extant reptiles, they differ with respect to the position of the root. Here, we apply a novel molecular dataset, the presence versus absence of specific microRNAs, to the problem of the phylogenetic position of turtles and the root of the reptilian tree, and find that this dataset unambiguously supports a turtle + lepidosaur group. We find that turtles and lizards share four unique miRNA gene families that are not found in any other organisms' genome or small RNA library, and no miRNAs are found in all diapsids but not turtles, or in turtles and archosaurs but not in lizards. The concordance between our result and some morphological analyses suggests that there have been numerous morphological convergences and reversals in reptile phylogeny, including the loss of temporal fenestrae.

A natural deletion of the HoxC cluster in elasmobranch fishes

Hox proteins are a metazoan-specific family of transcription factors that are required for developmental patterning. The genomic arrangement of Hox genes into four paralogous clusters is a primitive feature of jawed vertebrates. By using high-throughput sequencing, we demonstrate the absence of all HoxC transcripts from embryos of the shark Scyliorhinus canicula and the skate Leucoraja erinacea and the absence of all HoxC genes and two HoxC-associated microRNAs from the genome of L. erinacea. These data suggest a loss of the entire HoxC cluster in elasmobranch fishes and represent evidence for the natural deletion of an entire Hox cluster in vertebrates.

Datgan, a reusable software system for facile interrogation and visualization of complex transcription profiling data

BACKGROUND

We introduce Glaucoma Discovery Platform (GDP), an online environment for facile visualization and interrogation of complex transcription profiling datasets for glaucoma. We also report the availability of Datgan, the suite of scripts that was developed to construct GDP. This reusable software system complements existing repositories such as NCBI GEO or EBI ArrayExpress as it allows the construction of searchable databases to maximize understanding of user-selected transcription profiling datasets.

DESCRIPTION

Datgan scripts were used to construct both the underlying data tables and the web interface that form GDP. GDP is populated using data from a mouse model of glaucoma. The data was generated using the DBA/2J strain, a widely used mouse model of glaucoma. The DBA/2J-Gpnmb+ strain provided a genetically matched control strain that does not develop glaucoma. We separately assessed both the retina and the optic nerve head, important tissues in glaucoma. We used hierarchical clustering to identify early molecular stages of glaucoma that could not be identified using morphological assessment of disease. GDP has two components. First, an interactive search and retrieve component provides the ability to assess gene(s) of interest in all identified stages of disease in both the retina and optic nerve head. The output is returned in graphical and tabular format with statistically significant differences highlighted for easy visual analysis. Second, a bulk download component allows lists of differentially expressed genes to be retrieved as a series of files compatible with Excel. To facilitate access to additional information available for genes of interest, GDP is linked to selected external resources including Mouse Genome Informatics and Online Medelian Inheritance in Man (OMIM).

CONCLUSION

Datgan-constructed databases allow user-friendly access to datasets that involve temporally ordered stages of disease or developmental stages. Datgan and GDP are available from http://glaucomadb.jax.org/glaucoma.

Molecular clustering identifies complement and endothelin induction as early events in a mouse model of glaucoma

Glaucoma is one of the most common neurodegenerative diseases. Despite this, the earliest stages of this complex disease are still unclear. This study was specifically designed to identify early stages of glaucoma in DBA/2J mice. To do this, we used genome-wide expression profiling of optic nerve head and retina and a series of computational methods. Eyes with no detectable glaucoma by conventional assays were grouped into molecularly defined stages of disease using unbiased hierarchical clustering. These stages represent a temporally ordered sequence of glaucoma states. We then determined networks and biological processes that were altered at these early stages. Early-stage expression changes included upregulation of both the complement cascade and the endothelin system, and so we tested the therapeutic value of separately inhibiting them. Mice with a mutation in complement component 1a (C1qa) were protected from glaucoma. Similarly, inhibition of the endothelin system with bosentan, an endothelin receptor antagonist, was strongly protective against glaucomatous damage. Since endothelin 2 is potently vasoconstrictive and was produced by microglia/macrophages, our data provide what we believe to be a novel link between these cell types and vascular dysfunction in glaucoma. Targeting early molecular events, such as complement and endothelin induction, may provide effective new treatments for human glaucoma.

Prefoldin 5 is required for normal sensory and neuronal development in a murine model

Molecular chaperones and co-chaperones are crucial for cellular development and maintenance as they assist in protein folding and stabilization of unfolded or misfolded proteins. Prefoldin (PFDN), a ubiquitously expressed heterohexameric co-chaperone, is necessary for proper folding of nascent proteins, in particular, tubulin and actin. Here we show that a genetic disruption in the murine Pfdn5 gene, a subunit of prefoldin, causes a syndrome characterized by photoreceptor degeneration, central nervous system abnormalities, and male infertility. Our data indicate that a missense mutation in Pfdn5, may cause these phenotypes through a reduction in formation of microtubules and microfilaments, which are necessary for the development of cilia and cytoskeletal structures, respectively. The diversity of phenotypes demonstrated by models carrying mutations in different PFDN subunits suggests that each PFDN subunit must confer a distinct substrate specificity to the prefoldin holocomplex.