SCFFBXW11 Complex Targets Interleukin-17 Receptor A for Ubiquitin-Proteasome-Mediated Degradation.

Interleukin-17 (IL-17) is a pro-inflammatory cytokine that participates in innate and adaptive immune responses and plays an important role in host defense, autoimmune diseases, tissue regeneration, metabolic regulation, and tumor progression. Post-translational modifications (PTMs) are crucial for protein function, stability, cellular localization, cellular transduction, and cell death. However, PTMs of IL-17 receptor A (IL-17RA) have not been investigated. Here, we show that human IL-17RA was targeted by F-box and WD repeat domain-containing 11 (FBXW11) for ubiquitination, followed by proteasome-mediated degradation. We used bioinformatics tools and biochemical techniques to determine that FBXW11 ubiquitinated IL-17RA through a lysine 27-linked polyubiquitin chain, targeting IL-17RA for proteasomal degradation. Domain 665-804 of IL-17RA was critical for interaction with FBXW11 and subsequent ubiquitination. Our study demonstrates that FBXW11 regulates IL-17 signaling pathways at the IL-17RA level.

Five Critical Gene-Based Biomarkers With Optimal Performance for Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is one of the most fatal cancers in the world. There is an urgent need to understand the molecular background of HCC to facilitate the identification of biomarkers and discover effective therapeutic targets. Published transcriptomic studies have reported a large number of genes that are individually significant for HCC. However, reliable biomarkers remain to be determined. In this study, built on max-linear competing risk factor models, we developed a machine learning analytical framework to analyze transcriptomic data to identify the most miniature set of differentially expressed genes (DEGs). By analyzing 9 public whole-transcriptome datasets (containing 1184 HCC samples and 672 nontumor controls), we identified 5 critical differentially expressed genes (DEGs) (ie, CCDC107, CXCL12, GIGYF1, GMNN, and IFFO1) between HCC and control samples. The classifiers built on these 5 DEGs reached nearly perfect performance in identification of HCC. The performance of the 5 DEGs was further validated in a US Caucasian cohort that we collected (containing 17 HCC with paired nontumor tissue). The conceptual advance of our work lies in modeling gene-gene interactions and correcting batch effect in the analytic framework. The classifiers built on the 5 DEGs demonstrated clear signature patterns for HCC. The results are interpretable, robust, and reproducible across diverse cohorts/populations with various disease etiologies, indicating the 5 DEGs are intrinsic variables that can describe the overall features of HCC at the genomic level. The analytical framework applied in this study may pave a new way for improving transcriptome profiling analysis of human cancers.

Subtyping of hepatocellular adenoma: a machine learning-based approach.

Subtyping of hepatocellular adenoma (HCA) is an important task in practice as different subtypes may have different clinical outcomes and management algorithms. Definitive subtyping is currently dependent on immunohistochemical and molecular testing. The association between some morphologic/clinical features and HCA subtypes has been reported; however, the predictive performance of these features has been controversial. In this study, we attempted machine learning based methods to select an efficient and parsimonious set of morphologic/clinical features for differentiating a HCA subtype from the others, and then assessed the performance of the selected features in identifying the correct subtypes. We first examined 50 liver HCA resection specimens collected at the University of Washington and Kobe University/Kings College London, including HNF1α-mutated HCA (H-HCA) (n = 16), inflammatory HCA (I-HCA) (n = 20), beta-catenin activated HCA (ß-HCA) (n = 8), and unclassified HCA (U-HCA) (n = 6). Twenty-six morphologic/clinical features were assessed. We used LASSO (least absolute shrinkage and selection operator) to select key features that could differentiate a subtype from the others. We further performed SVM (support vector machine) analysis to assess the performance (sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and accuracy) of the selected features in HCA subtyping in an independent cohort of liver resection samples (n = 20) collected at the University of Wisconsin-Madison. With some overlap, different combinations of morphologic/clinical features were selected for each subtype. Based on SVM analysis, the selected features classified HCA into correct subtypes with an overall accuracy of at least 80%. Our findings are useful for initial diagnosis and subtyping of HCA, especially in clinical settings without access to immunohistochemical and molecular assays.

RNA sequencing analysis of hepatocellular carcinoma identified oxidative phosphorylation as a major pathologic feature.

Dysregulation of expression of functional genes and pathways plays critical roles in the etiology and progression of hepatocellular carcinoma (HCC). Next generation-based RNA sequencing (RNA-seq) offers unparalleled power to comprehensively characterize HCC at the whole transcriptome level. In this study, 17 fresh-frozen HCC samples with paired non-neoplastic liver tissue from Caucasian patients undergoing liver resection or transplantation were used for RNA-seq analysis. Pairwise differential expression analysis of the RNA-seq data was performed to identify genes, pathways, and functional terms differentially regulated in HCC versus normal tissues. At a false discovery rate (FDR) of 0.10, 13% (n = 4335) of transcripts were up-regulated and 19% (n = 6454) of transcripts were down-regulated in HCC versus non-neoplastic tissue. Eighty-five Kyoto Encyclopedia of Genes and Genomes pathways were differentially regulated (FDR, <0.10), with almost all pathways (n = 83) being up-regulated in HCC versus non-neoplastic tissue. Among the top up-regulated pathways was oxidative phosphorylation (hsa00190; FDR, 1.12E-15), which was confirmed by Database for Annotation, Visualization, and Integrated Discovery (DAVID) gene set enrichment analysis. Consistent with potential oxidative stress due to activated oxidative phosphorylation, DNA damage-related signals (e.g., the up-regulated hsa03420 nucleotide excision repair [FDR, 1.14E-04] and hsa03410 base excision repair [FDR, 2.71E-04] pathways) were observed. Among down-regulated genes (FDR, <0.10), functional terms related to cellular structures (e.g., cell membrane [FDR, 3.05E-21] and cell junction [FDR, 2.41E-07], were highly enriched, suggesting compromised formation of cellular structure in HCC at the transcriptome level. Interestingly, the olfactory transduction (hsa04740; FDR, 1.53E-07) pathway was observed to be down-regulated in HCC versus non-neoplastic tissue, suggesting impaired liver chemosensory functions in HCC. Our findings suggest oxidative phosphorylation and the associated DNA damage may be the major driving pathologic feature in HCC.

Dcx expression defines a subpopulation of Gdf5 + cells with chondrogenic potentials in E12.5 mouse embryonic limbs.

Growth differentiation factor 5 (Gdf5) and doublecortin (Dcx) genes are both expressed in joint interzone cells during synovial joint development. In this study, we re-analyzed the single cell RNA-sequencing data (Gene Expression Omnibus GSE151985) generated from Gdf5 + cells of mouse knee joints at embryonic stages of E12.5, E13.5, E14.5, and E15.5, with a new focus on Dcx. We found that Dcx expression was enriched in clusters of Gdf5 + cells, with high expression levels of pro-chondrogenic genes including sex determining region Y-box transcription factor 5 (Sox5), Sox6, Sox9, Gdf5, versican, matrilin 4, collagen type II α 1 chain (Col2a1), Col9a1, Col9a2, and Col9a3 at E12.5. Dcx + and Dcx - cells had differential gene expression profiles. The up-regulated genes in Dcx + vs. Dcx - cells at E12.5 and E13.5 were enriched in chondrocyte differentiation and cartilage development, whereas those genes up-regulated at E14.5 and E15.5 were enriched in RNA splicing, protein stability, cell proliferation, and cell growth. Gene expression profiles in Dcx + cells showed rapid daily changes from E12.5 to E15.5, with limited number of genes shared across the time period. Expression of Gdf5, Sox5, Sox6, melanoma inhibitory activity, noggin, odd-skipped related transcription factor 2, matrilin 4, and versican was positively correlated with Dcx expression. Our results demonstrate that Dcx expression defines a subpopulation of Gdf5 + cells with chondrogenic potentials in E12.5 mouse embryonic limbs.

Role of DNA Flow Cytometry in the Diagnosis of Malignancy in Bile Duct Biopsies Using Paraffin-Embedded Tissue.

OBJECTIVES: Histopathologic evaluation of bile biopsies for biliary strictures is frequently challenging and is affected by interobserver disagreement. Reliable ancillary tests that can help differentiate benign from malignant are not available. This study aimed to evaluate whether DNA content abnormalities detected by flow cytometry on formalin-fixed, paraffin-embedded (FFPE) tissue can help differentiate benign/reactive, dysplastic from malignant cell populations in bile duct biopsies. METHODS: We performed DNA flow cytometry on 30 FFPE bile duct biopsies in 5 well-defined diagnostic categories: (1) negative for dysplasia (NED), (2) low-grade dysplasia (LGD), (3) high-grade dysplasia (HGD), (4) carcinoma (CA), and (5) indefinite for dysplasia (IND). RESULTS: Abnormal DNA content was detected in 0 NED, 5 LGD (62.5%), 2 HGD (33.3%), 3 CA (60%), and 4 IND (80%) samples. As a diagnostic marker, the estimated sensitivity, specificity, positive predictive value, and negative predictive value were 63%, 100%, 100%, and 50%, respectively, for diagnosing HGD or CA. CONCLUSIONS: DNA flow cytometry analysis is a useful ancillary test for the interpretation of bile duct biopsies. DNA content abnormalities, when correlated with histologic findings, will not only help confirm the morphologic impression but also identify patients who are at a higher risk of developing malignancy.

Increased Proviral DNA in Circulating Cells Correlates with Plasma Viral Rebound in Simian Immunodeficiency Virus-Infected Rhesus Macaques after Antiretroviral Therapy Interruption.

The human immunodeficiency virus (HIV) reservoir is responsible for persistent viral infection, and a small number of mosaic latent cellular reservoirs promote viral rebound upon antiretroviral therapy interruption, which is the major obstacle to a cure. However, markers that determine effective therapy and viral rebound posttreatment interruption remain unclear. In this study, we comprehensively and longitudinally tracked dynamic decay of cell-associated viral RNA/DNA in systemic and lymphoid tissues in simian immunodeficiency virus (SIV)-infected rhesus macaques on prolonged combined antiretroviral therapy (cART) and evaluated predictors of viral rebound after treatment cessation. The results showed that suppressive ART substantially reduced plasma SIV RNA, cell-associated unspliced, and multiply spliced SIV RNA to undetectable levels, yet viral DNA remained detectable in systemic tissues and lymphoid compartments throughout cART. Intriguingly, a rapid increase of integrated proviral DNA in peripheral mononuclear cells was detected once treatment was withdrawn, accompanied by the emergence of detectable plasma viral load. Notably, the increase of peripheral proviral DNA after treatment interruption correlated with the emergence and degree of viral rebound. These findings suggest that measuring total viral DNA in SIV infection may be a relatively simple surrogate marker of reservoir size and may predict viral rebound after treatment interruption and inform treatment strategies.IMPORTANCE Viral reservoirs are involved in persistent HIV infection, and a small number of mosaic latent cellular reservoirs promote viral rebound upon analytical treatment interruption, which is the major obstacle to a cure. However, early indicators that can predict resurgence of viremia after treatment interruption may aid treatment decisions in people living with HIV. Utilizing the rhesus macaque model, we demonstrated that increased proviral DNA in peripheral cells after treatment interruption, rather than levels of proviral DNA, was a useful marker to predict the emergence and degree of viral rebound after treatment interruption, providing a rapid approach for monitoring HIV rebound and informing decisions.

Interobserver agreement in pathologic evaluation of bile duct biopsies.

Intraductal biopsy is commonly used for preoperative evaluation of the etiology of biliary strictures. Interpretation of intraductal biopsies is frequently challenging. The diagnosis often suffers from interobserver disagreement, which has not been studied in the literature. We sought to assess interobserver concordance in the interpretation of intraductal biopsies. Eighty-five biopsies were retrieved, falling into five diagnostic categories: negative for dysplasia (NED), indefinite for dysplasia (IND), low-grade dysplasia (LGD), high-grade dysplasia (HGD), and carcinoma (CA). Eight gastrointestinal pathologists blindly reviewed all the slides. Agreement among pathologists was analyzed using Fleiss κ and weighted concordance coefficient S∗. A face-to-face consensus/training session was held to discuss the classification criteria, followed by a second round review. The overall interobserver agreement was fair in the first round review (κ = 0.39; S∗ = 0.56) and improved to moderate in the second round review (κ = 0.48; S∗ = 0.69). The agreement before and after consensus meeting was substantial to nearly perfect for CA (κ = 0.65, S∗ = 0.83; and κ = 0.80, S∗ = 0.91), fair for HGD (κ = 0.28, S∗ = 0.69; and κ = 0.40, S∗ = 0.63), and moderate for NED (κ = 0.47, S∗ = 0.50; and κ = 0.47, S∗ = 0.53). Agreement improved from fair to moderate for LGD (κ = 0.36, S∗ = 0.61; and κ = 0.49, S∗ = 0.71) and slight to fair for IND (κ = 0.16, S∗ = 0.51; and κ = 0.33, S∗ = 0.50). Compared with Hollande's fixed specimens, the agreement was higher in almost all diagnostic categories in formalin-fixed biopsies. Overall, interobserver concordance was improved after a consensus/training session. Interobserver reproducibility was high at the end of the diagnostic spectrum (CA) but fair to moderate for other diagnostic categories.

Genome and transcriptome profiling of FBXW family in human prostate cancer.

F-box and WD repeat domain containing (FBXW) family of E3 ligases has 10 members that ubiquitinate substrate proteins for proteasome-mediated degradation. Publicly archived datasets from The Cancer Genome Atlas (TCGA), Prostate Cancer Transcriptome Atlas (PCTA), and cBioPortal were analyzed for mRNA expression and genetic alterations of 10 FBXW genes. We found that FBXW7 mRNA expression was significantly decreased in primary prostate cancers compared to normal prostate tissues, whereas mRNA expression of FBXW8-10 was significantly increased in primary prostate cancers compared to normal prostate tissues. FBXW7 mRNA expression was also significantly decreased in breast invasive carcinoma, glioblastoma multiforme, head and neck squamous cell carcinoma, lung squamous cell carcinoma, and uterine corpus endometrial carcinoma. In contrast, FBXW7 mRNA expression was significantly increased in cholangiocarcinoma, colon adenocarcinoma, kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, liver hepatocellular carcinoma, lung adenocarcinoma, pheochromocytoma and paraganglioma, and thyroid carcinoma. Compared to normal tissues, FBXW5 mRNA expression was significantly increased in breast invasive carcinoma, cholangiocarcinoma, kidney chromophobe, kidney renal clear cell carcinoma, liver hepatocellular carcinoma, lung adenocarcinoma, lung squamous cell carcinoma, prostate adenocarcinoma, thyroid carcinoma, and uterine corpus endometrial carcinoma, whereas FBXW5 mRNA expression was only significantly decreased in colon adenocarcinoma. There were not any significant differences in gene copy number gains, losses, or gene simple somatic mutations between primary prostate cancers and normal prostate tissues. The mRNA expression levels of FBXW5, 7, 8, 9, and 12 were significantly higher in metastatic castration-resistant prostate cancers (mCRPCs) than primary prostate cancers, whereas mRNA expression levels of FBXW1 and 4 were significantly lower in mCRPCs than primary prostate cancers. All 10 FBXW genes had significantly more overall gene alterations including gene amplifications in mCRPCs than primary prostate cancers. FBXW5 and 7 had significantly more gene deep deletions in mCRPCs than primary prostate cancers and FBXW7 had significantly more gene missense mutations in mCRPCs than primary prostate cancers. Our findings suggest that different FBXW genes have differential mRNA expression in prostate cancer and other cancer types and their gene amplifications are significantly more in mCRPCs than primary prostate cancers. FBXW7 mRNA expression is consistently decreased in primary prostate cancers compared to normal prostate tissues.

The Impact of the Deepwater Horizon Oil Spill upon Lung Health-Mouse Model-Based RNA-Seq Analyses.

We used a transcriptomic approach to interrogate the effects of a saline-accommodated fraction from the Macondo 252 well (MC252) oil and Corexit dispersants on lung tissue. Wild-type C57BL/6 male and female mice were exposed on days 0, 7 and 13 by oropharyngeal aspiration to saline accommodated fractions (SAF) of crude oil from the Macondo (MC252) well, Corexit 9500, Corexit 9527, 9500+oil and 9527+oil or a saline solution as the vehicle control. These treatments did not cause overt toxicity, with the exception of the Corexit exposures which caused brief weight loss after the first exposure. On day 14, total RNA was isolated from the left lung for RNA-seq analyses. KEGG-pathway-based differential expression revealed that Corexit 9527 elicited the strongest changes involving the upregulation of 19 KEGG pathways (FDR < 0.10), followed by Corexit 9500 with the upregulation of seven pathways (FDR < 0.10). As an important signature, pathways related to a response to DNA damage (e.g., p53 signaling and mismatch repair) dominate those upregulated by Corexit 9527 and Corexit 9500. In addition, pro-inflammatory pathways (e.g., cytokine-cytokine receptor interaction, IL-17 signaling pathway and TNF signaling pathways) were upregulated selectively in oil-treated male mice. Surprisingly, oil + dispersant combinations caused lesser effects than the individual treatments at the transcriptomic level. Overall, these findings support potential genotoxicity, inflammation and cell death due to dispersant or oil exposures. Similar exposures to lung tumor bearing K-RasLA1 mice provided evidence for tumor promotion by oil and Corexit dispersant treatments. Our mouse RNA-seq analyses may be relevant to the pulmonary health hazards of MC252 oil and dispersants experienced in exposed populations.

Rheumatoid arthritis-associated DNA methylation sites in peripheral blood mononuclear cells.

OBJECTIVES: To identify novel DNA methylation sites significant for rheumatoid arthritis (RA) and comprehensively understand their underlying pathological mechanism. METHODS: We performed (1) genome-wide DNA methylation and mRNA expression profiling in peripheral blood mononuclear cells from RA patients and health controls; (2) correlation analysis and causal inference tests for DNA methylation and mRNA expression data; (3) differential methylation genes regulatory network construction; (4) validation tests of 10 differential methylation positions (DMPs) of interest and corresponding gene expressions; (5) correlation between PARP9 methylation and its mRNA expression level in Jurkat cells and T cells from patients with RA; (6) testing the pathological functions of PARP9 in Jurkat cells. RESULTS: A total of 1046 DNA methylation positions were associated with RA. The identified DMPs have regulatory effects on mRNA expressions. Causal inference tests identified six DNA methylation-mRNA-RA regulatory chains (eg, cg00959259-PARP9-RA). The identified DMPs and genes formed an interferon-inducible gene interaction network (eg, MX1, IFI44L, DTX3L and PARP9). Key DMPs and corresponding genes were validated their differences in additional samples. Methylation of PARP9 was correlated with mRNA level in Jurkat cells and T lymphocytes isolated from patients with RA. The PARP9 gene exerted significant effects on Jurkat cells (eg, cell cycle, cell proliferation, cell activation and expression of inflammatory factor IL-2). CONCLUSIONS: This multistage study identified an interferon-inducible gene interaction network associated with RA and highlighted the importance of PARP9 gene in RA pathogenesis. The results enhanced our understanding of the important role of DNA methylation in pathology of RA.

Proportions of resting memory T cells and monocytes in blood have prognostic significance in idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a chronic lung disease characterized by progressive decline of lung function. Here, we tested the importance of differential proportions of blood immune cells to IPF clinical outcomes. We used Cibersort to deconvolute immune cell components based on PBMCs or whole blood IPF genomics datasets. We found that a higher proportion of resting memory (RM) T cells was associated with a better survival and a higher DLco (diffusing capacity for carbon monoxide) in IPF patients. The association was also found in opposite direction for monocytes. Additionally, in IPF patients as compared to healthy controls, proportions of monocytes were observed to be higher, yet RM T cells were observed to be lower. Taken together, our result suggests a beneficial effect of RM T cells and a detrimental effect of monocytes for IPF. Future genomics studies of IPF should be more focused on these two types of cells.

Interaction of Epstein-Barr virus genes with human gastric carcinoma transcriptome.

Gastric carcinoma (GC) is a leading cause of mortality. 10% of GC cases are related with EBV (Epstein-Barr virus) infection. The detailed mechanistic roles EBV genes play and especially the interaction between the viral genes and human genes in GC remain unclear. In this study, raw fastq data from 285 GC samples were downloaded from TCGA (The Cancer Genome Atlas), including 25 EBV positive (EBV+) GC samples and 260 EBV negative (EBV-) GC samples. RNA-seq based expression data were generated for both human genes (among all the samples) and for the EBV genes (among the 25 EBV+ samples). Bioinformatics analyses were performed to identify differentially expressed (DEx) human genes and DEx KEGG pathways in EBV+ vs. EBV- samples and co-expressed human gene modules and hub genes among the DEx genes. Within the EBV+ samples, analyses were conducted to find correlation between EBV gene expression and the human gene expression modules, between EBV gene expression and the human hub genes, and between EBV gene expression and the DEx human pathways. EBV genes LMP-1, BALF1 and BALF2 were found to have significant correlation with human hub genes, CNTD2 and VANGL2. EBV genes BALF4 and BALF5 were found to correlate with human pathways, including Jak-STAT signaling and Phosphatidylinositol Signaling System. Our study has revealed the coordinated expression patterns between EBV and human GC transcriptome and identified several key EBV genes that may play an important role in EBV+ GC pathogenesis through their interactions with human genes and pathways.

Carcinogenic effects of oil dispersants: A KEGG pathway-based RNA-seq study of human airway epithelial cells.

The health impacts of the BP oil spill are yet to be further revealed as the toxicological effects of oil products and dispersants on human respiratory system may be latent and complex, and hence difficult to study and follow up. Here we performed RNA-seq analyses of a system of human airway epithelial cells treated with the BP crude oil and/or dispersants Corexit 9500 and Corexit 9527 that were used to help break up the oil spill. Based on the RNA-seq data, we then systemically analyzed the transcriptomic perturbations of the cells at the KEGG pathway level using two pathway-based analysis tools, GAGE (generally applicable gene set enrichment) and GSNCA (Gene Sets Net Correlations Analysis). Our results suggested a pattern of change towards carcinogenesis for the treated cells marked by upregulation of ribosomal biosynthesis (hsa03008) (p=1.97E-13), protein processing (hsa04141) (p=4.09E-7), Wnt signaling (hsa04310) (p=6.76E-3), neurotrophin signaling (hsa04722) (p=7.73E-3) and insulin signaling (hsa04910) (p=1.16E-2) pathways under the dispersant Corexit 9527 treatment, as identified by GAGE analysis. Furthermore, through GSNCA analysis, we identified gene co-expression changes for several KEGG cancer pathways, including small cell lung cancer pathway (hsa05222, p=9.99E-5), under various treatments of oil/dispersant, especially the mixture of oil and Corexit 9527. Overall, our results suggested carcinogenic effects of dispersants (in particular Corexit 9527) and their mixtures with the BP crude oil, and provided further support for more stringent safety precautions and regulations for operations involving long-term respiratory exposure to oil and dispersants.

RNA-sequencing study of peripheral blood monocytes in chronic periodontitis.

BACKGROUND: Monocytes are an important cell type in chronic periodontitis (CP) by interacting with oral bacteria and mediating host immune response. The aim of this study was to reveal new functional genes and pathways for CP at monocyte transcriptomic level. METHODS: We performed an RNA-sequencing (RNA-seq) study of peripheral blood monocytes (PBMs) in 5 non-smoking moderate to severe CP (case) individuals vs. 5 controls. We took advantage of a microarray study of periodontitis to support our findings. We also performed pathway-based analysis on the identified differentially expressed (DEx) transcripts/isoforms using DAVID (Database for Annotation, Visualization and Integrated Discovery). RESULTS: Through differential expression analyses at both whole gene (or whole non-coding RNA) and isoform levels, we identified 380 DEx transcripts and 5955 DEx isoforms with a PPEE (posterior probability of equal expression) of <0.05. Pervasive up-regulation of transcripts at isoform level in CP vs. control individuals was observed, suggesting a more functionally active monocyte transcriptome for CP. By comparing with the microarray dataset, we identified several CP-associated novel genes (e.g., FACR and CUX1) that have functions to interact with invading microorganisms or enhance TNF production on lipopolysaccharide stimulation. DAVID analysis of both the RNA-seq and the microarray datasets leads to converging evidence supporting "endocytosis", "cytokine production" and "apoptosis" as significant biological processes in CP. CONCLUSIONS: As the first RNA-seq study of PBMs for CP, this study provided novel findings at both gene (e.g., FCAR and CUX1) and biological process level. The findings will contribute to better understanding of CP disease mechanisms.

The impact of oil spill to lung health--Insights from an RNA-seq study of human airway epithelial cells.

The Deepwater Horizon oil spill (BP oil spill) in the Gulf of Mexico was a unique disaster event, where a huge amount of oil spilled from the sea bed and a large volume of dispersants were applied to clean the spill. The operation lasted for almost 3 months and involved >50,000 workers. The potential health hazards to these workers may be significant as previous research suggested an association of persistent respiratory symptoms with exposure to oil and oil dispersants. To reveal the potential effects of oil and oil dispersants on the respiratory system at the molecular level, we evaluated the transcriptomic profile of human airway epithelial cells grown under treatment of crude oil, the dispersants Corexit 9500 and Corexit 9527, and oil-dispersant mixtures. We identified a very strong effect of Corexit 9500 treatment, with 84 genes (response genes) differentially expressed in treatment vs. control samples. We discovered an interactive effect of oil-dispersant mixtures; while no response gene was found for Corexit 9527 treatment alone, cells treated with Corexit 9527+oil mixture showed an increased number of response genes (46 response genes), suggesting a synergic effect of 9527 with oil on airway epithelial cells. Through GO (gene ontology) functional term and pathway-based analysis, we identified upregulation of gene sets involved in angiogenesis and immune responses and downregulation of gene sets involved in cell junctions and steroid synthesis as the prevailing transcriptomic signatures in the cells treated with Corexit 9500, oil, or Corexit 9500+oil mixture. Interestingly, these key molecular signatures coincide with important pathological features observed in common lung diseases, such as asthma, cystic fibrosis and chronic obstructive pulmonary disease. Our study provides mechanistic insights into the detrimental effects of oil and oil dispersants to the respiratory system and suggests significant health impacts of the recent BP oil spill to those people involved in the cleaning operation.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density.

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10(-8)) level: 14q24.2 (rs227425, P-value 3.98 × 10(-13), SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10(-9), CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n = 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis.

Pathway-based association analyses identified TRAIL pathway for osteoporotic fractures.

INTRODUCTION: Hip OF carries the highest morbidity and mortality. Previous studies revealed that individual genes/loci in the Tumor Necrosis Factor (TNF)-Related Apoptosis-Inducing Ligand (TRAIL) pathway were associated with bone metabolism. This study aims to verify the potential association between hip OF and TRAIL pathway. METHODS: Using genome-wide genotype data from Affymetrix 500 K SNP arrays, we performed novel pathway-based association analyses for hip OF in 700 elderly Chinese Han subjects (350 with hip OF and 350 healthy matched controls). RESULTS: The TRAIL pathway achieved a significant p value (p = 0.01) for association with hip OF. Among the 38 genes in the TRAIL pathway, seven genes achieved nominally significant association with hip OF (p<0.05); the TNFSF10 (TRAIL) gene obtained the most significant p value (p = 1.70×10(-4)). SNPs (rs719126, rs6533015, rs9594738, rs1805034, rs11160706) from five genes (CFLAR, NFKB1, TNFSF11, TNFRSF11A, TRAF3) of the pathway had minor alleles that appear to be protective to hip OF. SNPs (rs6445063 and rs4259415) from two genes (TNFSF10 and TNFRSF10B) of the pathway had minor alleles (A) that are associated with an increased risk of hip OF, with the ORs (odds ratios) of 16.51 (95%CI:3.83-71.24) and 1.37 (95%CI:1.08-1.74), respectively. CONCLUSIONS: Our study supports the potential role of the TRAIL pathway in the pathogenesis of hip OF in Chinese Han population. Further functional study of this pathway will be pursued to determine the mechanism by which it confers risk to hip OF.

Impact of female cigarette smoking on circulating B cells in vivo: the suppressed ICOSLG, TCF3, and VCAM1 gene functional network may inhibit normal cell function.

As pivotal immune guardians, B cells were found to be directly associated with the onset and development of many smoking-induced diseases. However, the in vivo molecular response of B cells underlying the female cigarette smoking remains unknown. Using the genome-wide Affymetrix HG-133A GeneChip microarray, we firstly compared the gene expression profiles of peripheral circulating B cells between 39 smoking and 40 non-smoking healthy US white women. A total of 125 differential expressed genes were identified in our study, and 75.2% of them were down-regulated in smokers. We further obtained genotypes of 702 single nucleotide polymorphisms in those promising genes and assessed their associations with smoking status. Using a novel multicriteria evaluation model integrating information from microarray and the association studies, several genes were further revealed to play important roles in the response of smoking, including ICOSLG (CD275, inducible T-cell co-stimulator ligand), TCF3 (E2A immunoglobulin enhancer binding factors E12/E47), VCAM1 (CD106, vascular cell adhesion molecule 1), CCR1 (CD191, chemokine C-C motif receptor 1) and IL13 (interleukin 13). The differential expression of ICOSLG (p = 0.0130) and TCF3 (p = 0.0125) genes between the two groups were confirmed by real-time reverse transcription PCR experiment. Our findings support the functional importance of the identified genes in response to the smoking stimulus. This is the first in vivo genome-wide expression study on B cells at today's context of high prevalence rate of smoking for women. Our results highlight the potential usage of integrated analyses for unveiling the novel pathogenesis mechanism and emphasized the significance of B cells in the etiology of smoking-induced disease.