Molecular strategies for genetic diversity analysis and development of markers linked to resistance traits in apple

The goal of many plant scientists’ research is to explain natural phenotypic variation in term of simple changes in DNA sequence. DNA-based molecular markers are extensively used for the construction of genome-wide molecular maps and to perform genetic analysis for simple and complex traits. The PhD thesis was divided into two main research lines according to the different approaches adopted. The first research line is to analyze the genetic diversity in an Italian apple germplasm collection for the identification of markers tightly linked to targeted genes by an association genetic method. This made it possible to identify synomym and homonym accessions and triploids. The fruit red skin color trait has been used to test the reliability of the genetic approaches in this species. The second line is related to the development of molecular markers closely linked to the Rvi13 and Rvi5 scab resistance genes, previously mapped on apple’s chromosome 10 and 17 respectively by using the traditional linkage mapping method. Both region have been fine-mapped with various type of markers that could be used for marker-assisted selection in future breeding programs and to isolate the two resistance genes.

Gwas analysis for the identification of molecular basis responsible for yield related traits and disease resistance in durum wheat

The PhD thesis was developed in the framework of Innovar H2020 project. This project aimed at using genomics, transcriptomics and phenotyping techniques to update varietal registration procedure used in Europe for Value of Cultivation and Use (VCU) and Distinctiness Uniformity and Stability (DUS) protocols. The phenotypic and genotypic diversity of a durum wheat panel were assessed for different agronomic traits, connected with wheat development, disease resistance and spike fertility. A panel of 253 durum wheat varieties was characterized for VCU and DUS traits and genotyped with Illumina 90K SNP Chip array (Wang et al., 2014). GWAS analysis was performed, detecting strong QTLs confirmed also by literature review. Candidate genes were identified for each trait and molecular markers will be developed to be used for marker assisted selection in breeding programs. As for disease resistance, the panel was evaluated for resistance to Soil-Borne-Cereal-Mosaic-Virus (SBCMV). A major QTL, sbm2, was detected on chromosome 2B responsible for durum wheat resistance (Maccaferri et al., 2011). The sbm2 interval was explored by fine mapping on segregant population using KASP markers and by RNASeq analysis, detecting candidate genes involved in plant-pathogen reaction. As regards yield related traits, detailed analysis was performed on the GNI-2A QTL (Milner et al., 2016), responsible for increased number spike fertility. Fine mapping analysis was performed on durum panel identifying hox2 a strong candidate gene, codifying for transcription factor protein. The gene is paralogue of GNI-1 (Sakuma et al., 2019), and it has a 4 kbp deletion responsible for increased number of florets per spikelet. To conclude, the herein reported thesis shows a complete characterization of agronomic and disease resistance traits in modern durum wheat varieties. The results obtained will augment available information for each variety, identifying informative molecular markers for breeding purposes and QTLs/candidate genes responsible for different agronomic traits.

Identificazione di un QTL principale per resistenza a ruggine bruna sul cromosoma 7B di frumento duro

Leaf rust caused by Puccinia triticina is a serious disease of durum wheat (Triticum durum) worldwide. However, genetic and molecular mapping studies aimed at characterizing leaf rust resistance genes in durum wheat have been only recently undertaken. The Italian durum wheat cv. Creso shows a high level of resistance to P. triticina that has been considered durable and that appears to be due to a combination of a single dominant gene and one or more additional factors conferring partial resistance. In this study, the genetic basis of leaf rust resistance carried by Creso was investigated using 176 recombinant inbred lines (RILs) from the cross between the cv. Colosseo (C, leaf rust resistance donor) and Lloyd (L, susceptible parent). Colosseo is a cv. directly related to Creso with the leaf rust resistance phenotype inherited from Creso, and was considered as resistance donor because of its better adaptation to local (Emilia Romagna, Italy) cultivation environment. RILs have been artificially inoculated with a mixture of 16 Italian P. triticina isolates that were characterized for virulence to seedlings of 22 common wheat cv. Thatcher isolines each carrying a different leaf rust resistance gene, and for molecular genotypes at 15 simple sequence repeat (SSR) loci, in order to determine their specialization with regard to the host species. The characterization of the leaf rust isolates was conducted at the Cereal Disease Laboratory of the University of Minnesota (St. Paul, USA) (Chapter 2). A genetic linkage map was constructed using segregation data from the population of 176 RILs from the cross CL. A total of 662 loci, including 162 simple sequence repeats (SSRs) and 500 Diversity Arrays Technology markers (DArTs), were analyzed by means of the package EasyMap 0.1. The integrated SSR-DArT linkage map consisted of 554 loci (162 SSR and 392 DArT markers) grouped into 19 linkage blocks with an average marker density of 5.7 cM/marker. The final map spanned a total of 2022 cM, which correspond to a tetraploid genome (AABB) coverage of ca. 77% (Chapter 3). The RIL population was phenotyped for their resistance to leaf rust under artificial inoculation in 2006; the percentage of infected leaf area (LRS, leaf rust susceptibility) was evaluated at three stages through the disease developmental cycle and the area under disease progress curve (AUDPC) was then calculated. The response at the seedling stage (infection type, IT) was also investigated. QTL analysis was carried out by means of the Composite Interval Mapping method based on a selection of markers from the CL map. A major QTL (QLr.ubo-7B.2) for leaf rust resistance controlling both the seedling and the adult plant response, was mapped on the distal region of chromosome arm 7BL (deletion bin 7BL10-0.78-1.00), in a gene-dense region known to carry several genes/QTLs for resistance to rusts and other major cereal fungal diseases in wheat and barley. QLr.ubo-7B.2 was identified within a supporting interval of ca. 5 cM tightly associated with three SSR markers (Xbarc340.2, Xgwm146 e Xgwm344.2), and showed an R2 and an LOD peak value for the AUDPC equal to 72.9% an 44.5, respectively. Three additional minor QTLs were also detected (QLr.ubo-7B.1 on chr. 7BS; QLr.ubo-2A on chr. 2AL and QLr.ubo-3A on chr. 3AS) (Chapter 4). The presence of the major QTL (QLr.ubo-7B.2) was validated by a linkage disequilibrium (LD)-based test using field data from two different plant materials: i) a set of 62 advanced lines from multiple crosses involving Creso and his directly related resistance derivates Colosseo and Plinio, and ii) a panel of 164 elite durum wheat accessions representative of the major durum breeding program of the Mediterranean basin. Lines and accessions were phenotyped for leaf rust resistance under artificial inoculation in two different field trials carried out at Argelato (BO, Italy) in 2006 and 2007; the durum elite accessions were also evaluated in two additional field experiments in Obregon (Messico; 2007 and 2008) and in a green-house experiment (seedling resistance) at the Cereal Disease Laboratory (St. Paul, USA, 2008). The molecular characterization involved 14 SSR markers mapping on the 7BL chromosome region found to harbour the major QTL. Association analysis was then performed with a mixed-linear-model approach. Results confirmed the presence of a major QTL for leaf rust resistance, both at adult plant and at seedling stage, located between markers Xbarc340.2, Xgwm146 and Xgwm344.2, in an interval that coincides with the supporting interval (LOD-2) of QLr.ubo-7B.2 as resulted from the RIL QTL analysis. (Chapter 5). The identification and mapping of the major QTL associated to the durable leaf rust resistance carried by Creso, together with the identification of the associated SSR markers, will enhance the selection efficiency in durum wheat breeding programs (MAS, Marker Assisted Selection) and will accelerate the release of cvs. with durable resistance through marker-assisted pyramiding of the tagged resistance genes/QTLs most effective against wheat fungal pathogens.

Different approaches to identify markers for meat quality and production traits in pigs: analysis of the transcriptome in post mortem skeletal muscle and investigation of candidate genes

Due to the growing attention of consumers towards their food, improvement of quality of animal products has become one of the main focus of research. To this aim, the application of modern molecular genetics approaches has been proved extremely useful and effective. This innovative drive includes all livestock species productions, including pork. The Italian pig breeding industry is unique because needs heavy pigs slaughtered at about 160 kg for the production of high quality processed products. For this reason, it requires precise meat quality and carcass characteristics. Two aspects have been considered in this thesis: the application of the transcriptome analysis in post mortem pig muscles as a possible method to evaluate meat quality parameters related to the pre mortem status of the animals, including health, nutrition, welfare, and with potential applications for product traceability (chapters 3 and 4); the study of candidate genes for obesity related traits in order to identify markers associated with fatness in pigs that could be applied to improve carcass quality (chapters 5, 6, and 7). Chapter three addresses the first issue from a methodological point of view. When we considered this issue, it was not obvious that post mortem skeletal muscle could be useful for transcriptomic analysis. Therefore we demonstrated that the quality of RNA extracted from skeletal muscle of pigs sampled at different post mortem intervals (20 minutes, 2 hours, 6 hours, and 24 hours) is good for downstream applications. Degradation occurred starting from 48 h post mortem even if at this time it is still possible to use some RNA products. In the fourth chapter, in order to demonstrate the potential use of RNA obtained up to 24 hours post mortem, we present the results of RNA analysis with the Affymetrix microarray platform that made it possible to assess the level of expression of more of 24000 mRNAs. We did not identify any significant differences between the different post mortem times suggesting that this technique could be applied to retrieve information coming from the transcriptome of skeletal muscle samples not collected just after slaughtering. This study represents the first contribution of this kind applied to pork. In the fifth chapter, we investigated as candidate for fat deposition the TBC1D1 [TBC1 (tre-2/USP6, BUB2, cdc16) gene. This gene is involved in mechanisms regulating energy homeostasis in skeletal muscle and is associated with predisposition to obesity in humans. By resequencing a fragment of the TBC1D1 gene we identified three synonymous mutations localized in exon 2 (g.40A>G, g.151C>T, and g.172T>C) and 2 polymorphisms localized in intron 2 (g.219G>A and g.252G>A). One of these polymorphisms (g.219G>A) was genotyped by high resolution melting (HRM) analysis and PCR-RFLP. Moreover, this gene sequence was mapped by radiation hybrid analysis on porcine chromosome 8. The association study was conducted in 756 performance tested pigs of Italian Large White and Italian Duroc breeds. Significant results were obtained for lean meat content, back fat thickness, visible intermuscular fat and ham weight. In chapter six, a second candidate gene (tribbles homolog 3, TRIB3) is analyzed in a study of association with carcass and meat quality traits. The TRIB3 gene is involved in energy metabolism of skeletal muscle and plays a role as suppressor of adipocyte differentiation. We identified two polymorphisms in the first coding exon of the porcine TRIB3 gene, one is a synonymous SNP (c.132T> C), a second is a missense mutation (c.146C> T, p.P49L). The two polymorphisms appear to be in complete linkage disequilibrium between and within breeds. The in silico analysis of the p.P49L substitution suggests that it might have a functional effect. The association study in about 650 pigs indicates that this marker is associated with back fat thickness in Italian Large White and Italian Duroc breeds in two different experimental designs. This polymorphisms is also associated with lactate content of muscle semimembranosus in Italian Large White pigs. Expression analysis indicated that this gene is transcribed in skeletal muscle and adipose tissue as well as in other tissues. In the seventh chapter, we reported the genotyping results for of 677 SNPs in extreme divergent groups of pigs chosen according to the extreme estimated breeding values for back fat thickness. SNPs were identified by resequencing, literature mining and in silico database mining. analysis, data reported in the literature of 60 candidates genes for obesity. Genotyping was carried out using the GoldenGate (Illumina) platform. Of the analyzed SNPs more that 300 were polymorphic in the genotyped population and had minor allele frequency (MAF) >0.05. Of these SNPs, 65 were associated (P<0.10) with back fat thickness. One of the most significant gene marker was the same TBC1D1 SNPs reported in chapter 5, confirming the role of this gene in fat deposition in pig. These results could be important to better define the pig as a model for human obesity other than for marker assisted selection to improve carcass characteristics.

A candidate gene approach to identify DNA markers associated with meat quality and production traits: investigation of several cathepsin genes in Italian heavy pigs.

The cathepsin enzymes represent an important family of lysosomal proteinases with a broad spectrum of functions in many, if not in all, tissues and cell types. In addition to their primary role during the normal protein turnover, they possess highly specific proteolytic activities, including antigen processing in the immune response and a direct role in the development of obesity and tumours. In pigs, the involvement of cathepsin enzymes in proteolytic processes have important effects during the conversion of muscle to meat, due to their influence on meat texture and sensory characteristics, mainly in seasoned products. Their contribution is fundamental in flavour development of dry-curing hams. However, several authors have demonstrated that high cathepsin activity, in particular of cathepsin B, is correlated to defects of these products, such as an excessive meat softness together with abnormal free tyrosine content, astringent or metallic aftertastes and formation of a white film on the cut surface. Thus, investigation of their genetic variability could be useful to identify DNA markers associated with these dry cured hams parameters, but also with meat quality, production and carcass traits in Italian heavy pigs. Unfortunately, no association has been found between cathepsin markers and meat quality traits so far, in particular with cathepsin B activity, suggesting that other genes, besides these, affect meat quality parameters. Nevertheless, significant associations were observed with several carcass and production traits in pigs. A recent study has demonstrated that different single nucleotide polymorphisms (SNPs) localized in cathepsin D (CTSD), F (CTSF), H and Z genes were highly associated with growth, fat deposition and production traits in an Italian Large White pig population. The aim of this thesis was to confirm some of these results in other pig populations and identify new cathepsin markers in order to evaluate their effects on cathepsin activity and other production traits. Furthermore, starting from the data obtained in previous studies on CTSD gene, we also analyzed the known polymorphism located in the insulin-like growth factor 2 gene (IGF2 intron3-g.3072G>A). This marker is considered the causative mutation for the quantitative trait loci (QTL) affecting muscle mass and fat deposition in pigs. Since IGF2 maps very close to CTSD on porcine chromosome (SSC) 2, we wanted to clarify if the effects of the CTSD marker were due to linkage disequilibrium with the IGF2 intron3-g.3072G>A mutation or not. In the first chapter, we reported the results from these two SSC2 gene markers. First of all, we evaluated the effects of the IGF2 intron3-g.3072G>A polymorphism in the Italian Large White breed, for which no previous studies have analysed this marker. Highly significant associations were identified with all estimated breeding values for production and carcass traits (P<0.00001), while no effects were observed for meat quality traits. Instead, the IGF2 intron3-g.3072G>A mutation did not show any associations with the analyzed traits in the Italian Duroc pigs, probably due to the low level of variability at this polymorphic site for this breed. In the same Duroc pig population, significant associations were obtained for the CTSD marker for all production and carcass traits (P < 0.001), after excluding possible confounding effects of the IGF2 mutation. The effects of the CTSD g.70G>A polymorphism were also confirmed in a group of Italian Large White pigs homozygous for the IGF2 intron3-g.3072G allele G (IGF2 intron3-g.3072GG) and by haplotype analysis between the markers of the two considered genes. Taken together, all these data indicated that the IGF2 intron3-g.3072G>A mutation is not the only polymorphism affecting fatness and muscle deposition in pigs. In the second chapter, we reported the analysis of two new SNPs identified in cathepsin L (CTSL) and cathepsin S (CTSS) genes and the association results with meat quality parameters (including cathepsin B activity) and several production traits in an Italian Large White pig population. Allele frequencies of these two markers were evaluated in 7 different pig breeds. Furthermore, we mapped using a radiation hybrid panel the CTSS gene on SSC4. Association studies with several production traits, carried out in 268 Italian Large White pigs, indicated positive effects of the CTSL polymorphism on average daily gain, weight of lean cuts and backfat thickness (P<0.05). The results for these latter traits were also confirmed using a selective genotype approach in other Italian Large White pigs (P<0.01). In the 268 pig group, the CTSS polymorphism was associated with feed:gain ratio and average daily gain (P<0.05). Instead, no association was observed between the analysed markers and meat quality parameters. Finally, we wanted to verify if the positive results obtained for the cathepsin L and S markers and for other previous identified SNPs (cathepsin F, cathepsin Z and their inhibitor cystatin B) were confirmed in the Italian Duroc pig breed (third chapter). We analysed them in two groups of Duroc pigs: the first group was made of 218 performance-tested pigs not selected by any phenotypic criteria, the second group was made of 100 Italian Duroc pigs extreme and divergent for visible intermuscular fat trait. In the first group, the CTSL polymorphism was associated with weight of lean cuts (P<0.05), while suggestive associations were obtained for average daily gain and backfat thickness (P<0.10). Allele frequencies of the CTSL gene marker also differed positively among the visible intermuscular extreme tails. Instead, no positive effects were observed for the other DNA markers on the analysed traits. In conclusion, in agreement with the present data and for the biological role of these enzymes, the porcine CTSD and CTSL markers: a) may have a direct effect in the biological mechanisms involved in determining fat and lean meat content in pigs, or b) these markers could be very close to the putative functional mutation(s) present in other genes. These findings have important practical applications, in particular the CTSD and CTSL mutations could be applied in a marker assisted selection (MAS) both in the Italian Large White and Italian Duroc breeds. Marker assisted selection could also increase in efficiency by adding information from the cathepsin S genotype, but only in the Italian Large White breed.

Candidate genes affecting fat deposition, carcass composition and meat quality traits in pigs

Pig meat quality is determined by several parameters, such as lipid content, tenderness, water-holding capacity, pH, color and flavor, that affect consumers’ acceptance and technological properties of meat. Carcass quality parameters are important for the production of fresh and dry-cure high-quality products, in particular the fat deposition and the lean cut yield. The identification of genes and markers associated with meat and carcass quality traits is of prime interest, for the possibility of improving the traits by marker-assisted selection (MAS) schemes. Therefore, the aim of this thesis was to investigate seven candidate genes for meat and carcass quality traits in pigs. In particular, we focused on genes belonging to the family of the lipid droplet coat proteins perilipins (PLIN1 and PLIN2) and to the calpain/calpastatin system (CAST, CAPN1, CAPN3, CAPNS1) and on the gene encoding for PPARg-coactivator 1A (PPARGC1A). In general, the candidate genes investigation included the protein localization, the detection of polymorphisms, the association analysis with meat and carcass traits and the analysis of the expression level, in order to assess the involvement of the gene in pork quality. Some of the analyzed genes showed effects on various pork traits that are subject to selection in genetic improvement programs, suggesting a possible involvement of the genes in controlling the traits variability. In particular, significant association results have been obtained for PLIN2, CAST and PPARGC1A genes, that are worthwhile of further validation. The obtained results contribute to a better understanding of biological mechanisms important for pig production as well as for a possible use of pig as animal model for studies regarding obesity in humans.

Dissecting the Heat-Stress Tolerance QTLome in Durum Wheat

Heat stress negatively affects wheat performance during its entire cycle, particularly during the reproductive stage. In view of the climate change and the prediction of a continued increase in temperature in the new future, it is urgent to concentrate efforts to discover novel genetic sources able to improve the resilience of wheat to heat stress. In this direction, this study addressed two different experiments in durum wheat to identify novel QTLs suitable to be applied in marker-assisted selection for heat tolerance. Chlorophyll fluorescence (ChlF) is a valuable indicator of plant response to environmental changes allowing a detailed assessment of PSII activity in view of its non-invasive measurement and high-throughput phenotyping. In the first study (Chapter 2), the Light-Induced Fluorescence Transient (LIFT) method was used to access ChlF data to map QTLs for ChlF-related traits during the vegetative growth stage in durum wheat under heat stress condition. Our results provide evidence that LIFT consistently measures ChlF at the level of high-throughput phenotyping combined with high accuracy which is required for Genome-Wide Association Study (GWAS) aimed at identifying genomic regions affecting PSII activity. The 50 QTLs identified for ChlF-related traits under heat stress mostly clustered into five chromosomes hotspots unrelated to phenology, a feature that makes these QTLs a valuable asset for marker-assisted breeding programs across different latitudes. In the second study (Chapter 3), a set of 183 accessions suitable for GWAS, was exposed to optimal and high temperature during two crop seasons under field conditions. Important agronomic traits were evaluated in order to identify valuable QTLs for GY and its components. The GWAS analysis identified several QTLs in the single years as well as in the joint analysis. From the total QTLs identified, 13 QTL clusters can be highlighted to be affecting heat tolerance across different years and/or different traits.

Exploring mediterranean durum wheat wild relatives and elite diversity panel for future breeding

Durum wheat (Triticum durum) is an important crop that has been used for millennia for human consumption, and modern breeding can take advantage of the wide variability useful for the adaptation to new challenges. Novel beneficial alleles can be found in wild relatives and landraces thus enhancing crop adaptation to many biotic and abiotic stresses. This dissertation considers the source of variability from both before and after wheat domestication, that caused a loss of potentially useful alleles. Chapter 1. is the thesis introduction, which outlines the importance of wheat in the world, providing an historical overview of the domestication, the evolution mechanisms that led to the current forms of durum wheat and the use of wild relatives as a source of germplasm for future breeding programs is crucial. Moreover, the emergence of Z. tritici has been considered as the main pathogen of wheat since it contains extremely high levels of genetic variability and is thus difficult to control. Chapter 2. Considers the contribution of the phenotypic diversity of 242 accessions of Aegilops tauschii from the Open Wild Wheat Consortium, involved in wheat domestication, provided with whole-genome resequencing. The accessions were phenotyped both in the field and in controlled conditions and A k-mer-based GWAS was performed to identify genomic regions involved in useful traits. Chapter 3. Describes the genetic basis of resistance to Z. tritici in a durum wheat elite diversity panel representative of the germplasm bred in Mediterranean. Quantitative trait loci (QTL) analysis results revealed several loci involved in the STB response that were found in several chromosome regions with a high infection rate. The genomic regions associated with STB resistance identified in this study could be of interest for marker assisted selection (MAS) in durum wheat breeding programs.

Association mapping of stem rust resistance in durum wheat at the seedling and adult plant stages

In wheat, stem rust is known to rapidly evolve new virulence to resistance genes. While more than 50 stem rust resistance (Sr) loci have been identified in wheat, only a few remain effective, particularly against the highly virulent race Ug99 (TTKSK race) and a mixture of durum-specific races. An association mapping (AM) study based on 183 durum wheat accessions was utilized to identify resistance loci for stem rust response in Ethiopia over four seasons and artificial inoculation with Ug99 (TTKSK race) and a mixture of durum-specific races under field conditions as well as in greenhouse test at seedling stage under controlled conditions for resistance to four highly virulent stem rust races: TRTTF, TTTTF, (TTKSK (Ug99) and JRCQC. The panel was profiled with 1,253 SSR and DArT markers. Twelve QTL-tagging markers were significant (P < 0.05) across three to four seasons. The role of Sr13, Sr9, Sr14, Sr17, and Sr28 was confirmed. Thirteen significant markers were in regions with no Sr genes/QTLs. The results under controlled conditions showed that 15, 20, 19 and 19 chromosome regions harbored markers that showed significant effects for races TRTTF, TTTTF, TTKSK and JRCQC, respectively. These genomic regions showed marker R2 values ranging from 1.13 to 8.34, 1.92 to 17.64, 1.75 to 23.12 and 1.51 to 15.33% for races TRTTF, TTTTF, TTKSK and JRCQC, respectively. The study demonstrates that stem rust resistance in durum wheat is governed in part by shared loci and in part by race-specific ones. The QTLs identified in this study through AM will be useful in the marker-assisted development of durum wheat cultivars with durable stem rust resistance.

Fine Mapping of qroot-yield-1.06, a QTL for Root, Plant Vigor and Yield in Maize

Root-yield-1.06 is a major QTL affecting root system architecture (RSA) and other agronomic traits in maize. The effect of this QTL has been evaluated with the development of near isogenic lines (NILs) differing at the QTL position. The objective of this study was to fine map qroot-yield-1.06 by marker-assisted searching for chromosome recombinants in the QTL interval and concurrent root phenotyping in both controlled and field conditions, through successive generations. Complementary approaches such as QTL meta-analysis and RNA-seq were deployed in order to help prioritizing candidate genes within the QTL target region. Using a selected group of genotypes, field based root analysis by ‘shovelomics’ enabled to accurately collect RSA information of adult maize plants. Shovelomics combined with software-assisted root imaging analysis proved to be an informative and relatively highly automated phenotyping protocol. A QTL interval mapping was conducted using a segregating population at the seedling stage grown in controlled environment. Results enabled to narrow down the QTL interval and to identify new polymorphic markers for MAS in field experiments. A collection of homozygous recombinant NILs was developed by screening segregating populations with markers flanking qroot-yield-1.06. A first set of lines from this collection was phenotyped based on the adapted shovelomics protocol. QTL analysis based on these data highlighted an interval of 1.3 Mb as completely linked with the target QTL but, a larger safer interval of 4.1 Mb was selected for further investigations. QTL meta-analysis allows to synthetize information on root QTLs and two mQTLs were identified in the qroot-yield-1.06 interval. Trascriptomics analysis based on RNA-seq data of the two contrasting QTL-NILs, confirmed alternative haplotypes at chromosome bin 1.06. qroot-yield-1.06 has now been delimited to a 4.1-Mb interval, and thanks to the availability of additional untested homozygous recombinant NILs, the potentially achievable mapping resolution at qroot-yield-1.06 is c. 50 kb.

Assessment of molecular sequences and enzyme activity of amylase trypsin inhibitors (ATIs) from a selection of ancient and modern grains

Wheat amylase-trypsin inhibitors (ATIs) are a family of wheat proteins, which play an important role in plant defence against pest attacks. ATIs are also of great interest for their impact on human health and recently ATIs have been identified as major stimulators of innate immune cells. In this study, ten selected wheat samples with different ploidy level and year of release were used for the agronomic trial, for in vitro enzymatic assays and for ATIs gene sequencing. Wheat samples were grown under organic farming management during three consecutive cropping years at two growing areas (Italy and USA). The PCA analysis performed on the deduced amino acid sequences of four representative ATIs genes (WMAI, WDAI, WTAI-CM3, CMx) evidenced that the ten wheat varieties can be differentiated on the basis of their ploidy level, but not with respect to ancient or recently developed wheat genotypes. The results from in vitro alpha-amylase and trypsin inhibitory activities showed high variability among the ten wheat genotypes and the contribution of the genotype and the cropping year was significant for both inhibitory activities. The hexaploid wheat genotypes showed the highest inhibitory activities. Einkorn showed a very low or even absent alpha-amylase inhibitory activity and the highest trypsin inhibitory activity. It was not possible to differentiate ancient and recently developed wheat genotypes on the basis of their ATIs activity. The weather conditions differently affected the two inhibitory activities. In both cultivation areas, higher precipitation and lower high mean temperatures correlated with lower alpha-amylase inhibitory activities, while there were different correlations considering trypsin inhibitory activity for the two growing areas. The protein content negatively correlated with both inhibitory activities in USA and Italy. This information can be important in the understanding of plant defence mechanisms in relation to the effect of both genotype and abiotic and biotic stress.

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.

Mass spectrometry-based protein profiling strategies for biomarker discovery in liver and inflammatory bowel diseases

The study of protein expression profiles for biomarker discovery in serum and in mammalian cell populations needs the continuous improvement and combination of proteins/peptides separation techniques, mass spectrometry, statistical and bioinformatic approaches. In this thesis work two different mass spectrometry-based protein profiling strategies have been developed and applied to liver and inflammatory bowel diseases (IBDs) for the discovery of new biomarkers. The first of them, based on bulk solid-phase extraction combined with matrix-assisted laser desorption/ionization - Time of Flight mass spectrometry (MALDI-TOF MS) and chemometric analysis of serum samples, was applied to the study of serum protein expression profiles both in IBDs (Crohn’s disease and ulcerative colitis) and in liver diseases (cirrhosis, hepatocellular carcinoma, viral hepatitis). The approach allowed the enrichment of serum proteins/peptides due to the high interaction surface between analytes and solid phase and the high recovery due to the elution step performed directly on the MALDI-target plate. Furthermore the use of chemometric algorithm for the selection of the variables with higher discriminant power permitted to evaluate patterns of 20-30 proteins involved in the differentiation and classification of serum samples from healthy donors and diseased patients. These proteins profiles permit to discriminate among the pathologies with an optimum classification and prediction abilities. In particular in the study of inflammatory bowel diseases, after the analysis using C18 of 129 serum samples from healthy donors and Crohn’s disease, ulcerative colitis and inflammatory controls patients, a 90.7% of classification ability and a 72.9% prediction ability were obtained. In the study of liver diseases (hepatocellular carcinoma, viral hepatitis and cirrhosis) a 80.6% of prediction ability was achieved using IDA-Cu(II) as extraction procedure. The identification of the selected proteins by MALDITOF/ TOF MS analysis or by their selective enrichment followed by enzymatic digestion and MS/MS analysis may give useful information in order to identify new biomarkers involved in the diseases. The second mass spectrometry-based protein profiling strategy developed was based on a label-free liquid chromatography electrospray ionization quadrupole - time of flight differential analysis approach (LC ESI-QTOF MS), combined with targeted MS/MS analysis of only identified differences. The strategy was used for biomarker discovery in IBDs, and in particular of Crohn’s disease. The enriched serum peptidome and the subcellular fractions of intestinal epithelial cells (IECs) from healthy donors and Crohn’s disease patients were analysed. The combining of the low molecular weight serum proteins enrichment step and the LCMS approach allowed to evaluate a pattern of peptides derived from specific exoprotease activity in the coagulation and complement activation pathways. Among these peptides, particularly interesting was the discovery of clusters of peptides from fibrinopeptide A, Apolipoprotein E and A4, and complement C3 and C4. Further studies need to be performed to evaluate the specificity of these clusters and validate the results, in order to develop a rapid serum diagnostic test. The analysis by label-free LC ESI-QTOF MS differential analysis of the subcellular fractions of IECs from Crohn’s disease patients and healthy donors permitted to find many proteins that could be involved in the inflammation process. Among them heat shock protein 70, tryptase alpha-1 precursor and proteins whose upregulation can be explained by the increased activity of IECs in Crohn’s disease were identified. Follow-up studies for the validation of the results and the in-depth investigation of the inflammation pathways involved in the disease will be performed. Both the developed mass spectrometry-based protein profiling strategies have been proved to be useful tools for the discovery of disease biomarkers that need to be validated in further studies.

Reverse genetic studies of Benyvirus - Polymyxa betae molecular interaction: Role of the RNA4-encoded protein in virus transmission

Beet necrotic yellow vein virus (BNYVV), the leading infectious agent that affects sugar beet, is included within viruses transmitted through the soil from plasmodiophorid as Polymyxa betae. BNYVV is the causal agent of Rhizomania, which induces abnormal rootlet proliferation and is widespread in the sugar beet growing areas in Europe, Asia and America; for review see (Peltier et al., 2008). In this latter continent, Beet soil-borne mosaic virus (BSBMV) has been identified (Lee et al., 2001) and belongs to the benyvirus genus together with BNYVV, both vectored by P. betae. BSBMV is widely distributed only in the United States and it has not been reported yet in others countries. It was first identified in Texas as a sugar beet virus morphologically similar but serologically distinct to BNYVV. Subsequent sequence analysis of BSBMV RNAs evidenced similar genomic organization to that of BNYVV but sufficient molecular differences to distinct BSBMV and BNYVV in two different species (Rush et al., 2003). Benyviruses field isolates usually consist of four RNA species but some BNYVV isolates contain a fifth RNA. RNAs -1 contains a single long ORF encoding polypeptide that shares amino acid homology with known viral RNA-dependent RNA polymerases (RdRp) and helicases. RNAs -2 contains six ORFs: capsid protein (CP), one readthrough protein, triple gene block proteins (TGB) that are required for cell-to-cell virus movement and the sixth 14 kDa ORF is a post-translation gene silencing suppressor. RNAs -3 is involved on disease symptoms and is essential for virus systemic movement. BSBMV RNA-3 can be trans-replicated, trans-encapsidated by the BNYVV helper strain (RNA-1 and -2) (Ratti et al., 2009). BNYVV RNA-4 encoded one 31 kDa protein and is essential for vector interactions and virus transmission by P. betae (Rahim et al., 2007). BNYVV RNA-5 encoded 26 kDa protein that improve virus infections and accumulation in the hosts. We are interest on BSBMV effect on Rhizomania studies using powerful tools as full-length infectious cDNA clones. B-type full-length infectious cDNA clones are available (Quillet et al., 1989) as well as A/P-type RNA-3, -4 and -5 from BNYVV (unpublished). A-type BNYVV full-length clones are also available, but RNA-1 cDNA clone still need to be modified. During the PhD program, we start production of BSBMV full-length cDNA clones and we investigate molecular interactions between plant and Benyviruses exploiting biological, epidemiological and molecular similarities/divergences between BSBMV and BNYVV. During my PhD researchrs we obtained full length infectious cDNA clones of BSBMV RNA-1 and -2 and we demonstrate that they transcripts are replicated and packaged in planta and able to substitute BNYVV RNA-1 or RNA-2 in a chimeric viral progeny (BSBMV RNA-1 + BNYVV RNA-2 or BNYVV RNA-1 + BSBMV RNA-2). During BSBMV full-length cDNA clones production, unexpected 1,730 nts long form of BSBMV RNA-4 has been detected from sugar beet roots grown on BSBMV infected soil. Sequence analysis of the new BSBMV RNA-4 form revealed high identity (~100%) with published version of BSBMV RNA-4 sequence (NC_003508) between nucleotides 1-608 and 1,138-1,730, however the new form shows 528 additionally nucleotides between positions 608-1,138 (FJ424610). Two putative ORFs has been identified, the first one (nucleotides 383 to 1,234), encode a protein with predicted mass of 32 kDa (p32) and the second one (nucleotides 885 to 1,244) express an expected product of 13 kDa (p13). As for BSBMV RNA-3 (Ratti et al., 2009), full-length BSBMV RNA-4 cDNA clone permitted to obtain infectious transcripts that BNYVV viral machinery (Stras12) is able to replicate and to encapsidate in planta. Moreover, we demonstrated that BSBMV RNA-4 can substitute BNYVV RNA-4 for an efficient transmission through the vector P. betae in Beta vulgaris plants, demonstrating a very high correlation between BNYVV and BSBMV. At the same time, using BNYVV helper strain, we studied BSBMV RNA-4’s protein expression in planta. We associated a local necrotic lesions phenotype to the p32 protein expression onto mechanically inoculated C. quinoa. Flag or GFP-tagged sequences of p32 and p13 have been expressed in viral context, using Rep3 replicons, based on BNYVV RNA-3. Western blot analyses of local lesions contents, using FLAG-specific antibody, revealed a high molecular weight protein, which suggest either a strong interaction of BSBMV RNA4’s protein with host protein(s) or post translational modifications. GFP-fusion sequences permitted the subcellular localization of BSBMV RNA4’s proteins. Moreover we demonstrated the absence of self-activation domains on p32 by yeast two hybrid system approaches. We also confirmed that p32 protein is essential for virus transmission by P. betae using BNYVV helper strain and BNYVV RNA-3 and we investigated its role by the use of different deleted forms of p32 protein. Serial mechanical inoculation of wild-type BSBMV on C. quinoa plants were performed every 7 days. Deleted form of BSBMV RNA-4 (1298 bp) appeared after 14 passages and its sequence analysis shows deletion of 433 nucleotides between positions 611 and 1044 of RNA-4 new form. We demonstrated that this deleted form can’t support transmission by P. betae using BNYVV helper strain and BNYVV RNA-3, moreover we confirmed our hypothesis that BSBMV RNA-4 described by Lee et al. (2001) is a deleted form. Interesting after 21 passages we identifed one chimeric form of BSBMV RNA-4 and BSBMV RNA-3 (1146 bp). Two putative ORFs has been identified on its sequence, the first one (nucleotides 383 to 562), encode a protein with predicted mass of 7 kDa (p7), corresponding to the N-terminal of p32 protein encoded by BSBMV RNA-4; the second one (nucleotides 562 to 789) express an expected product of 9 kDa (p9) corresponding to the C-terminal of p29 encoded by BSBMV RNA-3. Results obtained by our research in this topic opened new research lines that our laboratories will develop in a closely future. In particular BSBMV p32 and its mutated forms will be used to identify factors, as host or vector protein(s), involved in the virus transmission through P. betae. The new results could allow selection or production of sugar beet plants able to prevent virus transmission then able to reduce viral inoculum in the soil.

Immunohistochemical and Molecular Prognostic/Predictive Markers in Neoplastic Diseases

Traditional morphological examinations are not anymore sufficient for a complete evaluation of tumoral tissue and the use of neoplastic markers is of utmost importance. Neoplastic markers can be classified in: diagnostic, prognostic and predictive markers. Three markers were analyzed. 1) Insulin-like growth factor binding protein 2 (IGFBP2) was immunohistochemically examined in prostatic tissues: 40 radical prostatectomies from hormonally untreated patients with their preoperative biopsies, 10 radical prostatectomies from patients under complete androgen ablation before surgery and 10 simple prostatectomies from patients with bladder outlet obstruction. Results were compared with α-methylacyl-CoA racemase (AMACR). IGFBP2 was expressed in the cytoplasm of untreated adenocarcinomas and, to a lesser extent, in HG-PIN; the expression was markedly lower in patients after complete androgen ablation. AMACR was similarly expressed in both adenocarcinoma and HG-PIN, the level being similar in both lesions; the expression was slightly lower in patients after complete androgen ablation. IGFBP2 may be used a diagnostic marker of prostatic adenocarcinomas. 2) Heparan surface proteoglycan immunohistochemical expression was examined in 150 oral squamous cell carcinomas. Follow up information was available in 93 patients (range: 6-34 months, mean: 19±7). After surgery, chemotherapy was performed in 8 patients and radiotherapy in 61 patients. Multivariate and univariate overall survival analyses showed that high expression of syndecan-1 (SYN-1) was associated with a poor prognosis. In patients treated with radiotherapy, such association was higher. SYN-1 is a prognostic marker in oral squamous cell carcinomas; it may also represent a predictive factor for responsiveness to radiotherapy. 3) EGFR was studied in 33 pulmonary adenocarcinomas with traditional DNA sequencing methods and with two mutation-specific antibodies. Overall, the two antibodies had 61.1% sensitivity and 100% specificity in detecting EGFR mutations. EGFR mutation-specific antibodies may represent a predictive marker to identify patients candidate to tyrosine kinase inhibitors therapy.