Multiple Human Tissues Research Articles

Novel replication-competent circular DNA molecules from healthy cattle serum and milk and multiple sclerosis-affected human brain tissue.

Epidemiological data point to the involvement of a cow milk factor in the etiology of multiple sclerosis (MS). Eleven circular DNA molecules closely related to transmissible spongiform encephalopathy (TSE)-associated isolate Sphinx 1.76 were isolated from healthy cattle serum, cow milk, and serum and brain tissue from MS patients.

CMTM1_v17 is a novel potential therapeutic target in breast cancer.

Chemokine-like factor (CKLF)-like MARVEL transmembrane domain-containing1 (CMTM1) consists of at least 23 alternatively spliced isoforms designated CMTM1_v1-v23. In the present study, we detected CMTM1_v17 expression in multiple human normal and tumor tissues and found that CMTM1_v17 was highly expressed in testis and many tumor tissues including breast tumor. The overexpression of CMTM1_v17 in the breast cancer cell line MDA-MB-231 promoted cell proliferation and resistance to tumor necrosis factor-α (TNF-α)-induced apoptosis. Moreover, siRNA-mediated silencing of CMTM1_v17 sensitized MDA-MB-231 cells to TNF-α-induced apoptosis. We propose that CMTM1_v17 may be a novel potential target for therapy in breast cancer patients. The present study provides insight into a novel mechanism by which CMTM1_v17 enhances cellular proliferation and abrogates TNF-α-induced apoptosis. These findings also have implications for clinical practice as they highlight the potential for therapeutic targeting of CMTM1_v17 for the treatment of breast and other cancers in which CMTM1_v17 impacts cellular proliferation and survival.

Stem potentialities of the human iris – An in situ immunohistochemical study

According to recent findings multiple human tissues harbor stem cells which, in turn, have different levels of stemness. We performed an immunohistochemical study on paraffin-embedded samples to test if the in situ stromal cells of the iris of the human eye (EI) have immune stem/progenitor phenotypes. Eviscerated post-traumatic eyes from eight patients were studied. These irises were found to contain fibroblastoid stromal cells with a CD34+/CD45+/CD105+/CD117+/DOG1+/PDGFR-α+/vimentin+/nestin−/collagen III− phenotype. These were assumed to be possible stem/progenitor cells involved in physiological processes of iridial stromal maintenance. All the vascular endothelia were CD34+/CD105+/vimentin+. Newly formed nestin+ endothelia were also found; this finding was supported by evidence of filopodia-projecting CD34+ endothelial tip cells, which demonstrated active processes of sprouting angiogenesis. The phenotype of the stromal cells also suggests a role of the circulating fibrocytes in iridial regenerative processes.

An animal model with a cardiomyocyte-specific deletion of estrogen receptor alpha: functional, metabolic, and differential network analysis.

Estrogen exerts diverse biological effects in multiple tissues in both animals and humans. Much of the accumulated knowledge on the role of estrogen receptor (ER) in the heart has been obtained from studies using ovariectomized mice, whole body ER gene knock-out animal models, ex vivo heart studies, or from isolated cardiac myocytes. In light of the wide systemic influence of ER signaling in regulating a host of biological functions in multiple tissues, it is difficult to infer the direct role of ER on the heart. Therefore, we developed a mouse model with a cardiomyocyte-specific deletion of the ERα allele (cs-ERα−/−). Male and female cs-ERα−/− mice with age/sex-matched wild type controls were examined for differences in cardiac structure and function by echocardiogram and differential gene expression microarray analysis. Our study revealed sex-differences in structural parameters in the hearts of cs-ERα−/− mice, with minimal functional differences. Analysis of microarray data revealed differential variations in the expression of 208 genes affecting multiple transcriptional networks. Furthermore, we report sex-specific differences in the expression of 56 genes. Overall, we developed a mouse model with cardiac-specific deletion of ERα to characterize the role of ERα in the heart independent of systemic effects. Our results suggest that ERα is involved in controlling the expression of diverse genes and networks in the cardiomyocyte in a sex-dependent manner.

Ubiquitin Specific Protease 26 (USP26) expression analysis in human testicular and extragonadal tissues indicates diverse action of USP26 in cell differentiation and tumorigenesis.

Ubiquitin specific protease 26 (USP26), a deubiquitinating enzyme, is highly expressed early during murine spermatogenesis, in round spermatids, and at the blood-testis barrier. USP26 has also been recognized as a regulator of androgen receptor (AR) hormone-induced action involved in spermatogenesis and steroid production in in vitro studies. Prior mutation screening of USP26 demonstrated an association with human male infertility and low testosterone production, but protein localization and expression in the human testis has not been characterized previously. USP26 expression analysis of mRNA and protein was completed using murine and human testis tissue and human tissue arrays. USP26 and AR mRNA levels in human testis were quantitated using multiplex qRT-PCR. Immunofluorescence colocalization studies were performed with formalin-fixed/paraffin-embedded and frozen tissues using primary and secondary antibodies to detect USP26 and AR protein expression. Human microarray dot blots were used to identify protein expression in extra-gonadal tissues. For the first time, expression of USP26 and colocalization of USP26 with androgen receptor in human testis has been confirmed predominantly in Leydig cell nuclei, with less in Leydig cell cytoplasm, spermatogonia, primary spermatocytes, round spermatids, and Sertoli cells. USP26 likely affects regulatory proteins of early spermatogenesis, including androgen receptor with additional activity in round spermatids. This X-linked gene is not testis-specific, with USP26 mRNA and protein expression identified in multiple other human organ tissues (benign and malignant) including androgen-dependent tissues such as breast (myoepithelial cells and secretory luminal cells) and thyroid tissue (follicular cells). USP26/AR expression and interaction in spermatogenesis and androgen-dependent cancer warrants additional study and may prove useful in diagnosis and management of male infertility.

Programmed Translational Readthrough Generates Antiangiogenic VEGF-Ax

Translational readthrough, observed primarily in less complex organisms from viruses to Drosophila, expands the proteome by translating select transcripts beyond the canonical stop codon. Here, we show that vascular endothelial growth factor A (VEGFA) mRNA in mammalian endothelial cells undergoes programmed translational readthrough (PTR) generating VEGF-Ax, an isoform containing a unique 22-amino-acid C terminus extension. A cis-acting element in the VEGFA 3' UTR serves a dual function, not only encoding the appended peptide but also directing the PTR by decoding the UGA stop codon as serine. Heterogeneous nuclear ribonucleoprotein (hnRNP) A2/B1 binds this element and promotes readthrough. Remarkably, VEGF-Ax exhibits antiangiogenic activity in contrast to the proangiogenic activity of VEGF-A. Pathophysiological significance of VEGF-Ax is indicated by robust expression in multiple human tissues but depletion in colon adenocarcinoma. Furthermore, genome-wide analysis revealed AGO1 and MTCH2 as authentic readthrough targets. Overall, our studies reveal a novel protein-regulated PTR event in a vertebrate system.

A global profiling of alternative splicing exon expression patterns across multiple human tissues with RNA-Seq data

A global profiling of alternative splicing exon expression patterns across multiple human tissues with RNA-Seq data

CDKN2B expression and subcutaneous adipose tissue expandability: Possible influence of the 9p21 atherosclerosis locus

Risk alleles within a gene desert at the 9p21 locus constitute the most prevalent genetic determinant of cardiovascular disease. Previous research has demonstrated that 9p21 risk variants influence gene expression in vascular tissues, yet the biological mechanisms by which this would mediate atherosclerosis merits further investigation. To investigate possible influences of this locus on other tissues, we explored expression patterns of 9p21-regulated genes in a panel of multiple human tissues and found that the tumor suppressor CDKN2B was highly expressed in subcutaneous adipose tissue (SAT). CDKN2B expression was regulated by obesity status, and this effect was stronger in carriers of 9p21 risk alleles. Covariation between expression of CDKN2B and genes implemented in adipogenesis was consistent with an inhibitory effect of CDKN2B on SAT proliferation. Moreover, studies of postprandial triacylglycerol clearance indicated that CDKN2B is involved in down-regulation of SAT fatty acid trafficking. CDKN2B expression in SAT correlated with indicators of ectopic fat accumulation, including markers of hepatic steatosis. Among genes regulated by 9p21 risk variants, CDKN2B appears to play a significant role in the regulation of SAT expandability, which is a strong determinant of lipotoxicity and therefore might contribute to the development of atherosclerosis.

Ahnak functions as a tumor suppressor via modulation of TGFβ/Smad signaling pathway

We provide detailed mechanisms of Ahnak-mediated potentiation of transforming growth factor β (TGFβ) signaling, which leads to a negative regulation of cell growth. We show that Smad3 interacts with Ahnak through MH2 domain and that Ahnak stimulates Smad3 localization into nucleus leading to potentiating TGFβ-induced transcriptional activity of R-Smad. Moreover, overexpression of Ahnak resulted in growth retardation and cell cycle arrest through downregulation of c-Myc and cyclin D1/D2. We describe results from analyses of Ahnak−/− mouse model expressing middle T antigen in a mammary gland-specific manner (MMTVTg/+Ahnak−/−), which showed significantly progressed hyperplasia of mammary glands compared with MMTVTg/+Ahnak+/+. Finally, we screened multiple human breast cancer tissues and showed that the expression of Ahnak in cancer tissues is lower than that in control tissues by 50%. Taken together, these data indicate that Ahnak mediates a negative regulation of cell growth and acts as novel tumor suppressor through potentiation of TGFβ signaling.

A tissue-specific gene expression template portrays heart development and pathology

Congenital heart defects (CHD) are the most common cause of death in children under the age of 1. Tetralogy of Fallot (TOF) is a severe CHD that results from developmental defects in the conotruncal outflow tract. Recently, a tissue-specific gene expression template (GET) was derived from microarray data that accurately characterized multiple normal human tissues. We used the GET to examine spatial, temporal, and a pathological condition (TOF) within a single organ, the heart. The GET, as previously defined, generally identified temporal and spatial differences in the cardiac tissue. Differences in the stoichiometry of the GET reflected the severe developmental disturbance associated with TOF. Our analysis suggests that the homoeostatic equilibrium assessed by the GET at the inter-organ level is generally maintained at the intra-organ level as well.

Gene Family Level Comparative Analysis of Gene Expression in Mammals Validates the Ortholog Conjecture

The ortholog conjecture (OC), which is central to functional annotation of genomes, posits that orthologous genes are functionally more similar than paralogous genes at the same level of sequence divergence. However, a recent study challenged the OC by reporting a greater functional similarity, in terms of Gene Ontology (GO) annotations and expression profiles, among within-species paralogs compared with orthologs. These findings were taken to indicate that functional similarity of homologous genes is primarily determined by the cellular context of the genes, rather than evolutionary history. However, several subsequent studies suggest that GO annotations and microarray data could artificially inflate functional similarity between paralogs from the same organism. We sought to test the OC using approaches distinct from those used in previous studies. Analysis of a large RNAseq data set from multiple human and mouse tissues shows that expression similarity (correlations coefficients, rank‘s, or Z-scores) between orthologs is substantially greater than that for between-species paralogs with the same sequence divergence, in agreement with the OC and the results of recent detailed analyses. These findings are further corroborated by a fine-grain analysis in which expression profiles of orthologs and paralogs were compared separately for individual gene families. Expression profiles of within-species paralogs are more strongly correlated than profiles of orthologs but it is shown that this is caused by high background noise, that is, correlation between profiles of unrelated genes in the same organism. Z-scores and rank scores show a nonmonotonic dependence of expression profile similarity on sequence divergence. This complexity of gene expression evolution after duplication might be at least partially caused by selection for protein dosage rebalancing following gene duplication.

Understanding the phenotypic impact of non-coding sequence variation

Molecular phenotypes are important phenotypes that informs about genetic and environmental effects on cellular state. The elucidation of the genetics of gene expression and other cellular phenotypes are highly informative of the impact of genetic variants in the cell and the subsequent consequences in the organism. In this talk, I will discuss recent advances in three key areas of the analysis of the genomics of gene expression and cellular phenotypes in human populations and multiple tissues and how this assists in the interpretation of regulatory networks and human disease variants. I will also discuss how the recent advances in next generation sequencing and functional genomics are bringing closer our hopes for personalized medicine.

Tumor Cell Expression of Vascular Endothelial Growth Factor Receptor 2 Is an Adverse Prognostic Factor in Patients with Squamous Cell Carcinoma of the Lung

A robust immunohistochemical (IHC) assay for VEGFR2 was developed to investigate its utility for patient tailoring in clinical trials. The sensitivity, specificity, and selectivity of the IHC assay were established by siRNA knockdown, immunoblotting, mass spectrometry, and pre-absorption experiments. Characterization of the assay included screening a panel of multiple human cancer tissues and an independent cohort of non-small cell lung carcinoma (NSCLC, n = 118) characterized by TTF-1, p63, CK5/6, and CK7 IHC. VEGFR2 immunoreactivity was interpreted qualitatively (VEGFR2 positive/negative) in blood vessels and by semi-quantitative evaluation using H-scores in tumor cells (0–300). Associations were determined among combinations of VEGFR2 expression in blood vessels and tumor cells, and clinico-pathologic characteristics (age, sex, race, histologic subtype, disease stage) and overall survival using Kaplan-Meier analyses and appropriate statistical models. VEGFR2 expression both in blood vessels and in tumor cells in carcinomas of the lung, cervix, larynx, breast, and others was demonstrated. In the validation cohort, 99/118 (83.9%) NSCLC tissues expressed VEGFR2 in the blood vessels and 46/118 (39.0%) showed high tumor cell positivity (H-score ≥10). Vascular and tumor cell expression were inversely correlated (p = 0.0175). High tumor cell expression of VEGFR2 was associated with a 3.7-fold reduction in median overall survival in lung squamous-cell carcinoma (SCC, n = 25, p = 0.0134). The inverse correlation between vascular and tumor cell expression of VEGFR2 and the adverse prognosis associated with high VEGFR2 expression in immunohistochemically characterized pulmonary SCC are new findings with potential therapeutic implications. The robustness of this novel IHC assay will support further evaluation of its utility for patient tailoring in clinical trials of antiangiogenic agents.

MC13-0057 A comprehensive IT platform to support GTEx operation

Background: Designed to build upon Genome Wide Association Study (GWAS) findings, the NIH Common Fund’s Genotype-Tissue Expression (GTEx) project aims to study gene expression and regulation across multiple human tissues (30+ tissue types) from approximately 1000 healthy normal donors. It is expected to provide valuable insights into gene regulation and its tissue specificity, identify correlation between genetic variations and variations in gene expression levels as expression quantitative trait loci (eQTLs), and help to understand inherited susceptibility to diseases. Purpose/Objective: To meet the challenge of GTEx requirements for collecting and tracking high quality biospecimen samples, a custom-built software system named Comprehensive Data Resources (CDR) was developed to support sample collection work flow, clinical data entry, case management, and review and curation of study data. Materials and Methods: CDR is built with combination of technologies from Grails, Oracle, Groovy, jQuery, Apache Solr. Results: The CDR provides secure user access to case and sample data based on pre-defined roles and privileges. Personally Identifiable Information (PII) and Protected Health Information (PHI) are restricted to a limited data set (LDS) and to authorized users through dynamic content redaction. Intuitive graphic user interfaces for the Biopecimen Source Sites (BSS) streamline data entry workflow by strictly following SOPs for sample collection and processing. Contextual automated data checks and business rule validations confirm data integrity and SOP adherence simultaneously. Web services APIs allow the Pathology Resource Center to access digital imaging data from tissue slides housed remotely at the Comprehensive Biospecimen Resource (CBR). API’s connect to CBR’s LIMS systems for real-time sample inventory data. De-identified GTEx data is provided via a private API with the Broad Institute (LDACC) before the final release into dbGaP. The reporting and analytics module supports data analysis and aggregation, report generation and real-time operational data snapshots. Conclusions: CDR is a distributed web-based system designed to support GTEx operation from pilot phase to full scale-up stage. It manages and maintains multi-dimensional data models around each donor case (average 500+ data elements/case). As an efficient case management tool capable of connecting to various remote informatics systems, CDR could be adapted to the broader biobanking community with the flexibility of building user-defined work flows in the system.

Quantitative expression analysis of two NAGPA isoforms in multiple human cDNA tissue panels

Uncovering enzyme (UCE), encoded by the human NAGPA, is a trans-Golgi enzyme that adds the mannose-6-phosphate recognition tag on lysosomal enzymes destined for the lysosome. Mutations in NAGPA are known to cause stuttering, a common speech disorder with unknown etiology. The human NAGPA gene is transcribed into two different forms, probably due to alternative splicing. One of them, known as a brain isoform, is lacking exon 8 (102-bp). We performed quantitative real-time PCR for the NAGPA brain and non-brain isoforms in a cDNA panel originating from 16 human tissues and 24 sub-brain regions. According to our findings, the relative quantity of the NAGPA brain isoform in the brain was 4.7 times more than that in the control cDNA, a pooled mixture of equal amounts of cDNAs from the 16 different tissues. Further analysis using the cDNA panel originating from 24 different sub-brain regions revealed that the cerebral cortex contained the largest amount of NAGPA brain isoform. Relative quantity in the cerebral cortex was 8.6 times more than that in the control cDNA (P=0.00004). The lowest quantity of this isoform was detected in cDNA from the pituitary gland. In conclusion, findings of the current study suggest that the cerebral cortex, expressing the highest quantity of the NAGPA brain isoform, might be the region associated with speech function.

The Immunoglobulin Superfamily Protein Differentiation of Embryonic Stem Cells 1 (Dies1) Has a Regulatory Role in Preadipocyte to Adipocyte Conversion

Differentiation of Embryonic Stem Cells 1 (Dies1) was recently identified as a novel type I immunoglobulin (IgG) domain-containing plasma membrane protein important for effective differentiation of a murine pluripotent embryonic stem cell line. In this setting, Dies1 enhances bone morphogenetic protein 4 (BMP4) signaling. Here we show Dies1 transcript expression is induced ∼225-fold during in vitro adipogenesis of 3T3-L1 murine preadipocytes. Immunocytochemical imaging using ectopic expression of Flag-tagged Dies1 in 3T3-L1 adipocytes revealed localization to the adipocyte plasma membrane. Modulation of adipocyte phenotype with with tumor necrosis factor-α (TNFα) treatment or by siRNA knockdown of the master pro-adipogenic transcription factor peroxisome proliferator activated receptor gamma (PPARγ) resulted in a 90% and 60% reduction of Dies1 transcript levels, respectively. Moreover, siRNA-mediated Dies1 knockdown in 3T3-L1 preadipocytes inhibited adipogenic conversion. Such cultures had a 35% decrease in lipid content and a 45%–65% reduction in expression of key adipocyte transcripts, including that for PPARγ. The standard protocol for full in vitro adipogenic conversion of committed preadipocytes, such as 3T3-L1, does not include BMP4 treatment. Thus we posit the positive role of Dies1 in adipogenesis, unlike that for Dies1 in differentiation of embryonic stem cells, does not include its pro-BMP4 effects. In support of this idea, 3T3-L1 adipocytes knocked down for Dies1 did not evidence decreased phospho-Smad1 levels upon BMP4 exposure. qPCR analysis of Dies1 transcript in multiple murine and human tissues reveals high enrichment in white adipose tissue (WAT). Interestingly, we observed a 10-fold induction of Dies1 transcript in WAT of fasted vs. fed mice, suggesting a role for Dies1 in nutritional response of mature fat cells in vivo. Together our data identify Dies1 as a new differentiation-dependent adipocyte plasma membrane protein whose expression is required for effective adipogenesis and that may also play a role in regard to nutritional status in WAT.

Quantitative assessment of Tet-induced oxidation products of 5-methylcytosine in cellular and tissue DNA

Recent studies showed that Ten-eleven translocation (Tet) family dioxygenases can oxidize 5-methyl-2’-deoxycytidine (5-mdC) in DNA to yield the 5-hydroxymethyl, 5-formyl and 5-carboxyl derivatives of 2’-deoxycytidine (5-HmdC, 5-FodC and 5-CadC). 5-HmdC in DNA may be enzymatically deaminated to yield 5-hydroxymethyl-2’-deoxyuridine (5-HmdU). After their formation at CpG dinucleotide sites, these oxidized pyrimidine nucleosides, particularly 5-FodC, 5-CadC, and 5-HmdU, may be cleaved from DNA by thymine DNA glycosylase, and subsequent action of base-excision repair machinery restores unmethylated cytosine. These processes are proposed to be important in active DNA cytosine demethylation in mammals. Here we used a reversed-phase HPLC coupled with tandem mass spectrometry (LC-MS/MS/MS) method, along with the use of stable isotope-labeled standards, for accurate measurements of 5-HmdC, 5-FodC, 5-CadC and 5-HmdU in genomic DNA of cultured human cells and multiple mammalian tissues. We found that overexpression of the catalytic domain of human Tet1 led to marked increases in the levels of 5-HmdC, 5-FodC and 5-CadC, but only a modest increase in 5-HmdU, in genomic DNA of HEK293T cells. Moreover, 5-HmdC is present at a level that is approximately 2–3 and 3–4 orders of magnitude greater than 5-FodC and 5-CadC, respectively, and 35–400 times greater than 5-HmdU in the mouse brain and skin, and human brain. The robust analytical method built a solid foundation for dissecting the molecular mechanisms of active cytosine demethylation, for measuring these 5-mdC derivatives and assessing their involvement in epigenetic regulation in other organisms and for examining whether these 5-mdC derivatives can be used as biomarkers for human diseases.

Abstract 1801: ERK-mediated regulation of NFAT3 enhances CXCR4 expression in HeyA8 ovarian cell line.

Abstract The G-protein coupled chemokine (C-X-C motif) receptor CXCR4 is linked to cancer, HIV, and WHIM (Warts, Hypogammaglobulinemia, Infections, and Myelokathexis) syndrome. While CXCR4 is widely studied and characterized extensively to be over-expressed in multiple human tumor tissues, we have seen poor surface CXCR4 expression in many solid tumor cell lines in vitro. To better understand how CXCR4 is regulated, cells grown under 3-D spheroid conditions were compared to normal adherent culturing conditions. When cultured as 3D spheroids, HeyA8 cells showed a dramatic increase in surface CXCR4 protein levels as well as mRNA transcripts. Furthermore, HeyA8 3D spheroids showed a decrease in pERK levels when compared to adherent cells. Using an inhibitor of the MEK-ERK pathway, the treatment of adherent HeyA8 cells with U0126 resulted in a significant increase in surface CXCR4 expression similar to 3D spheroids. Additional investigation using PCR array showed a wide range of transcription factors being up-regulated but notably a > 20 fold increase in NFAT3 transcription factor. Using a calcineurin inhibitor, cyclosporin A diminished the U0126 treated up-regulation of CXCR4 surface expression. Finally, Chromatin Immunoprecipitation analysis showed direct binding of NFAT3 on the CXCR4 promoter to increase CXCR4 expression in HeyA8 ovarian cell line. Taken together, our results suggest that phospho-ERK levels regulate CXCR4 expression through the NFAT3 signaling pathway. Citation Format: Keven Huang, Christine Kiefer, Adeela Kamal. ERK-mediated regulation of NFAT3 enhances CXCR4 expression in HeyA8 ovarian cell line. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 1801. doi:10.1158/1538-7445.AM2013-1801 Note: This abstract was not presented at the AACR Annual Meeting 2013 because the presenter was unable to attend.

Expression profile and mitochondrial colocalization of Tdp1 in peripheral human tissues

Tyrosyl-DNA phosphodiesterase 1 (Tdp1) is a DNA repair enzyme that processes blocked 3' ends of DNA breaks. Functional loss of Tdp1 causes spinocerebellar ataxia with axonal neuropathy type 1 (SCAN1). Based on the prominent cytoplasmic expression of Tdp1 in the neurons presumably affected in SCAN1, we hypothesized that Tdp1 participates in the repair of mitochondrial DNA. As a step toward testing this hypothesis, we profiled Tdp1 expression in different human tissues by immunohistochemistry and immunofluorescence respectively and determined whether Tdp1 was expressed in the cytoplasm of tissues other than the neurons. We found that Tdp1 was ubiquitously expressed and present in the cytoplasm of many cell types. Within human skeletal muscle and multiple mouse tissues, Tdp1 partially colocalized with the mitochondria. In cultured human dermal fibroblasts, Tdp1 redistributed to the cytoplasm and partially colocalized with mitochondria following oxidative stress. These studies suggest that one role of cytoplasmic Tdp1 is the repair of mitochondrial DNA lesions arising from oxidative stress.

Expression Profile of Ectopic Olfactory Receptors Determined by Deep Sequencing

Olfactory receptors (ORs) provide the molecular basis for the detection of volatile odorant molecules by olfactory sensory neurons. The OR supergene family encodes G-protein coupled proteins that belong to the seven-transmembrane-domain receptor family. It was initially postulated that ORs are exclusively expressed in the olfactory epithelium. However, recent studies have demonstrated ectopic expression of some ORs in a variety of other tissues. In the present study, we conducted a comprehensive expression analysis of ORs using an extended panel of human tissues. This analysis made use of recent dramatic technical developments of the so-called Next Generation Sequencing (NGS) technique, which encouraged us to use open access data for the first comprehensive RNA-Seq expression analysis of ectopically expressed ORs in multiple human tissues. We analyzed mRNA-Seq data obtained by Illumina sequencing of 16 human tissues available from Illumina Body Map project 2.0 and from an additional study of OR expression in testis. At least some ORs were expressed in all the tissues analyzed. In several tissues, we could detect broadly expressed ORs such as OR2W3 and OR51E1. We also identified ORs that showed exclusive expression in one investigated tissue, such as OR4N4 in testis. For some ORs, the coding exon was found to be part of a transcript of upstream genes. In total, 111 of 400 OR genes were expressed with an FPKM (fragments per kilobase of exon per million fragments mapped) higher than 0.1 in at least one tissue. For several ORs, mRNA expression was verified by RT-PCR. Our results support the idea that ORs are broadly expressed in a variety of tissues and provide the basis for further functional studies.