Pan-cancer transcriptomic analysis dissects immune and proliferative functions of APOBEC3 cytidine deaminases.

APOBEC3 cytidine deaminases are largely known for their innate immune protection from viral infections. Recently, members of the family have been associated with a distinct mutational activity in some cancer types. We report a pan-tissue, pan-cancer analysis of RNA-seq data specific to the APOBEC3 genes in 8,951 tumours, 786 cancer cell lines and 6,119 normal tissues. By deconvolution of levels of different cell types in tumour admixtures, we demonstrate that APOBEC3B (A3B), the primary candidate as a cancer mutagen, shows little association with immune cell types compared to its paralogues. We present a pipeline called RESPECTEx (REconstituting SPecific Cell-Type Expression) and use it to deconvolute cell-type specific expression levels in a given cohort of tumour samples. We functionally annotate APOBEC3 co-expressing genes, and create an interactive visualization tool which 'barcodes' the functional enrichment (http://fraternalilab.kcl.ac.uk/apobec-barcodes/). These analyses reveal that A3B expression correlates with cell cycle and DNA repair genes, whereas the other APOBEC3 members display specificity for immune processes and immune cell populations. We offer molecular insights into the functions of individual APOBEC3 proteins in antiviral and proliferative contexts, and demonstrate the diversification this family of enzymes displays at the transcriptomic level, despite their high similarity in protein sequences and structures.

Functional and prognostic significance of the genomic amplification of frizzled 6 (FZD6) in breast cancer.

Frizzled receptors mediate Wnt ligand signalling, which is crucially involved in regulating tissue development and differentiation, and is often deregulated in cancer. In this study, we found that the gene encoding the Wnt receptor frizzled 6 (FZD6) is frequently amplified in breast cancer, with an increased incidence in the triple-negative breast cancer (TNBC) subtype. Ablation of FZD6 expression in mammary cancer cell lines: (1) inhibited motility and invasion; (2) induced a more symmetrical shape of organoid three-dimensional cultures; and (3) inhibited bone and liver metastasis in vivo. Mechanistically, FZD6 signalling is required for the assembly of the fibronectin matrix, interfering with the organization of the actin cytoskeleton. Ectopic delivery of fibronectin in FZD6-depleted, triple-negative MDA-MB-231 cells rearranged the actin cytoskeleton and restored epidermal growth factor-mediated invasion. In patients with localized, lymph node-negative (early) breast cancer, positivity of tumour cells for FZD6 protein identified patients with reduced distant relapse-free survival. Multivariate analysis indicated an independent prognostic significance of FZD6 expression in TNBC tumours, predicting distant, but not local, relapse. We conclude that the FZD6-fibronectin actin axis identified in our study could be exploited for drug development in highly metastatic forms of breast cancer, such as TNBC. © 2016 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Mouse mammary stem cells express prognostic markers for triple-negative breast cancer.

INTRODUCTION: Triple-negative breast cancer (TNBC) is a heterogeneous group of tumours in which chemotherapy, the current mainstay of systemic treatment, is often initially beneficial but with a high risk of relapse and metastasis. There is currently no means of predicting which TNBC will relapse. We tested the hypothesis that the biological properties of normal stem cells are re-activated in tumour metastasis and that, therefore, the activation of normal mammary stem cell-associated gene sets in primary TNBC would be highly prognostic for relapse and metastasis. METHODS: Mammary basal stem and myoepithelial cells were isolated by flow cytometry and tested in low-dose transplant assays. Gene expression microarrays were used to establish expression profiles of the stem and myoepithelial populations; these were compared to each other and to our previously established mammary epithelial gene expression profiles. Stem cell genes were classified by Gene Ontology (GO) analysis and the expression of a subset analysed in the stem cell population at single cell resolution. Activation of stem cell genes was interrogated across different breast cancer cohorts and within specific subtypes and tested for clinical prognostic power. RESULTS: A set of 323 genes was identified that was expressed significantly more highly in the purified basal stem cells compared to all other cells of the mammary epithelium. A total of 109 out of 323 genes had been associated with stem cell features in at least one other study in addition to our own, providing further support for their involvement in the biology of this cell type. GO analysis demonstrated an enrichment of these genes for an association with cell migration, cytoskeletal regulation and tissue morphogenesis, consistent with a role in invasion and metastasis. Single cell resolution analysis showed that individual cells co-expressed both epithelial- and mesenchymal-associated genes/proteins. Most strikingly, we demonstrated that strong activity of this stem cell gene set in TNBCs identified those tumours most likely to rapidly progress to metastasis. CONCLUSIONS: Our findings support the hypothesis that the biological properties of normal stem cells are drivers of metastasis and that these properties can be used to stratify patients with a highly heterogeneous disease such as TNBC.

Natural History of Germline BRCA1 Mutated and BRCA Wild-type Triple-negative Breast Cancer.

We report a deep next-generation sequencing analysis of 13 sequentially obtained tumor samples, eight sequentially obtained circulating tumor DNA (ctDNA) samples and three germline DNA samples over the life history of 3 patients with triple-negative breast cancer (TNBC), 2 of whom had germline pathogenic BRCA1 mutation, to unravel tumor evolution. Tumor tissue from all timepoints and germline DNA was subjected to whole-exome sequencing (WES), custom amplicon deep sequencing (30,000X) of a WES-derived somatic mutation panel, and SNP arrays for copy-number variation (CNV), while whole transcriptome sequencing (RNA-seq) was performed only on somatic tumor.There was enrichment of homologous recombination deficiency signature in all tumors and widespread CNV, which remained largely stable over time. Somatic tumor mutation numbers varied between patients and within each patient (range: 70-216, one outlier). There was minimal mutational overlap between patients with TP53 being the sole commonly mutated gene, but there was substantial overlap in sequential samples in each patient. Each patient's tumor contained a founding ("stem") clone at diagnosis, which persisted over time, from which all other clones ("subclone") were derived ("branching evolution"), which contained mutations in well-characterized cancer-related genes like PDGFRB, ARID2, TP53 (Patient_02), TP53, BRAF, BRIP1, CSF3R (Patient_04), and TP53, APC, EZH2 (Patient_07). Including stem and subclones, tumors from all patients were polyclonal at diagnosis and during disease progression. ctDNA recapitulated most tissue-derived stem clonal and subclonal mutations while detecting some additional subclonal mutations. RNA-seq revealed a stable basal-like pattern, with most highly expressed variants belonging to stem clone. SIGNIFICANCE: In germline BRCA1 mutated and BRCA wild-type patients, TNBC shows a branching evolutionary pattern of mutations with a single founding clone, are polyclonal throughout their disease course, and have widespread copy-number aberrations. This evolutionary pattern may be associated with treatment resistance or sensitivity and could be therapeutically exploited.

PAK6 acts downstream of IQGAP3 to promote contractility in triple negative breast cancer cells.

Breast cancer is a heterogeneous disease that remains the most common malignancy among women worldwide. During genomic analysis of breast tumours, mRNA levels of IQGAP3 were found to be upregulated in triple negative tumours. IQGAP3 was subsequently found to be expressed across a panel of triple negative breast cancer (TNBC) cell lines. Depleting expression levels of IQGAP3 delivered elongated cells, disrupted cell migration, and inhibited the ability of cells to form specialised invasive adhesion structures, termed invadopodia. The morphological changes induced by IQGAP3 depletion were found to be dependent on RhoA. Indeed, reduced expression of IQGAP3 disrupted RhoA activity and actomyosin contractility. Interestingly, IQGAP3 was also found to interact with p-21 activated kinase 6 (PAK6); a protein already associated with the regulation of cell morphology. Moreover, PAK6 depletion phenocopied IQGAP3 depletion in these cells. Whereas PAK6 overexpression rescued the IQGAP3 depletion phenotype. Our work points to an important PAK6-IQGAP3-RhoA pathway that drives the cellular contractility of breast cancer cells promoting both cell migration and adhesive invasion of these cells. As this phenotype is relevant to the process of metastasis and re-seeding of metastasis, the pharmacological targeting of PAK6 could lead to clinical benefit in TNBC patients.

Adipose-enriched peri-tumoral stroma, in contrast to myofibroblast-enriched stroma, prognosticates poorer survival in breast cancers.

Despite our understanding of the genetic basis of intra-tumoral heterogeneity, the role of stromal heterogeneity arising from an altered tumor microenvironment in affecting tumorigenesis is poorly understood. In particular, extensive study on the peri-tumoral stroma in the morphologically normal tissues surrounding the tumor is lacking. Here, we examine the heterogeneity in tumors and peri-tumoral stroma from 8 ER+/PR+/HER2- invasive breast carcinomas, through multi-region transcriptomic profiling by microarray. We describe the regional heterogeneity observed at the intrinsic molecular subtype, pathway enrichment, and cell type composition levels within each tumor and its peri-tumoral region, up to 7 cm from the tumor margins. Moreover, we identify a pro-inflammatory adipose-enriched peri-tumoral subtype which was significantly associated with poorer overall survival in breast cancer patients, in contrast to an adaptive immune cell- and myofibroblast-enriched subtype. These data together suggest that peri-tumoral heterogeneity may be an important determinant of the evolution and treatment of breast cancers.

Integrated Multimodal Analyses of DNA Damage Response and Immune Markers as Predictors of Response in Metastatic Triple-Negative Breast Cancer in the TNT Trial (NCT00532727).

PURPOSE: The TNT trial (NCT00532727) showed no evidence of carboplatin superiority over docetaxel in metastatic triple-negative breast cancer (mTNBC), but carboplatin benefit was observed in the germline BRCA1/2 mutation subgroup. Broader response-predictive biomarkers are needed. We explored the predictive ability of DNA damage response (DDR) and immune markers. EXPERIMENTAL DESIGN: Tumor-infiltrating lymphocytes were evaluated for 222 of 376 patients. Primary tumors (PT) from 186 TNT participants (13 matched recurrences) were profiled using total RNA sequencing. Four transcriptional DDR-related and 25 immune-related signatures were evaluated. We assessed their association with objective response rate (ORR) and progression-free survival (PFS). Conditional inference forest clustering was applied to integrate multimodal data. The biology of subgroups was characterized by 693 gene expression modules and other markers. RESULTS: Transcriptional DDR-related biomarkers were not predictive of ORR to either treatment overall. Changes from PT to recurrence were demonstrated; in chemotherapy-naïve patients, transcriptional DDR markers separated carboplatin responders from nonresponders (P values = 0.017; 0.046). High immune infiltration was associated with docetaxel ORR (interaction P values < 0.05). Six subgroups were identified; the immune-enriched cluster had preferential docetaxel response [62.5% (D) vs. 29.4% (C); P = 0.016]. The immune-depleted cluster had preferential carboplatin response [8.0% (D) vs. 40.0% (C); P = 0.011]. DDR-related subgroups were too small to assess ORR. CONCLUSIONS: High immune features predict docetaxel response, and high DDR signature scores predict carboplatin response in treatment-naïve mTNBC. Integrating multimodal DDR and immune-related markers identifies subgroups with differential treatment sensitivity. Treatment options for patients with immune-low and DDR-proficient tumors remains an outstanding need. Caution is needed using PT-derived transcriptional signatures to direct treatment in mTNBC, particularly DDR-related markers following prior chemotherapy.

CDK Inhibition Primes for Anti-PD-L1 Treatment in Triple-Negative Breast Cancer Models.

Triple-negative breast cancers (TNBC) expressing PD-L1 qualify for checkpoint inhibitor immunotherapy. Cyclin E/CDK2 is a potential target axis in TNBC; however, small-molecule drugs at efficacious doses may be associated with toxicity, and treatment alongside immunotherapy requires investigation. We evaluated CDK inhibition at suboptimal levels and its anti-tumor and immunomodulatory effects. Transcriptomic analyses of primary breast cancers confirmed higher cyclin E/CDK2 expression in TNBC compared with non-TNBC. Out of the three CDK2-targeting inhibitors tested, the CDK 2, 7 and 9 inhibitor SNS-032 was the most potent in reducing TNBC cell viability and exerted cytotoxicity against all eight TNBC cell lines evaluated in vitro. Suboptimal SNS-032 dosing elevated cell surface PD-L1 expression in surviving TNBC cells. In mice engrafted with human immune cells and challenged with human MDA-MB-231 TNBC xenografts in mammary fat pads, suboptimal SNS-032 dosing partially restricted tumor growth, enhanced the tumor infiltration of human CD45+ immune cells and elevated cell surface PD-L1 expression in surviving cancer cells. In tumor-bearing mice engrafted with human immune cells, the anti-PD-L1 antibody avelumab, given sequentially following suboptimal SNS-032 dosing, reduced tumor growth compared with SNS-032 alone or with avelumab without prior SNS-032 priming. CDK inhibition at suboptimal doses promotes immune cell recruitment to tumors, PD-L1 expression by surviving TNBC cells and may complement immunotherapy.

Genomic analysis defines clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancer.

Ductal carcinoma in situ (DCIS) is the most common form of preinvasive breast cancer and, despite treatment, a small fraction (5-10%) of DCIS patients develop subsequent invasive disease. A fundamental biologic question is whether the invasive disease arises from tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial DCIS lesion and paired invasive recurrent tumors in 95 patients together with single-cell DNA sequencing in a subset of cases. Our data show that in 75% of cases the invasive recurrence was clonally related to the initial DCIS, suggesting that tumor cells were not eliminated during the initial treatment. Surprisingly, however, 18% were clonally unrelated to the DCIS, representing new independent lineages and 7% of cases were ambiguous. This knowledge is essential for accurate risk evaluation of DCIS, treatment de-escalation strategies and the identification of predictive biomarkers.

Random Forest Modelling of High-Dimensional Mixed-Type Data for Breast Cancer Classification.

Advances in high-throughput technologies encourage the generation of large amounts of multiomics data to investigate complex diseases, including breast cancer. Given that the aetiologies of such diseases extend beyond a single biological entity, and that essential biological information can be carried by all data regardless of data type, integrative analyses are needed to identify clinically relevant patterns. To facilitate such analyses, we present a permutation-based framework for random forest methods which simultaneously allows the unbiased integration of mixed-type data and assessment of relative feature importance. Through simulation studies and machine learning datasets, the performance of the approach was evaluated. The results showed minimal multicollinearity and limited overfitting. To further assess the performance, the permutation-based framework was applied to high-dimensional mixed-type data from two independent breast cancer cohorts. Reproducibility and robustness of our approach was demonstrated by the concordance in relative feature importance between the cohorts, along with consistencies in clustering profiles. One of the identified clusters was shown to be prognostic for clinical outcome after standard-of-care adjuvant chemotherapy and outperformed current intrinsic molecular breast cancer classifications.

Systemic immune reaction in axillary lymph nodes adds to tumor-infiltrating lymphocytes in triple-negative breast cancer prognostication.

The level of stromal tumor-infiltrating lymphocytes (sTILs) in triple-negative (TNBC) and HER2-positive breast cancers convey prognostic information. The importance of systemic immunity to local immunity is unknown in breast cancer. We previously demonstrated that histological alterations in axillary lymph nodes (LNs) carry clinical relevance. Here, we capture local immune responses by scoring TILs at the primary tumor and systemic immune responses by recording the formation of secondary follicles, also known as germinal centers, in 2,857 cancer-free and involved axillary LNs on haematoxylin and eosin (H&E) stained sections from a retrospective cohort of 161 LN-positive triple-negative and HER2-positive breast cancer patients. Our data demonstrate that the number of germinal center formations across all cancer-free LNs, similar to high levels of TILs, is associated with a good prognosis in low TILs TNBC. This highlights the importance of assessing both primary and LN immune responses for prognostication and for future breast cancer research.

Tumor-Infiltrating B Lymphocyte Profiling Identifies IgG-Biased, Clonally Expanded Prognostic Phenotypes in Triple-Negative Breast Cancer.

In breast cancer, humoral immune responses may contribute to clinical outcomes, especially in more immunogenic subtypes. Here, we investigated B lymphocyte subsets, immunoglobulin expression, and clonal features in breast tumors, focusing on aggressive triple-negative breast cancers (TNBC). In samples from patients with TNBC and healthy volunteers, circulating and tumor-infiltrating B lymphocytes (TIL-B) were evaluated. CD20+CD27+IgD- isotype-switched B lymphocytes were increased in tumors, compared with matched blood. TIL-B frequently formed stromal clusters with T lymphocytes and engaged in bidirectional functional cross-talk, consistent with gene signatures associated with lymphoid assembly, costimulation, cytokine-cytokine receptor interactions, cytotoxic T-cell activation, and T-cell-dependent B-cell activation. TIL-B-upregulated B-cell receptor (BCR) pathway molecules FOS and JUN, germinal center chemokine regulator RGS1, activation marker CD69, and TNFα signal transduction via NFκB, suggesting BCR-immune complex formation. Expression of genes associated with B lymphocyte recruitment and lymphoid assembly, including CXCL13, CXCR4, and DC-LAMP, was elevated in TNBC compared with other subtypes and normal breast. TIL-B-rich tumors showed expansion of IgG but not IgA isotypes, and IgG isotype switching positively associated with survival outcomes in TNBC. Clonal expansion was biased toward IgG, showing expansive clonal families with specific variable region gene combinations and narrow repertoires. Stronger positive selection pressure was present in the complementarity determining regions of IgG compared with their clonally related IgA in tumor samples. Overall, class-switched B lymphocyte lineage traits were conspicuous in TNBC, associated with improved clinical outcomes, and conferred IgG-biased, clonally expanded, and likely antigen-driven humoral responses. SIGNIFICANCE: Tumor-infiltrating B lymphocytes assemble in clusters, undergoing B-cell receptor-driven activation, proliferation, and isotype switching. Clonally expanded, IgG isotype-biased humoral immunity associates with favorable prognosis primarily in triple-negative breast cancers.

TGF-ß1 potentiates Vγ9Vδ2 T cell adoptive immunotherapy of cancer.

Despite its role in cancer surveillance, adoptive immunotherapy using γδ T cells has achieved limited efficacy. To enhance trafficking to bone marrow, circulating Vγ9Vδ2 T cells are expanded in serum-free medium containing TGF-ß1 and IL-2 (γδ[T2] cells) or medium containing IL-2 alone (γδ[2] cells, as the control). Unexpectedly, the yield and viability of γδ[T2] cells are also increased by TGF-ß1, when compared to γδ[2] controls. γδ[T2] cells are less differentiated and yet display increased cytolytic activity, cytokine release, and antitumor activity in several leukemic and solid tumor models. Efficacy is further enhanced by cancer cell sensitization using aminobisphosphonates or Ara-C. A number of contributory effects of TGF-ß are described, including prostaglandin E2 receptor downmodulation, TGF-ß insensitivity, and upregulated integrin activity. Biological relevance is supported by the identification of a favorable γδ[T2] signature in acute myeloid leukemia (AML). Given their enhanced therapeutic activity and compatibility with allogeneic use, γδ[T2] cells warrant evaluation in cancer immunotherapy.

A Four-gene Decision Tree Signature Classification of Triple-negative Breast Cancer: Implications for Targeted Therapeutics.

The molecular complexity of triple-negative breast cancers (TNBCs) provides a challenge for patient management. We set out to characterize this heterogeneous disease by combining transcriptomics and genomics data, with the aim of revealing convergent pathway dependencies with the potential for treatment intervention. A Bayesian algorithm was used to integrate molecular profiles in two TNBC cohorts, followed by validation using five independent cohorts (n = 1,168), including three clinical trials. A four-gene decision tree signature was identified, which robustly classified TNBCs into six subtypes. All four genes in the signature (EXO1, TP53BP2, FOXM1, and RSU1) are associated with either genomic instability, malignant growth, or treatment response. One of the six subtypes, MC6, encompassed the largest proportion of tumors (∼50%) in early diagnosed TNBCs. In TNBC patients with metastatic disease, the MC6 proportion was reduced to 25%, and was independently associated with a higher response rate to platinum-based chemotherapy. In TNBC cell line data, platinum sensitivity was recapitulated, and a sensitivity to the inhibition of the phosphatase PPM1D was revealed. Molecularly, MC6-TNBCs displayed high levels of telomeric allelic imbalances, enrichment of CD4+ and CD8+ immune signatures, and reduced expression of genes negatively regulating the MAPK signaling pathway. These observations suggest that our integrative classification approach may identify TNBC patients with discernible and theoretically pharmacologically tractable features that merit further studies in prospective trials.

Integrated genomics and functional validation identifies malignant cell specific dependencies in triple negative breast cancer.

Triple negative breast cancers (TNBCs) lack recurrent targetable driver mutations but demonstrate frequent copy number aberrations (CNAs). Here, we describe an integrative genomic and RNAi-based approach that identifies and validates gene addictions in TNBCs. CNAs and gene expression alterations are integrated and genes scored for pre-specified target features revealing 130 candidate genes. We test functional dependence on each of these genes using RNAi in breast cancer and non-malignant cells, validating malignant cell selective dependence upon 37 of 130 genes. Further analysis reveals a cluster of 13 TNBC addiction genes frequently co-upregulated that includes genes regulating cell cycle checkpoints, DNA damage response, and malignant cell selective mitotic genes. We validate the mechanism of addiction to a potential drug target: the mitotic kinesin family member C1 (KIFC1/HSET), essential for successful bipolar division of centrosome-amplified malignant cells and develop a potential selection biomarker to identify patients with tumors exhibiting centrosome amplification.

Mutation Processes in 293-Based Clones Overexpressing the DNA Cytosine Deaminase APOBEC3B.

Molecular, cellular, and clinical studies have combined to demonstrate a contribution from the DNA cytosine deaminase APOBEC3B (A3B) to the overall mutation load in breast, head/neck, lung, bladder, cervical, ovarian, and other cancer types. However, the complete landscape of mutations attributable to this enzyme has yet to be determined in a controlled human cell system. We report a conditional and isogenic system for A3B induction, genomic DNA deamination, and mutagenesis. Human 293-derived cells were engineered to express doxycycline-inducible A3B-eGFP or eGFP constructs. Cells were subjected to 10 rounds of A3B-eGFP exposure that each caused 80-90% cell death. Control pools were subjected to parallel rounds of non-toxic eGFP exposure, and dilutions were done each round to mimic A3B-eGFP induced population fluctuations. Targeted sequencing of portions of TP53 and MYC demonstrated greater mutation accumulation in the A3B-eGFP exposed pools. Clones were generated and microarray analyses were used to identify those with the greatest number of SNP alterations for whole genome sequencing. A3B-eGFP exposed clones showed global increases in C-to-T transition mutations, enrichments for cytosine mutations within A3B-preferred trinucleotide motifs, and more copy number aberrations. Surprisingly, both control and A3B-eGFP clones also elicited strong mutator phenotypes characteristic of defective mismatch repair. Despite this additional mutational process, the 293-based system characterized here still yielded a genome-wide view of A3B-catalyzed mutagenesis in human cells and a system for additional studies on the compounded effects of simultaneous mutation mechanisms in cancer cells.

PIM1 kinase regulates cell death, tumor growth and chemotherapy response in triple-negative breast cancer.

Triple-negative breast cancers (TNBCs) have poor prognosis and lack targeted therapies. Here we identified increased copy number and expression of the PIM1 proto-oncogene in genomic data sets of patients with TNBC. TNBC cells, but not nonmalignant mammary epithelial cells, were dependent on PIM1 for proliferation and protection from apoptosis. PIM1 knockdown reduced expression of the anti-apoptotic factor BCL2, and dynamic BH3 profiling of apoptotic priming revealed that PIM1 prevents mitochondrial-mediated apoptosis in TNBC cell lines. In TNBC tumors and their cellular models, PIM1 expression was associated with several transcriptional signatures involving the transcription factor MYC, and PIM1 depletion in TNBC cell lines decreased, in a MYC-dependent manner, cell population growth and expression of the MYC target gene MCL1. Treatment with the pan-PIM kinase inhibitor AZD1208 impaired the growth of both cell line and patient-derived xenografts and sensitized them to standard-of-care chemotherapy. This work identifies PIM1 as a malignant-cell-selective target in TNBC and the potential use of PIM1 inhibitors for sensitizing TNBC to chemotherapy-induced apoptotic cell death.

Do Genetic Susceptibility Variants Associate with Disease Severity in Early Active Rheumatoid Arthritis?

OBJECTIVE: Genetic variants affect both the development and severity of rheumatoid arthritis (RA). Recent studies have expanded the number of RA susceptibility variants. We tested the hypothesis that these associated with disease severity in a clinical trial cohort of patients with early, active RA. METHODS: We evaluated 524 patients with RA enrolled in the Combination Anti-Rheumatic Drugs in Early RA (CARDERA) trials. We tested validated susceptibility variants - 69 single-nucleotide polymorphisms (SNP), 15 HLA-DRB1 alleles, and amino acid polymorphisms in 6 HLA molecule positions - for their associations with progression in Larsen scoring, 28-joint Disease Activity Scores, and Health Assessment Questionnaire (HAQ) scores over 2 years using linear mixed-effects and latent growth curve models. RESULTS: HLA variants were associated with joint destruction. The *04:01 SNP (rs660895, p = 0.0003), *04:01 allele (p = 0.0002), and HLA-DRß1 amino acids histidine at position 13 (p = 0.0005) and valine at position 11 (p = 0.0012) significantly associated with radiological progression. This association was only significant in anticitrullinated protein antibody (ACPA)-positive patients, suggesting that while their effects were not mediated by ACPA, they only predicted joint damage in ACPA-positive RA. Non-HLA variants did not associate with radiograph damage (assessed individually and cumulatively as a weighted genetic risk score). Two SNP - rs11889341 (STAT4, p = 0.0001) and rs653178 (SH2B3-PTPN11, p = 0.0004) - associated with HAQ scores over 6-24 months. CONCLUSION: HLA susceptibility variants play an important role in determining radiological progression in early, active ACPA-positive RA. Genome-wide and HLA-wide analyses across large populations are required to better characterize the genetic architecture of radiological progression in RA.