Using Prior Information from the Medical Literature in GWAS of Oral Cancer Identifies Novel Susceptibility Variant on Chromosome 4-the AdAPT Method

Background: Genome-wide association studies (GWAS) require large sample sizes to obtain adequate statistical power, but it may be possible to increase the power by incorporating complementary data. In this study we investigated the feasibility of automatically retrieving information from the medical literature and leveraging this information in GWAS. Methods: We developed a method that searches through PubMed abstracts for pre-assigned keywords and key concepts, and uses this information to assign prior probabilities of association for each single nucleotide polymorphism (SNP) with the phenotype of interest - the Adjusting Association Priors with Text (AdAPT) method. Association results from a GWAS can subsequently be ranked in the context of these priors using the Bayes False Discovery Probability (BFDP) framework. We initially tested AdAPT by comparing rankings of known susceptibility alleles in a previous lung cancer GWAS, and subsequently applied it in a two-phase GWAS of oral cancer. Results: Known lung cancer susceptibility SNPs were consistently ranked higher by AdAPT BFDPs than by p-values. In the oral cancer GWAS, we sought to replicate the top five SNPs as ranked by AdAPT BFDPs, of which rs991316, located in the ADH gene region of 4q23, displayed a statistically significant association with oral cancer risk in the replication phase (per-rare-allele log additive p-value [p(trend)] = 2.5 x 10(-3)). The combined OR for having one additional rare allele was 0.83 (95% CI: 0.76-0.90), and this association was independent of previously identified susceptibility SNPs that are associated with overall UADT cancer in this gene region. We also investigated if rs991316 was associated with other cancers of the upper aerodigestive tract (UADT), but no additional association signal was found. Conclusion: This study highlights the potential utility of systematically incorporating prior knowledge from the medical literature in genome-wide analyses using the AdAPT methodology. AdAPT is available online (url: http://services.gate.ac.uk/lld/gwas/service/config).

Higher incidence of death in multi-vessel coronary artery disease patients associated with polymorphisms in chromosome 9p21

Background: We investigated whether 9p21 polymorphisms are associated with cardiovascular events in a group of 611 patients enrolled in the Medical, Angioplasty or Surgery Study II (MASS II), a randomized trial comparing treatments for patients with coronary artery disease (CAD) and preserved left ventricular function. Methods: The participants of the MASS II were genotyped for 9p21 polymorphisms (rs10757274, rs2383206, rs10757278 and rs1333049). Survival curves were calculated with the Kaplan-Meier method and compared with the log-rank statistic. We assessed the relationship between baseline variables and the composite end-point of death, death from cardiac causes and myocardial infarction using a Cox proportional hazards survival model. Results: We observed significant differences between patients within each polymorphism genotype group for baseline characteristics. The frequency of diabetes was lower in patients carrying GG genotype for rs10757274, rs2383206 and rs10757278 (29.4%, 32.8%, 32.0%) compared to patients carrying AA or AG genotypes (49.1% and 39.2%, p = 0.01; 52.4% and 40.1%, p = 0.01; 47.8% and 37.9%, p = 0.04; respectively). Significant differences in genotype frequencies between double and triple vessel disease patients were observed for the rs10757274, rs10757278 and rs1333049. Finally, there was a higher incidence of overall mortality in patients with the GG genotype for rs2383206 compared to patients with AA and AG genotypes (19.5%, 11.9%, 11.0%, respectively; p = 0.04). Moreover, the rs2383206 was still significantly associated with a 1.75-fold increased risk of overall mortality (p = 0.02) even after adjustment of a Cox multivariate model for age, previous myocardial infarction, diabetes, smoking and type of coronary anatomy. Conclusions: Our data are in accordance to previous evidence that chromosome 9p21 genetic variation may constitute a genetic modulator in the cardiovascular system in different scenarios. In patients with established CAD, we observed an association between the rs2383206 and higher incidence of overall mortality and death from cardiac causes in patients with multi-vessel CAD.

Post-receptor IGF1 insensitivity restricted to the MAPK pathway in a Silver-Russell syndrome patient with hypomethylation at the imprinting control region on chromosome 11

Background: Hypomethylation of the paternal imprinting center region 1 (ICR1) is the most frequent molecular cause of Silver-Russell syndrome (SRS). Clinical evidence suggests that patients with this epimutation have mild IGF1 insensitivity. Objective: To assess in vitro IGF1 action in fibroblast culture from a patient with SRS and IGF1 insensitivity. Methods: Fibroblast cultures from one patient with SRS due to ICR1 demethylation and controls were established. The SRS patient has severe growth failure, elevated IGF1 level, and poor growth rate during human recombinant GH treatment. IGF1 action was assessed by cell proliferation, AKT, and p42/44-MAPK phosphorylation. Gene expression was determined by real-time PCR. Results: Despite normal IGF1R sequence and expression, fibroblast proliferation induced by IGF1 was 50% lower in SRS fibroblasts in comparison with controls. IGF1 and insulin promoted a p42/44-MAPK activation in SRS fibroblasts 40 and 36%, respectively, lower than that in control fibroblasts. On the other hand, p42/44-MAPK activation induced by EGF stimulation was only slightly reduced (75% in SRS fibroblasts in comparison with control), suggesting a general impairment in MAPK pathway with a greater impairment of the stimulation induced by insulin and IGF1 than by EGF. A PCR array analysis disclosed a defect in MAPK pathway characterized by an increase in DUSP4 and MEF2C gene expressions in patient fibroblasts. Conclusion: A post-receptor IGF1 insensitivity was characterized in one patient with SRS and ICR1 hypomethylation. Although based on one unique severely affected patient, these results raise an intriguing mechanism to explain the postnatal growth impairment observed in SRS patients that needs confirmation in larger cohorts.

Prognostication of Soft Tissue Sarcomas Based on Chromosome 17q Gene and Protein Status: Evaluation of TOP2A, HER-2/neu, and Survivin

Topoisomerase 2 alpha (), HER-2/ and are genes that lie on chromosome 17 and correlate with the prognosis and prediction of target-driven therapy against tumors. In a previous study, we showed that TOP2A transcripts levels were significantly higher in soft tissue sarcomas (STS) than in benign tumors and desmoid-type fibromatoses (FM). Because these genes have been insufficiently examined in STS, we aimed to identify alterations in TOP2A and HER-2 expression by fluorescent in situ hybridization and immunohistochemistry, as well as that of survivin, and correlate them with clinicopathologic findings to assess their prognostic value. Eighteen FM and 244 STS were included. Fluorescent in situ hybridization and immunohistochemistry were performed on a tissue microarray. TOP2A and survivin were more highly expressed in sarcomas than in FM. TOP2A was an independent predictor of an unfavorable prognosis; it was combined with formerly established prognostic factors (primarily histologic grade and tumor size at diagnosis) to create a prognostic index that evaluated overall survival. Gene amplification/polysomy (13%) did not correlate with protein overexpression. Survivin and HER-2 expression were not associated with patient outcomes. These findings might become valuable in the management of patients with STS and possibly in the prospective evaluation of responses to new target-driven therapies.

Interphase Chromosome Flow-FISH

A 2-day method using flow cytometry and FISH for interphase cells was developed to detect monosomy 7 cells in myelodysplastic syndrome patients. The method, Interphase Chromosome Flow-FISH (IC Flow-FISH), involves fixation of leukocytes from blood, membrane permeabilization, hybridization of cellular DNA with peptide nucleic acid probes with cells intact, and analysis by flow cytometry. Hundreds to thousands of monosomy 7 cells were consistently detected from 10-20 mL of blood in patients with monosomy 7. Proportions of monosomy 7 cells detected in IC Flow-FISH were compared with results from conventional cytogenetics; identification of monosomy 7 populations was verified with FACS; and patient and donor cells were mixed to test for sensitivity. IC Flow-FISH allows for detecting monosomy 7 without requiring bone marrow procurement or the necessity of metaphase spreads, and wider applications to other chromosomal abnormalities are in development. (Blood. 2012; 120(15): e54-e59)

Gene by environment QTL mapping through multiple trait analyses in blood pressure salt-sensitivity: identification of a novel QTL in rat chromosome 5

Background The genetic mechanisms underlying interindividual blood pressure variation reflect the complex interplay of both genetic and environmental variables. The current standard statistical methods for detecting genes involved in the regulation mechanisms of complex traits are based on univariate analysis. Few studies have focused on the search for and understanding of quantitative trait loci responsible for gene × environmental interactions or multiple trait analysis. Composite interval mapping has been extended to multiple traits and may be an interesting approach to such a problem. Methods We used multiple-trait analysis for quantitative trait locus mapping of loci having different effects on systolic blood pressure with NaCl exposure. Animals studied were 188 rats, the progenies of an F2 rat intercross between the hypertensive and normotensive strain, genotyped in 179 polymorphic markers across the rat genome. To accommodate the correlational structure from measurements taken in the same animals, we applied univariate and multivariate strategies for analyzing the data. Results We detected a new quantitative train locus on a region close to marker R589 in chromosome 5 of the rat genome, not previously identified through serial analysis of individual traits. In addition, we were able to justify analytically the parametric restrictions in terms of regression coefficients responsible for the gain in precision with the adopted analytical approach. Conclusion Future work should focus on fine mapping and the identification of the causative variant responsible for this quantitative trait locus signal. The multivariable strategy might be valuable in the study of genetic determinants of interindividual variation of antihypertensive drug effectiveness.

Genetic mapping of a new heart rate QTL on chromosome 8 of spontaneously hypertensive rats

Abstract Background Tachycardia is commonly observed in hypertensive patients, predominantly mediated by regulatory mechanisms integrated within the autonomic nervous system. The genetic loci and genes associated with increased heart rate in hypertension, however, have not yet been identified. Methods An F2 intercross of Spontaneously Hypertensive Rats (SHR) × Brown Norway (BN) linkage analysis of quantitative trait loci mapping was utilized to identify candidate genes associated with an increased heart rate in arterial hypertension. Results Basal heart rate in SHR was higher compared to that of normotensive BN rats (365 ± 3 vs. 314 ± 6 bpm, p < 0.05 for SHR and BN, respectively). A total genome scan identified one quantitative trait locus in a 6.78 cM interval on rat chromosome 8 (8q22–q24) that was responsible for elevated heart rate. This interval contained 241 genes, of which 65 are known genes. Conclusion Our data suggest that an influential genetic region located on the rat chromosome 8 contributes to the regulation of heart rate. Candidate genes that have previously been associated with tachycardia and/or hypertension were found within this QTL, strengthening our hypothesis that these genes are, potentially, associated with the increase in heart rate in a hypertension rat model.

Mechanisms of ring chromosome formation, ring instability and clinical consequences

Background The breakpoints and mechanisms of ring chromosome formation were studied and mapped in 14 patients. Methods Several techniques were performed such as genome-wide array, MLPA (Multiplex Ligation-Dependent Probe Amplification) and FISH (Fluorescent in situ Hybridization). Results The ring chromosomes of patients I to XIV were determined to be, respectively: r(3)(p26.1q29), r(4)(p16.3q35.2), r(10)(p15.3q26.2), r(10)(p15.3q26.13), r(13)(p13q31.1), r(13)(p13q34), r(14)(p13q32.33), r(15)(p13q26.2), r(18)(p11.32q22.2), r(18)(p11.32q21.33), r(18)(p11.21q23), r(22)(p13q13.33), r(22)(p13q13.2), and r(22)(p13q13.2). These rings were found to have been formed by different mechanisms, such as: breaks in both chromosome arms followed by end-to-end reunion (patients IV, VIII, IX, XI, XIII and XIV); a break in one chromosome arm followed by fusion with the subtelomeric region of the other (patients I and II); a break in one chromosome arm followed by fusion with the opposite telomeric region (patients III and X); fusion of two subtelomeric regions (patient VII); and telomere-telomere fusion (patient XII). Thus, the r(14) and one r(22) can be considered complete rings, since there was no loss of relevant genetic material. Two patients (V and VI) with r(13) showed duplication along with terminal deletion of 13q, one of them proved to be inverted, a mechanism known as inv-dup-del. Ring instability was detected by ring loss and secondary aberrations in all but three patients, who presented stable ring chromosomes (II, XIII and XIV). Conclusions We concluded that the clinical phenotype of patients with ring chromosomes may be related with different factors, including gene haploinsufficiency, gene duplications and ring instability. Epigenetic factors due to the circular architecture of ring chromosomes must also be considered, since even complete ring chromosomes can result in phenotypic alterations, as observed in our patients with complete r(14) and r(22).

Comparative genomic sequencing analysis of a region in human chromosome 11p15.3, mouse chromosome 7 and analysis of the novel STK33, Stk33 gene

The comparative genomic sequence analysis of a region in human chromosome 11p15.3 and its homologous segment in mouse chromosome 7 between ST5 and LMO1 genes has been performed. 158,201 bases were sequenced in the mouse and compared with the syntenic region in human, partially available in the public databases. The analysed region exhibits the typical eukaryotic genomic structure and compared with the close neighbouring regions, strikingly reflexes the mosaic pattern distribution of (G+C) and repeats content despites its relative short size. Within this region the novel gene STK33 was discovered (Stk33 in the mouse), that codes for a serine/threonine kinase. The finding of this gene constitutes an excellent example of the strength of the comparative sequencing approach. Poor gene-predictions in the mouse genomic sequence were corrected and improved by the comparison with the unordered data from the human genomic sequence publicly available. Phylogenetical analysis suggests that STK33 belongs to the calcium/calmodulin-dependent protein kinases group and seems to be a novelty in the chordate lineage. The gene, as a whole, seems to evolve under purifying selection whereas some regions appear to be under strong positive selection. Both human and mouse versions of serine/threonine kinase 33, consists of seventeen exons highly conserved in the coding regions, particularly in those coding for the core protein kinase domain. Also the exon/intron structure in the coding regions of the gene is conserved between human and mouse. The existence and functionality of the gene is supported by the presence of entries in the EST databases and was in vivo fully confirmed by isolating specific transcripts from human uterus total RNA and from several mouse tissues. Strong evidence for alternative splicing was found, which may result in tissue-specific starting points of transcription and in some extent, different protein N-termini. RT-PCR and hybridisation experiments suggest that STK33/Stk33 is differentially expressed in a few tissues and in relative low levels. STK33 has been shown to be reproducibly down-regulated in tumor tissues, particularly in ovarian tumors. RNA in-situ hybridisation experiments using mouse Stk33-specific probes showed expression in dividing cells from lung and germinal epithelium and possibly also in macrophages from kidney and lungs. Preliminary experimentation with antibodies designed in this work, performed in parallel to the preparation of this manuscript, seems to confirm this expression pattern. The fact that the chromosomal region 11p15 in which STK33 is located may be associated with several human diseases including tumor development, suggest further investigation is necessary to establish the role of STK33 in human health.

Breakpoint analysis of human chromosome 3 inversions during hominoid evolution

Comparative fluorescence in situ hybridization (FISH) mapping revealed four large DNA segments which have been conserved in their entirety between human chromosome 3 and Bornean orangutan chromosome 2 as well as three evolutionary breakpoints which distinguish between the human and Bornean orangutan chromosome forms. Examination of the structural and functional features of evolutionary breakpoints provides new insights into the possible effects of evolutionary rearrangements on genome function and the relationship between human chromosome pathology and evolution. FISH of human BAC clones which were assesssed in human genomic sequence to primate chromosomes, combined with precise breakpoint localizations by polymerase chain reaction (PCR) analysis of flow-sorted chromosomes and in silico analysis, were used to characterize the evolutionary breakpoints. None of the three breakpoints studied disrupts a validated gene(s), however they are all associated with segmental duplications. At least eleven DNA segments (&a

Genome characterization through a mathematical model of the genetic code: an analysis of the whole chromosome 1 of A. thaliana

The objective of this work is to characterize the genome of the chromosome 1 of A.thaliana, a small flowering plants used as a model organism in studies of biology and genetics, on the basis of a recent mathematical model of the genetic code. I analyze and compare different portions of the genome: genes, exons, coding sequences (CDS), introns, long introns, intergenes, untranslated regions (UTR) and regulatory sequences. In order to accomplish the task, I transformed nucleotide sequences into binary sequences based on the definition of the three different dichotomic classes. The descriptive analysis of binary strings indicate the presence of regularities in each portion of the genome considered. In particular, there are remarkable differences between coding sequences (CDS and exons) and non-coding sequences, suggesting that the frame is important only for coding sequences and that dichotomic classes can be useful to recognize them. Then, I assessed the existence of short-range dependence between binary sequences computed on the basis of the different dichotomic classes. I used three different measures of dependence: the well-known chi-squared test and two indices derived from the concept of entropy i.e. Mutual Information (MI) and Sρ, a normalized version of the “Bhattacharya Hellinger Matusita distance”. The results show that there is a significant short-range dependence structure only for the coding sequences whose existence is a clue of an underlying error detection and correction mechanism. No doubt, further studies are needed in order to assess how the information carried by dichotomic classes could discriminate between coding and noncoding sequence and, therefore, contribute to unveil the role of the mathematical structure in error detection and correction mechanisms. Still, I have shown the potential of the approach presented for understanding the management of genetic information.

Caratterizzazione Citogenetico-Molecolare delle alterazioni cromosomiche 3q26 e 1p36 nelle Sindromi Mieloproliferative

L’overespressione dei geni EVI1(3q26) e PRDM16(1p36), è descritta sia in presenza che in assenza di riarrangiamenti 3q26 e 1p36 in specifici sottogruppi citogenetici di LAM, ed è associata ad una prognosi sfavorevole. Lo scopo principale del nostro studio è stato identificare e caratterizzare tramite FISH e RQ-PCR, alterazioni di EVI1 e PRDM16 in pazienti con alterazioni cromosomiche 3q e 1p.Riarrangiamenti di EVI1 si associavano ad alterazioni cromosomiche 3q26, ma, in 6 casi (6/35;17,1%) erano presenti in assenza di coinvolgimenti, in citogenetica convenzionale, della regione 3q26, a causa di meccanismi complessi e/o alterazioni ‘criptiche’. Inoltre, abbiamo identificato quattro nuovi riarrangiamenti di EVI1, tra cui due nuove traslocazioni simili presenti in due fratelli. Riarrangiamenti e/o amplificazioni di PRDM16 erano spesso associate ad alterazioni 1p36 (7/14;50%). L’analisi di EVI1 e PRDM16 è stata estesa ad altri casi con alterazioni -7/7q-, con cariotipo normale, con alterazioni 3q per PRDM16 e con alterazioni 1p per EVI1. L’overespressione di EVI1 era presente solo nel gruppo -7/7q- (10/58;17.2%) ed in un caso si associava ad amplificazione genica, mentre PRDM16 era overespresso in casi di tutti i gruppi analizzati,sia con cariotipi complessi, dove si associava in alcuni casi ad amplificazione genica, sia con cariotipi normali o con singole alterazioni. Il nostro studio dimostra come la FISH permetta di identificare alterazioni dei geni EVI1 e PRDM16, anche in assenza di coinvolgimenti delle regioni 3q26 e 1p36. Riarrangiamenti complessi e/o una scarsa qualità dei preparati citogenetici sono le cause principali per la mancata identificazione di queste alterazioni. La RQ-PCR permette di identificare l’overespressione anche nei casi in cui non sia dovuta ad alterazioni citogenetiche. È importante confermare con FISH e/o RQ-PCR il coinvolgimento di questi due geni, per individuare alla diagnosi pazienti con prognosi sfavorevole e che potranno beneficiare di terapie maggiormente aggressive e/o di trapianto allogenico di cellule staminali.

Next-Generation Sequencing-Based Mutations Scanning Strategy of the BCR-ABL Kinase Domain in Patients with PhiladelPhia-Chromosome Positive Leukemias Treated with Tyrosine Kinase Inhibitors

In chronic myeloid leukemia and Philadelphia-positive acute lymphoblastic leukemia patients resistant to tyrosine kinase inhibitors (TKIs), BCR-ABL kinase domain mutation status is an essential component of the therapeutic decision algorithm. The recent development of Ultra-Deep Sequencing approach (UDS) has opened the way to a more accurate characterization of the mutant clones surviving TKIs conjugating assay sensitivity and throughput. We decided to set-up and validated an UDS-based for BCR-ABL KD mutation screening in order to i) resolve qualitatively and quantitatively the complexity and the clonal structure of mutated populations surviving TKIs, ii) study the dynamic of expansion of mutated clones in relation to TKIs therapy, iii) assess whether UDS may allow more sensitive detection of emerging clones, harboring critical 2GTKIs-resistant mutations predicting for an impending relapse, earlier than SS. UDS was performed on a Roche GS Junior instrument, according to an amplicon sequencing design and protocol set up and validated in the framework of the IRON-II (Interlaboratory Robustness of Next-Generation Sequencing) International consortium.Samples from CML and Ph+ ALL patients who had developed resistance to one or multiple TKIs and collected at regular time-points during treatment were selected for this study. Our results indicate the technical feasibility, accuracy and robustness of our UDS-based BCR-ABL KD mutation screening approach. UDS was found to provide a more accurate picture of BCR-ABL KD mutation status, both in terms of presence/absence of mutations and in terms of clonal complexity and showed that BCR-ABL KD mutations detected by SS are only the “tip of iceberg”. In addition UDS may reliably pick 2GTKIs-resistant mutations earlier than SS in a significantly greater proportion of patients.The enhanced sensitivity as well as the possibility to identify low level mutations point the UDS-based approach as an ideal alternative to conventional sequencing for BCR-ABL KD mutation screening in TKIs-resistant Ph+ leukemia patients

Examination of FGFRL1 as a candidate gene for diaphragmatic defects at chromosome 4p16.3 shows that Fgfrl1 null mice have reduced expression of Tpm3, sarcomere genes and Lrtm1 in the diaphragm

Fgfrl1 (also known as Fgfr5; OMIM 605830) homozygous null mice have thin, amuscular diaphragms and die at birth because of diaphragm hypoplasia. FGFRL1 is located at 4p16.3, and this chromosome region can be deleted in patients with congenital diaphragmatic hernia (CDH). We examined FGFRL1 as a candidate gene for the diaphragmatic defects associated with 4p16.3 deletions and re-sequenced this gene in 54 patients with CDH. We confirmed six known coding single nucleotide polymorphisms (SNPs): c.209G > A (p.Pro20Pro), c.977G > A (p.Pro276Pro), c.1040T > C (p.Asp297Asp), c.1234C > A (p.Pro362Gln), c.1420G > T (p.Arg424Leu), and c.1540C > T (p.Pro464Leu), but we did not identify any gene mutations. We genotyped additional CDH patients for four of these six SNPs, including the three non-synonymous SNPs, to make a total of 200 chromosomes, and found that the allele frequency for the four SNPs, did not differ significantly between patients and normal controls (p > or = 0.05). We then used Affymetrix Genechip Mouse Gene 1.0 ST arrays and found eight genes with significantly reduced expression levels in the diaphragms of Fgfrl1 homozygous null mice when compared with wildtype mice-Tpm3, Fgfrl1 (p = 0.004), Myl2, Lrtm1, Myh4, Myl3, Myh7 and Hephl1. Lrtm1 is closely related to Slit3, a protein associated with herniation of the central tendon of the diaphragm in mice. The Slit proteins are known to regulate axon branching and cell migration, and inhibition of Slit3 reduces cell motility and decreases the expression of Rac and Cdc42, two genes that are essential for myoblast fusion. Further studies to determine if Lrtm1 has a similar function to Slit3 and if reduced Fgfrl1 expression can cause diaphragm hypoplasia through a mechanism involving decreased myoblast motility and/or myoblast fusion, seem indicated.