Evaluation of AlphaFold structure-based protein stability prediction on missense variations in cancer.

Identifying pathogenic missense variants in hereditary cancer is critical to the efforts of patient surveillance and risk-reduction strategies. For this purpose, many different gene panels consisting of different number and/or set of genes are available and we are particularly interested in a panel of 26 genes with a varying degree of hereditary cancer risk consisting of ABRAXAS1, ATM, BARD1, BLM, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, EPCAM, MEN1, MLH1, MRE11, MSH2, MSH6, MUTYH, NBN, PALB2, PMS2, PTEN, RAD50, RAD51C, RAD51D, STK11, TP53, and XRCC2. In this study, we have compiled a collection of the missense variations reported in any of these 26 genes. More than a thousand missense variants were collected from ClinVar and the targeted screen of a breast cancer cohort of 355 patients which contributed to this set with 160 novel missense variations. We analyzed the impact of the missense variations on protein stability by five different predictors including both sequence- (SAAF2EC and MUpro) and structure-based (Maestro, mCSM, CUPSAT) predictors. For the structure-based tools, we have utilized the AlphaFold (AF2) protein structures which comprise the first structural analysis of this hereditary cancer proteins. Our results agreed with the recent benchmarks that computed the power of stability predictors in discriminating the pathogenic variants. Overall, we reported a low-to-medium-level performance for the stability predictors in discriminating pathogenic variants, except MUpro which had an AUROC of 0.534 (95% CI [0.499-0.570]). The AUROC values ranged between 0.614-0.719 for the total set and 0.596-0.682 for the set with high AF2 confidence regions. Furthermore, our findings revealed that the confidence score for a given variant in the AF2 structure could alone predict pathogenicity more robustly than any of the tested stability predictors with an AUROC of 0.852. Altogether, this study represents the first structural analysis of the 26 hereditary cancer genes underscoring 1) the thermodynamic stability predicted from AF2 structures as a moderate and 2) the confidence score of AF2 as a strong descriptor for variant pathogenicity.

Analysis of Mitochondrial DNA Control Region D-Loop in Gliomas: Result of 52 Patients.

AIM: To investigate the effect of mitochondrial DNA (mtDNA) variants mainly in D-loop on glioma biology. MATERIAL AND METHODS: Sanger sequencing of D-loop (15971?16451 bp) for 52 glioma patients was performed and the variations were statistically analyzed for gender, WHO classification, morphological grade, IDH/TERT status. RESULTS: Total of 122 variations (51 unique) were identified in 52 patients. C16223T, T16189C, T16311C and T16126C variants were frequently detected. The total variation number was statistically non-significant among the analyzed categories. When individual variants were considered, T16311C and T16224C were statistically significant for WHO classification (p=0.033), morphological grade (p=0.036) and gender (p=0.039), respectively. CONCLUSION: Total variation number in D-loop was not found to be related with clinical variables. Our data suggests that individual variants may play a critical role in glioma biology.

Mutations and Copy Number Alterations in IDH Wild-Type Glioblastomas Are Shaped by Different Oncogenic Mechanisms.

Little is known about the mutational processes that shape the genetic landscape of gliomas. Numerous mutational processes leave marks on the genome in the form of mutations, copy number alterations, rearrangements or their combinations. To explore gliomagenesis, we hypothesized that gliomas with different underlying oncogenic mechanisms would have differences in the burden of various forms of these genomic alterations. This was an analysis on adult diffuse gliomas, but IDH-mutant gliomas as well as diffuse midline gliomas H3-K27M were excluded to search for the possible presence of new entities among the very heterogenous group of IDH-WT glioblastomas. The cohort was divided into two molecular subsets: (1) Molecularly-defined GBM (mGBM) as those that carried molecular features of glioblastomas (including TERT promoter mutations, 7/10 pattern, or EGFR-amplification), and (2) those who did not (others). Whole exome sequencing was performed for 37 primary tumors and matched blood samples as well as 8 recurrences. Single nucleotide variations (SNV), short insertion or deletions (indels) and copy number alterations (CNA) were quantified using 5 quantitative metrics (SNV burden, indel burden, copy number alteration frequency-wGII, chromosomal arm event ratio-CAER, copy number amplitude) as well as 4 parameters that explored underlying oncogenic mechanisms (chromothripsis, double minutes, microsatellite instability and mutational signatures). Findings were validated in the TCGA pan-glioma cohort. mGBM and "Others" differed significantly in their SNV (only in the TCGA cohort) and CNA metrics but not indel burden. SNV burden increased with increasing age at diagnosis and at recurrences and was driven by mismatch repair deficiency. On the contrary, indel and CNA metrics remained stable over increasing age at diagnosis and with recurrences. Copy number alteration frequency (wGII) correlated significantly with chromothripsis while CAER and CN amplitude correlated significantly with the presence of double minutes, suggesting separate underlying mechanisms for different forms of CNA.

Toward In Vitro Epigenetic Drug Design for Thyroid Cancer: The Promise of PF-03814735, an Aurora Kinase Inhibitor.

Thyroid cancer (TC) is a very common malignancy worldwide. Chief among the innovative molecular drug targets for TC are epigenetic modifications. Increased telomerase activity in cancer cells makes telomerase a novel target for epigenetic anticancer drug innovation. Recently, telomerase reverse transcriptase (TERT) gene promoter (TERTp) mutations (C228T and C250T) were reported at high frequency in TC cell lines and tumor biopsies. In this study, three representative TC cell lines, mutant TERTp (TPC1), mutant BRAF/TERTp (KTC2), and wild-type TERTp (WRO), were screened with a drug library composed of 51 epigenetic drugs: 14 Aurora kinase inhibitors; 23 histone deacetylase inhibitors; 5 sirtuin modifiers; 3 hypoxia-inducible factor inhibitors; 2 DNA methyltransferase inhibitors; 2 histone methyltransferase inhibitors, a histone demethylase inhibitor, and a bromodomain inhibitor. Effects of the drugs on cell growth at 48 and 72 h were compared. PF-03814735, a small-molecule inhibitor of Aurora kinase A (IC50 = 0.8 nM) and B (IC50 = 5 nM), was the most potent on KTC2 cells, whereas CUDC-101, a multitarget inhibitor, was effective on both WRO and KTC2 cells. Notably, PF-03814735 was found to be the most effective epigenetic drug on cell lines harboring the C228T mutation. In conclusion, these new findings offer specific guidance on dose and time course selection to design novel therapeutic interventions against TC using PF-03814735, and specifically target cells carrying the TERTpC228T mutation. In a larger context of drug discovery science, these findings inform new strategies to forecast optimal treatment regimens for TC, particularly with Aurora kinase inhibitors and in ways guided by epigenetic drug design.

Whole exome sequencing-based analysis to identify DNA damage repair deficiency as a major contributor to gliomagenesis in adult diffuse gliomas.

OBJECTIVE: Processes that cause or contribute to cancer, such as aging, exposure to carcinogens, or DNA damage repair deficiency (DDRd), create predictable and traceable nucleotide alterations in one's genetic code (termed "mutational signatures"). Large studies have previously identified various such mutational signatures across cancers that can be attributed to the specific causative processes. To gain further insight into the processes in glioma development, the authors analyzed mutational signatures in adult diffuse gliomas (DGs). METHODS: Twenty-five DGs and paired blood samples were whole exome sequenced. Somatic mutational signatures were identified using 2 different methods. Associations of the signatures with age at diagnosis, molecular subset, and mutational load were investigated. As DDRd-related signatures were frequently observed, germline and somatic DDR gene mutations as well as microsatellite instability (MSI) status were determined for all samples. For validation of signature prevalence, publicly available data from The Cancer Genome Atlas (TCGA) were used. RESULTS: Each tumor had a unique combination of signatures. The most common signatures were signature 1 (88%, aging related), signature 3 (52%, homologous recombination related), and signature 15 (56%, mismatch repair related). Eighty-four percent of the tumors contained at least 1 DDRd signature. The findings were validated using public TCGA data. The weight of signature 1 positively correlated with age (r = 0.43) while cumulative weight of DDRd signatures negatively correlated with age (r = -0.16). Each subject had at least 1 germline/somatic alteration in a DDR gene, the most common being the risk single nucleotide polymorphism rs1800734 in MLH1. The rs1800734-AA genotype had a higher cumulative DDRd weight as well as higher mutational load; TP53 was the most common somatically altered DDR gene. MSI was observed in 24% of the tumors. No significant associations of MSI status with mutational load, rs1800734, or the cumulative weight of DDRd signatures were identified. CONCLUSIONS: Current findings suggest that DDRd may act as a fundamental mechanism in gliomagenesis rather than being a random, secondary event.

Use of telomerase promoter mutations to mark specific molecular subsets with reciprocal clinical behavior in IDH mutant and IDH wild-type diffuse gliomas.

OBJECTIVE Recent studies have established that hemispheric diffuse gliomas may be grouped into subsets on the basis of molecular markers; these subsets are loosely correlated with the histopathological diagnosis but are strong predictors of clinical tumor behavior. Based on an analysis of molecular and clinical parameters, the authors hypothesized that mutations of the telomerase promoter (TERTp-mut) mark separate oncogenic programs among isocitrate dehydrogenase 1 and/or 2 (IDH) mutant (IDH-mut) and IDH wild-type (IDH-wt) diffuse gliomas independent of histopathology or WHO grade. METHODS Four molecular subsets of the combined statuses of IDH and TERT-promoter mutations (double mutant, IDH only, TERT only, and double negative) were defined. Differences in age, anatomical location, molecular genetics, and survival rates in a surgical cohort of 299 patients with a total of 356 hemispheric diffuse gliomas (WHO Grade II, III, or IV) were analyzed. RESULTS TERTp-mut were present in 38.8% of IDH-mut and 70.2% of IDH-wt gliomas. The mutational status was stable in each patient at 57 recurrence events over a 2645-month cumulative follow-up period. Among patients with IDH-mut gliomas, those in the double-mutant subset had better survival and a lower incidence of malignant degeneration than those in the IDH-only subset. Of patients in the double-mutant subset, 96.3% were also positive for 1p/19q codeletions. All patients with 1p/19q codeletions had TERTp-mut. In patients with IDH-mut glioma, epidermal growth factor receptor or phosphatase and tensin homolog mutations were not observed, and copy-number variations were uncommon. Among IDH-wt gliomas, the TERT-only subset was associated with significantly higher age, higher Ki-67 labeling index, primary glioblastoma-specific oncogenic changes, and poor survival. The double-negative subset was genetically and biologically heterogeneous. Survival analyses (Kaplan-Meier, multivariate, and regression-tree analyses) confirmed that patients in the 4 molecular subsets had distinct prognoses. CONCLUSIONS Molecular subsets result in different tumor biology and clinical behaviors in hemispheric diffuse gliomas.

Evidence from autoimmune thyroiditis of skewed X-chromosome inactivation in female predisposition to autoimmunity.

The etiologic factors in the development of autoimmune thyroid diseases (AITDs) are not fully understood. We investigated the role of skewed X-chromosome inactivation (XCI) mosaicism in female predisposition to AITDs. One hundred and ten female AITDs patients (81 Hashimoto's thyroiditis (HT), 29 Graves' disease (GD)), and 160 female controls were analyzed for the androgen receptor locus by the HpaII/polymerase chain reaction assay to assess XCI patterns in DNA extracted from peripheral blood cells. In addition, thyroid biopsy, buccal mucosa, and hair follicle specimens were obtained from five patients whose blood revealed an extremely skewed pattern of XCI, and the analysis was repeated. Skewed XCI was observed in DNA from peripheral blood cells in 28 of 83 informative patients (34%) as compared with 10 of 124 informative controls (8%, P<0.0001). Extreme skewing was present in 16 patients (19%), but only in three controls (2.4%, P<0.0001). The buccal mucosa, and although less marked, the thyroid specimens also showed skewing. Analysis of two familial cases showed that only the affected individuals demonstrate skewed XCI patterns. Based on these results, skewed XCI mosaicism may play a significant role in the pathogenesis of AITDs.

Endothelial progenitor cells display clonal restriction in multiple myeloma.

BACKGROUND: In multiple myeloma (MM), increased neoangiogenesis contributes to tumor growth and disease progression. Increased levels of endothelial progenitor cells (EPCs) contribute to neoangiogenesis in MM, and, importantly, covary with disease activity and response to treatment. In order to understand the mechanisms responsible for increased EPC levels and neoangiogenic function in MM, we investigated whether these cells were clonal by determining X-chromosome inactivation (XCI) patterns in female patients by a human androgen receptor assay (HUMARA). In addition, EPCs and bone marrow cells were studied for the presence of clonotypic immunoglobulin heavy-chain (IGH) gene rearrangement, which indicates clonality in B cells; thus, its presence in EPCs would indicate a close genetic link between tumor cells in MM and endothelial cells that provide tumor neovascularization. METHODS: A total of twenty-three consecutive patients who had not received chemotherapy were studied. Screening in 18 patients found that 11 displayed allelic AR in peripheral blood mononuclear cells, and these patients were further studied for XCI patterns in EPCs and hair root cells by HUMARA. In 2 patients whose EPCs were clonal by HUMARA, and in an additional 5 new patients, EPCs were studied for IGH gene rearrangement using PCR with family-specific primers for IGH variable genes (VH). RESULTS: In 11 patients, analysis of EPCs by HUMARA revealed significant skewing (> or = 77% expression of a single allele) in 64% (n = 7). In 4 of these patients, XCI skewing was extreme (> or = 90% expression of a single allele). In contrast, XCI in hair root cells was random. Furthermore, PCR amplification with VH primers resulted in amplification of the same product in EPCs and bone marrow cells in 71% (n = 5) of 7 patients, while no IGH rearrangement was found in EPCs from healthy controls. In addition, in patients with XCI skewing in EPCs, advanced age was associated with poorer clinical status, unlike patients whose EPCs had random XCI. CONCLUSION: Our results suggest that EPCs in at least a substantial subpopulation of MM patients are related to the neoplastic clone and that this is an important mechanism for upregulation of tumor neovascularization in MM.

IDH-mutant glioma specific association of rs55705857 located at 8q24.21 involves MYC deregulation.

The single nucleotide polymorphism rs55705857, located in a non-coding but evolutionarily conserved region at 8q24.21, is strongly associated with IDH-mutant glioma development and was suggested to be a causal variant. However, the molecular mechanism underlying this association has remained unknown. With a case control study in 285 gliomas, 316 healthy controls, 380 systemic cancers, 31 other CNS-tumors, and 120 IDH-mutant cartilaginous tumors, we identified that the association was specific to IDH-mutant gliomas. Odds-ratios were 9.25 (5.17-16.52; 95% CI) for IDH-mutated gliomas and 12.85 (5.94-27.83; 95% CI) for IDH-mutated, 1p/19q co-deleted gliomas. Decreasing strength with increasing anaplasia implied a modulatory effect. No somatic mutations were noted at this locus in 114 blood-tumor pairs, nor was there a copy number difference between risk-allele and only-ancestral allele carriers. CCDC26 RNA-expression was rare and not different between the two groups. There were only minor subtype-specific differences in common glioma driver genes. RNA sequencing and LC-MS/MS comparisons pointed to significantly altered MYC-signaling. Baseline enhancer activity of the conserved region specifically on the MYC promoter and its further positive modulation by the SNP risk-allele was shown in vitro. Our findings implicate MYC deregulation as the underlying cause of the observed association.

Expression of IFITM1 in chronic myeloid leukemia patients.

We investigated the peripheral blood gene expression profile of interferon induced transmembrane protein 1 (IFITM1) in sixty chronic myeloid leukemia (CML) patients classified according to new prognostic score (NPS). IFITM1 is a component of a multimeric complex involved in the trunsduction of antiproliferative and cell adhesion signals. Expression level of IFITM1 was found significantly different between the high- and low-risk groups (P = 9.7976 x 10(-11)) by real-time reverse transcription polymerase chain reaction (RT-PCR). Higher IFITM1 expression correlated with improved survival (P = 0.01). These results indicate that IFITM1 expression profiling could be used for molecular classification of CML, which may also predict survival.

Somatic mosaicism for a MECP2 mutation associated with classic Rett syndrome in a boy.

Rett syndrome is a severe neurodevelopmental disorder that arises from mutations in the X-linked MECP2 gene. It is almost exclusively seen in girls due to the predominant occurrence of the mutations on the paternal X-chromosome, and also the early postnatal lethal effect of the disease causing mutations in hemizygous boys. We identified a boy with features of classic Rett syndrome who is mosaic for the truncating MECP2 mutation R270X. Chromosome analysis showed normal karyotype. These results indicate that a MECP2 mutation associated with Rett syndrome in females could lead to a similar phenotype in males as a result of somatic mosaicism.

Lack of association between RNASEL Arg462Gln variant and the risk of breast cancer.

BACKGROUND: The RNASEL G1385A variant was recently found to be implicated in the development of prostate cancer. Considering the function of RNase L and the pleiotropic effects of mutations associated with cancer, we sought to investigate whether the RNASEL G1385A variant is a risk factor for breast cancer. PATIENTS AND METHODS: A total of 453 breast cancer patients and 382 age- and sex-matched controls from Greece and Turkey were analyzed. Genotyping for the RNASEL G1385A variant was performed using an Amplification Refractory Mutation System (ARMS). RESULTS: Statistical evaluation of the RNASEL G1385A genotype distribution among breast cancer patients and controls revealed no significant association between the presence of the risk genotype and the occurrence of breast cancer. CONCLUSION: Although an increasing number of studies report an association between the RNASEL G1385A variant and prostate cancer risk; this variant does not appear to be implicated in the development of breast cancer.

The SOCS-1 gene methylation in chronic myeloid leukemia patients.

SOCS-1, an important protein in the JAK/STAT pathway, has a role in the down stream of BCR-ABL protein kinase. We investigated 56 CML patients and 16 controls for the methylation status of SOCS-1 gene promoter and Exon 2 regions. Exon 2 was found to be methylated in 58.9% of the patients and 93.8% of the controls [P = 0.020, OR = 0.121(0.015-0.957)%95CI]. The promoter region was found unmethylated in all patient samples and controls. Although previous studies revealed a relation between SOCS1 gene Exon-2 hypermethylation and CML development or progression, the results of this study showed no such correlation. On the contrary, our results might be indicating hypomethylation in CML patients, this hypothesis need to be studied in larger study population.

Extremely skewed X-chromosome inactivation patterns in women with recurrent spontaneous abortion.

BACKGROUND: The role of extremely skewed X-chromosome inactivation (XCI) has been questioned in the pathogenesis of recurrent spontaneous abortion (RSA) but the results obtained were conflicting. AIMS: We therefore investigated the XCI patterns in peripheral blood DNA obtained from 80 patients who had RSA and 160 age-matched controls. METHODS: Pregnancy history, age, karyotype, and disease information was collected from all subjects. The methylation status of a highly polymorphic cytosine-adenine-guanine repeat in the androgen-receptor (AR) gene was determined by use of methylation-sensitive restriction enzyme HpaII and polymerase chain reaction. RESULTS: Skewed XCI (> 85% skewing) was observed in 13 of the 62 patients informative for the AR polymorphism (20.9%), and eight of the 124 informative controls (6.4%) (P = 0.0069; chi2 test). More importantly, extremely skewed XCI, defined as > 90% inactivation of one allele, was present in 11 (17.7%) patients, and in only two controls (P = 0.0002; chi2 test). CONCLUSIONS: These results support the interpretation that disturbances in XCI mosaicism may be involved in the pathogenesis of RSA.

Skewed X chromosome inactivation in blood cells of women with scleroderma.

OBJECTIVE: Scleroderma (SSc) is an autoimmune disease of unknown etiology. The disease is 3-8 times more frequent in women than in men. The role of X chromosome inactivation (XCI) in the predisposition of women to autoimmunity has been questioned. Until now this has not been illustrated experimentally. This study was undertaken to test the hypothesis that disturbances in XCI mosaicism may be involved in the pathogenesis of the disease in female patients with SSc. METHODS: Seventy female SSc patients and 160 female controls were analyzed for the androgen receptor locus by the Hpa II/polymerase chain reaction assay to assess XCI patterns in DNA extracted from peripheral blood cells. Furthermore, skin biopsy samples were obtained from 5 patients whose blood revealed an extremely skewed pattern of XCI, and the analysis repeated. Since microchimerism in SSc was reported, Y chromosome sequences were investigated in all samples. RESULTS: Skewed XCI was observed in DNA from peripheral blood cells in 35 of 55 informative patients (64%), as compared with 10 of 124 informative controls (8%) (P < 0.0001). Extreme skewing was present in 27 of the patient group (49%), as compared with only 3 of the controls (2.4%) (P < 0.0001). However, XCI was random in all skin biopsy samples. The potential contribution of microchimerism to the random XCI pattern is highly unlikely based on the medical histories of the patients. CONCLUSION: Skewed XCI mosaicism may play a significant role in the pathogenesis of SSc.