Influence of genetic variants and post transcriptional factors on imatinib transporters as determinants of the pharmacological response

Introduction – Although imatinib (IM) is a recognized gold standard in chronic myeloid leukemia (CML) therapy, resistance has emerged in a significant proportion of patients. Aim – The aim of this study was: (1) to investigate the role of genetic variants in genes encoding for IM transporters, as candidate of IM responsiveness and (2) to test the influence of miRNAs on IM response, focusing on efflux transporters. Methods – As a first step, a panel of polymorphisms (SNPs) was genotyped in a subgroup population of 189 patients enrolled in the Tyrosine Kinase Inhibitor Optimization and Selectivity (TOPS) trial. The association with cytogenetic response and molecular response (MR) was assessed for each SNP. As a second step, an in vitro IM-resistant model (K-562 CML cell line) was established. miRNAs profiles were analyzed using Taqman arrays and in silico search was performed for miRNAs deregulated after IM treatment. mRNA and protein expression were quantified using TaqMan realtime PCR and Western blotting, respectively. Results – (1) Among Caucasian patients, ABCB1 rs60023214 significantly correlated with complete MR (P = 0.005). Concerning SNPs combination in IM uptake transporters, the associations with treatment outcomes were statistically significant for both major and complete MR (P = 0.005 and P = 0.01, respectively). (2) ABCB1 protein was not expressed under any conditions of treatment, differently from ABCG2. Two deregulated miRNAs, namely miR-212 and miR-328, were identified to be inversely correlated with ABCG2 (r2= 0.57; p=0.03 and r2=0.47; p=0.06, respectively). Experiments of loss and gain of function confirmed the functional influence of these miRNAs on ABCG2. Conclusion – The multiple candidate gene approach identified single and combination of SNPs that can be proposed as predictor of IM response. The in vitro study suggested that IM resistance could be mediated by miRNA-dependent mechanism. Further studies are needed to validate these preliminary findings.

Identification of Drosophila heart-specific Cis-Regulatory Modules under Hox control

Cardiac morphogenesis is a complex process governed by evolutionarily conserved transcription factors and signaling molecules. The Drosophila cardiac tube is linear, made of 52 pairs of cardiomyocytes (CMs), which express specific transcription factor genes that have human homologues implicated in Congenital Heart Diseases (CHDs) (NKX2-5, GATA4 and TBX5). The Drosophila cardiac tube is linear and composed of a rostral portion named aorta and a caudal one called heart, distinguished by morphological and functional differences controlled by Hox genes, key regulators of axial patterning. Overexpression and inactivation of the Hox gene abdominal-A (abd-A), which is expressed exclusively in the heart, revealed that abd-A controls heart identity. The aim of our work is to isolate the heart-specific cisregulatory sequences of abd-A direct target genes, the realizator genes granting heart identity. In each segment of the heart, four pairs of cardiomyocytes (CMs) express tinman (tin), homologous to NKX2-5, and acquire strong contractile and automatic rhythmic activities. By tyramide amplified FISH, we found that seven genes, encoding ion channels, pumps or transporters, are specifically expressed in the Tin-CMs of the heart. We initially used online available tools to identify their heart-specific cisregutatory modules by looking for Conserved Non-coding Sequences containing clusters of binding sites for various cardiac transcription factors, including Hox proteins. Based on these data we generated several reporter gene constructs and transgenic embryos, but none of them showed reporter gene expression in the heart. In order to identify additional abd-A target genes, we performed microarray experiments comparing the transcriptomes of aorta versus heart and identified 144 genes overexpressed in the heart. In order to find the heart-specific cis-regulatory regions of these target genes we developed a new bioinformatic approach where prediction is based on pattern matching and ordered statistics. We first retrieved Conserved Noncoding Sequences from the alignment between the D.melanogaster and D.pseudobscura genomes. We scored for combinations of conserved occurrences of ABD-A, ABD-B, TIN, PNR, dMEF2, MADS box, T-box and E-box sites and we ranked these results based on two independent strategies. On one hand we ranked the putative cis-regulatory sequences according to best scored ABD-A biding sites, on the other hand we scored according to conservation of binding sites. We integrated and ranked again the two lists obtained independently to produce a final rank. We generated nGFP reporter construct flies for in vivo validation. We identified three 1kblong heart-specific enhancers. By in vivo and in vitro experiments we are determining whether they are direct abd-A targets, demonstrating the role of a Hox gene in the realization of heart identity. The identified abd-A direct target genes may be targets also of the NKX2-5, GATA4 and/or TBX5 homologues tin, pannier and Doc genes, respectively. The identification of sequences coregulated by a Hox protein and the homologues of transcription factors causing CHDs, will provide a mean to test whether these factors function as Hox cofactors granting cardiac specificity to Hox proteins, increasing our knowledge on the molecular mechanisms underlying CHDs. Finally, it may be investigated whether these Hox targets are involved in CHDs.

Salt stress responses in pear and quince: physiological and molecular aspects

The productivity of agricultural crops is seriously limited by salinity. This problem is rapidly increasing, particularly in irrigated lands. Like almost all the fruit tree species, Pyrus communis is generally considered a salt sensitive species, but only little information is available on its behavior under saline conditions. Previous studies, carried out in the Department of Fruit Tree and Woody Plant Science (University of Bologna), focused their attention on pear and quince salt stress responses to understand which rootstock would be the most suitable for pear in order to tolerate a salt stress condition. It has been reported that pear and quince have different ability in the uptake, translocation and accumulation of chloride (Cl-) and sodium (Na+) ions, when plants were irrigated for one season with saline water (5 dS/m). The aim of the present work was to deepen these aspects and investigate salt stress responses in pear and quince. Two different experiments have been performed: a “short-term” trial in a growth chamber and a “long-term” experiment in the open field. In the short-term experiment, three different genotypes usually adopted as pear rootstocks (MC, BA29 and Farold®40) and the pear variety Abbé Fétel own rooted have been compared under salt stress conditions. The trial was performed in a hydroponic culture system, applying a 90 mM NaCl stress to half of the plants, after five weeks of normal growth in Hoagland’s solution. During the three-weeks of salt stress treatment, physiological, mineral and molecular analyses were performed in order to monitor, for each genotype, the development of the salt stress responses in comparison with the corresponding “unstressed” plants. Farold®40 and Abbé Fétel own rooted showed the onset of leaf necrosis, due to salt toxicity, one week before quinces. Moreover, quinces displayed a significant delay in premature senescence of old leaves, while pears emerged for their ability to regenerate new leaves from apparently dead foliage with the salt stress still running. Physiological measurements, such as shoots length, chlorophyll (Chl) content, and photosynthesis, have been carried out and revealed that pears exhibited a significant reduction in water content and a wilting aspect, while for quinces a decrease in Chl content and a growth slowdown were observed. At the end of the trial, all plants were collected and organs separated for dry weight estimation and mineral analyses (Cu, Fe, Mn, Zn Mg, Ca, K, Na and Cl). Mineral contents have been affected by salinity; same macro/micro nutrients were altered in some organs or relocated within the plant. This plant response could have partially contributed to face the salt stress. Leaves and roots have been harvested for molecular analyses at four different times during stress conditions. Molecular analyses consisted of the gene expression study of three main ion transporters, well known in Arabidopsis thaliana as salt-tolerance determinants in the “SOS” pathway: NHX1 (tonoplast Na+/H+ antiporter), SOS1 (plasmalemma Na+/H+ antiporter) and HKT1 (K+ high-affinity and Na+ low-affinity transporter). These studies showed that two quince rootstocks adopted different responsive mechanisms to NaCl stress. BA29 increased its Na+ sequestration activity into leaf vacuoles, while MC enhanced temporarily the same ability, but in roots. Farold®40, instead, exhibited increases in SOS1 and HKT1 expression mainly at leaf level in the attempt to retrieve Na+ from xylem, while Abbé Fétel differently altered the expression of these genes in roots. Finally, each genotype showed a peculiar response to salt stress that was the sum of its ability in Na+ exclusion, osmotic tolerance and tissue tolerance. In the long-term experiment, potted trees of the pear variety Abbé Fétel grafted on different rootstocks (MC, BA29 and Farold®40), or own rooted and also rootstocks only were subjected to a salt stress through saline water irrigation with an electrical conductivity of 5 dS/m for two years. The purposes of this study were to evaluate salinity effects on physiological (shoot length, number of buds, photosynthesis, etc.) and yield parameters of cultivar Abbé Fétel in the different combinations and to determine the salt amount that pear is able to tolerate over the years. With this work, we confirmed the previous hypothesis that pear, despite being classified as a salt-sensitive fruit tree, can be cultivated for two years under saline water irrigation, without showing any salt toxicity symptoms or severe drawbacks on plant development and production. Among different combinations, Abbé Fétel grafted on MC resulted interesting for its peculiar behaviors under salt stress conditions. In the near future, further investigations on physiological and molecular aspects will be necessary to enrich and broaden the knowledge of salt stress responses in pear.

Detection and localization of GLUTs in spermatozoa from different domestic species

Sperm cells need hexoses as a substrate for their function, for both the maintenance of membrane homeostasis and the movement of the tail. These cells have a peculiar metabolism that has not yet been fully understood, but it is clear that they obtain energy from hexoses through glycolisis and/or oxidative phosphorylation. Spermatozoa are in contact with different external environments, beginning from the testicular and epididymal fluid, passing to the seminal plasma and finally to the female genital tract fluids; in addition, with the spread of reproductive biotechnologies, sperm cells are diluted and stored in various media, containing different energetic substrates. To utilize these energetic sources, sperm cells, as other eukaryotic cells, have a well-constructed protein system, that is mainly represented by the GLUT family proteins. These transporters have a membrane-spanning α-helix structure and work as an enzymatic pump that permit a fast gradient dependent passage of sugar molecules through the lipidic bilayer of sperm membrane. Many GLUTs have been studied in man, bull and rat spermatozoa; the presence of some GLUTs has been also demonstrated in boar and dog spermatozoa. The aims of the present study were - to determine the presence of GLUTs 1, 2, 3, 4 and 5 in boar, horse, dog and donkey spermatozoa and to describe their localization; - to study eventual changes in GLUTs location after capacitation and acrosome reaction in boar, stallion and dog spermatozoa; - to determine possible changes in GLUTs localization after capacitation induced by insulin and IGF stimulation in boar spermatozoa; - to evaluate changes in GLUTs localization after flow-cytometric sex sorting in boar sperm cells. GLUTs 1, 2, 3 and 5 presence and localization have been demonstrated in boar, stallion, dog and donkey spermatozoa by western blotting and immunofluorescence analysis; a relocation in GLUTs after capacitation has been observed only in dog sperm cells, while no changes have been observed in the other species examined. As for boar, the stimulation of the capacitation with insulin and IGF didn’t cause any change in GLUTs localization, as well as for the flow cytometric sorting procedure. In conclusion, this study confirms the presence of GLUTs 1, 2 ,3 and 5 in boar, dog, stallion and donkey spermatozoa, while GLUT 4 seems to be absent, as a confirmation of other studies. Only in dog sperm cells capacitating conditions induce a change in GLUTs distribution, even if the physiological role of these changes should be deepened.

The impact of polymorphisms in P-gp, DNA repair and folic acid metabolism genes in newly diagnosed multiple myeloma patients treated with thalidomide plus dexamethasone, with or without bortezomib

The principle aim of this study was to investigate biological predictors of response and resistance to multiple myeloma treatment. Two hypothesis had been proposed as responsible of responsiveness: SNPs in DNA repair and Folate pathway, and P-gp dependent efflux. As a first objective, panel of SNPs in DNA repair and Folate pathway genes, were analyzed. It was a retrospective study in a group of 454, previously untreated, MM patients enrolled in a randomized phase III open-label study. Results show that some SNPs in Folate pathway are correlated with response to MM treatment. MTR genotype was associated with favorable response in the overall population of MM patients. However, this relation, disappear after adjustment for treatment response. When poor responder includes very good partial response, partial response and stable/progressive disease MTFHR rs1801131 genotype was associated with poor response to therapy. This relation - unlike in MTR – was still significant after adjustment for treatment response. Identification of this genetic variant in MM patients could be used as an independent prognostic factor for therapeutic outcome in the clinical practice. In the second objective, basic disposition characteristics of bortezomib was investigated. We demonstrated that bortezomib is a P-gp substrate in a bi-directional transport study. We obtain apparent permeability rate values that together with solubility values can have a crucial implication in better understanding of bortezomib pharmacokinetics with respect to the importance of membrane transporters. Subsequently, in view of the importance of P-gp for bortezomib responsiveness a panel of SNPs in ABCB1 gene - coding for P-gp - were analyzed. In particular we analyzed five SNPs, none of them however correlated with treatment responsiveness. However, we found a significant association between ABCB1 variants and cytogenetic abnormalities. In particular, deletion of chromosome 17 and t(4;14) translocation were present in patients harboring rs60023214 and rs2038502 variants respectively.

A clustering method for robust and reliable large scale functional and structural protein sequence annotation

Bioinformatics, in the last few decades, has played a fundamental role to give sense to the huge amount of data produced. Obtained the complete sequence of a genome, the major problem of knowing as much as possible of its coding regions, is crucial. Protein sequence annotation is challenging and, due to the size of the problem, only computational approaches can provide a feasible solution. As it has been recently pointed out by the Critical Assessment of Function Annotations (CAFA), most accurate methods are those based on the transfer-by-homology approach and the most incisive contribution is given by cross-genome comparisons. In the present thesis it is described a non-hierarchical sequence clustering method for protein automatic large-scale annotation, called “The Bologna Annotation Resource Plus” (BAR+). The method is based on an all-against-all alignment of more than 13 millions protein sequences characterized by a very stringent metric. BAR+ can safely transfer functional features (Gene Ontology and Pfam terms) inside clusters by means of a statistical validation, even in the case of multi-domain proteins. Within BAR+ clusters it is also possible to transfer the three dimensional structure (when a template is available). This is possible by the way of cluster-specific HMM profiles that can be used to calculate reliable template-to-target alignments even in the case of distantly related proteins (sequence identity < 30%). Other BAR+ based applications have been developed during my doctorate including the prediction of Magnesium binding sites in human proteins, the ABC transporters superfamily classification and the functional prediction (GO terms) of the CAFA targets. Remarkably, in the CAFA assessment, BAR+ placed among the ten most accurate methods. At present, as a web server for the functional and structural protein sequence annotation, BAR+ is freely available at http://bar.biocomp.unibo.it/bar2.0.

SMO inhibitor specifically targets the Hedgehog Pathway and reverts the drug-resistance of Leukemic Stem Cells

Abnormal Hedgehog signaling is associated with human malignancies. Smo, a key player of that signaling, is the most suitable target to inhibit this pathway. To this aim several molecules, antagonists of Smo, have been synthesized, and some of them have started the phase I in clinical trials. Our hospital participated to one of these studies which investigated the oral administration of a new selective inhibitor of Smo (SMOi). To evaluate ex vivo SMOi efficacy and to identify new potential clinical biomarkers of responsiveness, we separated bone marrow CD34+ cells from 5 acute myeloid leukemia (AML), 1 myelofibrosis (MF), 2 blastic phases chronic myeloid leukemia (CML) patients treated with SMOi by immunomagnetic separation, and we analysed their gene expression profile using Affimetrix HG-U133 Plus 2.0 platform. This analysis, showed differential expression after 28 days start of therapy (p-value ≤ 0.05) of 1,197 genes in CML patients and 589 genes in AML patients. This differential expression is related to Hedgehog pathway with a p-value = 0.003 in CML patients and with a p-value = 0.0002 in AML patients, suggesting that SMOi targets specifically this pathway. Among the genes differentially expressed we observed strong up-regulation of Gas1 and Kif27 genes, which may work as biomarkers of responsiveness of SMOi treatment in CML CD34+ cells whereas Hedgehog target genes (such as Smo, Gli1, Gli2, Gli3), Bcl2 and Abca2 were down-regulated, in both AML and CML CD34+ cells. It has been reported that Bcl-2 expression could be correlated with cancer therapy resistance and that Hedgehog signaling modulate ATP-binding (ABC) cassette transporters, whose expression has been correlated with chemoresistance. Moreover we confirmed that in vitro SMOi treatment targets Hedgehog pathway, down-regulate ABC transporters, Abcg2 and Abcb1 genes, and in combination with tyrosine kinase inhibitors (TKIs) could revert the chemoresistance mechanism in K562 TKIs-resistant cell line.

Valutazione del coinvolgimento dell'oncogene c-myc e dei trasportatori abc nella farmacoresistenza dell'osteosarcoma ad alto grado

L’insorgenza di fenomeni coinvolti nello sviluppo della farmacoresistenza costituisce al momento la principale causa di mancata risposta al trattamento chemioterapico nell’osteosarcoma. Questo è in parte dovuto ad una sovraespressione di diversi trasportatori ABC nelle cellule tumorali che causano un aumento dell’efflusso extracellulare del chemioterapico e pertanto una ridotta risposta al trattamento farmacologico. L'oncogene C-MYC è coinvolto nella resistenza al metothrexate, alla doxorubicina e al cisplatino ed è un fattore prognostico avverso, se sovraespresso al momento della diagnosi, in pazienti affetti da osteosarcoma. C-MYC è in grado di regolare l'espressione di diversi trasportatori ABC, probabilmente coinvolti nella resistenza ai farmaci nell’osteosarcoma, e questo potrebbe spiegare l’impatto prognostico avverso dell’oncogene in questo tumore. L’espressione genica di C-MYC e di 16 trasportatori ABC, regolati da C-MYC e / o responsabili dell'efflusso di diversi chemioterapici, è stata valutata su due diverse casistiche cliniche e su un pannello di linee cellulari di osteosarcoma umano mediante real-time PCR. L'espressione della proteina è stata valutata per i 9 trasportatori ABC risultati più rilevanti.Infine l'efficacia in vitro di un inibitore, specifico per ABCB1 e ABCC1, è stata valutata su linee cellulari di osteosarcoma. ABCB1 e ABCC1 sono i trasportatori più espressi nelle linee cellulari di osteosarcoma. ABCB1 è sovraespresso al momento della diagnosi in circa il 40-45% dei pazienti affetti da osteosarcoma e si conferma essere un fattore prognostico avverso se sovraespresso al momento della diagnosi. Pertanto ABCB1 diventa il bersaglio di elezione per lo sviluppo di strategie terapeutiche alternative, nel trattamento dell’osteosarcoma, atte al superamento della farmacoresistenza. L’inibizione dell'attività di tale trasportatore causa un aumento della sensibilità al trattamento chemioterapico nelle linee cellulari di osteosarcoma farmacoresistenti, indicando questo approccio come una possibile strategia per superare il problema della mancata risposta al trattamento farmacologico nei pazienti con osteosarcoma che sovraesprimono ABCB1.

Mechanisms Contributing to Tyrosin Kinase Inhibitor Resistance in GISTs: Toward a Personalized Therapy

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal tumors in the gastrointestinal tract. This work considers the pharmacological response in GIST patients treated with imatinib by two different angles: the genetic and somatic point of view. We analyzed polymorphisms influence on treatment outcome, keeping in consideration SNPs in genes involved in drug transport and folate pathway. Naturally, all these intriguing results cannot be considered as the only main mechanism in imatinib response. GIST mainly depends by oncogenic gain of function mutations in tyrosin kinase receptor genes, KIT or PDGFRA, and the mutational status of these two genes or acquisition of secondary mutation is considered the main player in GIST development and progression. To this purpose we analyzed the secondary mutations to better understand how these are involved in imatinib resistance. In our analysis we considered both imatinib and the second line treatment, sunitinib, in a subset of progressive patients. KIT/PDGFRA mutation analysis is an important tool for physicians, as specific mutations may guide therapeutic choices. Currently, the only adaptations in treatment strategy include imatinib starting dose of 800 mg/daily in KIT exon-9-mutated GISTs. In the attempt to individualize treatment, genetic polymorphisms represent a novelty in the definition of biomarkers of imatinib response in addition to the use of tumor genotype. Accumulating data indicate a contributing role of pharmacokinetics in imatinib efficacy, as well as initial response, time to progression and acquired resistance. At the same time it is becoming evident that genetic host factors may contribute to the observed pharmacokinetic inter-patient variability. Genetic polymorphisms in transporters and metabolism may affect the activity or stability of the encoded enzymes. Thus, integrating pharmacogenetic data of imatinib transporters and metabolizing genes, whose interplay has yet to be fully unraveled, has the potential to provide further insight into imatinib response/resistance mechanisms.