Association of three SNPs in TOX3 and breast cancer risk: Evidence from 97275 cases and 128686 controls

Expression- and splicing-based multi-tissue transcriptome-wide association studies identified multiple genes for breast cancer by estrogen-receptor status

BACKGROUND

Although several transcriptome-wide association studies (TWASs) have been performed to identify genes associated with overall breast cancer (BC) risk, only a few TWAS have explored the differences in estrogen receptor-positive (ER+) and estrogen receptor-negative (ER-) breast cancer. Additionally, these studies were based on gene expression prediction models trained primarily in breast tissue, and they did not account for alternative splicing of genes.

METHODS

In this study, we utilized two approaches to perform multi-tissue TWASs of breast cancer by ER subtype: (1) an expression-based TWAS that combined TWAS signals for each gene across multiple tissues and (2) a splicing-based TWAS that combined TWAS signals of all excised introns for each gene across tissues. To perform this TWAS, we utilized summary statistics for ER + BC from the Breast Cancer Association Consortium (BCAC) and for ER- BC from a meta-analysis of BCAC and the Consortium of Investigators of Modifiers of BRCA1 and BRCA2 (CIMBA).

RESULTS

In total, we identified 230 genes in 86 loci that were associated with ER + BC and 66 genes in 29 loci that were associated with ER- BC at a Bonferroni threshold of significance. Of these genes, 2 genes associated with ER + BC at the 1q21.1 locus were located at least 1 Mb from published GWAS hits. For several well-studied tumor suppressor genes such as TP53 and CHEK2 which have historically been thought to impact BC risk through rare, penetrant mutations, we discovered that common variants, which modulate gene expression, may additionally contribute to ER + or ER- etiology.

CONCLUSIONS

Our study comprehensively examined how differences in common variation contribute to molecular differences between ER + and ER- BC and introduces a novel, splicing-based framework that can be used in future TWAS studies.

Risk Association of TOX3 and MMP7 Gene Polymorphisms with Sporadic Breast Cancer in Mexican Women

Breast cancer (BC) has one of the highest incidences and mortality worldwide. Single nucleotide polymorphisms (SNPs) in TOX3 rs3803662 and MMP7 rs1943779 have been associated with susceptibility to BC. In this case-control study, we evaluated the association of rs3803662 (TOX3)/rs1943779 (MMP7) SNPs with clinical features, immunohistochemical reactivity, and risk association with BC in women from northeastern Mexico. We compared 212 BC cases and 212 controls. DNA was isolated from peripheral blood to perform the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay. We calculated genotype frequencies, odds ratios, and 95% confidence intervals. We found that CT (Cytocine-Thymine) and TT (Thymine -Thymine) genotypes, and T alleles of TOX3 rs3803662, were associated with BC risk (p = 0.034, p = 0.011, respectively). SNP TOX3 rs3803662 was associated with positive progesterone receptors (PR) and triple-negative BC (TNBC) but not with estrogen receptor (ER) or HER2 reactivity. CT and TT genotypes (p = 0.006) and T alleles (p = 0.002) of SNP MMP7 rs1943779 were associated with risk of BC. We found that T alleles of TOX3 rs3803662 and MMP7 rs1943779 SNPs are associated with BC risk. These findings contribute to personalized medicine in Mexican women.

A review of databases predicting the effects of SNPs in miRNA genes or miRNA-binding sites

Modern precision medicine comprises the knowledge and understanding of individual differences in the genomic sequence of patients to provide tailor-made treatments. Regularly, such variants are considered in coding regions only, and their effects are predicted based on their impact on the amino acid sequence of expressed proteins. However, assessing the effects of variants in noncoding elements, in particular microRNAs (miRNAs) and their binding sites, is important as well, as a single miRNA can influence the expression patterns of many genes at the same time. To analyze the effects of variants in miRNAs and their target sites, several databases storing variant impact predictions have been published. In this review, we will compare the core functionalities and features of these databases and discuss the importance of up-to-date data resources in the context of web applications. Finally, we will outline some recommendations for future developments in the field.

An in silico approach to characterize nonsynonymous SNPs and regulatory SNPs in human TOX3 gene

Cancer is one of the deadliest complex diseases having multigene nature where the role of single-nucleotide polymorphism (SNP) has been well explored in multiple genes. TOX high mobility group box family member 3 (TOX3) is one such gene, in which SNPs have been found to be associated with breast cancer. In this study, we have examined the potentially damaging nonsynonymous SNPs(nsSNPs) in TOX3 gene using in silico tools, namely PolyPhen2, SNP&GO, PhD-SNP and PROVEAN, which were further confirmed by I-Mutant, MutPred1.2 and ConSurf for their stability, functional and structural effects. nsSNPs rs368713418 (A266D), rs751141352 (P273S, P273T), rs200878352 (A275T) have been found to be the most deleterious that may have a vital role in breast cancer. Premature stop codon producing SNPs (Q527STOP), rs1259790811 (G495STOP), rs1294465822 (S395STOP) and rs1335372738 (G8STOP) were also found having prime importance in truncated and malfunctional protein formation. We also characterized regulatory SNPs for its potential effect on TOX3 gene regulation and found nine SNPs that may affect the gene regulation. Further, we have also designed 3D models using I-TASSER for the wild type and four mutant TOX3 proteins. Our study concludes that these SNPs can be of prime importance while studying breast cancer and other associated diseases as well. They are required to be studied in model organisms and cell cultures, and may have potential importance in personalized medicines and gene therapy.

Significant association of TOX3/LOC643714 locus-rs3803662 and breast cancer risk in a cohort of Iranian population

Genome-wide association studies normally focus on low penetrance and moderate to high-frequency single nucleotide polymorphisms (SNPs), which lead to genetic susceptibility to breast cancer. In this regard, the T allele of rs3803662 has been associated with breast cancer risk and with lower expression level of TOX3. We aimed to assess the risk of breast cancer associated with this polymorphism in an Iranian population. Using Tetra Primer ARMS PCR, rs3803662 was analyzed in a total of 943 individuals (430 cases and 513 healthy controls form North East of Iran). Allele frequencies and genotype distribution were analyzed in case and control samples to find out any association using the Chi-squared test and Logistic regression. All cases were pathologically confirmed; all controls were mainly healthy individuals. Genotype frequencies were found to be in agreement with HWE in controls and cases. TOX3-rs3803662 SNP was associated with breast cancer risk in our study (T vs. C allele contrast model: OR 1.36, 95% CI 1.12-1.64, P_value = 0.002; TT vs. CT + TT dominant model: OR 0.67, 95% CI 0.51-0.87, P_value = 0.003; TT vs. CT + CC recessive model: OR 1.54, 95% CI 1.02-2.30, P_vlue = 0.036). Moreover, after adjusting for age, BMI, history of previous cancer and also family history of cancer, all results, except for the recessive model, were remained significant. TOX3-rs3803662, may confer some degrees of risk of breast cancer in Iranian population. This finding is in line with similar results in other populations. It highlights the importance of TOX3 pathway in tumorigenesis.

Germline mutations in 40 cancer susceptibility genes among Chinese patients with high hereditary risk breast cancer

Multigene panel testing of breast cancer predisposition genes have been extensively conducted in Europe and America, which is relatively rare in Asia however. In this study, we assessed the frequency of germline mutations in 40 cancer predisposition genes, including BRCA1 and BRCA2, among a large cohort of Chinese patients with high hereditary risk of BC. From 2015 to 2016, consecutive BC patients from 26 centers of China with high hereditary risk were recruited (n = 937). Clinical information was collected and next-generation sequencing (NGS) was performed using blood samples of participants to identify germline mutations. In total, we acquired 223 patients with putative germline mutations, including 159 in BRCA1/2, 61 in 15 other BC susceptibility genes and 3 in both BRCA1/2 and non-BRCA1/2 gene. Major mutant non-BRCA1/2 genes were TP53 (n = 18), PALB2 (n = 11), CHEK2 (n = 6), ATM (n = 6) and BARD1 (n = 5). No factors predicted pathologic mutations in non-BRCA1/2 genes when treated as a whole. TP53 mutations were associated with HER-2 positive BC and younger age at diagnosis; and CHEK2 and PALB2 mutations were enriched in patients with luminal BC. Among high hereditary risk Chinese BC patients, 23.8% contained germline mutations, including 6.8% in non-BRCA1/2 genes. TP53 and PALB2 had a relatively high mutation rate (1.9 and 1.2%). Although no factors predicted for detrimental mutations in non-BRCA1/2 genes, some clinical features were associated with mutations of several particular genes.

The breast cancer susceptibility-related polymorphisms at the TOX3/LOC643714 locus associated with lung cancer risk in a Han Chinese population

It has been well established that besides environmental factors, genetic factors are also associated with lung cancer risk. However, to date, the prior identified genetic variants and loci only explain a small fraction of the familial risk of lung cancer. Hence it is vital to investigate the remaining missing heritability to understand the development and process of lung cancer. In the study, to test our hypothesis that the previously identified breast cancer risk-associated genetic polymorphisms at the TOX3/LOC643714 locus might contribute to lung cancer risk, 16 SNPs at the TOX3/LOC643714 locus were evaluated in a Han Chinese population based on a case-control study. Pearson's chi-square test or Fisher's exact test revealed that rs9933638, rs12443621, and rs3104746 were significantly associated with lung cancer risk (P < 0.001, P < 0.001, and P = 0.005, respectively). Logistic regression analyses displayed that lung cancer risk of individuals with rs9933638(GG+GA) were 1.89 times higher than that of rs9933638AA carriers (OR = 1.893, 95% CI = 1.308-2.741, P = 0.001). Similar findings were manifested for rs12443621 (OR = 1.824, 95% CI = 1.272-2.616, P = 0.001, rs12443621(GG+GA) carriers vs. rs12443621AA carriers) and rs3104746 (OR = 1.665, 95% CI = 1.243-2.230, P = 0.001, rs3104746TT carriers vs. rs3104746(TA+AA) carriers). The study discovered for the first time that three SNPs (rs9933638, rs12443621, and rs3104746) at the TOX3/LOC643714 locus contributed to lung cancer risk, providing new evidences that lung cancer and breast cancer are linked at the molecular and genetic level to a certain extent.

Polymorphisms in TIM-3 and breast cancer susceptibility in Chinese women: A case-control study

Previous studies have found associations between polymorphisms in T cell immunoglobulin and mucin domain 3 (TIM-3) and increased risks of various cancers. However, the association between TIM-3 polymorphisms and breast cancer (BC) remains uncertain. In this study, a total of 560 BC patients and 583 age, sex, and ethnicity-matched healthy controls from Northwest China were included. The polymorphisms were genotyped using Sequenom MassARRAY. The expression level of TIM-3 protein was detected by immunohistochemistry. We observed rs10053538 had a significantly increased risk of BC, comparing with the wild-type genotype even after Bonferroni correction. In addition, the rs4704853 G>A variants were more frequent among BC patients than the controls (GA + AA vs. GG: OR = 1.32, 95% CI = 1.03-1.69, P = 0.026); However, the significance was lost after Bonferroni correction (P = 0.078). Furthermore, rs10053538 was associated with lymph node metastasis. Age stratification revealed that among patients aged <49 years, those with the rs4704853 GA/AA genotype had a higher risk of BC; But there was no difference when Bonferroni correction was conducted. Immunohistochemical analysis showed that the expression of TIM-3 protein in the breast cancer tissues was higher in patients carrying the rs10053538 GT+TT genotype than those with GG genotype (P = 0.012). However, we failed to find any difference between BC patients and controls in any rs1036199 genetic model. These findings suggested that rs10053538 in TIM-3 might increase susceptibility to BC and promote the progression of BC in Chinese women.

Breast cancer: The translation of big genomic data to cancer precision medicine

Cancer is a complex genetic disease that develops from the accumulation of genomic alterations in which germline variations predispose individuals to cancer and somatic alterations initiate and trigger the progression of cancer. For the past 2 decades, genomic research has advanced remarkably, evolving from single-gene to whole-genome screening by using genome-wide association study and next-generation sequencing that contributes to big genomic data. International collaborative efforts have contributed to curating these data to identify clinically significant alterations that could be used in clinical settings. Focusing on breast cancer, the present review summarizes the identification of genomic alterations with high-throughput screening as well as the use of genomic information in clinical trials that match cancer patients to therapies, which further leads to cancer precision medicine. Furthermore, cancer screening and monitoring were enhanced greatly by the use of liquid biopsies. With the growing data complexity and size, there is much anticipation in exploiting deep machine learning and artificial intelligence to curate integrative "-omics" data to refine the current medical practice to be applied in the near future.