Landscape ecology and genetics of the wolf in Italy

This PhD Thesis includes five main parts on diverse topics. The first two parts deal with the trophic ecology of wolves in Italy consequently to a recent increase of wild ungulates abundance. Data on wolf diet across time highlighted how wild ungulates are important food resource for wolves in Italy. Increasing wolf population, increasing numbers of wild ungulates and decreasing livestock consume are mitigating wolf-man conflicts in Italy in the near future. In the third part, non-invasive genetic sampling techniques were used to obtain genotypes and genders of about 400 wolves. Thus, wolf packs were genetically reconstructed using diverse population genetic and parentage software. Combining the results on pack structure and genetic relatedness with sampling locations, home ranges of wolf packs and dispersal patterns were identified. These results, particularly important for the conservation management of wolves in Italy, illustrated detailed information that can be retrieved from genetic identification of individuals. In the fourth part, wolf locations were combined with environmental information obtained as GIS-layers. Modern species distribution models (niche models) were applied to infer potential wolf distribution and predation risk. From the resulting distribution maps, information pastures with the highest risk of depredation were derived. This is particularly relevant as it allows identifying those areas under danger of carnivore attack on livestock. Finally, in the fifth part, habitat suitability models were combined with landscape genetic analysis. On one side landscape genetic analyses on the Italian wolves provided new information on the dynamics and connectivity of the population and, on the other side, a profound analysis of the effects that habitat suitability methods had on the parameterization of landscape genetic analyses was carried out to contributed significantly to landscape genetic theory.

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.

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.

Genetic characterization, population history and evolutionary medicine perspective in two native south american populations: Yanesha and Wichi

Two Amerindian populations from the Peruvian Amazon (Yanesha) and from rural lowlands of the Argentinean Gran Chaco (Wichi) were analyzed. They represent two case study of the South American genetic variability. The Yanesha represent a model of population isolated for long-time in the Amazon rainforest, characterized by environmental and altitudinal stratifications. The Wichi represent a model of population living in an area recently colonized by European populations (the Criollos are the population of the admixed descendents), whose aim is to depict the native ancestral gene pool and the degree of admixture, in relation to the very high prevalence of Chagas disease. The methods used for the genotyping are common, concerning the Y chromosome markers (male lineage) and the mitochondrial markers (maternal lineage). The determination of the phylogeographic diagnostic polymorphisms was carried out by the classical techniques of PCR, restriction enzymes, sequencing and specific mini-sequencing. New method for the detection of the protozoa Trypanosoma cruzi was developed by means of the nested PCR. The main results show patterns of genetic stratification in Yanesha forest communities, referable to different migrations at different times, estimated by Bayesian analyses. In particular Yanesha were considered as a population of transition between the Amazon basin and the Andean Cordillera, evaluating the potential migration routes and the separation of clusters of community in relation to different genetic bio-ancestry. As the Wichi, the gene pool analyzed appears clearly differentiated by the admixed sympatric Criollos, due to strict social practices (deeply analyzed with the support of cultural anthropological tools) that have preserved the native identity at a diachronic level. A pattern of distribution of the seropositivity in relation to the different phylogenetic lineages (the adaptation in evolutionary terms) does not appear, neither Amerindian nor European, but in relation to environmental and living conditions of the two distinct subpopulations.

Identification and genetic diversity in phytoplasmas associated with diseases of cassava and other agronomic relevant crops in south-east Asia and Latin America

Identification and genetic diversity of phytoplasmas infecting tropical plant species, selected among those most agronomically relevant in South-east Asia and Latin America were studied. Correlation between evolutionary divergence of relevant phytoplasma strains and their geographic distribution by comparison on homologous genes of phytoplasma strains detected in the same or related plant species in other geographical areas worldwide was achieved. Molecular diversity was studied on genes coding ribosomal proteins, groEL, tuf and amp besides phytoplasma 16S rRNA. Selected samples infected by phytoplasmas belonging to diverse ribosomal groups were also studied by in silico RFLP followed by phylogenetic analyses. Moreover a partial genome annotation of a ‘Ca. P. brasiliense’ strain was done towards future application for epidemiological studies. Phytoplasma presence in cassava showing frog skin (CFSD) and witches’ broom (CWB) diseases in Costa Rica - Paraguay and in Vietnam – Thailand, respectively, was evaluated. In both cases, the diseases were associated with phytoplasmas related to aster yellows, apple proliferation and “stolbur” groups, while only phytoplasma related to X-disease group in CFSD, and to hibiscus witches’ broom, elm yellows and clover proliferation groups in CWB. Variability was found among strains belonging to the same ribosomal group but having different geographic origin and associated with different disease. Additionally, a dodder transmission assay to elucidate the role of phytoplasmas in CWB disease was carried out, and resulted in typical phytoplasma symptoms in periwinkle plants associated with the presence of aster yellows-related strains. Lethal wilt disease, a severe disease of oil palm in Colombia that is spreading throughout South America was also studied. Phytoplasmas were detected in symptomatic oil palm and identified as ‘Ca. P. asteris’, ribosomal subgroup 16SrI-B, and were distinguished from other aster yellows phytoplasmas used as reference strains; in particular, from an aster yellows strain infecting corn in the same country.

Identification of novel genetic alterations in pediatric cytogenetically normal acute myeloid leukemia by next-generation sequencing

Pediatric acute myeloid leukemia (AML) is a molecularly heterogeneous disease that arises from genetic alterations in pathways that regulate self-renewal and myeloid differentiation. While the majority of patients carry recurrent chromosomal translocations, almost 20% of childhood AML do not show any recognizable cytogenetic alteration and are defined as cytogenetically normal (CN)-AML. CN-AML patients have always showed a great variability in response to therapy and overall outcome, underlining the presence of unknown genetic changes, not detectable by conventional analyses, but relevant for pathogenesis, and outcome of AML. The development of novel genome-wide techniques such as next-generation sequencing, have tremendously improved our ability to interrogate the cancer genome. Based on this background, the aim of this research study was to investigate the mutational landscape of pediatric CN-AML patients negative for all the currently known somatic mutations reported in AML through whole-transcriptome sequencing (RNA-seq). RNA-seq performed on diagnostic leukemic blasts from 19 pediatric CN-AML cases revealed a considerable incidence of cryptic chromosomal rearrangements, with the identification of 21 putative fusion genes. Several of the fusion genes that were identified in this study are recurrent and might have a prognostic and/or therapeutic relevance. A paradigm of that is the CBFA2T3-GLIS2 fusion, which has been demonstrated to be a common alteration in pediatric CN-AML, predicting poor outcome. Important findings have been also obtained in the identification of novel therapeutic targets. On one side, the identification of NUP98-JARID1A fusion suggests the use of disulfiram; on the other, here we describe alteration-activating tyrosine kinases, providing functional data supporting the use of tyrosine kinase inhibitors to specifically inhibit leukemia cells. This study provides new insights in the knowledge of genetic alterations underlying pediatric AML, defines novel prognostic markers and putative therapeutic targets, and prospectively ensures a correct risk stratification and risk-adapted therapy also for the “all-neg” AML subgroup.