Mitochondrial DNA rearrangements during human aging: a study on liver, muscle and adipose tissue

Aging is characterized by a chronic, low-grade inflammatory state called “inflammaging”. Mitochondria are the main source of reactive oxygen species (ROS), which trigger the production of pro-inflammatory molecules. We are interested in studying the age-related modifications of the mitochondrial DNA (mtDNA), which can be affected by the lifelong exposure to ROS and are responsible of mitochondrial dysfunction. Moreover, increasing evidences show that telomere shortening, naturally occurring with aging, is involved in mtDNA damage processes and thus in the pathogenesis of age-related disorders. Thus the primary aim of this thesis was the analysis of mtDNA copy number, deletion level and integrity in different-age human biopsies from liver, vastus lateralis skeletal muscle of healthy subjects and patients with limited mobility of lower limbs (LMLL), as well as adipose tissue. The telomere length and the expression of nuclear genes related to mitobiogenesis, fusion and fission, mitophagy, mitochondrial protein quality control system, hypoxia, production and protection from ROS were also evaluated. In liver the decrease in mtDNA integrity with age is accompanied with an increase in mtDNA copy number, suggesting the existence of a “compensatory mechanism” able to maintain the functionality of this organ. Different is the case of vastus lateralis muscle, where any “compensatory pathway” is activated and mtDNA integrity and copy number decrease with age, both in healthy subjects and in patients. Interestingly, mtDNA rearrangements do not incur in adipose tissue with advancing age. Finally, in all tissues a marked gender difference appears, suggesting that aging and also gender diversely affect mtDNA rearrangements and telomere length in the three human tissues considered, likely depending on their different metabolic needs and inflammatory status.

Characterization of mitochondrial genomes in bivalve species with doubly uniparental inheritance of mitochondria

Many bivalve species possess two distinct mtDNA lineages, called F and M, respectively inherited maternally and paternally: this system is called doubly uniparental inheritance (DUI). The main experimental project of my PhD was the quantification of the two mtDNAs during the development of the DUI species Ruditapes philippinarum, from early embryos to sub-adults, using Real-Time qPCR. I identified the time interval in which M mtDNA is lost from female individuals, while it is retained in males (which are heteroplasmic through all of their life cycle). The results also suggested absence of mtDNA replication during early embryogenesis, a process constituting a bottleneck that highly reduces the copy number of mtDNA molecules in cells of developing larvae. In males this bottleneck may produce cells homoplasmic for M mtDNA, and could be considered as a first step of the segregation of M in the male germ line. Another finding was the characterization, in young clams approaching the first reproductive season, of a significant boost in copy number of F mtDNA in females and of M in males. Given the age of animals in which this mtDNA-specific growth was observed, the finding could probably be the outcome of the first round of gonads and gametes production. Other lines of research included the characterization of the unassigned regions in mt genomes of DUI bivalves. These regions can harbor signals involved in the control of replication and/or transcription of the mtDNA molecule, as well as additional open reading frames (ORFs) not related to oxidative phosphorylation. These features in DUI species could be associated to the maintenance of separate inheritance routes for the two mtDNAs. Additional ORFs are also found in other animal mt genomes: I summarized the presence of gene duplications as a co-author in a review focusing on animal mt genomes with unusual gene content.

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

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

Addressing tumor heterogeneity in esophageal adenocarcinoma through different molecular approaches

Several studies have shown epidemiologic, clinical, immune-histochemical and molecular differences among esophageal adenocarcinomas (EAC). Since pathogenesis and biology of this tumor are far to be well defined, our study aimed to examine intra- and inter-tumor heterogeneity and to solve crucial controversies through different molecular approaches. Target sequencing was performed for sorted cancer subpopulations from formalin embedded material obtained from 38 EACs, not treated with neoadjuvant therapy. 35 out 38 cases carried at least one somatic mutation, not present in the corresponding sorted stromal cells. 73.7% of cases carried mutations in TP53 and 10.5% in CDKN2A. Mutations in other genes occurred at lower frequency, including HNF1A, not previously associated with EAC. Sorting allowed us to isolate clones with different mutational loads and/or additional copy number amplifications, confirming the high intra-tumor heterogeneity of these cancers. In our cohort TP53 gene abnormalities correlated with a better survival (P = 0.028); conversely, loss of SMAD4 protein expression was associated with a higher recurrence rate (P = 0.015). Shifting the focus on the epigenetic characterization of EAC, miR-221 and miR-483-3p resulted upregulated from the MicroRNA Array card analysis and confirmed with further testing. The up-regulation of both miRNAs correlated with clinical outcomes, in particular with a reduced cancer-specific survival (miR483-3p P=0.0293; miR221 P=0.0059). In vitro analyses demonstrated an increase for miR-483-3p (fold-change=2.7) that appear to be inversely correlated with SMAD4 expression in FLO-1 cell-line. In conclusion, selective sorting allowed to define the real mutation status and to isolate different cancer subclones. MiRNA expression analysis revealed a significant up-regulation of miR-221 and miR-483-3p, which correlated with worst prognosis, implying that they can be considered oncogenic factors in EAC. Therefore, cell sorting technologies, coupled with next generation sequencing, and the analysis of microRNA profiles seem to be promising strategies to guide treatment and help classify cancer prognosis.

Harnessing the power of Drosophila model to obtain new insights on MYC and ABCC1 cellular functions: from neuroblastoma to autism spectrum disorder

MYCN amplification is a genetic hallmark of the childhood tumour neuroblastoma. MYCN-MAX dimers activate the expression of genes promoting cell proliferation. Moreover, MYCN seems to transcriptionally repress cell differentiation even in absence of MAX. We adopted the Drosophila eye as model to investigate the effect of high MYC to MAX expression ratio on cells. We found that dMyc overexpression in eye cell precursors inhibits cell differentiation and induces the ectopic expression of Antennapedia (the wing Hox gene). The further increase of MYC/MAX ratio results in an eye-to-wing homeotic transformation. Notably, dMyc overexpression phenotype is suppressed by low levels of transcriptional co-repressors and MYCN associates to the promoter of Deformed (the eye Hox gene) in proximity to repressive sites. Hence, we envisage that, in presence of high MYC/MAX ratio, the “free MYC” might inhibit Deformed expression, leading in turn to the ectopic expression of Antennapedia. This suggests that MYCN might reinforce its oncogenic role by affecting the physiological homeotic program. Furthermore, poor neuroblastoma outcome associates with a high level of the MRP1 protein, encoded by the ABCC1 gene and known to promote drug efflux in cancer cells. Intriguingly, this correlation persists regardless of chemotherapy and ABCC1 overexpression enhances neuroblastoma cell motility. We found that Drosophila dMRP contributes to the adhesion between the dorsal and ventral epithelia of the wing by inhibiting the function of integrin receptors, well known regulators of cell adhesion and migration. Besides, integrins play a crucial role during synaptogenesis and ABCC1 locus is included in a copy number variable region of the human genome (16p13.11) involved in neuropsychiatric diseases. Interestingly, we found that the altered dMRP/MRP1 level affects nervous system development in Drosophila embryos. These preliminary findings point out novel ABCC1 functions possibly defining ABCC1 contribution to neuroblastoma and to the pathogenicity of 16p13.11 deletion/duplication

Tumor heterogeneity in Esophageal Adenocarcinoma: a genomic approach

INTRODUCTION: Esophageal adenocarcinoma (EAC) is a severe malignancy in terms of prognosis and mortality rate. Because its great genetic heterogeneity, disputes regarding classification, prevention and treatments are still unsolved. AIM: We investigated intra- and inter-EAC heterogeneity by defining EAC’s somatic mutational profile and the role of candidate microRNAs, to correlate the molecular profile of tumors to clinical outcomes and to identify biomarkers for classification. METHODS: 38 EAC cases were analyzed via high-throughput cell sorting technology combined with targeted sequencing and whole genome low-pass sequencing. Targeted sequencing of further 169 cases was performed to widen the study. miR221 and miR483-3p expression was profiled via qPCR in 112 EACs and correlation with clinical outcomes was investigated. RESULTS: 35/38 EACs carried at least one somatic mutation absent in stromal cells. TP53 was found mutated in 73.7% of cases. Selective sorting revealed tumor subclones with different mutational loads and copy number alterations, confirming the high intra-tumor heterogeneity of EAC. Mutations were in most cases at homozygous state, and we identified alterations that were missed with the whole-tumor analysis. Mutations in HNF1A gene, not previously associated with EAC, were identified in both cohorts. Higher expression of miR483-3p and miR221 was associated with poorer cancer specific survival (P=0.0293 and P=0.0059), and recurrence in the Lauren intestinal subtype (P=0.0459 and P=0.0002). Median expression levels of miRNAs were higher in patients with advanced tumor stages. The loss of SMAD4 immunoreactivity was significantly associated with poorer cancer specific survival and recurrence (P=0.0452; P=0.022 respectively). CONCLUSION: Combining selective sorting technology and next generation sequencing allowed to better define EAC inter- and intra-tumor heterogeneity. We identified HNF1A as a new mutated gene associated to EAC that could be involved in tumor progression and promising biomarkers such as SMAD4, miR221 and miR483-3p to identify patients at higher risk for more aggressive tumors.

Phenotypic and molecular investigation of circulating tumor cells (CTCs) isolated from breast cancer patients at single cell resolution

Despite the paramount advances in cancer research, breast cancer (BC) still ranks one of the leading causes of cancer-related death worldwide. Thanks to the screening campaign started in developed countries, BC is often diagnosed at early stages (non-metastatic BC, nmBC), but disease relapse occurrence even after decades and at distant sites is not an uncommon phenomenon. Conversely, metastatic BC (mBC) is considered an incurable disease. The major perpetrators of tumor spread to secondary organs are circulating tumor cells (CTCs), a rare population of cells detectable in the peripheral blood of oncologic patients. In this study, CTCs from patients diagnosed with luminal nmBC and mBC (hormone receptor positive, Human Epidermal Growth Factor Receptor 2 (HER2) negative) were characterized at both phenotypic and molecular levels. To better understand the molecular mechanisms underlying their biology and their metastatic potential, next-generation sequencing (NGS) analyses were performed at single-cell resolution to assess copy number aberrations (CNAs), single nucleotide variants (SNVs) and gene expression profiling. The findings of this study arise hints in CTC detection, and pave the way to new application in CTC research.