Survey of the Heritability and Sparse Architecture of Gene Expression Traits across Human Tissues.

Understanding the genetic architecture of gene expression traits is key to elucidating the underlying mechanisms of complex traits. Here, for the first time, we perform a systematic survey of the heritability and the distribution of effect sizes across all representative tissues in the human body. We find that local h2 can be relatively well characterized with 59% of expressed genes showing significant h2 (FDR < 0.1) in the DGN whole blood cohort. However, current sample sizes (n ≤ 922) do not allow us to compute distal h2. Bayesian Sparse Linear Mixed Model (BSLMM) analysis provides strong evidence that the genetic contribution to local expression traits is dominated by a handful of genetic variants rather than by the collective contribution of a large number of variants each of modest size. In other words, the local architecture of gene expression traits is sparse rather than polygenic across all 40 tissues (from DGN and GTEx) examined. This result is confirmed by the sparsity of optimal performing gene expression predictors via elastic net modeling. To further explore the tissue context specificity, we decompose the expression traits into cross-tissue and tissue-specific components using a novel Orthogonal Tissue Decomposition (OTD) approach. Through a series of simulations we show that the cross-tissue and tissue-specific components are identifiable via OTD. Heritability and sparsity estimates of these derived expression phenotypes show similar characteristics to the original traits. Consistent properties relative to prior GTEx multi-tissue analysis results suggest that these traits reflect the expected biology. Finally, we apply this knowledge to develop prediction models of gene expression traits for all tissues. The prediction models, heritability, and prediction performance R2 for original and decomposed expression phenotypes are made publicly available (https://github.com/hakyimlab/PrediXcan).

X. couchianus and X. hellerii genome models provide genomic variation insight among Xiphophorus species.

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.

Genomic and physiological footprint of the Deepwater Horizon oil spill on resident marsh fishes.

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.

RNA-Seq reveals complex genetic response to Deepwater Horizon oil release in Fundulus grandis.

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.

Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates.

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.

Neonatal thyroxine activation modifies epigenetic programming of the liver.

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 < 0.00001). These data explain how the systemic hormone T3 acts locally during development to define future expression of hepatic genes.

Multi-step screening of neoantigens' HLA- and TCR-interfaces improves prediction of survival.

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.

Characteristic motifs for families of allergenic proteins.

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 MotifMate (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.

Mechanical stability and differentially conserved physical-chemical properties of titin Ig-domains.

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.

Structural analysis of linear and conformational epitopes of allergens.

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.

Single-molecule force spectroscopy reveals a stepwise unfolding of Caenorhabditis elegans giant protein kinase domains.

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 beta-sandwich domains (30-150 pN). 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.

Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics.

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.

Aorta macrophage inflammatory and epigenetic changes in a murine model of obstructive sleep apnea: Potential role of CD36.

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.

GR and ER Coactivation Alters the Expression of Differentiation Genes and Associates with Improved ER+ Breast Cancer Outcome.

UNLABELLED: 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.

UVB-induced gene expression in the skin of Xiphophorus maculatus Jp 163 B.

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.

Novel method for analysis of allele specific expression in triploid Oryzias latipes reveals consistent pattern of allele exclusion.

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.

Alternative strategies for development of a reference transcriptome for quantification of allele specific expression in organisms having sparse genomic resources.

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.