Deregulation of apoptosis-related genes is associated with PRV1 overexpression and JAK2 V617F allele burden in Essential Thrombocythemia and Myelofibrosis

Background: Essential Thrombocythemia (ET) and Primary Myelofibrosis (PMF) are Chronic Myeloproliferative Neoplasms (MPN) characterized by clonal myeloproliferation/myeloaccumulation without cell maturation impairment. The JAK2 V617F mutation and PRV1 gene overexpression may contribute to MPN physiopathology. We hypothesized that deregulation of the apoptotic machinery may also play a role in the pathogenesis of ET and PMF. In this study we evaluated the apoptosis-related gene and protein expression of BCL2 family members in bone marrow CD34(+) hematopoietic stem cells (HSC) and peripheral blood leukocytes from ET and PMF patients. We also tested whether the gene expression results were correlated with JAK2 V617F allele burden percentage, PRV1 overexpression, and clinical and laboratory parameters. Results: By real time PCR assay, we observed that A1, MCL1, BIK and BID, as well as A1, BCLW and BAK gene expression were increased in ET and PMF CD34(+) cells respectively, while pro-apoptotic BAX and anti-apoptotic BCL2 mRNA levels were found to be lower in ET and PMF CD34(+) cells respectively, in relation to controls. In patients' leukocytes, we detected an upregulation of anti-apoptotic genes A1, BCL2, BCL-XL and BCLW. In contrast, pro-apoptotic BID and BIMEL expression were downregulated in ET leukocytes. Increased BCL-XL protein expression in PMF leukocytes and decreased BID protein expression in ET leukocytes were observed by Western Blot. In ET leukocytes, we found a correlation between JAK2 V617F allele burden and BAX, BIK and BAD gene expression and between A1, BAX and BIK and PRV1 gene expression. A negative correlation between PRV1 gene expression and platelet count was observed, as well as a positive correlation between PRV1 gene expression and splenomegaly. Conclusions: Our results suggest the participation of intrinsic apoptosis pathway in the MPN physiopathology. In addition, PRV1 and JAK2 V617F allele burden were linked to deregulation of the apoptotic machinery.

Expression of Clock Genes in Human Subcutaneous and Visceral Adipose Tissues

Disrupted circadian rhythms are associated with obesity and metabolic alterations, but little is known about the participation of peripheral circadian clock machinery in these processes. The aim of the present study was to analyze RNA expression of clock genes in subcutaneous (SAT) and visceral (VAT) adipose tissues of male and female subjects in AM (morning) and PM (afternoon) periods, and its interactions with body mass index (BMI). Ninety-one subjects (41 +/- 11 yrs of age) presenting a wide range of BMI (21.4 to 48.6 kg/m(2)) were included. SAT and VAT biopsies were obtained from patients undergoing abdominal surgeries. Clock genes expressions were evaluated by qRT-PCR. The only clock gene that showed higher expression (p < .0001) in SAT in comparison to VAT was PER1 of female (372%) and male (326%) subjects. Different patterns of expression between the AM and PM periods were observed, in particular REV-ERBa, which was reduced (p < .05) at the PM period in SAT and VAT of both women and men (women: similar to 53% lower; men: similar to 78% lower), whereas CLOCK expression was not altered. Relationships between clock genes were different in SAT vs. VAT. BMI was negatively correlated with SATPER1 (r = -.549; p = .001) and SATPER2 (r = -.613; p = .0001) and positively with VATCLOCK (r = .541; p = .001) and VATBMAL1 (r = .468; p = .007) only in women. These data suggest that the circadian clock machinery of adipose tissue depots differs between female and male subjects, with a sex-specific effect observed for some genes. BMI correlated with clock genes, but at this moment it is not possible to establish the cause-effect relationship. (Author correspondence: mzanquetta@gmail.com)

Bosentan, an endothelin receptor antagonist, ameliorates collagen-induced arthritis: the role of TNF-alpha in the induction of endothelin system genes

Endothelins (ETs) are involved in several inflammatory events. The present study investigated the efficacy of bosentan, a dual ETA/ETB receptor antagonist, in collagen-induced arthritis (CIA) in mice. CIA was induced in DBA/1J mice. Arthritic mice were treated with bosentan (100 mg/kg) once a day, starting from the day when arthritis was clinically detectable. CIA progression was assessed by measurements of visual clinical score, paw swelling and hypernociception. Histological changes, neutrophil infiltration and pro-inflammatory cytokines were evaluated in the joints. Gene expression in the lymph nodes of arthritic mice was evaluated by microarray technology. PreproET-1 mRNA expression in the lymph nodes of mice and in peripheral blood mononuclear cells (PBMCs) was evaluated by real-time PCR. The differences were evaluated by one-way ANOVA or Student's t test. Oral treatment with bosentan markedly ameliorated the clinical aspects of CIA (visual clinical score, paw swelling and hyperalgesia). Bosentan treatment also reduced joint damage, leukocyte infiltration and pro-inflammatory cytokine levels (IL-1 beta, TNF alpha and IL-17) in the joint tissues. Changes in gene expression in the lymph nodes of arthritic mice returned to the levels of the control mice after bosentan treatment. PreproET mRNA expression increased in PBMCs from rheumatoid arthritis (RA) patients but returned to basal level in PBMCs from patients under anti-TNF therapy. In-vitro treatment of PBMCs with TNF alpha upregulated ET system genes. These findings indicate that ET receptor antagonists, such as bosentan, might be useful in controlling RA. Moreover, it seems that ET mediation of arthritis is triggered by TNF alpha.

Exploring the Schistosoma mansoni adult male transcriptome using RNA-seq

Schistosoma mansoni is one of the agents of schistosomiasis, a chronic and debilitating disease. Here we, present a transcriptome-wide characterization of adult S. mansoni males by high-throughput RNA-sequencing. We obtained 1,620,432 high-quality ESTs from a directional strand-specific cDNA library, resulting in a 26% higher coverage of genome bases than that of the public ESTs available at NCBI. With a 15 x-deep coverage of transcribed genomic regions, our data were able to (i) confirm for the first time 990 predictions without previous evidence of transcription; (ii) correct gene predictions; (iii) discover 989 and 1196 RNA-seq contigs that map to intergenic and intronic genomic regions, respectively, where no gene had been predicted before. These contigs could represent new protein-coding genes or non-coding RNAs (ncRNAs). Interestingly, we identified 11 novel Micro-exon genes (MEGs). These data reveal new features of the S. mansoni transcriptional landscape and significantly advance our understanding of the parasite transcriptome. (c) 2011 Elsevier Inc. All rights reserved.

Validation of reference genes for RT-qPCR normalization in common bean during biotic and abiotic stresses

Selection of reference genes is an essential consideration to increase the precision and quality of relative expression analysis by the quantitative RT-PCR method. The stability of eight expressed sequence tags was evaluated to define potential reference genes to study the differential expression of common bean target genes under biotic (incompatible interaction between common bean and fungus Colletotrichum lindemuthianum) and abiotic (drought; salinity; cold temperature) stresses. The efficiency of amplification curves and quantification cycle (C (q)) were determined using LinRegPCR software. The stability of the candidate reference genes was obtained using geNorm and NormFinder software, whereas the normalization of differential expression of target genes [beta-1,3-glucanase 1 (BG1) gene for biotic stress and dehydration responsive element binding (DREB) gene for abiotic stress] was defined by REST software. High stability was obtained for insulin degrading enzyme (IDE), actin-11 (Act11), unknown 1 (Ukn1) and unknown 2 (Ukn2) genes during biotic stress, and for SKP1/ASK-interacting protein 16 (Skip16), Act11, Tubulin beta-8 (beta-Tub8) and Unk1 genes under abiotic stresses. However, IDE and Act11 were indicated as the best combination of reference genes for biotic stress analysis, whereas the Skip16 and Act11 genes were the best combination to study abiotic stress. These genes should be useful in the normalization of gene expression by RT-PCR analysis in common bean, the most important edible legume.

Different patterns of evolution for duplicated DNA repair genes in bacteria of the Xanthomonadales group

Abstract Background DNA repair genes encode proteins that protect organisms against genetic damage generated by environmental agents and by-products of cell metabolism. The importance of these genes in life maintenance is supported by their high conservation, and the presence of duplications of such genes may be easily traced, especially in prokaryotic genomes. Results The genome sequences of two Xanthomonas species were used as the basis for phylogenetic analyses of genes related to DNA repair that were found duplicated. Although 16S rRNA phylogenetic analyses confirm their classification at the basis of the gamma proteobacteria subdivision, differences were found in the origin of the various genes investigated. Except for lexA, detected as a recent duplication, most of the genes in more than one copy are represented by two highly divergent orthologs. Basically, one of such duplications is frequently positioned close to other gamma proteobacteria, but the second is often positioned close to unrelated bacteria. These orthologs may have occurred from old duplication events, followed by extensive gene loss, or were originated from lateral gene transfer (LGT), as is the case of the uvrD homolog. Conclusions Duplications of DNA repair related genes may result in redundancy and also improve the organisms' responses to environmental challenges. Most of such duplications, in Xanthomonas, seem to have arisen from old events and possibly enlarge both functional and evolutionary genome potentiality.

Genome analysis of DNA repair genes in the alpha proteobacterium Caulobacter crescentus

Abstract Background The integrity of DNA molecules is fundamental for maintaining life. The DNA repair proteins protect organisms against genetic damage, by removal of DNA lesions or helping to tolerate them. DNA repair genes are best known from the gamma-proteobacterium Escherichia coli, which is the most understood bacterial model. However, genome sequencing raises questions regarding uniformity and ubiquity of these DNA repair genes and pathways, reinforcing the need for identifying genes and proteins, which may respond to DNA damage in other bacteria. Results In this study, we employed a bioinformatic approach, to analyse and describe the open reading frames potentially related to DNA repair from the genome of the alpha-proteobacterium Caulobacter crescentus. This was performed by comparison with known DNA repair related genes found in public databases. As expected, although C. crescentus and E. coli bacteria belong to separate phylogenetic groups, many of their DNA repair genes are very similar. However, some important DNA repair genes are absent in the C. crescentus genome and other interesting functionally related gene duplications are present, which do not occur in E. coli. These include DNA ligases, exonuclease III (xthA), endonuclease III (nth), O6-methylguanine-DNA methyltransferase (ada gene), photolyase-like genes, and uracil-DNA-glycosylases. On the other hand, the genes imuA and imuB, which are involved in DNA damage induced mutagenesis, have recently been described in C. crescentus, but are absent in E. coli. Particularly interesting are the potential atypical phylogeny of one of the photolyase genes in alpha-proteobacteria, indicating an origin by horizontal transfer, and the duplication of the Ada orthologs, which have diverse structural configurations, including one that is still unique for C. crescentus. Conclusion The absence and the presence of certain genes are discussed and predictions are made considering the particular aspects of the C. crescentus among other known DNA repair pathways. The observed differences enlarge what is known for DNA repair in the Bacterial world, and provide a useful framework for further experimental studies in this organism.

Expression profiling and cross-species RNA interference (RNAi) of desiccation-induced transcripts in the anhydrobiotic nematode Aphelenchus avenae

Abstract Background Some organisms can survive extreme desiccation by entering a state of suspended animation known as anhydrobiosis. The free-living mycophagous nematode Aphelenchus avenae can be induced to enter anhydrobiosis by pre-exposure to moderate reductions in relative humidity (RH) prior to extreme desiccation. This preconditioning phase is thought to allow modification of the transcriptome by activation of genes required for desiccation tolerance. Results To identify such genes, a panel of expressed sequence tags (ESTs) enriched for sequences upregulated in A. avenae during preconditioning was created. A subset of 30 genes with significant matches in databases, together with a number of apparently novel sequences, were chosen for further study. Several of the recognisable genes are associated with water stress, encoding, for example, two new hydrophilic proteins related to the late embryogenesis abundant (LEA) protein family. Expression studies confirmed EST panel members to be upregulated by evaporative water loss, and the majority of genes was also induced by osmotic stress and cold, but rather fewer by heat. We attempted to use RNA interference (RNAi) to demonstrate the importance of this gene set for anhydrobiosis, but found A. avenae to be recalcitrant with the techniques used. Instead, therefore, we developed a cross-species RNAi procedure using A. avenae sequences in another anhydrobiotic nematode, Panagrolaimus superbus, which is amenable to gene silencing. Of 20 A. avenae ESTs screened, a significant reduction in survival of desiccation in treated P. superbus populations was observed with two sequences, one of which was novel, while the other encoded a glutathione peroxidase. To confirm a role for glutathione peroxidases in anhydrobiosis, RNAi with cognate sequences from P. superbus was performed and was also shown to reduce desiccation tolerance in this species. Conclusions This study has identified and characterised the expression profiles of members of the anhydrobiotic gene set in A. avenae. It also demonstrates the potential of RNAi for the analysis of anhydrobiosis and provides the first genetic data to underline the importance of effective antioxidant systems in metazoan desiccation tolerance.

Disruption of vitellogenin gene function in adult honeybees by intra-abdominal injection of double-stranded RNA

Abstract Background The ability to manipulate the genetic networks underlying the physiological and behavioural repertoires of the adult honeybee worker (Apis mellifera) is likely to deepen our understanding of issues such as learning and memory generation, ageing, and the regulatory anatomy of social systems in proximate as well as evolutionary terms. Here we assess two methods for probing gene function by RNA interference (RNAi) in adult honeybees. Results The vitellogenin gene was chosen as target because its expression is unlikely to have a phenotypic effect until the adult stage in bees. This allowed us to introduce dsRNA in preblastoderm eggs without affecting gene function during development. Of workers reared from eggs injected with dsRNA derived from a 504 bp stretch of the vitellogenin coding sequence, 15% had strongly reduced levels of vitellogenin mRNA. When dsRNA was introduced by intra-abdominal injection in newly emerged bees, almost all individuals (96 %) showed the mutant phenotype. An RNA-fragment with an apparent size similar to the template dsRNA was still present in this group after 15 days. Conclusion Injection of dsRNA in eggs at the preblastoderm stage seems to allow disruption of gene function in all developmental stages. To dissect gene function in the adult stage, the intra-abdominal injection technique seems superior to egg injection as it gives a much higher penetrance, it is much simpler, and it makes it possible to address genes that are also expressed in the embryonic, larval or pupal stages.

Retroviral transfer of the p16INK4a cDNA inhibits C6 glioma formation in Wistar rats

Abstract Background The p16INK4A gene product halts cell proliferation by preventing phosphorylation of the Rb protein. The p16INK4a gene is often deleted in human glioblastoma multiforme, contributing to unchecked Rb phosphorylation and rapid cell division. We show here that transduction of the human p16INK4a cDNA using the pCL retroviral system is an efficient means of stopping the proliferation of the rat-derrived glioma cell line, C6, both in tissue culture and in an animal model. C6 cells were transduced with pCL retrovirus encoding the p16INK4a, p53, or Rb genes. These cells were analyzed by a colony formation assay. Expression of p16INK4a was confirmed by immunohistochemistry and Western blot analysis. The altered morphology of the p16-expressing cells was further characterized by the senescence-associated β-galactosidase assay. C6 cells infected ex vivo were implanted by stereotaxic injection in order to assess tumor formation. Results The p16INK4a gene arrested C6 cells more efficiently than either p53 or Rb. Continued studies with the p16INK4a gene revealed that a large portion of infected cells expressed the p16INK4a protein and the morphology of these cells was altered. The enlarged, flat, and bi-polar shape indicated a senescence-like state, confirmed by the senescence-associated β-galactosidase assay. The animal model revealed that cells infected with the pCLp16 virus did not form tumors. Conclusion Our results show that retrovirus mediated transfer of p16INK4a halts glioma formation in a rat model. These results corroborate the idea that retrovirus-mediated transfer of the p16INK4a gene may be an effective means to arrest human glioma and glioblastoma.

Identification of protein-coding and non-coding RNA expression profiles in CD34+ and in stromal cells in refractory anemia with ringed sideroblasts

Abstract Background Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods In this study, gene expression profiles of CD34+ cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results In CD34+ cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value ≤ 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value ≤ 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34+ cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

Abstract Background RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known. Results A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages. Conclusion We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

Storage of bovine reproductive tissues and RNA extracts on ice for 24h or repeated freeze–thaw cycles do not affect RNA integrity

The aims of this study were to test (i) the effect of time of tissue and RNA extracts storage on ice and (ii) the effect of repeated freeze–thaw cycles on RNA integrity and gene expression of bovine reproductive tissues. Fragments of endometrium (ENDO), corpus luteum (CL) and ampulla (AMP) were subdivided and incubated for 0, 1, 3, 6, 12 or 24 h on ice. RNA extracts were incubated on ice for 0, 3, 12 or 24 h, or exposed to 1, 2, 4 or 6 freeze–thaw cycles. RNA integrity number (RIN) was estimated. Expression of progesterone receptor (PGR) and cyclophilin genes from RNA extracts stored on ice for 0 or 24 h, and 1 or 6 freeze–thaw cycles was measured by qPCR. Tissue and RNA extract incubation on ice, and repeated freeze–thaw cycles did not affect RIN values of RNA from ENDO, CL or AMP. Storage on ice or exposure to freeze–thaw cycles did not affect Cq values for PGR or cyclophilin genes. In conclusion, neither generalized RNA degradation nor specific RNA degradation was affected by storage of tissue or RNA extracts on ice for up to 24 h, or by up to 6 freeze–thaw cycles of RNA extracts obtained from bovine ENDO, CL and AMP.

The intronic long noncoding RNA ANRASSF1 recruits PRC2 to the RASSF1A promoter, reducing the expression of RASSF1A and increasing cell proliferation

The down-regulation of the tumor-suppressor gene RASSF1A has been shown to increase cell proliferation in several tumors. RASSF1A expression is regulated through epigenetic events involving the polycomb repressive complex 2 (PRC2); however, the molecular mechanisms modulating the recruitment of this epigenetic modifier to the RASSF1 locus remain largely unknown. Here, we identify and characterize ANRASSF1, an endogenous unspliced long noncoding RNA (lncRNA) that is transcribed from the opposite strand on the RASSF1 gene locus in several cell lines and tissues and binds PRC2. ANRASSF1 is transcribed through RNA polymerase II and is 5'-capped and polyadenylated; it exhibits nuclear localization and has a shorter half-life compared with other lncRNAs that bind PRC2. ANRASSF1 endogenous expression is higher in breast and prostate tumor cell lines compared with non-tumor, and an opposite pattern is observed for RASSF1A. ANRASSF1 ectopic overexpression reduces RASSF1A abundance and increases the proliferation of HeLa cells, whereas ANRASSF1 silencing causes the opposite effects. These changes in ANRASSF1 levels do not affect the RASSF1C isoform abundance. ANRASSF1 overexpression causes a marked increase in both PRC2 occupancy and histone H3K27me3 repressive marks, specifically at the RASSF1A promoter region. No effect of ANRASSF1 overexpression was detected on PRC2 occupancy and histone H3K27me3 at the promoter regions of RASSF1C and the four other neighboring genes, including two well-characterized tumor suppressor genes. Additionally, we demonstrated that ANRASSF1 forms an RNA/DNA hybrid and recruits PRC2 to the RASSF1A promoter. Together, these results demonstrate a novel mechanism of epigenetic repression of the RASSF1A tumor suppressor gene involving antisense unspliced lncRNA, in which ANRASSF1 selectively represses the expression of the RASSF1 isoform overlapping the antisense transcript in a location-specific manner. In a broader perspective, our findings suggest that other non-characterized unspliced intronic lncRNAs transcribed in the human genome might contribute to a location-specific epigenetic modulation of genes.

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.