Evaluating Reproducibility in Computational Biology Research

For my Honors Senior Project, I read five research papers in the field of computational biology and attempted to reproduce the results. However, for the most part, this proved a challenge, as many details vital to utilizing relevant software and data had been excluded. Using Geir Kjetil Sandve\u27s paper Ten Simple Rules for Reproducible Computational Research as a guide, I discuss how authors of these five papers did and did not obey these rules of reproducibility and how this affected my ability to reproduce their results

Author Correction:GWAS of thyroid stimulating hormone highlights the pleiotropic effects and inverse association with thyroid cancer

The original version of this article contained an error in the results, in the second paragraph of the subsection entitled “Fine-mapping for potentially causal variants among TSH loci”, in which effect sizes for two variants were incorrectly reported

GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer

Correction: Volume12, Issue1 Article Number7354 DOI10.1038/s41467-021-27675-w PublishedDEC 16 2021Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors. Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. Here, the authors conduct a GWAS and suggest protective effect of higher TSH on risk of thyroid cancer and goitre.Peer reviewe

GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors

MHC Class II Expression Influences the Composition and Distribution of Immune Cells in the Metastatic Colorectal Cancer Microenvironment

Despite advances in therapy over the past decades, metastatic colorectal cancer (mCRC) remains a highly morbid disease. While the impact of MHC-I on immune infiltration in mCRC has been well studied, data on the consequences of MHC-II loss are lacking. Multiplex fluorescent immunohistochemistry (mfIHC) was performed on 149 patients undergoing curative intent resection for mCRC and stratified into high and low human leukocyte antigen isotype DR (HLA-DR) expressing tumors. Intratumoral HLA-DR expression was found in stromal bands, and its expression level was associated with different infiltrating immune cell makeup and distribution. Low HLA-DR expression was associated with increased intercellular distances and decreased population mixing of T helper cells and antigen-presenting cells (APC), suggestive of decreased interactions. This was associated with less co-localization of tumor cells and cytotoxic T lymphocytes (CTLs), which tended to be in a less activated state as determined by Ki67 and granzyme B expression. These findings suggest that low HLA-DR in the tumor microenvironment of mCRC may reflect a state of poor helper T-cell interactions with APCs and CTL-mediated anti-tumor activity. Efforts to restore/enhance MHC-II presentation may be a useful strategy to enhance checkpoint inhibition therapy in the future

Heterogeneity in intrahepatic macrophage populations and druggable target expression in patients with steatotic liver disease-related fibrosis

Background & Aims: Clinical trials for reducing fibrosis in steatotic liver disease (SLD) have targeted macrophages with variable results. We evaluated intrahepatic macrophages in patients with SLD to determine if activity scores or fibrosis stages influenced phenotypes and expression of druggable targets, such as CCR2 and galectin-3. Methods: Liver biopsies from controls or patients with minimal or advanced fibrosis were subject to gene expression analysis using nCounter to determine differences in macrophage-related genes (n = 30). To investigate variability among individual patients, we compared additional biopsies by staining them with multiplex antibody panels (CD68/CD14/CD16/CD163/Mac387 or CD163/CCR2/galectin-3/Mac387) followed by spectral imaging and spatial analysis. Algorithms that utilize deep learning/artificial intelligence were applied to create cell cluster plots, phenotype profile maps, and to determine levels of protein expression (n = 34). Results: Several genes known to be pro-fibrotic (e.g. CD206, TREM2, CD163, and ARG1) showed either no significant differences or significantly decreased with advanced fibrosis. Although marked variability in gene expression was observed in individual patients with cirrhosis, several druggable targets and their ligands (e.g. CCR2, CCR5, CCL2, CCL5, and LGALS3) were significantly increased when compared to patients with minimal fibrosis. Antibody panels identified populations that were significantly increased (e.g. Mac387+), decreased (e.g. CD14+), or enriched (e.g. interactions of Mac387) in patients that had progression of disease or advanced fibrosis. Despite heterogeneity in patients with SLD, several macrophage phenotypes and druggable targets showed a positive correlation with increasing NAFLD activity scores and fibrosis stages. Conclusions: Patients with SLD have markedly varied macrophage- and druggable target-related gene and protein expression in their livers. Several patients had relatively high expression, while others were like controls. Overall, patients with more advanced disease had significantly higher expression of CCR2 and galectin-3 at both the gene and protein levels. Impact and implications: Appreciating individual differences within the hepatic microenvironment of patients with SLD may be paramount to developing effective treatments. These results may explain why such a small percentage of patients have responded to macrophage-targeting therapies and provide additional support for precision medicine-guided treatment of chronic liver diseases

Subclonal evolution and expansion of spatially distinct THY1-positive cells is associated with recurrence in glioblastoma

Purpose: Glioblastoma(GBM) is a lethal disease characterized by inevitable recurrence. Here we investigate the molecular pathways mediating resistance, with the goal of identifying novel therapeutic opportunities. Experimental design: We developed a longitudinal in vivo recurrence model utilizing patient-derived explants to produce paired specimens(pre- and post-recurrence) following temozolomide(TMZ) and radiation(IR). These specimens were evaluated for treatment response and to identify gene expression pathways driving treatment resistance. Findings were clinically validated using spatial transcriptomics of human GBMs. Results: These studies reveal in replicate cohorts, a gene expression profile characterized by upregulation of mesenchymal and stem-like genes at recurrence. Analyses of clinical databases revealed significant association of this transcriptional profile with worse overall survival and upregulation at recurrence. Notably, gene expression analyses identified upregulation of TGFβ signaling, and more than one-hundred-fold increase in THY1 levels at recurrence. Furthermore, THY1-positive cells represented <10% of cells in treatment-naïve tumors, compared to 75-96% in recurrent tumors. We then isolated THY1-positive cells from treatment-naïve patient samples and determined that they were inherently resistant to chemoradiation in orthotopic models. Additionally, using image-guided biopsies from treatment-naïve human GBM, we conducted spatial transcriptomic analyses. This revealed rare THY1+ regions characterized by mesenchymal/stem-like gene expression, analogous to our recurrent mouse model, which co-localized with macrophages within the perivascular niche. We then inhibited TGFBRI activity in vivo which decreased mesenchymal/stem-like protein levels, including THY1, and restored sensitivity to TMZ/IR in recurrent tumors. Conclusions: These findings reveal that GBM recurrence may result from tumor repopulation by pre-existing, therapy-resistant, THY1-positive, mesenchymal cells within the perivascular niche

Author Correction: GWAS of thyroid stimulating hormone highlights the pleiotropic effects and inverse association with thyroid cancer

To access publisher's full text version of this article click on the hyperlink belo