Identification of immune-related hub genes and construction of an immune-related gene prognostic index for low-grade glioma

Low-grade glioma (LGG) poses significant management challenges and has a dismal prognosis. While immunotherapy has shown significant promise in cancer treatment, its progress in glioma has confronted with challenges. In our study, we aimed to develop an immune-related gene prognostic index (IRGPI) which could be used to evaluate the response and efficacy of LGG patients with immunotherapy. We included a total of 529 LGG samples from TCGA database and 1152 normal brain tissue samples from the GTEx database. Immune-related differentially expressed genes (DEGs) were screened. Then, we used weighted gene co-expression network analysis (WGCNA) to identify immune-related hub genes in LGG patients and performed Cox regression analysis to construct an IRGPI. The median IRGPI was used as the cut-off value to categorize LGG patients into IRGPI-high and low subgroups, and the molecular and immune mechanism in IRGPI-defined subgroups were analysed. Finally, we explored the relationship between IRGPI-defined subgroups and immunotherapy related indicators in patients after immunotherapy. Three genes (RHOA, NFKBIA and CCL3) were selected to construct the IRGPI. In a survival analysis using TCGA cohort as a training set, patients in the IRGPI-low subgroup had a better OS than those in IRGPI-high subgroup, consistent with the results in CGGA cohort. The comprehensive results showed that IRGPI-low subgroup had a more abundant activated immune cell population and lower TIDE score, higher MSI, higher TMB score, lower T cell dysfunction score, more likely benefit from ICIs therapy. IRGPI is a promising biomarker in the field of LGG ICIs therapy to distinguish the prognosis, the molecular and immunological characteristics of patients.     

Identification of immune microenvironment subtypes that predicted the prognosis of patients with ovarian cancer

Ovarian cancer (OC) is associated with high mortality rate. However, the correlation between immune microenvironment and prognosis of OC remains unclear. This study aimed to explore prognostic significance of OC tumour microenvironment. The OC data set was selected from the cancer genome atlas (TCGA), and 307 samples were collected. Hierarchical clustering was performed according to the expression of 756 genes. The immune and matrix scores of all immune subtypes were determined, and Kruskal-Wallis test was used to analyse the differences in the immune and matrix scores between OC samples with different immune subtypes. The model for predicting prognosis was constructed based on the expression of immune-related genes. TIDE platform was applied to predict the effect of immunotherapy on patients with OC of different immune subtypes. The 307 OC samples were classified into three immune subtypes A-C. Patients in subtype B had poorer prognosis and lower survival rate. The infiltration of helper T cells and macrophages in microenvironment indicated significant differences between immune subtypes. Enrichment analyses of immune cell molecular pathways showed that JAK–STAT3 pathway changed significantly in subtype B. Furthermore, predictive response to immunotherapy in subtype B was significantly higher than that in subtype A and C. Immune subtyping can be used as an independent predictor of the prognosis of OC patients, which may be related to the infiltration patterns of immune cells in tumour microenvironment. In addition, patients in immune subtype B have superior response to immunotherapy, suggesting that patients in subtype B are suitable for immunotherapy.     

Integrated Multiple “-omics” Data Reveal Subtypes of Hepatocellular Carcinoma

Hepatocellular carcinoma is one of the most heterogeneous cancers, as reflected by its multiple grades and difficulty to subtype. In this study, we integrated copy number variation, DNA methylation, mRNA, and miRNA data with the developed “cluster of cluster” method and classified 256 HCC samples from TCGA (The Cancer Genome Atlas) into five major subgroups (S1-S5). We observed that this classification was associated with specific mutations and protein expression, and we detected that each subgroup had distinct molecular signatures. The subclasses were associated not only with survival but also with clinical observations. S1 was characterized by bulk amplification on 8q24, TP53 mutation, low lipid metabolism, highly expressed onco-proteins, attenuated tumor suppressor proteins and a worse survival rate. S2 and S3 were characterized by telomere hypomethylation and a low expression of TERT and DNMT1/3B. Compared to S2, S3 was associated with less copy number variation and some good prognosis biomarkers, including CRP and CYP2E1. In contrast, the mutation rate of CTNNB1 was higher in S3. S4 was associated with bulk amplification and various molecular characteristics at different biological levels. In summary, we classified the HCC samples into five subgroups using multiple “-omics” data. Each subgroup had a distinct survival rate and molecular signature, which may provide information about the pathogenesis of subtypes in HCC.     

Molecular subtype classification is a determinant of non-sentinel lymph node metastasis in breast cancer patients with positive sentinel lymph nodes

Background

Previous studies suggested that the molecular subtypes were strongly associated with sentinel lymph node (SLN) status. The purpose of this study was to determine whether molecular subtype classification was associated with non-sentinel lymph nodes (NSLN) metastasis in patients with a positive SLN.

Methodology and Principal Findings

Between January 2001 and March 2011, a total of 130 patients with a positive SLN were recruited. All these patients underwent a complete axillary lymph node dissection. The univariate and multivariate analyses of NSLN metastasis were performed. In univariate and multivariate analyses, large tumor size, macrometastasis and high tumor grade were all significant risk factors of NSLN metastasis in patients with a positive SLN. In univariate analysis, luminal B subgroup showed higher rate of NSLN metastasis than other subgroup (P = 0.027). When other variables were adjusted in multivariate analysis, the molecular subtype classification was a determinant of NSLN metastasis. Relative to triple negative subgroup, both luminal A (P = 0.047) and luminal B (P = 0.010) subgroups showed a higher risk of NSLN metastasis. Otherwise, HER2 over-expression subgroup did not have a higher risk than triple negative subgroup (P = 0.183). The area under the curve (AUC) value was 0.8095 for the Cambridge model. When molecular subtype classification was added to the Cambridge model, the AUC value was 0.8475.

Conclusions

Except for other factors, molecular subtype classification was a determinant of NSLN metastasis in patients with a positive SLN. The predictive accuracy of mathematical models including molecular subtype should be determined in the future.     

Gene Expression Profiling of Ampullary Carcinomas Classifies Ampullary Carcinomas into Biliary-Like and Intestinal-Like Subtypes That Are Prognostic of Outcome

Background

Adenocarcinomas of the ampulla of Vater are classified as biliary cancers, though the exact epithelium of origin for these cancers is not known. We sought to molecularly classify ampullary adenocarcinomas in comparison to known adenocarcinomas of the pancreas, bile duct, and duodenum by gene expression analysis.

Methods

We analyzed 32 fresh-frozen resected, untreated periampullary adenocarcinomas (8 pancreatic, 2 extrahepatic biliary, 8 duodenal, and 14 ampullary) using the Affymetrix U133 Plus 2.0 genome array. Unsupervised and supervised hierarchical clustering identified two subtypes of ampullary carcinomas that were molecularly and histologically characterized.

Results

Hierarchical clustering of periampullary carcinomas segregated ampullary carcinomas into two subgroups, which were distinctly different from pancreatic carcinomas. Non-pancreatic periampullary adenocarcinomas were segregated into two subgroups with differing prognoses: 5 year RFS (77% vs. 0%, p = 0.007) and 5 year OS (100% vs. 35%, p = 0.005). Unsupervised clustering analysis of the 14 ampullary samples also identified two subgroups: a good prognosis intestinal-like subgroup and a poor prognosis biliary-like subgroup with 5 year OS of 70% vs. 28%, P = 0.09. Expression of CK7+/CK20- but not CDX-2 correlated with these two subgroups. Activation of the AKT and MAPK pathways were both increased in the poor prognostic biliary-like subgroup. In an independent 80 patient ampullary validation dataset only histological subtype (intestinal vs. pancreaticobiliary) was significantly associated with OS in both univariate (p = 0.006) and multivariate analysis (P = 0.04).

Conclusions

Gene expression analysis discriminated pancreatic adenocarcinomas from other periampullary carcinomas and identified two prognostically relevant subgroups of ampullary adenocarcinomas. Histological subtype was an independent prognostic factor in ampullary adenocarcinomas.     

CTLA4 Variants and Haplotype Contribute Genetic Susceptibility to Myasthenia Gravis in Northern Chinese Population

Background

Cytotoxic T lymphocyte-associated antigen-4 (CTLA4), a critical negative regulator of the T-cell response, has been considered a candidate for many autoimmune diseases. Evidence from Caucasians supported a genetic predisposition of CTLA4 to myasthenia gravis (MG), but the contribution in East Asians has not been established.

Objectives

To investigate the role of CTLA4 variants in the susceptibility to MG and the contribution to subtypes of MG.

Methods

Six autoimmune disease-related risk alleles of CTLA4 (rs1863800, rs733618, rs4553808, rs5742909, rs231775, and rs3087243) were investigated for MG in northern Chinese. 168 patients with MG (mean age 37.1±20.5 years, 64 men and 104 women) and 233 healthy controls (mean age 53.3±8.7 years, 96 men and 137 women) were screened, and the contribution of CTLA4 to the general risk of MG and each subgroup was explored.

Results

rs1863800*C, rs733618*C, and rs231775*G were significantly associated with the whole cohort of patients with MG after permutation correction for multiple-testing adjustment (P = 0.027, 0.001, and 0.032, respectively). A risk haplotype (CCACG) [odds ratio (OR) = 1.535, range = 1.150–2.059, P = 0.004)] was also identified. The stratified subtype analysis indicated that the positive contribution was possibly derived from early onset MG (EOMG), seropositive MG (SPMG), female patients, and MG without thymoma. No association was observed in juvenile MG/LOMG, and MG coupled with thymoma.

Conclusion

A predisposing effect of rs1863800*C, rs733618*C, and rs231775*G of CTLA4 gene to general risk of MG in Chinese was demonstrated for the first time, which was likely derived from EOMG, SPMG, MG without thymoma and the female patients.     

Quantitative proteomic analysis of A549 cells infected with human respiratory syncytial virus subgroup B using SILAC coupled to LC-MS/MS

Human respiratory syncytial virus (HRSV) is a leading cause of serious lower respiratory tract infections in infants. The virus has two subgroups A and B, which differ in prevalence and (nucleotide) sequence. The interaction of subgroup A viruses with the host cell is relatively well characterized, whereas for subgroup B viruses it is not. Therefore quantitative proteomics was used to investigate the interaction of subgroup B viruses with A549 cells, a respiratory cell line. Changes in the cellular proteome and potential canonical pathways were determined using SILAC coupled to LC-MS/MS and Ingenuity Pathway Analysis. To reduce sample complexity and investigate potential trafficking both nuclear and cytoplasmic fractions were analyzed. A total of 904 cellular and six viral proteins were identified and quantified, of which 112 cellular proteins showed a twofold or more change in HRSV-infected cells. Data sets were validated using indirect immunofluorescence confocal microscopy on independent samples. Major changes were observed in constituents of mitochondria including components of the electron transport chain complexes and channels, as well as increases in the abundance of the products of interferon-stimulated genes. This is the first quantitative proteomic analysis of cells infected with HRSV-subgroup B.     

The Association of PTPN22 R620W Polymorphism Is Stronger with Late-Onset AChR-Myasthenia Gravis in Turkey

A functional single nucleotide polymorphism (SNP) of the PTPN22 gene encoding a protein tyrosine phosphatase has been associated with autoimmune disorders including myasthenia gravis (MG). As the PTPN22 R620W polymorphism has a wide variation of allele frequencies among different populations, this polymorphism was investigated in MG in Turkey. An emphasis is put on MG subgroups according to autoantibody (Abs) production and presence of thymoma. DNA samples from 416 patients with clinically diagnosed generalized MG (231 with Abs to acetylcholine receptor, AChR-MG), 53 with Abs to muscle-specific kinase (MuSK-MG), 55 patients with no detectable Abs (SN-MG), 77 patients with thymoma (TAMG) and 293 healthy controls (HC) were genotyped for the SNP (PTPN22 R620W, C1858T, rs2476601). The PTPN22 T allele was increased in AChR-MG patients (odds ratio [OR]: 2.5, 95%CI: 1.2–5.1). The association was stronger in late disease-onset AChR (LOMG, OR: 3.1, 95%CI: 1.2–8.2). MuSK-MG, SN-MG and TAMG groups did not carry the variant allele more frequently than the HC. In contrast to findings in other autoimmune diseases, the distribution of the PTPN22 polymorphism in this population provides a susceptibility marker for AChR-MG. The strongest association is detected in patients with LOMG.     

生物Fas和Bcl-2在胸腺瘤患者瘤组织表达中与并发重症肌无力的相关性研究

探讨胸腺瘤患者瘤组织中凋亡诱导基因Fas和凋亡抑制基因Bcl-2的表达情况及其与胸腺瘤患者并发重症肌无力（MG）的相关性。通过免疫组化S-P法,检测经手术治疗的胸腺瘤伴MG患者切除瘤组织中Fas、Bcl-2的表达水平,并以胸腺瘤不伴MG患者瘤组织中Fas、Bcl-2的表达水平作为对照。Fas在伴有MG的胸腺瘤中的表达高于对照组,差别有显著性意义（P〈0.01）;Bcl-2在伴有MG的胸腺瘤中的表达低于对照组,差别有显著性意义（P〈0.01）。胸腺瘤合并MG患者瘤组织中凋亡诱导基因Fas呈高表达和凋亡抑制基因Bcl-2呈低表达与胸腺瘤患者并发MG有相关性。     

Evolutionary Rate and Genetic Heterogeneity of Human T-Cell Lymphotropic Virus Type II (HTLV-II) Using Isolates from European Injecting Drug Users

Seven new Italian and two new British HTLV-II isolates were obtained from injecting drug users and the entire long terminal repeat (LTR) region was sequenced. Restriction analysis showed that all the Italian isolates are of the IIb subtype, whereas the British isolates are of the IIa subtype. To understand whether the further differentiation of each two principal HTLV-II subtypes in several subgroups could be statistically supported by phylogenetic analysis, the neighbor-joining, parsimony, and maximum likelihood methods were used. The separation between IIa and IIb is very well supported by all three methods. At least two phylogenetic subgroups exist within the HTLV-IIa and at least three within the HTLV-IIb subtype. In the present analysis, no statistical support was obtained for additional phylogroups. Two particular subgroups seem interesting because they include all European and North American injecting drug user strains within the IIa and IIb subtypes, respectively. These data confirm that European HTLV-II infection among drug users is probably derived from North America. They also suggest that though a certain differentiation by restriction analysis in different subgroups is possible, carefully interpreted phylogenetic analyses remain necessary. Using the likelihood ratio test, a molecular clock for the drug user strains was calibrated. A fixation rate between 1.08 × 10⁻⁴ and 2.7 × 10⁻⁵ nucleotide substitutions per site per year was calculated for the IIa and IIb injecting drug user strains. This is the lowest fixation rate so far reported for RNA viruses, including for HIV, which typically range between 10⁻² and 10⁻⁴.     

Non-Sentinel Lymph Node Metastasis Prediction in Breast Cancer with Metastatic Sentinel Lymph Node: Impact of Molecular Subtypes Classification

Introduction

To decipher the interaction between the molecular subtype classification and the probability of a non-sentinel node metastasis in breast cancer patients with a metastatic sentinel lymph-node, we applied two validated predictors (Tenon Score and MSKCC Nomogram) on two large independent datasets.

Materials and Methods

Our datasets consisted of 656 and 574 early-stage breast cancer patients with a metastatic sentinel lymph-node biopsy treated at first by surgery. We applied both predictors on the whole dataset and on each molecular immune-phenotype subgroups. The performances of the two predictors were analyzed in terms of discrimination and calibration. Probability of non-sentinel lymph node metastasis was detailed for each molecular subtype.

Results

Similar results were obtained with both predictors. We showed that the performance in terms of discrimination was as expected in ER Positive HER2 negative subgroup in both datasets (MSKCC AUC Dataset 1 = 0.73 [0.69–0.78], MSKCC AUC Dataset 2 = 0.71 (0.65–0.76), Tenon Score AUC Dataset 1 = 0.7 (0.65–0.75), Tenon Score AUC Dataset 2 = 0.72 (0.66–0.76)). Probability of non-sentinel node metastatic involvement was slightly under-estimated. Contradictory results were obtained in other subgroups (ER negative HER2 negative, HER2 positive subgroups) in both datasets probably due to a small sample size issue. We showed that merging the two datasets shifted the performance close to the ER positive HER2 negative subgroup.

Discussion

We showed that validated predictors like the Tenon Score or the MSKCC nomogram built on heterogeneous population of breast cancer performed equally on the different subgroups analyzed. Our present study re-enforce the idea that performing subgroup analysis of such predictors within less than 200 samples subgroup is at major risk of misleading conclusions.     

Serological characterization of human reassortant rotaviruses.

We analyzed the serological properties of two human wild-type cell culture-adapted rotaviruses (strains 308 and 46) and of 308 X 46 reassortants which were previously obtained and genetically characterized. Strain 308, exhibiting a so-called long RNA pattern, was found to belong to human rotavirus subgroup II, serotype 1, whereas strain 46, exhibiting a so-called short RNA pattern, represented subgroup I, serotype 2. Among the 308 X 46 reassortants we analyzed, two belonged to subgroup II, serotype 1, and exhibited short RNA patterns. This showed that the correlation observed between human subgroups I and II rotaviruses and the short and long electrophoretic patterns is not supported by any molecular basis (i.e., gene segment 10 or 11 was not involved in the subgroup specificity).     

果蝇Drosophila saltans种亚组（双翅目：果蝇科）系统发育关系：形态学能否解决分子冲突？

正确的系统发生重建对于理解进化事件至关重要.尽管分子系统学对于解决此类问题取得了极大的成功,由于一些诸如密码子使用偏性等的内在约束,来源于DNA的信息可能仍然存在着局限.因为发生在祖先的替代性转换,果蝇Drosophila saltans 5个种亚组由不同基因构建的分子系统树之间存在着冲突(在以往发表的分子系统学研究中,这些种组的每一个种亚组全少有一个代表).本文用40个形态学特征重新分析了这些种组.不同于以前发表的大多数假说,本研究支序分类学的结果表明,果蝇sturtevanti种亚组是一个较早的分支,而剩下的4个亚组形成一个支持度较高的类群;后者又可以再分为两个姐妹群:一个包含cordata和elliptica亚组,另一个包含parasaltans和saltans亚来组.本研究结果修正了果蝇saltans种组的分子进化(密码子使用偏性),并强调形态学对于系统发生重建和理解分子进化现象的重要作用.     

Detection of Patient Subgroups with Differential Expression in Omics Data: A Comprehensive Comparison of Univariate Measures

Detection of yet unknown subgroups showing differential gene or protein expression is a frequent goal in the analysis of modern molecular data. Applications range from cancer biology over developmental biology to toxicology. Often a control and an experimental group are compared, and subgroups can be characterized by differential expression for only a subgroup-specific set of genes or proteins. Finding such genes and corresponding patient subgroups can help in understanding pathological pathways, diagnosis and defining drug targets. The size of the subgroup and the type of differential expression determine the optimal strategy for subgroup identification. To date, commonly used software packages hardly provide statistical tests and methods for the detection of such subgroups. Different univariate methods for subgroup detection are characterized and compared, both on simulated and on real data. We present an advanced design for simulation studies: Data is simulated under different distributional assumptions for the expression of the subgroup, and performance results are compared against theoretical upper bounds. For each distribution, different degrees of deviation from the majority of observations are considered for the subgroup. We evaluate classical approaches as well as various new suggestions in the context of omics data, including outlier sum, PADGE, and kurtosis. We also propose the new FisherSum score. ROC curve analysis and AUC values are used to quantify the ability of the methods to distinguish between genes or proteins with and without certain subgroup patterns. In general, FisherSum for small subgroups and -test for large subgroups achieve best results. We apply each method to a case-control study on Parkinson''s disease and underline the biological benefit of the new method.     

Molecular genetic testing in the United States: comparison with international practice

OBJECTIVE: To compare data on the practices of molecular genetic testing (MGT) in laboratories in the United States with those in 18 other countries. METHODS: A Web-based survey of MGT laboratory directors (n = 827; response rate 63%) in 18 countries on three continents was carried out, and the response from U.S. laboratories compared to all others. Quality assurance and reporting indices were developed and calculated for each responding laboratory. RESULTS: A comparison of U.S. results with all other countries identified differences in laboratory setting, personnel qualifications, and the specific tests being offered, but similar rates of adherence to MGT quality standards and reporting practices were found. The survey also documented substantial transborder flow of specimens, most commonly due to the lack of availability of the test in the United States or because the test was available only through a research protocol, highlighting the need for common reporting and practice guidelines for the international MGT community. CONCLUSION: The findings presented here provide further support for the need to consider the application of the Organisation for Economic Cooperation and Development (OECD) Guidelines and the establishment of compatible accreditation programs or equivalent mechanisms across national borders to ensure the quality of laboratory services and the clinical usefulness of molecular genetic test reports for referred specimens.     

重症肌无力患者抗乙酰胆碱受体抗体的测定及意义

用ELISA法对56例重症肌无力(MG)患者治疗前后的血清乙酰胆碱受体(AchR)抗体进行了检测。检测结果为MG患者治疗前后的血清AchR抗体阳性率46.4%。而且发病年龄越大,患者体内AchR抗体阳性率越高。患者患病时间越长体内AchR抗体含量降低。伴发胸腺瘤的患者AchR抗体阳性率明显高于胸腺正常者,全身型MG患者AchR抗体检测阳性率高于眼肌型患者。     

Identification of two molecular subtypes of dysregulated immune lncRNAs in ovarian cancer

Long non-coding RNA (lncRNA) has increasingly been identified as a key regulator in pathologies such as cancer. Multiple platforms were used for comprehensive analysis of ovarian cancer to identify molecular subgroups. However, lncRNA and its role in mapping the ovarian cancer subpopulation are still largely unknown. RNA-sequencing and clinical characteristics of ovarian cancer were acquired from The Cancer Genome Atlas database (TCGA). A total of 52 lncRNAs were identified as aberrant immune lncRNAs specific to ovarian cancer. We redefined two different molecular subtypes, C1(188) and C2(184 samples), in “iClusterPlus” R package, among which C2 grouped ovarian cancer samples have higher survival probability and longer median survival time (P <0.05) with activated IFN-gamma response, Wound Healing and Cytotoxic lymphocytes signal; 456 differentially expressed genes were acquired in C1 and C2 subtypes using limma (3.40.6) package, among which 419 were up-regulated and 37 were down-regulated, in TCGA dataset. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis revealed that these genes were actively involved in ECM-receptor interaction, PI3K-Akt signaling pathway interaction KEGG pathway. Compared with the existing immune subtype, the Cluster2 sample showed a substantial increase in the proportion of the existing C2 immune subtype, accounting for 81.37%, which was associated with good prognosis. Our C1 subtype contains only 56.49% of the existing immune C1 and C4, which also explains the poor prognosis of C1. Furthermore, 52 immune-related lncRNAs were used to divide the TCGA-endometrial cancer and cervical cancer samples into two categories, and C2 had a good prognosis. The differentially expressed genes were highly correlated with immune-cell-related pathways. Based on lncRNA, two molecular subtypes of ovarian cancer were identified and had significant prognostic differences and immunological characteristics.     

Expression Profiling of Rectal Tumors Defines Response to Neoadjuvant Treatment Related Genes

To date, no effective method exists that predicts the response to preoperative chemoradiation (CRT) in locally advanced rectal cancer (LARC). Nevertheless, identification of patients who have a higher likelihood of responding to preoperative CRT could be crucial in decreasing treatment morbidity and avoiding expensive and time-consuming treatments. The aim of this study was to identify signatures or molecular markers related to response to pre-operative CRT in LARC. We analyzed the gene expression profiles of 26 pre-treatment biopsies of LARC (10 responders and 16 non-responders) without metastasis using Human WG CodeLink microarray platform. Two hundred and fifty seven genes were differentially over-expressed in the responder patient subgroup. Ingenuity Pathway Analysis revealed a significant ratio of differentially expressed genes related to cancer, cellular growth and proliferation pathways, and c-Myc network. We demonstrated that high Gng4, c-Myc, Pola1, and Rrm1 mRNA expression levels was a significant prognostic factor for response to treatment in LARC patients (p<0.05). Using this gene set, we were able to establish a new model for predicting the response to CRT in rectal cancer with a sensitivity of 60% and 100% specificity. Our results reflect the value of gene expression profiling to gain insight about the molecular pathways involved in the response to treatment of LARC patients. These findings could be clinically relevant and support the use of mRNA levels when aiming to identify patients who respond to CRT therapy.