Two HLA Class II Gene Variants Are Independently Associated with Pediatric Osteosarcoma Risk

Genetic variation near GRB10 associated with bone growth and osteosarcoma risk in canine and human populations

BACKGROUND

Canine and human osteosarcoma are similar in clinical presentation and tumor genomics. Giant breed dogs experience elevated osteosarcoma incidence, and taller stature remains a consistent risk factor for human osteosarcoma. Whether evolutionarily conserved genes contribute to both human and canine osteosarcoma predisposition merits evaluation.

METHODS

A multi-center sample of childhood osteosarcoma patients and controls underwent genome-wide genotyping and imputation. Ancestry-adjusted SNP associations were calculated within each dataset using logistic regression, then meta-analyzed across the three datasets, totaling 1091 patients and 3026 controls. Ten regions previously associated with canine osteosarcoma risk were mapped to the human genome, spanning ∼6 Mb. We prioritized association testing of 5985 human SNPs mapping to candidate osteosarcoma risk regions detected in Irish wolfhounds, the largest dog breed studied. Secondary analyses explored 6289 additional human SNPs mapping to candidate osteosarcoma risk regions identified in Rottweilers and greyhounds.

RESULTS

Fourteen SNPs were associated with human osteosarcoma risk after adjustment for multiple comparisons, all within a 42 kb region of human Chromosome 7p12.1. The lead variant was rs17454681 (OR=1.25, 95 %CI: 1.12-1.39; P=4.1×10-5), and independent risk variants were not observed in conditional analyses. While the associated region spanned 2.1 Mb and contained eight genes in Irish wolfhounds, associations were localized to a 50-fold smaller region of the human genome and strongly implicate GRB10 (growth factor receptor-bound protein 10) in canine and human osteosarcoma predisposition. PheWAS analysis in UK Biobank data identified noteworthy associations of the rs17454681 risk allele with varied measures of height and pubertal timing.

CONCLUSIONS

Our comparative oncology analysis identified a novel human osteosarcoma risk allele near GRB10, a growth inhibitor that suppresses activated receptor tyrosine kinases including IGF1R, PDGFRB, and EGFR. Epidemiologists may benefit from leveraging cross-species comparisons to identify haplotypes in highly susceptible but genetically homogenous populations of domesticated animals, then fine-mapping these associations in diverse human populations.

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers.

Identification of potential crucial genes and key pathways in osteosarcoma

Background

The aim of this study is to identify the potential pathogenic and metastasis-related differentially expressed genes (DEGs) in osteosarcoma through bioinformatic analysis of Gene Expression Omnibus (GEO) database.

Results

Gene expression profiles of GSE14359, GSE16088, and GSE33383, in total 112 osteosarcoma tissue samples and 7 osteoblasts, were analyzed. Seventy-four normal-primary DEGs (NPDEGs) and 764 primary-metastatic DEGs (PMDEGs) were screened. VAMP8, A2M, HLA-DRA, SPARCL1, HLA-DQA1, APOC1 and AQP1 were identified continuously upregulating during the oncogenesis and metastasis of osteosarcoma. The enriched functions and pathways of NPDEGs include procession and presentation of antigens, activation of MHC class II receptors and phagocytosis. The enriched functions and pathways of PMDEGs include mitotic nuclear division, cell adhesion molecules (CAMs) and focal adhesion. With protein-protein interaction (PPI) network analyzed by Molecular Complex Detection (MCODE) plug-in of Cytoscape software, one hub NPDEG (HLA-DRA) and 7 hub PMDEGs (CDK1, CDK20, CCNB1, MTIF2, MRPS7, VEGFA and EGF) were eventually selected, and the most significant pathways in NPDEGs module and PMDEGs module were enriched in the procession and presentation of exogenous peptide antigen via MHC class II and the nuclear division, respectively.

Conclusions

By integrated bioinformatic analysis, numerous DEGs related to osteosarcoma were screened, and the hub DEGs identified in this study are possibly part of the potential biomarkers for osteosarcoma. However, further experimental studies are still necessary to elucidate the biological function and mechanism of these genes.

Common genetic variation and risk of osteosarcoma in a multi-ethnic pediatric and adolescent population

Osteosarcoma, a malignant primary bone tumor most commonly diagnosed in children and adolescents, has a poorly understood genetic etiology. Genome-wide association studies (GWAS) and candidate-gene analyses have identified putative risk variants in subjects of European ancestry. However, despite higher incidence among African-American and Hispanic children, little is known regarding common heritable variation that contributes to osteosarcoma incidence and clinical presentation across racial/ethnic groups. In a multi-ethnic sample of non-Hispanic white, Hispanic, African-American and Asian/Pacific Islander children (537 cases, 2165 controls), we performed association analyses assessing previously-reported loci for osteosarcoma risk and metastasis, including meta-analysis across racial/ethnic groups. We also assessed a previously described association between genetic predisposition to longer leukocyte telomere length (LTL) and osteosarcoma risk in this independent multi-ethnic dataset. In our sample, we were unable to replicate previously-reported loci for osteosarcoma risk or metastasis detected in GWAS of European-ancestry individuals in either ethnicity-stratified analyses or meta-analysis across ethnic groups. Our analyses did confirm that genetic predisposition to longer LTL is a risk factor for osteosarcoma (ORmeta: 1.22; 95% CI: 1.09-1.36; P = 3.8 × 10-4), and the strongest effect was seen in Hispanic subjects (OR: 1.32; 95% CI: 1.12-1.54, P = 6.2 × 10-4). Our findings shed light on the replicability of osteosarcoma risk loci across ethnicities and motivate further characterization of these genetic factors in diverse clinical cohorts.

Genome-wide expression profiling-based copy number variations and colorectal cancer risk in Chinese

Genetic factors play important roles in colorectal carcinogenesis. This study was aimed to evaluate the effects of gene expression-related copy number variations (CNVs) on the risk of colorectal cancer in Chinese. Expression Quantitative Trait Locus (eQTL) mapping was conducted to explore the most regulatable gene expressions by CNVs among the whole genome based on publicly available data. Then a case-control study was performed to evaluate the associations between copy numbers of the most regulatable genes and colorectal cancer. The influence of the target CNVs on the expression of corresponding gene and protein was verified in colorectal tissue, and the biological effects of these CNVs on cell-cycle arrest and apoptosis of colon cancer cell lines were further detected. The eQTL revealed the most significant association between CNV of HM3_CNP_342 and gene expressions of human leukocyte antigen (HLA)-DQA1 and HLA-DQB1 among the whole genome. The later case-control study found that amplified HLA-DQB1 was inversely associated with colorectal cancer risk (odds ratio = 0.73; 95% confidence interval: 0.58-0.93), especially among those with a family history of cancer. The positive association between amplified HLA-DQB1 and upregulation of gene and protein was validated in colorectal tissue. In addition, overexpression of HLA-DQB1 in dendritic cells promoted cell-cycle arrest and apoptosis of cocultured SW480 and HCT116 cell lines, and vice versa. Our study suggests that the amplified copy number of HLA-DQB1 is associated with lower risk of colorectal cancer and able to induce the apoptosis of colon cancer cells, which implies the potential of HLA class II in cancer predisposition and immunotherapy.