Search

The type 2 deiodinase (D2) in the neonatal liver accelerates local thyroid hormone triiodothyronine (T3) production and expression of T3-responsive genes. Here we show that this surge in T3 permanently modifies hepatic gene expression. Liver-specific Dio2 inactivation (Alb-D2KO) transiently increases H3K9me3 levels during post-natal days 1–5 (P1–P5), and results in methylation of 1,508 DNA sites (H-sites) in the adult mouse liver. These sites are associated with 1,551 areas of reduced chromatin accessibility (RCA) within core promoters and 2,426 within intergenic regions, with reduction in the expression of 1,363 genes. There is strong spatial correlation between density of H-sites and RCA sites. Chromosome conformation capture (Hi-C) data reveals a set of 81 repressed genes with a promoter RCA in contact with an intergenic RCA ~300 Kbp apart, within the same topologically associating domain (χ2 = 777; p, Whether thyroid hormones affect gene expression via DNA methylation is not well known. Here the authors show that type 2 deiodinase (D2) converts T4 to produce T3, which prevents DNA methylation of discrete areas in the neonatal liver. In the absence of D2, DNA methylation occurs and is associated with reduced chromatin accessibility in promoters and enhancers and affects gene expression.

The National Cancer Institute (NCI) Genomic Data Commons (GDC) contains more than 2.9 petabytes of genomic and associated clinical data from more than 60 NCI-funded and other contributed cancer genomics research projects. The GDC consists of five applications over a common data model and a common application programming interface.

Improvement of risk stratification through prognostic biomarkers may enhance the personalization of cancer patient monitoring and treatment. We used Ancer, an immunoinformatic CD8, CD4, and regulatory T cell neoepitope screening system, to perform an advanced neoantigen analysis of genomic data derived from the urothelial cancer cohort of The Cancer Genome Atlas. Ancer demonstrated improved prognostic stratification and five-year survival prediction compared to standard analyses using tumor mutational burden or neoepitope identification using NetMHCpan and NetMHCIIpan. The superiority of Ancer, shown in both univariate and multivariate survival analyses, is attributed to the removal of neoepitopes that do not contribute to tumor immunogenicity based on their homology with self-epitopes. This analysis suggests that the presence of a higher number of unique, non-self CD8- and CD4-neoepitopes contributes to cancer survival, and that prospectively defining these neoepitopes using Ancer is a novel prognostic or predictive biomarker.

Background Improvement of current prognosis biomarkers will enhance our ability to identify cancer patients at higher risk of recurrence and will further advance the personalization of patient monitoring and treatment. We hypothesized that the presence of a mutation alone is not sufficient to generate an immunogenic neoepitope, but that significant differences must exist between the Human Leukocyte Antigen (HLA)- and/or T Cell Receptor (TCR)-interfaces of the neoepitope and its non-mutated form, or with other self-epitopes, in order to be recognized as non-self by the immune system. As such, cancer patient clinical outcomes may be better understood by neoepitope analyses that integrate these considerations. Methods We analyzed large scale (n=412) bladder cancer genomic data from The Cancer Genome Atlas (TCGA) using Ancer, an automated machine-learning-based pipeline we designed for neoantigen screening and vaccine design. Ancer shares components with other commercial-grade screening platforms used routinely in immunogenicity assessments of biologics and infectious disease antigens, such as the EpiMatrix algorithm for HLA-I and HLA-II neoepitope identification, and the JanusMatrix algorithm for tolerated, tolerogenic, and cross-reactive T cell epitope identification. Evaluation of patient survival with Ancer was compared to other analyses employing tumor mutational burden (TMB) or neoepitopes identified with the commonly used NetMHCpan-4.0 and NetMHCIIpan-3.1 T cell epitope prediction tools. Results We stratified bladder patients based on their Ancer HLA-I and HLA-II neoepitope burdens and observed significantly prolonged disease free and overall survival in patients whose tumor contained both high numbers of HLA-I and HLA-II neoepitopes compared to other individuals. Stratifications performed with Ancer were superior to comparative analyses performed with TMB or with NetMHCpan and NetMHCIIpan. In addition, we showed that Ancer’s precise filtering and characterization of mutated epitopes contributed to this increased association with survival, as showcased by gradual improvements in survival analyses performed after each of its filtering step. Multivariate survival analyses indicated that Ancer neoepitope content remained a significant factor in patient overall survival even when adjusted for TMB, and other clinical covariates such as age at diagnosis and disease stage, unlike analyses involving NetMHCpan and NetMHCIIpan neoepitopes. Conclusions Our analysis suggests that enhanced presence of CD8, CD4 T cell epitopes, and limited inclusion of Treg epitopes in the tumor genome plays a key role in cancer survival. Ancer scoring provides a predictive method for predicting patient outcomes, by defining the number of true neoepitopes and by identifying Treg epitopes that would interfere with T cell-based immune activation and response to the tumor.

In estrogen receptor (ER)–negative breast cancer, high tumor glucocorticoid receptor (GR) expression has been associated with a relatively poor outcome. In contrast, using a meta-analysis of several genomic datasets, here we find that tumor GR mRNA expression is associated with improved ER+ relapse-free survival (RFS; independently of progesterone receptor expression). To understand the mechanism by which GR expression is associated with a better ER+ breast cancer outcome, the global effect of GR-mediated transcriptional activation in ER+ breast cancer cells was studied. Analysis of GR chromatin immunoprecipitation followed by high-throughput sequencing in ER+/GR+ MCF-7 cells revealed that upon coactivation of GR and ER, GR chromatin association became enriched at proximal promoter regions. Furthermore, following ER activation, increased GR chromatin association was observed at ER, FOXO, and AP1 response elements. In addition, ER associated with GR response elements, suggesting that ER and GR interact in a complex. Coactivation of GR and ER resulted in increased expression (relative to ER activation alone) of transcripts that encode proteins promoting cellular differentiation (e.g., KDM4B, VDR) and inhibiting the Wnt signaling pathway (IGFBP4). Finally, expression of these individual prodifferentiation genes was associated with significantly improved RFS in ER+ breast cancer patients. Together, these data suggest that the coexpression and subsequent activity of tumor cell GR and ER contribute to the less aggressive natural history of early-stage breast cancer by coordinating the altered expression of genes favoring differentiation. Implications: The interaction between ER and GR activity highlights the importance of context-dependent nuclear receptor function in cancer. Mol Cancer Res; 14(8); 707–19. ©2016 AACR.

Hexavalent chromium Cr(VI) is known to produce cytotoxic effects in humans and is a highly toxic environmental contaminant. Interestingly, it has been shown that free ranging sperm whales (Phyester macrocephalus) may have exceedingly high levels of Cr in their skin. Also, it has been demonstrated that skin cells from whales appear more resistant to both cytotoxicity and clastogenicity upon Cr exposure compared to human cells. However, the molecular genetic mechanisms employed in whale skin cells that might lead to Cr tolerance are unknown. In an effort to understand the underlying mechanisms of Cr(VI) tolerance and to illuminate global gene expression patterns modulated by Cr, we exposed whale skin cells in culture to varying levels of Cr(VI) (i.e., 0.0, 0.5, 1.0 and 5.0 μg/cm²) followed by short read (100 bp) next generation RNA sequencing (RNA-seq). RNA-seq reads from all exposures (≈280 million reads) were pooled to generate a de novo reference transcriptome assembly. The resulting whale reference assembly had 11K contigs and an N50 of 2954 bp. Using the reads from each dose (0.0, 0.5, 1.0 and 5.0 μg/cm²) we performed RNA-seq based gene expression analysis that identified 35 up-regulated genes and 19 down-regulated genes. The experimental results suggest that low dose exposure to Cr (1.0 μg/cm²) serves to induce up-regulation of oxidative stress response genes, DNA repair genes and cell cycle regulator genes. However, at higher doses (5.0 μg/cm²) the DNA repair genes appeared down-regulated while other genes that were induced suggest the initiation of cytotoxicity. The set of genes identified that show regulatory modulation at different Cr doses provide specific candidates for further studies aimed at determination of how whales exhibit resistance to Cr toxicity and what role(s) reactive oxygen species (ROS) may play in this process.

Obstructive sleep apnea (OSA) affects 8–10% of the population, is characterized by chronic intermittent hypoxia (CIH), and causally associates with cardiovascular morbidities. In CIH-exposed mice, closely mimicking the chronicity of human OSA, increased accumulation and proliferation of pro-inflammatory metabolic M1-like macrophages highly expressing CD36, emerged in aorta. Transcriptomic and MeDIP-seq approaches identified activation of pro-atherogenic pathways involving a complex interplay of histone modifications in functionally-relevant biological pathways, such as inflammation and oxidative stress in aorta macrophages. Discontinuation of CIH did not elicit significant improvements in aorta wall macrophage phenotype. However, CIH-induced aorta changes were absent in CD36 knockout mice, Our results provide mechanistic insights showing that CIH exposures during sleep in absence of concurrent pro-atherogenic settings (i.e., genetic propensity or dietary manipulation) lead to the recruitment of CD36(+)high macrophages to the aortic wall and trigger atherogenesis. Furthermore, long-term CIH-induced changes may not be reversible with usual OSA treatment.

Relationship of GO categories that are enriched in genes with high impact variants. (PPTX 93 kb)

Dot plots of location of one-to-one orthologues in the chromosomes of X. hellerii and X. maculatus. (PDF 111 kb)

Background Xiphophorus fishes are represented by 26 live-bearing species of tropical fish that express many attributes (e.g., viviparity, genetic and phenotypic variation, ecological adaptation, varied sexual developmental mechanisms, ability to produce fertile interspecies hybrids) that have made attractive research models for over 85 years. Use of various interspecies hybrids to investigate the genetics underlying spontaneous and induced tumorigenesis has resulted in the development and maintenance of pedigreed Xiphophorus lines specifically bred for research. The recent availability of the X. maculatus reference genome assembly now provides unprecedented opportunities for novel and exciting comparative research studies among Xiphophorus species. Results We present sequencing, assembly and annotation of two new genomes representing Xiphophorus couchianus and Xiphophorus hellerii. The final X. couchianus and X. hellerii assemblies have total sizes of 708 Mb and 734 Mb and correspond to 98 % and 102 % of the X. maculatus Jp 163 A genome size, respectively. The rates of single nucleotide change range from 1 per 52 bp to 1 per 69 bp among the three genomes and the impact of putatively damaging variants are presented. In addition, a survey of transposable elements allowed us to deduce an ancestral TE landscape, uncovered potential active TEs and document a recent burst of TEs during evolution of this genus. Conclusions Two new Xiphophorus genomes and their corresponding transcriptomes were efficiently assembled, the former using a novel guided assembly approach. Three assembled genome sequences within this single vertebrate order of new world live-bearing fishes will accelerate our understanding of relationship between environmental adaptation and genome evolution. In addition, these genome resources provide capability to determine allele specific gene regulation among interspecies hybrids produced by crossing any of the three species that are known to produce progeny predisposed to tumor development. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2361-z) contains supplementary material, which is available to authorized users.

Dot plots of location of one-to-one orthologues in the 24 chromosomes of X. couchianus and X. maculatus. (PDF 110 kb)

Variations in gene expression are essential for the evolution of novel phenotypes and for speciation. Studying allelic specific gene expression (ASGE) within interspecies hybrids provides a unique opportunity to reveal underlying mechanisms of genetic variation. Using Xiphophorus interspecies hybrid fishes and high-throughput next generation sequencing technology, we were able to assess variations between two closely related vertebrate species, Xiphophorus maculatus and Xiphophorus couchianus, and their F(1) interspecies hybrids. We constructed transcriptome-wide SNP polymorphism sets between two highly inbred X. maculatus lines (JP 163 A and B), and between X. maculatus and a second species, X. couchianus. The X. maculatus JP 163 A and B parental lines have been separated in the laboratory for ≈70 years and we were able to identify SNPs at a resolution of 1 SNP per 49 kb of transcriptome. In contrast, SNP polymorphisms between X. couchianus and X. maculatus species, which diverged ≈5-10 million years ago, were identified about every 700 bp. Using 6524 transcripts with identified SNPs between the two parental species (X. maculatus and X. couchianus), we mapped RNA-seq reads to determine ASGE within F(1) interspecies hybrids. We developed an in silico X. couchianus transcriptome by replacing 90,788 SNP bases for X. maculatus transcriptome with the consensus X. couchianus SNP bases and provide evidence that this procedure overcomes read mapping biases. Employment of the in silico reference transcriptome and tolerating 5 mismatches during read mapping allow direct assessment of ASGE in the F(1) interspecies hybrids. Overall, these results show that Xiphophorus is a tractable vertebrate experimental model to investigate how genetic variations that occur during speciation may affect gene interactions and the regulation of gene expression.

The biological consequences of the Deepwater Horizon oil spill are unknown, especially for resident organisms. Here, we report results from a field study tracking the effects of contaminating oil across space and time in resident killifish during the first 4 mo of the spill event. Remote sensing and analytical chemistry identified exposures, which were linked to effects in fish characterized by genome expression and associated gill immunohistochemistry, despite very low concentrations of hydrocarbons remaining in water and tissues. Divergence in genome expression coincides with contaminating oil and is consistent with genome responses that are predictive of exposure to hydrocarbon-like chemicals and indicative of physiological and reproductive impairment. Oil-contaminated waters are also associated with aberrant protein expression in gill tissues of larval and adult fish. These data suggest that heavily weathered crude oil from the spill imparts significant biological impacts in sensitive Louisiana marshes, some of which remain for over 2 mo following initial exposures.

Background: Although organ-confined, moderately graded prostate cancer has encouraging cure rates, progression to castration-resistant prostate cancer (CRPC) from advanced, high-grade lesions represents the lethal form of this malignancy, for which current therapeutic modalities are only palliative; patients with CRPC will die from their disease. Importantly, African-American (AA) men exhibit both a higher incidence of prostate cancer diagnoses and are at higher risk of developing CRPC as compared to other racial groups. Currently, the mechanisms driving prostate cancer disparities in AA men are poorly understood, and new biomarkers and novel drugs to address this disparity are in high demand to improve personalized medicine and clinical management of prostate cancer in AA men. Work by our group and others has identified SOX2 (sex determining region Y-box 2), an essential transcription factor for maintaining the survival and pluripotency of undifferentiated embryonic stem cells (ESCs), to be associated with aggressive prostate cancer and that the constitutive overexpression of SOX2 in a hormone-sensitive, SOX2-negative prostate cancer cell line is sufficient to generate a castration-resistant phenotype both in vitro and in vivo. Further, our data demonstrate canonical SOX2 transcriptional cofactors OCT4 and NANOG are frequently not expressed in prostate cancer, implying that SOX2 has unique, non-stem cell gene targets and binding partners within malignant prostate malignancies. Methods and Results: To elucidate such SOX2 gene targets in CRPC, we performed SOX2 chromatin immunoprecipitation and sequencing (ChIP-Seq) in prostate cancer cells, revealing SOX2 binding of several interesting and functionally important gene target promoters, including FOXP1, an AR target gene and transcription factor known to play a role in disparate aggressiveness of prostate cancer in African-American men. Based on these findings, we hypothesize that upon AR ablation in prostate cancer treatment, FOXP1, a tumor suppressor demonstrated to be associated with disparate prostate cancer aggressiveness in AA men, is repressed by SOX2 to drive CRPC in AA men. Work to test this hypothesis is currently under way, and includes (1) assessment of differential FOXP1 gene expression using MDA-PCa-2a, MDA-PCa-2b, and EA0066-hT, prostate cancer cell lines derived from AA men, and CWRR1, CWR22Rv1, and LNCaP, prostate cancer cell lines derived from CA men overexpressing SOX2 in the context of AR ablation; and (2) FOXP1 ChIP-Seq paired with RNA-Seq in tumor tissues prospectively collected from AA and CA men with prostate cancer undergoing radical prostatectomy. Conclusion: The work proposed herein represents an in-depth exploration of basic SOX2 biology in the context of CRPC, is highly innovative and translational, and has transformative potential to improve clinical patient management and eradicate disparities in CRPC. Note: This abstract was not presented at the conference. Citation Format: Anthony Williams, Larischa de Wet, Marc Gillard, Steve Kregel, Tzintzuni Garcia, Russell Szmulewitz, Don Vander Griend. SOX2 regulation of FOXP1 as a mechanism driving disparate aggressiveness of prostate cancer in African-American men [abstract]. In: Proceedings of the Tenth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2017 Sep 25-28; Atlanta, GA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2018;27(7 Suppl):Abstract nr B73.

Prostate cancer is the most common type of cancer in men, with approximately 181,000 new cases diagnosed in 2016. Due to the central role of the androgen receptor (AR) in prostate development and more importantly prostate cancer cell survival and proliferation, strategies targeting AR have been the mainstay therapy for over 70 years. However, despite potent inhibition of AR pathway activation, many patients develop castration-resistant prostate cancer (CRPC). Second-line therapies, such as enzalutamide, has increased overall survival in CRPC, but resistance to these therapies inevitably emerges, suggesting that other pathways apart from AR signalling are contributing to the failure of treatments. We have previously demonstrated that SOX2 [SRY (sex determining region Y)-box 2] is an AR-repressed gene that is expressed in a large percentage of high Gleason grade prostate tumours, as well as in most metastases. Additionally, expression of SOX2 within a castration-sensitive cell line is sufficient to enable castration-resistant tumor formation in vivo, and enzalutamide resistance in vitro. In prostate cancer, SOX2 is not found with its normal stem cell co-factors, NANOG and OCT4, leading to the hypothesis that SOX2 is interacting with a novel co-factor in prostate cancer to regulate expression of genes promoting castration-resistance. A chromatin-immunoprecipitation experiment in a castration-resistant prostate cancer cell line was performed to determine if SOX2 binds the same genes and activates similar pathways in prostate cancer as in embryonic stem cells. Approximately half of the SOX2 bound genes in the CRPC cell line were shared with the embryonic stem cell line, and these genes are involved in signalling pathways and maintenance of stem cell pluripotency; the genes bound uniquely in the CRPC cell line are present in pathways involved with metabolic processes. To identify potential binding partners of SOX2, computational analysis of the sites bound by SOX2 in the promoter regions of target genes determined the FOXA1 motif is within close proximity to the SOX2 motif, and the physical interaction of these proteins was confirmed through co-immunoprecipitation. Further understanding of SOX2 target genes and the pathways that SOX2 activates in the presence of enzalutamide will be crucial to identify mechanisms of resistance and enable the development of novel therapies for castration-resistant prostate cancer. Citation Format: Larischa de Wet, Anthony Williams, Marc Gillard, Steve Kregel, Tzintzuni Garcia, Erin McAuley, Ryan Brown, Donald Vander Griend. The role of SOX2 in promoting resistance to AR-targeted therapies in prostate cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3348.

Chikungunya virus (CHIKV), a mosquito-borne alphavirus, has traditionally circulated in Africa and Asia, causing human febrile illness accompanied by severe, chronic joint pain. In Africa, epidemic emergence of CHIKV involves the transition from an enzootic, sylvatic cycle involving arboreal mosquito vectors and nonhuman primates, into an urban cycle where peridomestic mosquitoes transmit among humans. In Asia, however, CHIKV appears to circulate only in the endemic, urban cycle. Recently, CHIKV emerged into the Indian Ocean and the Indian subcontinent to cause major epidemics. To examine patterns of CHIKV evolution and the origins of these outbreaks, as well as to examine whether evolutionary rates that vary between enzootic and epidemic transmission, we sequenced the genomes of 40 CHIKV strains and performed a phylogenetic analysis representing the most comprehensive study of its kind to date. We inferred that extant CHIKV strains evolved from an ancestor that existed within the last 500 years and that some geographic overlap exists between two main enzootic lineages previously thought to be geographically separated within Africa. We estimated that CHIKV was introduced from Africa into Asia 70 to 90 years ago. The recent Indian Ocean and Indian subcontinent epidemics appear to have emerged independently from the mainland of East Africa. This finding underscores the importance of surveillance to rapidly detect and control African outbreaks before exportation can occur. Significantly higher rates of nucleotide substitution appear to occur during urban than during enzootic transmission. These results suggest fundamental differences in transmission modes and/or dynamics in these two transmission cycles.

In many countries regulatory agencies have adopted safety guidelines, based on bioinformatics rules from the WHO/FAO and EFSA recommendations, to prevent potentially allergenic novel foods or agricultural products from reaching consumers. We created the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP/) to combine data that had previously been available only as flat files on Web pages or in the literature. SDAP was designed to be user friendly, to be of maximum use to regulatory agencies, clinicians, as well as to scientists interested in assessing the potential allergenic risk of a protein. We developed methods, unique to SDAP, to compare the physicochemical properties of discrete areas of allergenic proteins to known IgE epitopes. We developed a new similarity measure, the property distance (PD) value that can be used to detect related segments in allergens with clinical observed cross-reactivity. We have now expanded this work to obtain experimental validation of the PD index as a quantitative predictor of IgE cross-reactivity, by designing peptide variants with predetermined PD scores relative to known IgE epitopes. In complementary work we show how sequence motifs characteristic of allergenic proteins in protein families can be used as fingerprints for allergenicity.

The mechanisms that determine mechanical stabilities of protein folds remain elusive. Our understanding of these mechanisms is vital to both bioengineering efforts and to the better understanding and eventual treatment of pathogenic mutations affecting mechanically important proteins such as titin. We present a new approach to analyze data from single-molecule force spectroscopy for different domains of the giant muscle protein titin. The region of titin found in the I-band of a sarcomere is composed of about 40 Ig-domains and is exposed to force under normal physiological conditions and connects the free-hanging ends of the myosin filaments to the Z-disc. Recent single-molecule force spectroscopy data show a mechanical hierarchy in the I-band domains. Domains near the C-terminus in this region unfold at forces two to three times greater than domains near the beginning of the I-band. Though all of these Ig-domains are thought to share a fold and topology common to members of the Ig-like fold family, the sequences of neighboring domains vary greatly with an average sequence identity of only 25%. We examine in this study the relation of these unique mechanical stabilities of each I-band Ig domain to specific, conserved physical-chemical properties of amino acid sequences in related Ig domains. We find that the sequences of each individual titin Ig domain are very highly conserved, with an average sequence identity of 79% across species that are divergent as humans, chickens, and zebra fish. This indicates that the mechanical properties of each domain are well conserved and tailored to its unique position in the titin molecule. We used the PCPMer software to determine the conservation of amino acid properties in titin Ig domains grouped by unfolding forces into "strong" and "weak" families. We found two motifs unique to each family that may have some role in determining the mechanical properties of these Ig domains. A detailed statistical analysis of properties of individual residues revealed several positions that displayed differentially conserved properties in strong and weak families. In contrast to previous studies, we find evidence that suggests that the mechanical stability of Ig domains is determined by several residues scattered across the beta-sandwich fold, and force sensitive residues are not only confined to the A'-G region.

Myofibril assembly and disassembly are complex processes that regulate overall muscle mass. Titin kinase has been implicated as an initiating catalyst in signaling pathways that ultimately result in myofibril growth. In titin, the kinase domain is in an ideal position to sense mechanical strain that occurs during muscle activity. The enzyme is negatively regulated by intramolecular interactions occurring between the kinase catalytic core and autoinhibitory/regulatory region. Molecular dynamics simulations suggest that human titin kinase acts as a force sensor. However, the precise mechanism(s) resulting in the conformational changes that relieve the kinase of this autoinhibition are unknown. Here we measured the mechanical properties of the kinase domain and flanking Ig/Fn domains of the Caenorhabditis elegans titin-like proteins twitchin and TTN-1 using single-molecule atomic force microscopy. Our results show that these kinase domains have significant mechanical resistance, unfolding at forces similar to those for Ig/Fn β-sandwich domains (30–150pN). Further, our atomic force microscopy data is consistent with molecular dynamic simulations, which show that these kinases unfold in a stepwise fashion, first an unwinding of the autoinhibitory region, followed by a two-step unfolding of the catalytic core. These data support the hypothesis that titin kinase may function as an effective force sensor.

Xiphophorus fish and interspecies hybrids represent long-standing models to study the genetics underlying spontaneous and induced tumorigenesis. The recent release of the Xiphophorus maculatus genome sequence will allow global genetic regulation studies of genes involved in the inherited susceptibility to UVB-induced melanoma within select backcross hybrids. As a first step toward this goal, we report results of an RNA-Seq approach to identify genes and pathways showing modulated transcription within the skin of X. maculatus Jp 163 B upon UVB exposure. X. maculatus Jp 163 B were exposed to various doses of UVB followed by RNA-Seq analysis at each dose to investigate overall gene expression in each sample. A total of 357 genes with a minimum expression change of 4-fold (p-adj < 0.05) were identified as responsive to UVB. The molecular genetic response of Xiphophorus skin to UVB exposure permitted assessment of; (1) the basal expression level of each transcript for each skin sample, (2) the changes in expression levels for each gene in the transcriptome upon exposure to increasing doses of UVB, and (3) clusters of genes that exhibit similar patterns of change in expression upon UVB exposure. These data provide a foundation for understanding the molecular genetic response of fish skin to UVB exposure.

The authors would like to thank the staff of the Xiphophorus Genetic Stock Center (XGSC) and the Biocenter Fish Facility for maintaining the pedigreed fish lines used in this study. We gratefully acknowledge the sequencing efforts of C. Fronick, K. Delehaunty and the production sequencing group at the Genome Institute. This work was supported by US National Institutes of Health, National Center for Research Resources (NCRR) and Office of Research Infrastructure Programs (ORIP), Division of Comparative Medicine grants R24 RR024790 and R24 OD011120 (R. B. W.), including an American Recovery and Reinvestment Act supplement to this award, and R24OD011199 (R. B. W.), R24 RR032658 and R24 OD011198 (W. C. W.), R01 RR020833 and R01 OD011116 (J.H.P.), by the Deutsche Forschungsgemeinschaft, TRR 58/A5 (K. P. L.) and TRR 17 (M. S.) VolkswagenStiftung, grant I/84 815 (I. B.) and by the Agence Nationale de Recherche (J.-N.V.).; International audience; Several attributes intuitively considered to be typical mammalian features, such as complex behavior, live birth and malignant disease such as cancer, also appeared several times independently in lower vertebrates. The genetic mechanisms underlying the evolution of these elaborate traits are poorly understood. The platyfish, X. maculatus, offers a unique model to better understand the molecular biology of such traits. We report here the sequencing of the platyfish genome. Integrating genome assembly with extensive genetic maps identified an unexpected evolutionary stability of chromosomes in fish, in contrast to in mammals. Genes associated with viviparity show signatures of positive selection, identifying new putative functional domains and rare cases of parallel evolution. We also find that genes implicated in cognition show an unexpectedly high rate of duplicate gene retention after the teleost genome duplication event, suggesting a hypothesis for the evolution of the behavioral complexity in fish, which exceeds that found in amphibians and reptiles.

In recent years RNA-Seq technology has been used not only to quantify differences in gene expression but also to understand the underlying mechanisms that lead to these differences. Nucleotide sequence variation arising through evolution may differentially affect the expression profiles of divergent species. RNA-Seq technology, combined with techniques to differentiate parental alleles and quantify their abundance, have recently become popular methods for allele specific gene expression (ASGE) analyses. However, analysis of gene expression within interspecies hybrids may be difficult when one of the two parental genomes represented in the hybrid does not have robust genomic resources or available transcriptome data. Herein, we compare two strategies for analyzing allele specific expression within interspecies hybrids produced from crossing two Xiphophorus fish species. The first strategy relies upon a robust reference transcriptome assembly from one species followed by identification of SNPs and creation of an in silico reference transcriptome for the second species. The second strategy employs de novo assembly of reference transcriptomes for both parental species followed by identification of homologous transcripts prior to mapping hybrid reads to a combined hybrid reference. Our results show that, although both methods are able to achieve balanced allelic distribution upon read mapping of F(1) hybrid fish transcriptomes, the second "de novo" assembly approach is superior for ASGE analyses and leads to results more consistent with those found from quantitative real time PCR assessment of gene expression. In addition, our analysis indicates that indels between the two parental alleles are the major cause of the differences in results observed when employing these two methods.

Our original article (1) linked exposure of resident killifish to contaminating oil from the Deepwater Horizon (DWH) oil spill with significant biological responses, including genome expression, protein expression, and tissue morphology. Given decades of laboratory studies on the effects of crude oil in many species, including fish, and after extensive field studies following the Exxon Valdez spill, some of the responses we captured are recognized as diagnostic of exposure to, and effects from, the toxic components of weathered crude oil (e.g., ref. 2). Jenkins et … [↵][1]1To whom correspondence should be addressed. E-mail: awhitehead{at}ucdavis.edu. [1]: #xref-corresp-1-1

Development of massively parallel “next generation” sequencing technology (NGS) has dramatically revolutionized biological studies. Among the many applications of NGS, RNASeq is one of the most important uses of this technology. RNA-Seq enables investigators to accurately probe the current state of a transcriptome and assess many biologically important issues, such as; gene expression levels, differential splicing events, and allele-specific gene expression. Compared with previous technologies (e.g., microarrays, etc.) NGS has the clear advantage of not being limited to experimental systems having well characterized genomes or transcript sequence libraries. This positions RNA-seq approaches as important and versatile techniques for experimental systems and species where specific genetic information may be limited or altogether lacking. A major goal of most transcriptomic studies is the identification and characterization of all transcripts within a developmental stage or specific tissue. NGS techniques have made the massive amount of data required to carry out such studies both inexpensive and available to an unprecedented extent. Clever computer algorithms have made the assembly of these massive data sets the work of one or two people with reasonably powerful workstations or a moderate analytical server. Once a reference transcriptome has been assembled, analyses can be carried out that involve several steps, such as; mapping short sequence reads to transcriptome, quantifying the abundance of genes or gene sets, and comparing differential expression patterns among all samples. Herein we outline the processes from obtaining raw short read data to advanced comparative gene expression analysis and we review bioinformatic programs currently available, such as Tophat, Cufflinks, DESeq, that are specifically designed to address each of the above steps. We will discuss both accuracy and ease of use of these tools by biologists beginning to pursue these types of analyses. In addition to individual programs, we will also discuss integration of multiple programs into pipelines for more rapid and complete expression analyses. Overall, the future applications of RNA-Seq will open new avenues for transcriptome analyses of less well-studied and/or wild caught species that could not have previously been approached. This will yield a wealth of new comparative data highlighting the many ways plants and animals have developed to survive in this rapidly changing environment.

Research investigating telomere lengths and telomerase expression in vertebrates has progressively become important due to the association of these two biological endpoints with cellular aging and cancer in humans. Studies that rely upon the traditional use of laboratory mice have been faced with limitations largely due to inbred mice possessing large telomeres and ubiquitous expression of telomerase. Recently, a number of small fish species have been shown to provide potentially informative models for examining the role of telomeres and telomerase within intact vertebrate animals. Xiphophorus fishes represent a new world live-bearing genus that has not previously been assessed for telomere length or telomerase expression. To add to the knowledge base of telomere and telomerase biology in vertebrates we assessed telomere length and telomerase expression among several species of Xiphophorus. The telomere lengths in several organs (gill, brain, eyes, testis, ovary and liver) in three species (Xiphophorus hellerii, Xiphophorus maculatus, Xiphophorus couchianus) and also in F(1) interspecies hybrids were approximately 2-6 kb. This size was consistent within the same organs of the same species, as well as between species and F(1) hybrids. Despite possessing relatively short telomere lengths compared to humans, the consistency of size among Xiphophorus species and organs may allow experimental detection of telomere shortening. The relative expression of telomerase reverse transcriptase (TERT) was determined by quantitative real-time PCR. Expression levels of TERT was measured in seven organs (ovary, testis, liver, gill, brain, heart, skin) from X. maculatus, X. hellerii and in control and ultraviolet light (UVB) exposed skin samples from X. maculatus, X. hellerii, and F(1) interspecies hybrids. TERT gene expression was significantly higher in ovary and testis, while all other organs showed low relative TERT expression. Detectable increases in TERT expression were found in skin samples upon UVB exposure. Our findings suggest that Xiphophorus may serve as a suitable model for future studies investigating the association of telomere length and telomerase expression in regard to aging and disease.

For many researchers, next generation sequencing data holds the key to answering a category of questions previously unassailable. One of the important and challenging steps in achieving these goals is accurately assembling the massive quantity of short sequencing reads into full nucleic acid sequences. For research groups working with non-model or wild systems, short read assembly can pose a significant challenge due to the lack of pre-existing EST or genome reference libraries. While many publications describe the overall process of sequencing and assembly, few address the topic of how many and what types of reads are best for assembly. The goal of this project was use real world data to explore the effects of read quantity and short read quality scores on the resulting de novo assemblies. Using several samples of short reads of various sizes and qualities we produced many assemblies in an automated manner. We observe how the properties of read length, read quality, and read quantity affect the resulting assemblies and provide some general recommendations based on our real-world data set.

Variations in gene expression are essential for the evolution of novel phenotypes and for speciation. Studying allelic specific gene expression (ASGE) within interspecies hybrids provides a unique opportunity to reveal underlying mechanisms of genetic variation. Using Xiphophorus interspecies hybrid fishes and high-throughput next generation sequencing technology, we were able to assess variations between two closely related vertebrate species, Xiphophorus maculatus and Xiphophorus couchianus, and their F(1) interspecies hybrids. We constructed transcriptome-wide SNP polymorphism sets between two highly inbred X. maculatus lines (JP 163 A and B), and between X. maculatus and a second species, X. couchianus. The X. maculatus JP 163 A and B parental lines have been separated in the laboratory for ≈70 years and we were able to identify SNPs at a resolution of 1 SNP per 49 kb of transcriptome. In contrast, SNP polymorphisms between X. couchianus and X. maculatus species, which diverged ≈5-10 million years ago, were identified about every 700 bp. Using 6524 transcripts with identified SNPs between the two parental species (X. maculatus and X. couchianus), we mapped RNA-seq reads to determine ASGE within F(1) interspecies hybrids. We developed an in silico X. couchianus transcriptome by replacing 90,788 SNP bases for X. maculatus transcriptome with the consensus X. couchianus SNP bases and provide evidence that this procedure overcomes read mapping biases. Employment of the in silico reference transcriptome and tolerating 5 mismatches during read mapping allow direct assessment of ASGE in the F(1) interspecies hybrids. Overall, these results show that Xiphophorus is a tractable vertebrate experimental model to investigate how genetic variations that occur during speciation may affect gene interactions and the regulation of gene expression.

The identification of potential allergenic proteins is usually done by scanning a database of allergenic proteins and locating known allergens with a high sequence similarity. However, there is no universally accepted cut-off value for sequence similarity to indicate potential IgE cross-reactivity. Further, overall sequence similarity may be less important than discrete areas of similarity in proteins with homologous structure. To identify such areas, we first classified all allergens and their subdomains in the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP/) to their closest protein families as defined in Pfam, and identified conserved physicochemical property motifs characteristic of each group of sequences. Allergens populate only a small subset of all known Pfam families, as all allergenic proteins in SDAP could be grouped to only 130 (of 9318 total) Pfams, and 31 families contain more than four allergens. Conserved physicochemical property motifs for the aligned sequences of the most populated Pfam families were identified with the PCPMer program suite and catalogued in the webserver Motif-Mate (http://born.utmb.edu/motifmate/summary.php). We also determined specific motifs for allergenic members of a family that could distinguish them from non-allergenic ones. These allergen specific motifs should be most useful in database searches for potential allergens. We found that sequence motifs unique to the allergens in three families (seed storage proteins, Bet v 1, and tropomyosin) overlap with known IgE epitopes, thus providing evidence that our motif based approach can be used to assess the potential allergenicity of novel proteins.

Projectin and kettin are titin-like proteins mainly responsible for the high passive stiffness of insect indirect flight muscles, which is needed to generate oscillatory work during flight. Here we report the mechanical properties of kettin and projectin by single-molecule force spectroscopy. Force–extension and force-clamp curves obtained from Lethocerus projectin and Drosophila recombinant projectin or kettin fragments revealed that fibronectin type III domains in projectin are mechanically weaker (unfolding force, F u ≈ 50–150 pN) than Ig-domains (F u ≈ 150–250 pN). Among Ig domains in Sls/kettin, the domains near the N terminus are less stable than those near the C terminus. Projectin domains refolded very fast [85% at 15 s −1 (25°C)] and even under high forces (15–30 pN). Temperature affected the unfolding forces with a Q 10 of 1.3, whereas the refolding speed had a Q 10 of 2–3, probably reflecting the cooperative nature of the folding mechanism. High bending rigidities of projectin and kettin indicated that straightening the proteins requires low forces. Our results suggest that titin-like proteins in indirect flight muscles could function according to a folding-based-spring mechanism.

Assessing allele-specific gene expression (ASE) on a large scale continues to be a technically challenging problem. Certain biological phenomena, such as X chromosome inactivation and parental imprinting, affect ASE most drastically by completely shutting down the expression of a whole set of alleles. Other more subtle effects on ASE are likely to be much more complex and dependent on the genetic environment and are perhaps more important to understand since they may be responsible for a significant amount of biological diversity. Tools to assess ASE in a diploid biological system are becoming more reliable. Non-diploid systems are, however, not uncommon. In humans full or partial polyploid states are regularly found in both healthy (meiotic cells, polynucleated cell types) and diseased tissues (trisomies, non-disjunction events, cancerous tissues). In this work we have studied ASE in the medaka fish model system. We have developed a method for determining ASE in polyploid organisms from RNAseq data and we have implemented this method in a software tool set. As a biological model system we have used nuclear transplantation to experimentally produce artificial triploid medaka composed of three different haplomes. We measured ASE in RNA isolated from the livers of two adult, triploid medaka fish that showed a high degree of similarity. The majority of genes examined (82%) shared expression more or less evenly among the three alleles in both triploids. The rest of the genes (18%) displayed a wide range of ASE levels. Interestingly the majority of genes (78%) displayed generally consistent ASE levels in both triploid individuals. A large contingent of these genes had the same allele entirely suppressed in both triploids. When viewed in a chromosomal context, it is revealed that these genes are from large sections of 4 chromosomes and may be indicative of some broad scale suppression of gene expression.

The first step in the encounter between a host and a pathogen is attachment to the host epithelium. For uropathogenic Escherichia coli, these interactions are mediated by type 1 and P adhesive pili, which are long (approximately 1 microm) rods composed of more than 1000 protein subunits arranged in a helical structure. Here we used single-molecule atomic force microscopy to study the mechanical properties of type 1 pili. We found that type 1 pili readily extend under an applied force and that this extensibility is the result of unwinding the pilus rod's helical quaternary structure. The forced unraveling is also reversible, with helical rewinding taking place under considerable forces (approximately 60 pN). These data are similar to those obtained on P pili using optical tweezers, indicating that these are conserved properties of uropathogenic E. coli pili. We also show that our data can readily be reproduced using Monte Carlo simulation techniques based on a two-state kinetic model. This model provides a simple way to extrapolate the mechanical behavior of pili under a wide range of forces. We propose that type 1 pilus unraveling is an essential mechanism for absorbing physiological shear forces encountered during urinary tract infections and probably essential for adhesion and colonization of the bladder epithelium.

Scalable, integrative methods to understand mechanisms that link genetic variants with phenotypes are needed. Here we derive a mathematical expression to compute PrediXcan (a gene mapping approach) results using summary data (S-PrediXcan) and show its accuracy and general robustness to misspecified reference sets. We apply this framework to 44 GTEx tissues and 100+ phenotypes from GWAS and meta-analysis studies, creating a growing public catalog of associations that seeks to capture the effects of gene expression variation on human phenotypes. Replication in an independent cohort is shown. Most of the associations are tissue specific, suggesting context specificity of the trait etiology. Colocalized significant associations in unexpected tissues underscore the need for an agnostic scanning of multiple contexts to improve our ability to detect causal regulatory mechanisms. Monogenic disease genes are enriched among significant associations for related traits, suggesting that smaller alterations of these genes may cause a spectrum of milder phenotypes., Phenotypic variation and diseases are influenced by factors such as genetic variants and gene expression. Here, Barbeira et al. develop S-PrediXcan to compute PrediXcan results using summary data, and investigate the effects of gene expression variation on human phenotypes in 44 GTEx tissues and >100 phenotypes.

Background The release of oil resulting from the blowout of the Deepwater Horizon (DH) drilling platform was one of the largest in history discharging more than 189 million gallons of oil and subject to widespread application of oil dispersants. This event impacted a wide range of ecological habitats with a complex mix of pollutants whose biological impact is still not yet fully understood. To better understand the effects on a vertebrate genome, we studied gene expression in the salt marsh minnow Fundulus grandis, which is local to the northern coast of the Gulf of Mexico and is a sister species of the ecotoxicological model Fundulus heteroclitus. To assess genomic changes, we quantified mRNA expression using high throughput sequencing technologies (RNA-Seq) in F. grandis populations in the marshes and estuaries impacted by DH oil release. This application of RNA-Seq to a non-model, wild, and ecologically significant organism is an important evaluation of the technology to quickly assess similar events in the future. Results Our de novo assembly of RNA-Seq data produced a large set of sequences which included many duplicates and fragments. In many cases several of these could be associated with a common reference sequence using blast to query a reference database. This reduced the set of significant genes to 1,070 down-regulated and 1,251 up-regulated genes. These genes indicate a broad and complex genomic response to DH oil exposure including the expected AHR-mediated response and CYP genes. In addition a response to hypoxic conditions and an immune response are also indicated. Several genes in the choriogenin family were down-regulated in the exposed group; a response that is consistent with AH exposure. These analyses are in agreement with oligonucleotide-based microarray analyses, and describe only a subset of significant genes with aberrant regulation in the exposed set. Conclusion RNA-Seq may be successfully applied to feral and extremely polymorphic organisms that do not have an underlying genome sequence assembly to address timely environmental problems. Additionally, the observed changes in a large set of transcript expression levels are indicative of a complex response to the varied petroleum components to which the fish were exposed.