The search for Multiple Myeloma Stem cells: molecular characterization and self-renewal mechanisms involved in the disease persistence

The existence of Multiple Myeloma Stem cells (MMSCs)is supposed to be one of the major causes of MM drug-resistance. However, very little is known about the molecular characteristics of MMSCs, even if some studies suggested that these cells resembles the memory B cells. In order to molecularly characterize MMSCs, we isolated the 138+138- population. For each cell fraction we performed a VDJ rearrangement analysis. The complete set of aberrations were performed by SNP Array 6.0 and HG-U133 Plus 2.0 microarray analyses (Affymetrix). The VDJ rearrangement analyses confirmed the clonal relationship between the 138+ clone and the immature clone. Both BM and PBL 138+ clones showed exactly the same genomic macroalterations. In the BM and PBL 138-19+27+ cell fractions several micro-alterations (range: 1-350 Kb) unique of the memory B cells clone were highlighted. Any micro-alterations detected were located out of any genomic variants region and are presumably associated to the MM pathogenesis, as confirmed by the presence of KRAS, WWOX and XIAP genes among the amplified regions. To get insight into the biology of the clonotypic B cell population, we compared the gene expression profile of 8 MM B cells samples 5 donor B cells vs, thus showing a differential expression of 11480 probes (p-value: <0,05). Among the self-renewal mechanisms, we observed the down-regulation of Hedgehog pathway and the iperactivation of Notch and Wnt signaling. Moreover, these immature cells showed a particular phenotype correlated to resistance to proteasome inhibitors (IRE1α-XBP1: -18.0; -19.96. P<0,05). Data suggested that the MM 138+ clone might resume the end of the complex process of myelomagenesis, whereas the memory B cells have some intriguing micro-alterations and a specific transcriptional program, supporting the idea that these post germinal center cells might be involved in the transforming event that originate and sustain the neoplastic clone.

Analisi di “Copy Number Variants” ed identificazione di nuovi geni candidati per l’Autismo e Ritardo Mentale

I disturbi dello spettro autistico (DSA) ed il ritardo mentale (RM) sono caratterizzati da un’eziologia genetica complessa ed eterogenea. Grazie ai recenti sviluppi nella ricerca genomica, è stato possibile dimostrare il ruolo di numerose copy number variants (CNVs) nella patogenesi di questi disturbi, anche se nella maggior parte dei casi l’eziologia rimane ancora sconosciuta. Questo lavoro riguarda l’identificazione e la caratterizzazione dei CNVs in famiglie con DSA e RM. E’ stata studiata una microdelezione in 7q31 che coinvolge i geni IMMP2L e DOCK4, trasmessa dalla madre con dislessia a due figli con autismo ed una figlia con dislessia. Nella stessa famiglia segrega una seconda microdelezione in 2q14 che inattiva il gene CNTNAP5 ed è trasmessa dal padre (con tratti autistici) ai due figli con autismo. Abbiamo quindi ipotizzato che i geni DOCK4 e CNTNAP5 potessero essere implicati, rispettivamente, nella suscettibilità a dislessia e DSA. Lo screening di numerosi individui affetti ha supportato la nostra ipotesi, con l’identificazione di una nuova microdelezione di DOCK4 che segrega con la dislessia, e 3 nuove varianti missenso in CNTNAP5 in individui con autismo. Dall’analisi genomica comparativa su array (aCGH) di individui con RM, è stata identificata una delezione nella regione 7q31.32, che coinvolge il gene CADPS2, in due fratelli con RM e tratti autistici, probabilmente ereditata dalla madre. Lo screening di mutazione di questo gene in individui con autismo o RM, ha portato all’identificazione di 3 varianti non sinonime, assenti nei controlli, ed ereditate per via materna. Poiché CADPS2 risiede in una regione genomica che contiene loci soggetti ad imprinting, abbiamo ipotizzato che il gene CADPS2 possa essere anch’esso caratterizzato da imprinting, con espressione monoallelica materna. Lo studio di espressione di CADPS2 in cellule del sangue ha avvalorato questa ipotesi, implicando perciò CADPS2 come un nuovo gene di suscettibilità per il RM e DSA.

Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

Contribution of rare damaging variants in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a heterogeneous and highly heritable neurodevelopmental disorder with a complex genetic architecture, consisting of a combination of common low-risk and more penetrant rare variants. This PhD project aimed to explore the contribution of rare variants in ASD susceptibility through NGS approaches in a cohort of 106 ASD families including 125 ASD individuals. Firstly, I explored the contribution of inherited rare variants towards the ASD phenotype in a girl with a maternally inherited pathogenic NRXN1 deletion. Whole exome sequencing of the trio family identified an increased burden of deleterious variants in the proband that could modulate the CNV penetrance and determine the disease development. In the second part of the project, I investigated the role of rare variants emerging from whole genome sequencing in ASD aetiology. To properly manage and analyse sequencing data, a robust and efficient variant filtering and prioritization pipeline was developed, and by its application a stringent set of rare recessive-acting and ultra-rare variants was obtained. As a first follow-up, I performed a preliminary analysis on de novo variants, identifying the most likely deleterious variants and highlighting candidate genes for further analyses. In the third part of the project, considering the well-established involvement of calcium signalling in the molecular bases of ASD, I investigated the role of rare variants in voltage-gated calcium channels genes, that mainly regulate intracellular calcium concentration, and whose alterations have been correlated with enhanced ASD risk. Specifically, I functionally tested the effect of rare damaging variants identified in CACNA1H, showing that CACNA1H variation may be involved in ASD development by additively combining with other high risk variants. This project highlights the challenges in the analysis and interpretation of variants from NGS analysis in ASD, and underlines the importance of a comprehensive assessment of the genomic landscape of ASD individuals.

Genomic and non genomic effects of elevated concentration of anabolic steroids in human neuronal cells

Nandrolone and other anabolic androgenic steroids (AAS) at elevated concentration can alter the expression and function of neurotransmitter systems and contribute to neuronal cell death. This effect can explain the behavioural changes, drug dependence and neuro degeneration observed in steroid abuser. Nandrolone treatment (10-8M–10-5M) caused a time- and concentration-dependent downregulation of mu opioid receptor (MOPr) transcripts in SH-SY5Y human neuroblastoma cells. This effect was prevented by the androgen receptor (AR) antagonist hydroxyflutamide. Receptor binding assays confirmed a decrease in MOPr of approximately 40% in nandrolonetreated cells. Treatment with actinomycin D (10-5M), a transcription inhibitor, revealed that nandrolone may regulate MOPr mRNA stability. In SH-SY5Y cells transfected with a human MOPr luciferase promoter/reporter construct, nandrolone did not alter the rate of gene transcription. These results suggest that nandrolone may regulate MOPr expression through post-transcriptional mechanisms requiring the AR. Cito-toxicity assays demonstrated a time- and concentration dependent decrease of cells viability in SH-SY5Y cells exposed to steroids (10-6M–10-4M). This toxic effects is independent of activation of AR and sigma-2 receptor. An increased of caspase-3 activity was observed in cells treated with Nandrolone 10-6M for 48h. Collectively, these data support the existence of two cellular mechanisms that might explain the neurological syndromes observed in steroids abuser.

Mechanism of resistance to tyrosine kinase inhibitors in philadelphia-positive acute lymphblastic leukaemia (all): from genetic alterations to impaired RNA editing

The Ph chromosome is the most frequent cytogenetic aberration associated with adult ALL and it represents the single most significant adverse prognostic marker. Despite imatinib has led to significant improvements in the treatment of patients with Ph+ ALL, in the majority of cases resistance developed quickly and disease progressed. Some mechanisms of resistance have been widely described but the full knowledge of contributing factors, driving both the disease and resistance, remains to be defined. The observation of rapid development of lymphoblastic leukemia in mice expressing altered Ikaros (Ik) isoforms represented the background of this study. Ikaros is a zinc finger transcription factor required for normal hemopoietic differentiation and proliferation, particularly in the lymphoid lineages. By means of alternative splicing, Ikaros encodes several proteins that differ in their abilities to bind to a consensus DNA-binding site. Shorter, DNA nonbinding isoforms exert a dominant negative effect, inhibiting the ability of longer heterodimer partners to bind DNA. The differential expression pattern of Ik isoforms in Ph+ ALL patients was analyzed in order to determine if molecular abnormalities involving the Ik gene could associate with resistance to imatinib and dasatinib. Bone marrow and peripheral blood samples from 46 adult patients (median age 55 yrs, 18-76) with Ph+ ALL at diagnosis and during treatment with imatinib (16 pts) or dasatinib (30 pts) were collected. We set up a fast, high-throughput method based on capillary electrophoresis technology to detect and quantify splice variants. 41% Ph+ ALL patients expressed high levels of the non DNA-binding dominant negative Ik6 isoform lacking critical N-terminal zinc-fingers which display abnormal subcellular compartmentalization pattern. Nuclear extracts from patients expressed Ik6 failed to bind DNA in mobility shift assay using a DNA probe containing an Ikaros-specific DNA binding sequence. In 59% Ph+ ALL patients there was the coexistence in the same PCR sample and at the same time of many splice variants corresponded to Ik1, Ik2, Ik4, Ik4A, Ik5A, Ik6, Ik6 and Ik8 isoforms. In these patients aberrant full-length Ikaros isoforms in Ph+ ALL characterized by a 60-bp insertion immediately downstream of exon 3 and a recurring 30-bp in-frame deletion at the end of exon 7 involving most frequently the Ik2, Ik4 isoforms were also identified. Both the insertion and deletion were due to the selection of alternative splice donor and acceptor sites. The molecular monitoring of minimal residual disease showed for the first time in vivo that the Ik6 expression strongly correlated with the BCR-ABL transcript levels suggesting that this alteration could depend on the Bcr-Abl activity. Patient-derived leukaemia cells expressed dominant-negative Ik6 at diagnosis and at the time of relapse, but never during remission. In order to mechanistically demonstrated whether in vitro the overexpression of Ik6 impairs the response to tyrosine kinase inhibitors (TKIs) and contributes to resistance, an imatinib-sensitive Ik6-negative Ph+ ALL cell line (SUP-B15) was transfected with the complete Ik6 DNA coding sequence. The expression of Ik6 strongly increased proliferation and inhibited apoptosis in TKI sensitive cells establishing a previously unknown link between specific molecular defects that involve the Ikaros gene and the resistance to TKIs in Ph+ ALL patients. Amplification and genomic sequence analysis of the exon splice junction regions showed the presence of 2 single nucleotide polymorphisms (SNPs): rs10251980 [A/G] in the exon2/3 splice junction and of rs10262731 [A/G] in the exon 7/8 splice junction in 50% and 36% of patients, respectively. A variant of the rs11329346 [-/C], in 16% of patients was also found. Other two different single nucleotide substitutions not recognized as SNP were observed. Some mutations were predicted by computational analyses (RESCUE approach) to alter cis-splicing elements. In conclusion, these findings demonstrated that the post-transcriptional regulation of alternative splicing of Ikaros gene is defective in the majority of Ph+ ALL patients treated with TKIs. The overexpression of Ik6 blocking B-cell differentiation could contribute to resistance opening a time frame, during which leukaemia cells acquire secondary transforming events that confer definitive resistance to imatinib and dasatinib.

TNFRSF13B Genetic variability an anthropological - evolutionary approach to Biomedical Research

In the recent years TNFRSF13B coding variants have been implicated by clinical genetics studies in Common Variable Immunodeficiency (CVID), the most common clinically relevant primary immunodeficiency in individuals of European ancestry, but their functional effects in relation to the development of the disease have not been entirely established. To examine the potential contribution of such variants to CVID, the more comprehensive perspective of an evolutionary approach was applied in this study, underling the belief that evolutionary genetics methods can play a role in dissecting the origin, causes and diffusion of human diseases, representing a powerful tool also in human health research. For this purpose, TNFRSF13B coding region was sequenced in 451 healthy individuals belonging to 26 worldwide populations, in addition to 96 control, 77 CVID and 38 Selective IgA Deficiency (IgAD) individuals from Italy, leading to the first achievement of a global picture of TNFRSF13B nucleotide diversity and haplotype structure and making suggestion of its evolutionary history possible. A slow rate of evolution, within our species and when compared to the chimpanzee, low levels of genetic diversity geographical structure and the absence of recent population specific selective pressures were observed for the examined genomic region, suggesting that geographical distribution of its variability is more plausibly related to its involvement also in innate immunity rather than in adaptive immunity only. This, together with the extremely subtle disease/healthy samples differences observed, suggests that CVID might be more likely related to still unknown environmental and genetic factors, rather than to the nature of TNFRSF13B variants only.

Genomic characterization of the italian wolf (Canis lupus): the genes involved in black coat colour determination and application of microarray technique for snps detection.

This study provides a comprehensive genetic overview on the endangered Italian wolf population. In particular, it focuses on two research lines. On one hand, we focalised on melanism in wolf in order to isolate a mutation related with black coat colour in canids. With several reported black individuals (an exception at European level), the Italian wolf population constituted a challenging research field posing many unanswered questions. As found in North American wolf, we reported that melanism in the Italian population is caused by a different melanocortin pathway component, the K locus, in which a beta-defensin protein acts as an alternative ligand for the Mc1r. This research project was conducted in collaboration with Prof. Gregory Barsh, Department of Genetics and Paediatrics, Stanford University. On the other hand, we performed analysis on a high number of SNPs thanks to a customized Canine microarray useful to integrate or substitute the STR markers for genotyping individuals and detecting wolf-dog hybrids. Thanks to DNA microchip technology, we obtained an impressive amount of genetic data which provides a solid base for future functional genomic studies. This study was undertaken in collaboration with Prof. Robert K. Wayne, Department of Ecology and Evolutionary Biology, University of California, Los Angeles (UCLA).

Genomic and transcriptional analysis of allergen genes in apple (Malus x domestica)

Apple consumption is highly recomended for a healthy diet and is the most important fruit produced in temperate climate regions. Unfortunately, it is also one of the fruit that most ofthen provoks allergy in atopic patients and the only treatment available up to date for these apple allergic patients is the avoidance. Apple allergy is due to the presence of four major classes of allergens: Mal d 1 (PR-10/Bet v 1-like proteins), Mal d 2 (Thaumatine-like proteins), Mal d 3 (Lipid transfer protein) and Mal d 4 (profilin). In this work new advances in the characterization of apple allergen gene families have been reached using a multidisciplinary approach. First of all, a genomic approach was used for the characterization of the allergen gene families of Mal d 1 (task of Chapter 1), Mal d 2 and Mal d 4 (task of Chapter 5). In particular, in Chapter 1 the study of two large contiguos blocks of DNA sequences containing the Mal d 1 gene cluster on LG16 allowed to acquire many new findings on number and orientation of genes in the cluster, their physical distances, their regulatory sequences and the presence of other genes or pseudogenes in this genomic region. Three new members were discovered co-localizing with the other Mal d 1 genes of LG16 suggesting that the complexity of the genetic base of allergenicity will increase with new advances. Many retrotranspon elements were also retrieved in this cluster. Due to the developement of molecular markers on the two sequences, the anchoring of the physical and the genetic map of the region has been successfully achieved. Moreover, in Chapter 5 the existence of other loci for the Thaumatine-like protein family in apple (Mal d 2.03 on LG4 and Mal d 2.02 on LG17) respect the one reported up to now was demonstred for the first time. Also one new locus for profilins (Mal d 4.04) was mapped on LG2, close to the Mal d 4.02 locus, suggesting a cluster organization for this gene family, as is well reported for Mal d 1 family. Secondly, a methodological approach was used to set up an highly specific tool to discriminate and quantify the expression of each Mal d 1 allergen gene (task of Chapter 2). In aprticular, a set of 20 Mal d 1 gene specific primer pairs for the quantitative Real time PCR technique was validated and optimized. As a first application, this tool was used on leaves and fruit tissues of the cultivar Florina in order to identify the Mal d 1 allergen genes that are expressed in different tissues. The differential expression retrieved in this study revealed a tissue-specificity for some Mal d 1 genes: 10/20 Mal d 1 genes were expressed in fruits and, indeed, probably more involved in the allergic reactions; while 17/20 Mal d 1 genes were expressed in leaves challenged with the fungus Venturia inaequalis and therefore probably interesting in the study of the plant defense mechanism. In Chapter 3 the specific expression levels of the 10 Mal d 1 isoallergen genes, found to be expressed in fruits, were studied for the first time in skin and flesh of apples of different genotypes. A complex gene expression profile was obtained due to the high gene-, tissue- and genotype-variability. Despite this, Mal d 1.06A and Mal d 1.07 expression patterns resulted particularly associated with the degree of allergenicity of the different cultivars. They were not the most expressed Mal d 1 genes in apple but here it was hypotized a relevant importance in the determination of allergenicity for both qualitative and quantitative aspects of the Mal d 1 gene expression levels. In Chapter 4 a clear modulation for all the 17 PR-10 genes tested in young leaves of Florina after challenging with the fungus V. inaequalis have been reported but with a peculiar expression profile for each gene. Interestingly, all the Mal d 1 genes resulted up-regulated except Mal d 1.10 that was down-regulated after the challenging with the fungus. The differences in direction, timing and magnitude of induction seem to confirm the hypothesis of a subfunctionalization inside the gene family despite an high sequencce and structure similarity. Moreover, a modulation of PR-10 genes was showed both in compatible (Gala-V. inaequalis) and incompatible (Florina-V. inaequalis) interactions contribute to validate the hypothesis of an indirect role for at least some of these proteins in the induced defense responses. Finally, a certain modulation of PR-10 transcripts retrieved also in leaves treated with water confirm their abilty to respond also to abiotic stress. To conclude, the genomic approach used here allowed to create a comprehensive inventory of all the genes of allergen families, especially in the case of extended gene families like Mal d 1. This knowledge can be considered a basal prerequisite for many further studies. On the other hand, the specific transcriptional approach make it possible to evaluate the Mal d 1 genes behavior on different samples and conditions and therefore, to speculate on their involvement on apple allergenicity process. Considering the double nature of Mal d 1 proteins, as apple allergens and as PR-10 proteins, the gene expression analysis upon the attack of the fungus created the base for unravel the Mal d 1 biological functions. In particular, the knowledge acquired in this work about the PR-10 genes putatively more involved in the specific Malus-V. inaequalis interaction will be helpful, in the future, to drive the apple breeding for hypo-allergenicity genotype without compromise the mechanism of response of the plants to stress conditions. For the future, the survey of the differences in allergenicity among cultivars has to be be thorough including other genotypes and allergic patients in the tests. After this, the allelic diversity analysis with the high and low allergenic cultivars on all the allergen genes, in particular on the ones with transcription levels correlated to allergencity, will provide the genetic background of the low ones. This step from genes to alleles will allow the develop of molecular markers for them that might be used to effectively addressed the apple breeding for hypo-allergenicity. Another important step forward for the study of apple allergens will be the use of a specific proteomic approach since apple allergy is a multifactor-determined disease and only an interdisciplinary and integrated approach can be effective for its prevention and treatment.

Epidemiological and functional characterization of Streptococcus pneumoniae pili

Streptococcus pneumoniae is an important life threatening human pathogen causing agent of invasive diseases such as otitis media, pneumonia, sepsis and meningitis, but is also a common inhabitant of the respiratory tract of children and healthy adults. Likewise most streptococci, S. pneumoniae decorates its surface with adhesive pili, composed of covalently linked subunits and involved in the attachment to epithelial cells and virulence. The pneumococcal pili are encoded by two genomic regions, pilus islet 1 (PI-1), and pilus islet-2 (PI-2), which are present in about 30% and 16% of the pneumococcal strains, respectively. PI-1 exists in three clonally related variants, whereas PI-2 is highly conserved. The presence of the islets does not correlate with the serotype of the strains, but with the genotype (as determined by Multi Locus Sequence Typing). The prevalence of PI-1 and PI-2 positive strains is similar in isolates from invasive disease and carriage. To better dissect a possible association between PIs presence and disease we evaluated the distribution of the two PIs in a panel of 113 acute otitis media (AOM) clinical isolates from Israel. PI-1 was present in 30.1% (N=34) of the isolates tested, and PI-2 in 7% (N=8). We found that 50% of the PI-1 positive isolates belonged to the international clones Spain9V-3 (ST156) and Taiwan19F-14 (ST236), and that PI-2 was not present in the absence of Pl-1. In conclusion, there was no correlation between PIs presence and AOM, and, in general, the observed differences in PIs prevalence are strictly dependent upon regional differences in the distribution of the clones. Finally, in the AOM collection the prevalence of PI-1 was higher among antibiotic resistant isolates, confirming previous indications obtained by the in silico analysis of the MLST database collection. Since the pilus-1 subunits were shown to confer protection in mouse models of infection both in active and passive immunization studies, and were regarded as potential candidates for a new generation of protein-based vaccines, the functional characterization was mainly focused on S. pneumoniae pilus -1 components. The pneumococcal pilus-1 is composed of three subunits, RrgA, RrgB and RrgC, each stabilized by intra-molecular isopeptide bonds and covalently polymerized by means of inter-molecular isopeptide bonds to form an extended fibre. The pilus shaft is a multimeric structure mainly composed by the RrgB backbone subunit. The minor ancillary proteins are located at the tip and at the base of the pilus, where they have been proposed to act as the major adhesin (RrgA) and as the pilus anchor (RrgC), respectively. RrgA is protective in in vivo mouse models, and exists in two variants (clades I and II). Mapping of the sequence variability onto the RrgA structure predicted from X-ray data showed that the diversity was restricted to the “head” of the protein, which contains the putative binding domains, whereas the elongated “stalk” was mostly conserved. To investigate whether this variability could influence the adhesive capacity of RrgA and to map the regions important for binding, two full-length protein variants and three recombinant RrgA portions were tested for adhesion to lung epithelial cells and to purified extracellular matrix (ECM) components. The two RrgA variants displayed similar binding abilities, whereas none of the recombinant fragments adhered at levels comparable to those of the full-length protein, suggesting that proper folding and structural arrangement are crucial to retain protein functionality. Furthermore, the two RrgA variants were shown to be cross-reactive in vitro and cross-protective in vivo in a murine model of passive immunization. Taken together, these data indicate that the region implicated in adhesion and the functional epitopes responsible for the protective ability of RrgA may be conserved and that the considerable level of variation found within the “head” domain of RrgA may have been generated by immunologic pressure without impairing the functional integrity of the pilus.

Bioinformatic methods in applied genomic research

Here I will focus on three main topics that best address and include the projects I have been working in during my three year PhD period that I have spent in different research laboratories addressing both computationally and practically important problems all related to modern molecular genomics. The first topic is the use of livestock species (pigs) as a model of obesity, a complex human dysfunction. My efforts here concern the detection and annotation of Single Nucleotide Polymorphisms. I developed a pipeline for mining human and porcine sequences. Starting from a set of human genes related with obesity the platform returns a list of annotated porcine SNPs extracted from a new set of potential obesity-genes. 565 of these SNPs were analyzed on an Illumina chip to test the involvement in obesity on a population composed by more than 500 pigs. Results will be discussed. All the computational analysis and experiments were done in collaboration with the Biocomputing group and Dr.Luca Fontanesi, respectively, under the direction of prof. Rita Casadio at the Bologna University, Italy. The second topic concerns developing a methodology, based on Factor Analysis, to simultaneously mine information from different levels of biological organization. With specific test cases we develop models of the complexity of the mRNA-miRNA molecular interaction in brain tumors measured indirectly by microarray and quantitative PCR. This work was done under the supervision of Prof. Christine Nardini, at the “CAS-MPG Partner Institute for Computational Biology” of Shangai, China (co-founded by the Max Planck Society and the Chinese Academy of Sciences jointly) The third topic concerns the development of a new method to overcome the variety of PCR technologies routinely adopted to characterize unknown flanking DNA regions of a viral integration locus of the human genome after clinical gene therapy. This new method is entirely based on next generation sequencing and it reduces the time required to detect insertion sites, decreasing the complexity of the procedure. This work was done in collaboration with the group of Dr. Manfred Schmidt at the Nationales Centrum für Tumorerkrankungen (Heidelberg, Germany) supervised by Dr. Annette Deichmann and Dr. Ali Nowrouzi. Furthermore I add as an Appendix the description of a R package for gene network reconstruction that I helped to develop for scientific usage (http://www.bioconductor.org/help/bioc-views/release/bioc/html/BUS.html).

Depicting the role of CDKN2A/ARF alterations in adult BCR-ABL1-positive acute lymphoblastic leukemia patients: from genomic deletions to prognostic impact

This 9p21 locus, encode for important proteins involved in cell cycle regulation and apoptosis containing the p16/CDKN2A (cyclin-dependent kinase inhibitor 2a) tumor suppressor gene and two other related genes, p14/ARF and p15/CDKN2B. This locus, is a major target of inactivation in the pathogenesis of a number of human tumors, both solid and haematologic, and is a frequent site of loss or deletion also in acute lymphoblastic leukemia (ALL) ranging from 18% to 45% 1. In order to explore, at high resolution, the frequency and size of alterations affecting this locus in adult BCR-ABL1-positive ALL and to investigate their prognostic value, 112 patients (101 de novo and 11 relapse cases) were analyzed by genome-wide single nucleotide polymorphisms arrays and gene candidate deep exon sequencing. Paired diagnosis-relapse samples were further available and analyzed for 19 (19%) cases. CDKN2A/ARF and CDKN2B genomic alterations were identified in 29% and 25% of newly diagnosed patients, respectively. Deletions were monoallelic in 72% of cases and in 43% the minimal overlapping region of the lost area spanned only the CDKN2A/2B gene locus. The analysis at the time of relapse showed an almost significant increase in the detection rate of CDKN2A/ARF loss (47%) compared to diagnosis (p = 0.06). Point mutations within the 9p21 locus were found at very low level with only a non-synonymous substition in the exon 2 of CDKN2A. Finally, correlation with clinical outcome showed that deletions of CDKN2A/B are significantly associated with poor outcome in terms of overall survival (p = 0.0206), disease free-survival (p = 0.0010) and cumulative incidence of relapse (p = 0.0014). The inactivation of 9p21 locus by genomic deletions is a frequent event in BCR-ABL1-positive ALL. Deletions are frequently acquired at the leukemia progression and work as a poor prognostic marker.

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.

Functional and Genetic characterization of new genomic islands from an E. coli strain associated with neonatal meningitis.

Enterobacteriaceae genomes evolve through mutations, rearrangements and horizontal gene transfer (HGT). The latter evolutionary pathway works through the acquisition DNA (GEI) modules of foreign origin that enhances fitness of the host to a given environment. The genome of E. coli IHE3034, a strain isolated from a case of neonatal meningitis, has recently been sequenced and its subsequent sequence analysis has predicted 18 possible GEIs, of which: 8 have not been previously described, 5 fully meet the pathogenic island definition and at least 10 that seem to be of prophagic origin. In order to study the GEI distribution of our reference strain, we screened for the presence 18 GEIs a panel of 132 strains, representative of E. coli diversity. Also, using an inverse nested PCR approach we identified 9 GEI that can form an extrachromosomal circular intermediate (CI) and their respective attachment sites (att). Further, we set up a qPCR approach that allowed us to determine the excision rates of 5 genomic islands in different growth conditions. Four islands, specific for strains appertaining to the sequence type complex 95 (STC95), have been deleted in order to assess their function in a Dictyostelium discoideum grazing assays. Overall, the distribution data presented here indicate that 16 IHE3034 GEIs are more associated to the STC95 strains. Also the functional and genetic characterization has uncovered that GEI 13, 17 and 19 are involved in the resistance to phagocitation by Dictyostelium d thus suggesting a possible role in the adaptation of the pathogen during certain stages of infection.