Genome wide screen identifies microsatellite markers associated with acute adverse effects following radiotherapy in cancer patients
BMC Medical Genetics volume 11, Article number: 123 (2010)
The response of normal tissues in cancer patients undergoing radiotherapy varies, possibly due to genetic differences underlying variation in radiosensitivity.
Cancer patients (n = 360) were selected retrospectively from the RadGenomics project. Adverse effects within 3 months of radiotherapy completion were graded using the National Cancer Institute Common Toxicity Criteria; high grade group were grade 3 or more (n = 180), low grade group were grade 1 or less (n = 180). Pooled genomic DNA (gDNA) (n = 90 from each group) was screened using 23,244 microsatellites. Markers with different inter-group frequencies (Fisher exact test P < 0.05) were analyzed using the remaining pooled gDNA. Silencing RNA treatment was performed in cultured normal human skin fibroblasts.
Forty-seven markers had positive association values; including one in the SEMA3A promoter region (P = 1.24 × 10-5). SEMA3A knockdown enhanced radiation resistance.
This study identified 47 putative radiosensitivity markers, and suggested a role for SEMA3A in radiosensitivity.
A principle determinant of the efficiency of tumor eradication following radiotherapy is the total radiation dose given to a patient. The radiation tolerance of important organs located at the margins of the radiotherapy target volume is a critical issue in determining dose thresholds. However, variation in the genetic background of individuals also contributes to the severity of radiation-related adverse events [1–7].
Andreassen et al. have summarized the results of many genetic association studies that used single nucleotide polymorphisms (SNPs) as genetic markers, and compared allele frequencies in radiosensitive and nonradiosensitive individuals . Most studies use a candidate gene approach; with genes selected based on ontology. In particular, these studies have focused on genes involved in processes including response to DNA damage, cell death, cell cycle control, oxidative stress, radiation-induced fibrogenesis, and endothelial cell damage.
Systematic microarray gene expression analyses [8–13] and in vitro functional screening using siRNA knockdown of gene expression , have been used to identify potential radiation susceptibility genes. Significant association has been found between the risk of early adverse skin reactions (EASRs) following radiotherapy, and SNP haplotypes associated with six of 137 candidate genes (CD44, MAD2L2, PTTG1, RAD9A, LIG3 and REV3L) . This has led to the development of a novel DNA chip-based technique to analyze haplotype markers in individual cancer patients [16–19].
Although positive associations between genetic markers and radiosensitivity have been found, the search for strongly associated genetic markers has been unrewarding , and this is partly due to inadequate understanding of the molecular pathology of adverse reactions induced by radiotherapy.
Microsatellites are useful mapping tools as they are abundant and interspersed throughout the human genome, similar to SNPs. Importantly though, microsatellite polymorphism generally exceeds that of single SNPs, even reaching the degree of polymorphism provided by SNP haplotypes . Thus, association analyses using a relatively small number of microsatellites should still have adequate statistical power relative to that provided by SNPs . This is illustrated by the identification of genes associated with rheumatoid arthritis , narcolepsy  and Behcet's disease  using genome-wide association studies based on microsatellites.
Hence, a genome-wide association study was performed to identify candidate genes that are strongly associated with radiosensitivity in humans. The screen analyzed data from 23,244 microsatellites in 360 cancer patients who had undergone radiotherapy and been graded for normal tissue adverse reactions. Forty-seven markers were identified as being of interest, with a role for the involvement of SEMA3A in radiosensitivity suggested.
Grading of patients with low and high grade radiosensitivity
Since 2001, the RadGenomics project has enrolled more than 3000 patients who have undergone radiotherapy. All patients provided written informed consent to participate in this study, which was approved by the Institutional Review Board at the National Institute of Radiological Sciences and by each collaborating institution. The acute adverse reactions of individual patients up to three months after completion of radiotherapy were graded according to the National Cancer Institute's Common Toxicity Criteria (NCICTC) version 2. We retrospectively selected 180 patients, who presented with a severe acute reaction of equal to or greater than grade 3, as a high-grade group (HGG). We also retrospectively selected 180 patients with less than or equal to grade 1 acute reaction on any end point, as a low-grade group (LGG). The assignment of patients to the LGG and HGG took into consideration their cancer type, age, gender, treatment type, and radiation dose (table 1).
Preparation of pooled DNA samples
Extraction of gDNA from whole blood was performed using an automatic nucleic acid isolator, NA3000S (Kurabo, Osaka, Japan) or with the QIAamp DNA blood kit (Qiagen, Hilden, Germany). The gDNA concentrations were measured in triplicate using a PicoGreen double-stranded DNA quantification kit (Invitrogen, Carlsbad, USA) and an SF600 microtiter plate reader (Corona Electric, Ibaraki, Japan). To reduce the amount of genotyping required, gDNA samples were pooled according to the method of Collins et al . Concentrations of individual gDNA samples were adjusted to 8 ng/μL. An equal volume of each of 90 gDNA samples from the HGG was combined to generate the first set of pooled gDNA and termed HGG-1. Similarly, 90 gDNA samples from the LGG were pooled and termed LGG-1. A second set of pooled gDNA samples was also prepared from 90 samples of the HGG and 90 samples of the LGG, and these were termed HGG-2 and LGG-2, respectively.
Analysis of microsatellite markers
All microsatellite markers and the methods for microsatellite analysis used in this study are described in Tamiya et al . The genomic location of the microsatellite markers was investigated using the UCSC Genome Browser http://genome.ucsc.edu/cgi-bin/hgGateway. PCR primers to amplify microsatellites were designed to anneal at 57°C, with forward primers having a 5' fluorescent label (6-FAM or HEX). PCR was performed using the GeneAmp PCR system 9700 (GE Healthcare, Amersham Place, UK) in 20 μL containing 48 ng pooled DNA, 0.5 U AmpliTaq DNA polymerase, reaction buffer containing 1.5 mM MgCl2 (GE Healthcare, Amersham Place, UK), 5 μM of each primer, and 0.25 mM of each dNTP in 96- or 384-well plates. PCR profile was as follows; 96°C for 5 min, 57°C for 1 min, 72°C for 1 min; 40 cycles of 96°C for 45 s, 57°C for 45 s, 72°C for 1 min. For microsatellite typing of individual samples, PCR was performed in 12 μL containing 2 ng DNA, 0.25 U AmpliTaq Gold DNA polymerase (GE Healthcare, Amersham Place, UK), reaction buffer containing 1.5 mM MgCl2, 5 μM of each primer, and 0.2 mM of each dNTP in 96- or 384-well plates and amplified as above. PCR products were denatured in Hi-Di formamide (GE Healthcare, Amersham Place, UK) at 95°C for 3 min and separated by capillary electrophoresis using an ABI Prism 3700 Genetic Analyzer and ROX size standards (GE Healthcare, Amersham Place, UK). Analysis of fragment size and electrophoretograms was performed using GeneScan and Genotyper software (GE Healthcare, Little Chalfont, UK).
Allele frequencies in pooled DNA were estimated from the height of peaks in the electrophoretogram . Association of microsatellites with radiosensitivity was assessed using Fisher exact test and 2 × 2 contingency tables for each allele. The lowest P value for any allele of a particular microsatellite was used in analysis for that microsatellite marker and significance was set at 0.05. To account for multiple testing across the microsatellite markers, the P values of the second screening were adjusted using the false discovery rate (FDR) controlling procedure of Benjamini and Hochberg . Association between particular alleles of a microsatellite and grade of radiosensitivity was performed using the Cochran-Armitage test for trend .
Cell culture conditions
Normal human skin fibroblast NB1RGB cells were obtained from Riken Cell Bank (Tsukuba, Japan) and maintained in Eagle's minimum essential medium (Nissui, Tokyo, Japan) supplemented with 10% fetal bovine serum (FBS) and nonessential amino acids under a humidified atmosphere of 5% CO2 at 37°C.
siRNA treatment of human fibroblast cultures
Two different siRNAs (1: 135598, 2: 135597, Applied Biosystems/Ambion, Austin, USA) designed for the human SEMA3A gene were used to treat cells by reverse transfection. Cultured cells were harvested by incubation with 0.05% trypsin, 0.53 mM EDTA in phosphate-buffered saline (trypsin-EDTA/PBS) for 5 min at room temperature, followed by inactivation of the trypsin by adding complete culture medium. The number of cells with a diameter of 12 μm was measured using the Z1 Coulter Particle Counter (Beckman Coulter, Brea, USA). Transfection complexes were prepared in 266.6 μl Opti-MEM serum free medium by mixing 5.4 μL of siPORT NeoFX Transfection Reagent (Applied Biosystems/Ambion, Austin, USA) for 10 min prior to adding 2.7 μl of 10 μM siRNA (Applied Biosystems/Ambion, Austin, USA) for 10 min. The cell suspension containing 27,000 cells in 725.3 μl was added to the transfection complexes and the mixture plated onto a 35 mm plastic dish. Cells were maintained under a humidified atmosphere of 5% CO2 at 37°C for 24 hours. The mock control cells were treated similarly, except the transfection complexes were prepared without siRNA. The transfection medium was then changed to complete culture medium, and cells were maintained for the indicated time.
The post-transfection complete culture medium was recovered at the indicated time, centrifuged at 12,000 × g for 5 min and the supernatant transferred to a microfuge tube and stored at -30°C until use. Fifteen μL aliquots of supernatant were mixed with 5 μL of loading dye, incubated at 80°C for 10 min and loaded onto NuPAGE 4-12% polyacrylamide Bis-Tris gradient gels (Invitrogen, Carlsbad, USA) under denaturing conditions. Proteins in the gels were then transferred onto PVDF membranes by the iBlot Gel Transfer Device (Invitrogen, Carlsbad, USA). The membranes were probed with primary rabbit polyclonal antibody SEMA (H-300) against amino acids 103-402 of the human protein from SEMA3A (Santa Cruz Biotechnology, Santa Cruz, USA). A secondary donkey anti-rabbit IgG conjugated to horseradish peroxidase (GE Healthcare, Amersham Place, UK) was used to detect protein bands using an ECL Advance Western Blotting Detection Kit (GE Healthcare, Little Chalfont, UK).
Radiosensitivity assay of siRNA-treated human skin fibroblast cultures
Cells cultured for 48 hours following replacement of the siRNA transfection medium, were harvested by incubation with trypsin-EDTA/PBS for 5 min at room temperature. Trypsin was then inactivated by adding complete culture medium, and the number of cells with a diameter of more than 12 μm was counted using a Z1 Coulter Particle Counter. The appropriate number of cells (1,000 cells for 0 Gy, 1,500 cells for 1 Gy, 3,000 cells for 2 Gy, 6,000 cells for 3 Gy and 12,000 cells for 4 Gy) were then irradiated at room temperature with 200 kV X-rays (20 mA) with 0.5 mm aluminum and 0.5 mm copper filters. Immediately after irradiation, cells were plated onto 100 mm (0 and 1 Gy) or 150 mm (2-4 Gy) plastic dishes, and cultured under a humidified atmosphere of 5% CO2 at 37°C for 2 weeks. Cells were washed with PBS, then fixed in 100% methanol for a few minutes at room temperature. After removing the methanol, the cells were dried for 30 minutes and then stained with 3% Giemsa solution for 2 hours. Colonies consisting of more than 50 cells were scored as survivors. Experiments were performed with triplicate plating of cells. Relative colony survival as a function of irradiated dose was tested using the linear-quadratic model .
Genome-wide association study
The entire human genome was screened using 23,244 microsatellite markers and one set of the pooled DNA from each group (HGG-1 and LGG-1). This analysis identified 3,052 markers with allele/s that showed significantly different estimated frequency (P < 0.05) between the two groups (see figure 1A). These markers were further analyzed in the second round screening with the remaining set of pooled DNA samples from each group (HGG-2 and LGG-2). A total of 101 markers had allele/s with significantly different estimated frequency (P < 0.05) between the two groups, and had similar peak association patterns in the data from the first and second DNA pools (see figure 1B). The FDR of the markers analyzed in the second round screening was then estimated to correct for the effect of multiple comparisons. After correction for the multiple comparisons performed, 47 autosomal markers had a FDR < 0.05, and so were significantly associated with radiosensitivity in this study. These markers are summarized in table 2.
Individual patient typing of selected markers
The genomic location of the 47 positive markers was determined using the UCSC Genome Browser. Four markers were within 15 kb upstream of the transcription start site for the nearest gene (table 3). To confirm the associations observed in the pooled DNA analysis, these four markers were individually typed in all 360 patients. The particular alleles of D1S0288i and D7S0338i found to be statistically significantly associated with radiosensitivity in the experiments using pooled DNA, were also found to be so when experiments using individual DNA samples were used (table 3). One of the positive markers, D7S0338i, was selected for further study as it had the strongest association with radiosensitivity.
Detailed association analysis of the D7S0338i marker
Allelic distribution of the marker D7S0338i is represented in additional file 1. Eleven alleles were identified in D7S0338i, with three major alleles (292 bp, 294 bp and 296 bp) accounting for more than 90% of the alleles in the study population (see table 4 and additional file 1). Among them, two alleles (292 bp and 296 bp) had an overall P value < 0.05. Association analysis revealed the 292/292 genotype was associated with risk reduction (odds ratio: 0.490, 95% CI: 0.271-0.885), the 294/294 genotype was associated with intermediate risk (odds ratio: 0.794, 95% CI: 0.321-1.966), and the 296/296 genotype was associated with an increased risk (odds ratio: 2.861, 95% CI: 1.417-5.774) (see table 5). Interestingly, the heterozygote individuals carrying the risk reducing and risk increasing alleles (292/296) had intermediate risk (odds ratio: 1.152, 95% CI: 0.727-1.824). This suggests the inter-allele relationship may not behave in the usual dominant-recessive manner. A linear trend was identified between six combinations of alleles for the marker and radiosensitivity (P = 7.228 × 10-4) (see table 5).
Impact on cellular radiosensitivity of the D7S0338i associated gene SEMA3A
The D7S0338i marker is 1500 bp upstream of the transcription start site of the SEMA3A gene (table 3). The linkage disequilibrium block around this marker did not extend to neighbouring genes (see additional file 2). To assess possible functional impacts of the SEMA3A gene on cellular radiosensitivity (see additional file 3), a model system using in vitro irradiation of cultured fibroblasts derived from human skin was employed. Specific SEMA3A siRNA treatment of these cells reduced expression of the SEMA3A gene in normal human skin fibroblasts (see additional file 4). When these cells were irradiated with X-rays (200 keV), the colony survival assay of the irradiated cells showed enhancement of cellular radiation resistance in the SEMA3A siRNA-treated cells (see figure 2).
A human genome-wide microsatellite association study was performed in cancer patients who showed radiation-induced adverse reactions. After correction for multiple comparisons, this study identified 47 autosomal markers with a FDR < 0.05. One of these markers is within the proximal promoter region of the SEMA3A gene on chromosome 7. Knockdown of SEMA3A expression in a normal human skin fibroblast culture caused a significant change in the radiosensitivity of these cells.
The SEMA3A gene has not been previously described as having a role in radiosensitivity. SEMA3A encodes a secreted protein (semaphorin-3A), which is involved in a wide range of functional processes including regulation of axon guidance, cell survival, motility, immune responses and angiogenesis [28–49] (see additional file 3). This diversity of these roles provides many possible mechanisms for its involvement in radiosensitivity. Semaphorin-3A is also a competitor of the angiogenic growth factor coded for by the VEGF gene, as both bind to the same transmembrane receptor [29, 31, 46]. As VEGF expression is directly correlated with radiosensitivity , its competitor, semaphorin-3A, may also be associated with radiosensitivity.
The potential role of the D7S0338i marker in radiosensitivity is interesting. The microsatellite is located 1500 bp upstream from the transcription start site of SEMA3A, and is a GA dinucleotide repeat (see additional file 5). Two other polymorphic repetitive sequences also occur between the D7S0338i marker and the SEMA3A transcription start site (additional file 6). These three repetitive sequences are in a low nucleosome occupancy region (data not shown). Since nucleosomes play a major role in generating the higher order structure of chromatin that regulates gene expression , these sequences may affect the activity of the SEMA3A promoter. A study into the activity of the SEMA3A promoter may provide information on the functional impacts of the D7S0338i marker polymorphism, especially mechanisms underlying the phenotypes associated with various alleles (table 5).
A major limitation of association studies on rare phenotypes such as the severe (equal to or greater than grade 3), acute, adverse reactions induced by radiotherapy, is the ability to enroll sufficient numbers of patients to provide adequate statistical power [7, 52]. The reproducibility of any association identified must also be replicated , further increasing the required patient number. Hence in this study, cancer patients were selected with differing severe, acute, adverse reaction endpoints and various cancer types. Identical numbers of control patients were selected who did not develop severe, adverse reactions on any endpoint. The clinical characteristics and therapeutic protocols used in the control and subject patients were also similar (table 1). While this suggests inherent or genetic differences between patients are the cause of the variations in severity in patient reactions, the involvement of SEMA3A requires further validation using large numbers of patients with a unique cancer type.
A possible implication for this study's finding that semaphorin-3A may be involved in radiosensitivity, is the identification of a potential new agent for the treatment of radiotherapy-induced damage. SM-216289 (xanthofulvin) was originally isolated from the fermentation broth of a fungal strain, Penicillium sp. SPF-3059, and is a natural inhibitor of Semaphorin-3A . SM-216289 abolished the growth cone collapse of dorsal root ganglion neurons that was induced by Semaphorin-3A in vitro and in vivo, possibly through direct interference of the receptor-ligand association . Local administration of SM-216289 in the adult rat model of spinal cord injury, substantially enhanced functional recovery of injured axons, with decreases in apoptotic cell number and marked enhancement of angiogenesis . Therefore, locally administered SM-216289 may aid functional recovery in radiotherapy-induced injury.
Increasingly, phenotypic differences have been shown to be caused by diverse genetic variations. While SNP and repetitive DNA polymorphisms have been shown to be disease-associated, other disease-associated changes include copy number variation  and transgenerational epigenetic modification of the genome . Thus, diagnostic testing should be considered to identify highly radiosensitive cancer patients through the detection of genetic variants in individual patients. The rapidly developing next-generation genome sequencing technology  may be the most suitable one for this purpose.
A total of 47 putative markers for individual radiosensitivity were identified using a genome-wide screen based on microsatellite markers. One of these markers is in the proximal promoter region of SEMA3A, with knockdown of this gene using siRNA supporting its potential role in radiosensitivity.
Bentzen SM, Overgaard J: Patient-to-patient variability in the expression of radiation-induced normal tissue injury. Semin Radiat Oncol. 1994, 4: 68-80. 10.1016/S1053-4296(05)80034-7.
Iwakawa M, Noda S, Yamada S, Yamamoto N, Miyazawa Y, Yamazaki H, Kawakami Y, Matsui Y, Tsujii H, Mizoe J, Oda E, Fukunaga Y, Imai T: Analysis of non-genetic risk factors for adverse skin reactions to radiotherapy among 284 breast cancer patients. Breast Cancer. 2006, 13: 300-307. 10.2325/jbcs.13.300.
Andreassen CN, Alsner J, Overgaard J: Does variability in normal tissue reactions after radiotherapy have a genetic basis - where and how to look for it?. Radiother Oncol. 2002, 64: 131-140. 10.1016/S0167-8140(02)00154-8.
Fernet M, Hall J: Genetic biomarkers of therapeutic radiation sensitivity. DNA Repair. 2004, 3: 1237-1243. 10.1016/j.dnarep.2004.03.019.
Gatti RA: The inherited basis of human radiosensitivity. Acta Oncol. 2001, 40: 702-711. 10.1080/02841860152619115.
Barnett GC, West CM, Dunning AM, Elliott RM, Coles CE, Pharoah PD, Burnet NG: Normal tissue reactions to radiotherapy: towards tailoring treatment dose by genotype. Nat Rev Cancer. 2009, 9: 134-142. 10.1038/nrc2587.
Andreassen CN, Alsner J: Genetic variants and normal tissue toxicity after radiotherapy: a systematic review. Radiother Oncol. 2009, 92: 299-309. 10.1016/j.radonc.2009.06.015.
Iwakawa M, Noda S, Ohta T, Ohira C, Lee R, Goto M, Wakabayashi M, Matsui Y, Harada Y, Imai T: Different radiation susceptibility among five strains of mice detected by a skin reaction. J Radiat Res. 2003, 44: 7-13. 10.1269/jrr.44.7.
Ohta T, Iwakawa M, Oohira C, Noda S, Minfu Y, Goto M, Tanaka H, Harada Y, Imai T: Fractionated irradiation augments inter-strain variation of skin reactions among three strains of mice. J Radiat Res. 2004, 45: 515-519. 10.1269/jrr.45.515.
Iwakawa M, Noda S, Ohta T, Tanaka H, Tsuji A, Ishikawa A, Imai T: Strain dependent differences in a histological study of CD44 and collagen fibers with expression analysis of inflammatory response-related genes in irradiated murine lung. J Radiat Res. 2004, 45: 423-433. 10.1269/jrr.45.423.
Noda S, Iwakawa M, Ohta T, Iwata M, Yang M, Goto M, Tanaka H, Harada Y, Imai T: Inter-strain variance in late phase of erythematous reaction or leg contracture after local irradiation among three strains of mice. Cancer Detect Prev. 2005, 29: 376-382. 10.1016/j.cdp.2005.06.005.
Ishikawa K, Koyama-Saegusa K, Otsuka Y, Ishikawa A, Kawai S, Yasuda K, Suga T, Michikawa Y, Suzuki M, Iwakawa M, Imai T: Gene expression profile changes correlating with radioresistance in human cell lines. Int J Radiat Oncol Biol Phys. 2006, 65: 234-245. 10.1016/j.ijrobp.2005.12.048.
Ban S, Ishikawa K, Kawai S, Koyama-Saegusa K, Ishikawa A, Shimada Y, Inazawa J, Imai T: Potential in a single cancer cell to produce heterogeneous morphology, radiosensitivity and gene expression. J Radiat Res. 2005, 46: 43-50. 10.1269/jrr.46.43.
Tsuji AB, Sudo H, Sugyo A, Otsuki M, Miyagishi M, Taira K, Imai T, Harada YN: A fast, simple method for screening radiation susceptibility genes by RNA interference. Biochem Biophys Res Commun. 2005, 333: 1370-1377. 10.1016/j.bbrc.2005.06.047.
Suga T, Ishikawa A, Kohda M, Otsuka Y, Yamada S, Yamamoto N, Shibamoto Y, Ogawa Y, Nomura K, Sho K, Omura M, Sekiguchi K, Kikuchi Y, Michikawa Y, Noda S, Sagara M, Ohashi J, Yoshinaga S, Mizoe J, Tsujii H, Iwakawa M, Imai T: Haplotype-based analysis of genes associated with risk of adverse skin reactions after radiotherapy in breast cancer patients. Int J Radiat Oncol Biol Phys. 2007, 69: 685-693. 10.1016/j.ijrobp.2007.06.021.
Michikawa Y, Fujimoto K, Kinoshita K, Kawai S, Sugahara K, Suga T, Otsuka Y, Fujiwara K, Iwakawa M, Imai T: Reliable and fast allele-specific extension of 3'-LNA modified oligonucleotides covalently immobilized on a plastic base, combined with biotin-dUTP mediated optical detection. Anal Sci. 2006, 22: 1537-1545. 10.2116/analsci.22.1537.
Michikawa Y, Suga T, Ohtsuka Y, Matsumoto I, Ishikawa A, Ishikawa K, Iwakawa M, Imai T: Visible genotype sensor array. Sensors. 2008, 8: 2722-2735. 10.3390/s8042722.
Michikawa Y, Suga T, Ishikawa A, Ohtsuka Y, Iwakawa M, Imai T: Visible haplotype-tag SNP typing array device for human radiation sensitivity-associated genes. Oligonucleotide Array Sequence Analysis. Edited by: Moretti MK, Rizzo LJ. 2008, New York: Nova Publishers, 3-14.
Michikawa Y, Sugahara K, Suga T, Ohtsuka Y, Ishikawa K, Ishikawa A, Shiomi N, Shiomi T, Iwakawa M, Imai T: In-gel multiple displacement amplification of long DNA fragments diluted to the single molecule level. Anal Biochem. 2008, 383: 151-158. 10.1016/j.ab.2008.08.011.
Payseur BA, Jing PA: Genomewide comparison of population structure at STRPs and nearby SNPs in humans. Mol Biol Evol. 2009, 26: 1369-1377. 10.1093/molbev/msp052.
Tamiya G, Shinya M, Imanishi T, Ikuta T, Makino S, Okamoto K, Furugaki K, Matsumoto T, Mano S, Ando S, Nozaki Y, Yukawa W, Nakashige R, Yamaguchi D, Ishibashi H, Yonekura M, Nakami Y, Takayama S, Endo T, Saruwatari T, Yagura M, Yoshikawa Y, Fujimoto K, Oka A, Chiku S, Linsen SE, Giphart MJ, Kulski JK, Fukazawa T, Hashimoto H, Kimura M, Hoshina Y, Suzuki Y, Hotta T, Mochida J, Minezaki T, Komai K, Shiozawa S, Taniguchi A, Yamanaka H, Kamatani N, Gojobori T, Bahram S, Inoko H: Whole genome association study of rheumatoid arthritis using 27039 microsatellites. Hum Mol Genet. 2005, 14: 2305-2321. 10.1093/hmg/ddi234.
Kawashima M, Tamiya G, Oka A, Hohjoh H, Juji T, Ebisawa T, Honda Y, Inoko H, Tokunaga K: Genomewide Association Analysis of Human Narcolepsy and a New Resistance Gene. Am J Hum Genet. 2006, 79: 252-263. 10.1086/505539.
Meguro A, Ota M, Katsuyama Y, Oka A, Ohno S, Inoko H, Mizuki N: Association of the toll-like receptor 4 gene polymorphisms with Behçet's disease. Ann Rheum Dis. 2008, 67: 725-727. 10.1136/ard.2007.079871.
Collins HE, Li H, Inda SE, Anderson J, Laiho K, Tuomilehto J, Seldin MF: A simple and accurate method for determination of microsatellite total allele content differences between DNA pools. Hum Genet. 2000, 106: 218-226. 10.1007/s004390051031.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B. 1995, 57: 289-300.
Armitage P: Tests for linear trends in proportions and frequencies. Biometrics. 1955, 11: 375-386. 10.2307/3001775.
Tucker SL: Tests for the fit of the linear-quadratic model to radiation isoeffect data. Int J Radiat Oncol Biol Phys. 1984, 10: 1933-1939. 10.1016/0360-3016(84)90274-8.
Vacca A, Scavelli C, Serini G, Di Pietro G, Cirulli T, Merchionne F, Ribatti D, Bussolino F, Guidolin D, Piaggio G, Bacigalupo A, Dammacco F: Loss of inhibitory semaphorin 3A (SEMA3A) autocrine loops in bone marrow endothelial cells of patients with multiple myeloma. Blood. 2006, 108: 1661-1667. 10.1182/blood-2006-04-014563.
Narazaki M, Segarra M, Tosato G: Sulfated polysaccharides identified as inducers of neuropilin-1 internalization and functional inhibition of VEGF165 and semaphorin3A. Blood. 2008, 111: 4126-4136. 10.1182/blood-2007-09-112474.
Guttmann-Raviv N, Shraga-Heled N, Varshavsky A, Guimaraes-Sternberg C, Kessler O, Neufeld G: Semaphorin-3A and semaphorin-3F work together to repel endothelial cells and to inhibit their survival by induction of apoptosis. J Biol Chem. 2007, 282: 26294-26305. 10.1074/jbc.M609711200.
Serini G, Valdembri D, Zanivan S, Morterra G, Burkhardt C, Caccavari F, Zammataro L, Primo L, Tamagnone L, Logan M, Tessier-Lavigne M, Taniguchi M, Püschel AW, Bussolino F: Class 3 semaphorins control vascular morphogenesis by inhibiting integrin function. Nature. 2003, 424: 391-397. 10.1038/nature01784.
Barresi V, Vitarelli E, Cerasoli S: Semaphorin3A immunohistochemical expression in human meningiomas: correlation with the microvessel density. Virchows Arch. 2009, 454: 563-571. 10.1007/s00428-009-0757-3.
Moretti S, Procopio A, Lazzarini R, Rippo MR, Testa R, Marra M, Tamagnone L, Catalano A: Semaphorin3A signaling controls Fas (CD95)-mediated apoptosis by promoting Fas translocation into lipid rafts. Blood. 2008, 111: 2290-2299. 10.1182/blood-2007-06-096529.
Schmidt EF, Strittmatter SM: The CRMP family of proteins and their role in Sema3A signaling. Adv Exp Med Biol. 2007, 600: 1-11. full_text.
Tannemaat MR, Korecka J, Ehlert EM, Mason MR, van Duinen SG, Boer GJ, Malessy MJ, Verhaagen J: Human neuroma contains increased levels of semaphorin 3A, which surrounds nerve fibers and reduces neurite extension in vitro. J Neurosci. 2007, 27: 14260-14264. 10.1523/JNEUROSCI.4571-07.2007.
Herman JG, Meadows GG: Increased class 3 semaphorin expression modulates the invasive and adhesive properties of prostate cancer cells. Int J Oncol. 2007, 30: 1231-1238.
Kurschat P, Bielenberg D, Rossignol-Tallandier M, Stahl A, Klagsbrun M: Neuron restrictive silencer factor NRSF/REST is a transcriptional repressor of neuropilin-1 and diminishes the ability of semaphorin 3A to inhibit keratinocyte migration. J Biol Chem. 2006, 281: 2721-2729. 10.1074/jbc.M507860200.
Lepelletier Y, Smaniotto S, Hadj-Slimane R, Villa-Verde DM, Nogueira AC, Dardenne M, Hermine O, Savino W: Control of human thymocyte migration by neuropilin-1/semaphorin-3A-mediated interactions. Proc Natl Acad Sci USA. 2007, 104: 5545-5550. 10.1073/pnas.0700705104.
Bachelder RE, Lipscomb EA, Lin X, Wendt MA, Chadborn NH, Eickholt BJ, Mercurio AM: Competing autocrine pathways involving alternative neuropilin-1 ligands regulate chemotaxis of carcinoma cells. Cancer Res. 2003, 63: 5230-5233.
Ko JA, Morishige N, Yanai R, Nishida T: Up-regulation of semaphorin 3A in human corneal fibroblasts by epidermal growth factor released from cocultured human corneal epithelial cells. Biochem Biophys Res Commun. 2008, 377: 104-108. 10.1016/j.bbrc.2008.09.085.
Kashiwagi H, Shiraga M, Kato H, Kamae T, Yamamoto N, Tadokoro S, Kurata Y, Tomiyama Y, Kanakura Y: Negative regulation of platelet function by a secreted cell repulsive protein, semaphorin 3A. Blood. 2005, 106: 913-921. 10.1182/blood-2004-10-4092.
Lepelletier Y, Moura IC, Hadj-Slimane R, Renand A, Fiorentino S, Baude C, Shirvan A, Barzilai A, Hermine O: Immunosuppressive role of semaphorin-3A on T cell proliferation is mediated by inhibition of actin cytoskeleton reorganization. Eur J Immunol. 2006, 36: 1782-1793. 10.1002/eji.200535601.
Marzioni D, Tamagnone L, Capparuccia L, Marchini C, Amici A, Todros T, Bischof P, Neidhart S, Grenningloh G, Castellucci M: Restricted innervation of uterus and placenta during pregnancy: evidence for a role of the repelling signal semaphorin 3A. Dev Dyn. 2004, 231: 839-848. 10.1002/dvdy.20178.
Eastwood SL, Law AJ, Everall IP, Harrison PJ: The axonal chemorepellant semaphorin 3A is increased in the cerebellum in schizophrenia and may contribute to its synaptic pathology. Mol Psychiatry. 2003, 8: 148-155. 10.1038/sj.mp.4001233.
Potiron VA, Roche J, Drabkin HA: Semaphorins and their receptors in lung cancer. Cancer Lett. 2009, 273: 1-14. 10.1016/j.canlet.2008.05.032.
Catalano A, Caprari P, Rodilossi S, Betta P, Castellucci M, Casazza A, Tamagnone L, Procopio A: Cross-talk between vascular endothelial growth factor and semaphorin-3A pathway in the regulation of normal and malignant mesothelial cell proliferation. FASEB J. 2004, 18: 358-363.
Rieger J, Wick W, Weller M: Human malignant glioma cells express semaphorins and their receptors, neuropilins and plexins. Glia. 2003, 42: 379-389. 10.1002/glia.10210.
Müller MW, Giese NA, Swiercz JM, Ceyhan GO, Esposito I, Hinz U, Büchler P, Giese T, Büchler MW, Offermanns S, Friess H: Association of axon guidance factor semaphorin 3A with poor outcome in pancreatic cancer. Int J Cancer. 2007, 121: 2421-2433. 10.1002/ijc.22949.
Kigel B, Varshavsky A, Kessler O, Neufeld G: Successful inhibition of tumor development by specific class-3 semaphorins is associated with expression of appropriate semaphorin receptors by tumor cells. PLoS One. 2008, 3: e3287-10.1371/journal.pone.0003287.
Brieger J, Kattwinkel J, Berres M, Gosepath J, Mann WJ: Impact of vascular endothelial growth factor release on radiation resistance. Oncol Rep. 2007, 18: 1597-1601.
Horowitz RA, Agard DA, Sedat JW, Woodcock C: The three-dimensional architecture of chromatin in situ: electron tomography reveals fibers composed of a continuously variable zig-zag nucleosomal ribbon. J Cell Biol. 1994, 125: 1-10. 10.1083/jcb.125.1.1.
Spencer CC, Su Z, Donnelly P, Marchini J: Designing genome-wide association studies: sample size, power, imputation, and the choice of genotyping chip. PLoS Genet. 2009, 5: e1000477-10.1371/journal.pgen.1000477.
NCI-NHGRI working group on replication in association studies: Replicating genotype-phenotype associations. Nature. 2007, 447: 655-660. 10.1038/447655a.
Kikuchi K, Kishino A, Konishi O, Kumagai K, Hosotani N, Saji I, Nakayama C, Kimura T: In vitro and in vivo characterization of a novel semaphorin 3A inhibitor, SM-216289 or xanthofulvin. J Biol Chem. 2003, 278: 42985-42991. 10.1074/jbc.M302395200.
Kaneko S, Iwanami A, Nakamura M, Kishino A, Kikuchi K, Shibata S, Okano HJ, Ikegami T, Moriya A, Konishi O, Nakayama C, Kumagai K, Kimura T, Sato Y, Goshima Y, Taniguchi M, Ito M, He Z, Toyama Y, Okano H: A selective Sema3A inhibitor enhances regenerative responses and functional recovery of the injured spinal cord. Nat Med. 2006, 12: 1380-1389. 10.1038/nm1505.
Diskin SJ, Hou C, Glessner JT, Attiyeh EF, Laudenslager M, Bosse K, Cole K, Mossé YP, Wood A, Lynch JE, Pecor K, Diamond M, Winter C, Wang K, Kim C, Geiger EA, McGrady PW, Blakemore AI, London WB, Shaikh TH, Bradfield J, Grant SF, Li H, Devoto M, Rappaport ER, Hakonarson H, Maris JM: Copy number variation at 1q21.1 associated with neuroblastoma. Nature. 2009, 459: 987-991. 10.1038/nature08035.
Johannes F, Porcher E, Teixeira FK, Saliba-Colombani V, Simon M, Agier N, Bulski A, Albuisson J, Heredia F, Audigier P, Bouchez D, Dillmann C, Guerche P, Hospital F, Colot V: Assessing the impact of transgenerational epigenetic variation on complex traits. PLoS Genet. 2009, 5: e1000530-10.1371/journal.pgen.1000530.
Ansorge WJ: Next-generation DNA sequencing techniques. Nat Biotechnol. 2009, 25: 195-203.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/11/123/prepub
The authors thank the cancer patients for participating in this study. The authors also thank Ms. Masayo Terada for assistance with manuscript preparation. YM thanks Mrs Katsuko Noshiro of NIRS and Ms Erika Matsushita of TU for technical assistance.
The authors declare that they have no competing interests.
YM participated in study conception and design, conducted experiments and statistical analyses, participated in data interpretation, and drafted the manuscript. TS conducted experiments and statistical analysis. AI conducted statistical analyses and participated in data interpretation. HH, AO and HI conducted statistical analyses, participated in data interpretation, and revised the article critically for important intellectual content. MI and TI participated in study conception and design, conducted statistical analyses, participated in data interpretation, and revised the article critically for important intellectual content. All authors read and approved the final manuscript
Electronic supplementary material
Additional file 1: Allele frequencies for the D7S0338i marker. (A) First screening samples. (B) Second screening samples. (C) Overall results. The number of patients with the indicated allele size is plotted. Reproducible allele frequency differences were observed between the groups. White bar: LGG. Black bar: HGG. Significant differences (P < 0.05 by Fisher's exact test based on 2 × 2 contingency tables) are indicated by the vertical lines. (TIFF 1023 KB)
Additional file 2: Linkage disequilibrium map around the D7S0338i marker. Diagram showing the region of linkage disequilibrium (D') surrounding the SEMA3A gene and neighbouring genes. An increase in the degree of linkage disequilibrium between the SEMA3A gene and its neighbours is represented by an increase in the intensity of red color. Thus, the D7S0338i marker is in the linkage disequilibrium block that spans the transcriptional regulatory region of the SEMA3A gene. The D' data was calculated from JPT and CHB HapMap data. Positions are NCBI build 36 coordinates. Gene region and direction of transcription are indicated by an arrow. (TIFF 3 MB)
Additional file 4: Reduction in the amount of semaphorin-3A protein following siRNA treatment. The amount of semaphorin-3A in the culture medium of cells treated by the siRNA was quantified by Western blotting. The protein amounts relative to that of mock control cells are plotted. (TIFF 2 MB)
Additional file 5: Sequence of the D7S0338i marker. This marker is the GA dinucleotide repeat sequence indicated by the red color. The primer sequences for PCR amplification are indicated by the arrows. (TIFF 3 MB)
Additional file 6: Location of the D7S0338i marker in the promoter region of the SEMA3Agene. The GA dinucleotide repeat sequence is indicated by an open box. Two other polymorphic repetitive sequences (closed box) exist in the region between the D7S0338i marker and the transcription start site (+1). (TIFF 3 MB)
About this article
Cite this article
Michikawa, Y., Suga, T., Ishikawa, A. et al. Genome wide screen identifies microsatellite markers associated with acute adverse effects following radiotherapy in cancer patients. BMC Med Genet 11, 123 (2010). https://doi.org/10.1186/1471-2350-11-123