PCR-Based Analysis of Mitochondrial DNA Copy Number, Mitochondrial DNA Damage, and Nuclear DNA Damage.

Because of the role that DNA damage and depletion play in human disease, it is important to develop and improve tools to assess these endpoints. This unit describes PCR-based methods to measure nuclear and mitochondrial DNA damage and copy number. Long amplicon quantitative polymerase chain reaction (LA-QPCR) is used to detect DNA damage by measuring the number of polymerase-inhibiting lesions present based on the amount of PCR amplification; real-time PCR (RT-PCR) is used to calculate genome content. In this unit, we provide step-by-step instructions to perform these assays in Homo sapiens, Mus musculus, Rattus norvegicus, Caenorhabditis elegans, Drosophila melanogaster, Danio rerio, Oryzias latipes, Fundulus grandis, and Fundulus heteroclitus, and discuss the advantages and disadvantages of these assays.

Maternal Copy Number Variants (CNV) transmission to their Autism Spectrum Disorder (ASD) sons correlates with phenotypic traits

Background: Common neurodevelopmental disorder, global prevalence ~1 %; Persistent deficits in social communication and social interaction; restricted and repetitive behavior, interests, or activities; Highly heterogeneous clinical presentation; Male to female ratio ~4:1.

Small-plasmid-mediated antibiotic resistance is enhanced by increases in plasmid copy number and bacterial fitness

Plasmids play a key role in the horizontal spread of antibiotic resistance determinants among bacterial pathogens. When an antibiotic resistance plasmid arrives in a new bacterial host, it produces a fitness cost, causing a competitive disadvantage for the plasmid-bearing bacterium in the absence of antibiotics. On the other hand, in the presence of antibiotics, the plasmid promotes the survival of the clone. The adaptations experienced by plasmid and bacterium in the presence of antibiotics during the first generations of coexistence will be crucial for the progress of the infection and the maintenance of plasmid-mediated resistance once the treatment is over. Here we developed a model system using the human pathogen Haemophilus influenzae carrying the small plasmid pB1000 conferring resistance to β-lactam antibiotics to investigate host and plasmid adaptations in the course of a simulated ampicillin therapy. Our results proved that plasmid-bearing clones compensated for the fitness disadvantage during the first 100 generations of plasmid-host adaptation. In addition, ampicillin treatment was associated with an increase in pB1000 copy number. The augmentation in both bacterial fitness and plasmid copy number gave rise to H. influenzae populations with higher ampicillin resistance levels. In conclusion, we show here that the modulations in bacterial fitness and plasmid copy number help a plasmid-bearing bacterium to adapt during antibiotic therapy, promoting both the survival of the host and the spread of the plasmid.

TiMMing: developing an innovative suite of bioinformatic tools to harmonize and track the origin of copy number alterations in the evolutive history of multiple myeloma

Multiple Myeloma (MM) is a hematologic cancer with heterogeneous and complex genomic landscape, where Copy Number Alterations (CNAs) play a key role in the disease's pathogenesis and prognosis. It is of biological and clinical interest to study the temporal occurrence of early alterations, as they play a disease "driver" function by deregulating key tumor pathways. This study presents an innovative bioinformatic tools suite created for harmonizing and tracing the origin of CNAs throughout the evolutionary history of MM. To this aim, large cohorts of newly-diagnosed MM (NDMM, N=1582) and Smoldering-MM (SMM, N=282) were aggregated. The tools developed in this study enable the harmonization of CNAs as obtained from different genomic platforms in such a way that a high statistical power can be obtained. By doing so, the high numerosity of those cohorts was harnessed for the identification of novel genes characterized as "driver" (NFKB2, NOTCH2, MAX, EVI5 and MYC-ME2-enhancer), and the generation of an innovative timing model, implemented with a statistical method to introduce confidence intervals in the CNAs-calls. By applying this model on both NDMM and SMM cohorts, it was possible to identify specific CNAs (1q(CKS1B)amp, 13q(RB1)del, 11q(CCND1)amp and 14q(MAX)del) and categorize them as "early"/ "driver" events. A high level of precision was guaranteed by the narrow confidence intervals in the timing estimates. These CNAs were proposed as critical MM alterations, which play a foundational role in the evolutionary history of both SMM and NDMM. Finally, a multivariate survival model was able to identify the independent genomic alterations with the greatest effect on patients’ survival, including RB1-del, CKS1B-amp, MYC-amp, NOTCH2-amp and TRAF3-del/mut. In conclusion, the alterations that were identified as both "early-drivers” and correlated with patients’ survival were proposed as biomarkers that, if included in wider survival models, could provide a better disease stratification and an improved prognosis definition.

Identification of outcome-oriented cut-offs for copy number alterations in multiple myeloma: predictive biomarkers with improved prognostic accuracy

Aim of the present study was to develop a statistical approach to define the best cut-off Copy number alterations (CNAs) calling from genomic data provided by high throughput experiments, able to predict a specific clinical end-point (early relapse, 18 months) in the context of Multiple Myeloma (MM). 743 newly diagnosed MM patients with SNPs array-derived genomic and clinical data were included in the study. CNAs were called both by a conventional (classic, CL) and an outcome-oriented (OO) method, and Progression Free Survival (PFS) hazard ratios of CNAs called by the two approaches were compared. The OO approach successfully identified patients at higher risk of relapse and the univariate survival analysis showed stronger prognostic effects for OO-defined high-risk alterations, as compared to that defined by CL approach, statistically significant for 12 CNAs. Overall, 155/743 patients relapsed within 18 months from the therapy start. A small number of OO-defined CNAs were significantly recurrent in early-relapsed patients (ER-CNAs) - amp1q, amp2p, del2p, del12p, del17p, del19p -. Two groups of patients were identified either carrying or not ≥1 ER-CNAs (249 vs. 494, respectively), the first one with significantly shorter PFS and overall survivals (OS) (PFS HR 2.15, p<0001; OS HR 2.37, p<0.0001). The risk of relapse defined by the presence of ≥1 ER-CNAs was independent from those conferred both by R-IIS 3 (HR=1.51; p=0.01) and by low quality (< stable disease) clinical response (HR=2.59 p=0.004). Notably, the type of induction therapy was not descriptive, suggesting that ER is strongly related to patients’ baseline genomic architecture. In conclusion, the OO- approach employed allowed to define CNAs-specific dynamic clonality cut-offs, improving the CNAs calls’ accuracy to identify MM patients with the highest probability to ER. As being outcome-dependent, the OO-approach is dynamic and might be adjusted according to the selected outcome variable of interest.

Down-regulation Of Slc8a1 As A Putative Apoptosis-evasion Mechanism By Modulation Of Calcium Levels In Penile Carcinoma.

The SLC8A1 gene, which encodes the Na(+)/Ca(2+) exchanger, plays a key role in calcium homeostasis. Our previous gene expression oligoarray data revealed SLC8A1 underexpression in penile carcinoma (PeCa). The aim of this study was to investigate whether the dysregulation of SLC8A1 expression is associated with apoptosis and cell proliferation in PeCa, via modulation of calcium concentration. The underlying mechanisms of SLC8A1 underexpression were also explored, focusing on copy number alteration and microRNA. Transcript levels of SLC8A1 gene and miR-223 were evaluated by quantitative PCR, comparing PeCa samples with normal glans tissues. SLC8A1 copy number was evaluated by microarray-based comparative genomic hybridization (array-CGH). Caspase-3 and Ki-67 immunostaining, as well as calcium distribution by Laser Ablation Imaging Inductively Coupled Plasma Mass Spectrometry [LA(i)-ICP-MS], were investigated in both normal and tumor samples. Confirming our previous data, SLC8A1 underexpression was detected in PeCa samples (P=0.001) and was not associated with gene copy number loss. In contrast, overexpression of miR-223 (P=0.002) was inversely correlated with SLC8A1 (P=0.015, r=-0.426), its putative repressor. In addition, SLC8A1 underexpression was associated with decreased calcium distribution, high Ki-67 and low caspase-3 immunoexpression in PeCa when compared with normal tissues. Down-regulation of the SLC8A1 gene, most likely mediated by its regulator miR-223, can lead to reduced calcium levels in PeCa and, consequently, to suppression of apoptosis and increased tumor cell proliferation. These data suggest that the miR-223-NCX1-calcium-signaling axis may represent a potential therapeutic approach in PeCa.

Mutation accumulation in Tetrahymena

Background - The rate and fitness effects of mutations are key in understanding the evolution of every species. Traditionally, these parameters are estimated in mutation accumulation experiments where replicate lines are propagated in conditions that allow mutations to randomly accumulate without the purging effect of natural selection. These experiments have been performed with many model organisms but we still lack empirical estimates of the rate and effects of mutation in the protists. Results - We performed a mutation accumulation (MA) experiment in Tetrahymena thermophila, a species that can reproduce sexually and asexually in nature, and measured both the mean decline and variance increase in fitness of 20 lines. The results obtained with T. thermophila were compared with T. pyriformis that is an obligate asexual species. We show that MA lines of T. thermophila go to extinction at a rate of 1.25 clonal extinctions per bottleneck. In contrast, populations of T. pyriformis show a much higher resistance to extinction. Variation in gene copy number is likely to be a key factor in explaining these results, and indeed we show that T. pyriformis has a higher mean copy number per cell than T. thermophila. From fitness measurements during the MA experiment, we infer a rate of mutation to copy number variation of 0.0333 per haploid MAC genome of T. thermophila and a mean effect against copy number variation of 0.16. A strong effect of population size in the rate of fitness decline was also found, consistent with the increased power of natural selection. Conclusions - The rate of clonal extinction measured for T. thermophila is characteristic of a mutational degradation and suggests that this species must undergo sexual reproduction to avoid the deleterious effects detected in the laboratory experiments. We also suggest that an increase in chromosomal copy number associated with the phenotypic assortment of amitotic divisions can provide an alternative mechanism to escape the deleterious effect of random chromosomal copy number variation in species like T. pyriformis that lack the resetting mechanism of sexual reproduction. Our results are relevant to the understanding of cell line longevity and senescence in ciliates.

Low-grade and high-grade mammary carcinomas in WAP-T transgenic mice are independent entities distinguished by Met expression.

Mammary carcinomas developing in SV40 transgenic WAP-T mice arise in two distinct histological phenotypes: as differentiated low-grade and undifferentiated high-grade tumors. We integrated different types of information such as histological grading, analysis of aCGH-based gene copy number and gene expression profiling to provide a comprehensive molecular description of mammary tumors in WAP-T mice. Applying a novel procedure for the correlation of gene copy number with gene expression on a global scale, we observed in tumor samples a global coherence between genotype and transcription. This coherence can be interpreted as a matched transcriptional regulation inherited from the cells of tumor origin and determined by the activity of cancer driver genes. Despite common recurrent genomic aberrations, e.g. gain of chr. 15 in most WAP-T tumors, loss of chr. 19 frequently occurs only in low-grade tumors. These tumors show features of "basal-like" epithelial differentiation, particularly expression of keratin 14. The high-grade tumors are clearly separated from the low-grade tumors by strong expression of the Met gene and by coexpression of epithelial (e.g. keratin 18) and mesenchymal (e.g. vimentin) markers. In high-grade tumors, the expression of the nonmutated Met protein is associated with Met-locus amplification and Met activity. The role of Met as a cancer driver gene is supported by the contribution of active Met signaling to motility and growth of mammary tumor-derived cells. Finally, we discuss the independent origin of low- and high-grade tumors from distinct cells of tumor origin, possibly luminal progenitors, distinguished by Met gene expression and Met signaling.

Molecular evolution of the Drosophila olfactory system

The olfactory system is an attractive model to study the genetic mechanisms underlying evolution of the nervous system. This sensory system mediates the detection and behavioural responses to an enormous diversity of volatile chemicals in the environment and displays rapid evolution, as species acquire, modify and discard olfactory receptors and circuits to adapt to new olfactory stimuli. Drosophilids provide an attractive model to study these processes. The availability of 12 sequenced genomes of Drosophila species occupying diverse ecological niches provides a rich resource for genomic analyses. Moreover, one of these species, Drosophila melanogaster, is amenable to a powerful combination of genetic and electrophysiological analyses. D. melanogaster has two distinct families of olfactory receptors to detect odours, the well-characterised Odorant Receptors (ORs) and the recently identified lonotropic Receptors (IRs). In my thesis, I have provided new insights into the genetic mechanisms underlying olfactory system evolution through three distinct, but interrelated projects. First, I performed a comparative genomic analysis of the IR repertoire in 12 sequenced Drosophila species, which has revealed that the olfactory IRs are highly conserved across species. By contrast, a large fraction of IRs that are not expressed in the olfactory system - and which may be gustatory receptors - are much more variable in sequence and gene copy number. Second, to identify ligands for IR expressing olfactory sensory neurons, I have performed an electrophysiological screen in D. melanogaster using a panel of over 160 odours. I found that the IRs respond to a number of amines, aldehydes and acids, contrasting with the chemical specificity of the OR repertoire, which is mainly tuned to esters, alcohols and ketones. Finally, the identification of ligands for IRs in this species allowed me to investigate in detail the molecular and functional evolution of a tandem array of IRs, IR75a/IR75b/IR75c, in D. sechellia. This species is endemic to the Seychelles archipelago and highly specialised to breed on the fruits of Morinda citrifolia, which is repulsive and toxic for other Drosophila species. These studies led me to discover that receptor loss, changes in receptor specificity and changes in receptor expression have likely played an important role during the evolution of these IRs in D. sechellia. These changes may explain, in part, the unique chemical ecology of this species. - Le système olfactif est un excellent modèle pour étudier les mécanismes génétiques impliqués dans l'étude de l'évolution du système nerveux. Ce système sensoriel permet la détection de nombreux composés volatils présents dans l'environnement et est à la base des réponses comportementales. Il est propre à chaque espèce et évolue rapidement en modifiant ou en éliminant des récepteurs et leurs circuits olfactifs correspondants pour s'adapter à de nouvelles odeurs. Pour étudier le système olfactif et son évolution, nous avons décidé d'utiliser la drosophile comme modèle. Le séquençage complet de 12 souches de drosophiles habitant différentes niches écologiques permet une analyse génomique conséquente. De plus, l'une de ces espèces Drosophila melanogaster permet la combinaison d'analyses génétiques et électrophysiologiques. En effet, D. melanogaster possède 2 familles distinctes de récepteurs olfactifs qui permettent la détection d'odeurs: les récepteurs olfactifs (ORs) étant les mieux caractérisés et les récepteurs ionotropiques (IRs), plus récemment identifiés. Au cours de ma thèse, j'ai apporté des nouvelles connaissances qui m'ont permis de mieux comprendre les mécanismes génétiques à la base de l'évolution du système olfactif au travers de trois projets différents, mais interdépendants. Premièrement, j'ai réalisé une analyse génomique comparative de l'ensemble des IRs dans les 12 souches de drosophiles séquencées jusqu'à présent. Ceci a montré que les récepteurs olfactifs IRs sont hautement conservés parmi l'ensemble de ces espèces. Au contraire, une grande partie des IRs qui ne sont pas exprimés dans le système olfactif, et qui semblent être des récepteurs gustatifs, sont beaucoup plus variables dans leur séquence et dans le nombre de copie de gènes. Deuxièmement, pour identifier les ligands des récepteurs IRs exprimés par les neurones sensoriels olfactifs, j'ai réalisé une étude électrophysiologique chez D. melanogaster e η testant l'effet de plus de 160 composés chimiques sur les IRs. J'ai trouvé que les IRs répondent à un nombre d'amines, d'aldéhydes et d'acides, contrairement aux récepteurs olfactifs ORs qui eux répondent principalement aux esthers, alcools et cétones. Finalement, l'identification de ligands pour les IRs dans ces espèces m'a permis d'étudier en détail l'évolution fonctionnelle et moléculaire des IR75a/IR75b/IR75c dans D. sechellia. Cette espèce est endémique de l'archipel des Seychelles et se nourrit spécifiquement du fruit Morinda citrifolia qui est répulsif et toxique pour d'autres souches de drosophiles. Ces études m'ont poussé à découvrir que, la perte de IR75a, le changement dans la spécificité de IR75b ainsi que le changement dans l'expression de IR75c ont probablement joué un rôle important dans l'évolution des IRs chez D. sechellia. Ces changements peuvent expliquer, en partie, l'écologie chimique propre à cette espèce. Résumé français large public Le système olfactif permet aux animaux de détecter des milliers de molécules odorantes, les aidant ainsi à trouver de la nourriture, à distinguer si elle est fraîche ou avariée, à trouver des partenaires sexuels, ainsi qu'à éviter les prédateurs. Selon l'environnement et le mode de vie des espèces, le système olfactif doit détecter des odeurs très diverses ; en effet, un moustique qui recherche du sang humain pour se nourrir doit détecter des odeurs bien différentes d'une abeille qui recherche des fleurs. Dans ma thèse, j'ai essayé de comprendre comment les systèmes olfactifs d'une espèce évoluent pour s'adapter aux exigences induites par son environnement. Un très bon modèle pour étudier cela est la drosophile dont les différentes espèces se nichent dans des habitats très divers. Pour ce faire, j'ai étudié les récepteurs olfactifs de différentes espèces de la drosophile. Ces récepteurs sont des protéines qui se lient à des odeurs spécifiques. Lorsqu'ils se lient, ils activent un neurone qui envoie un signal électrique au cerveau. Ce signal est ensuite traité par ce dernier qui indique à la mouche si l'odeur est attractive ou répulsive. J'ai identifié les récepteurs olfactifs de plusieurs espèces de drosophile et étudié s'il y avait des différences entre elles. La plupart des récepteurs sont similaires entre les espèces, cependant dans l'une d'entre elles, certains récepteurs sont différents. Ce fait est particulièrement intéressant car cette espèce de drosophile se nourrit de fruits que les autres espèces n'apprécient pas. Comme nous ne savons pas quels récepteurs se lient à quelles odeurs, j'ai testé un grand nombre de composants odorants. Ceci m'a permis de constater que, effectivement, certains changements produits dans ces récepteurs expliquent pourquoi cette espèce aime particulièrement ces fruits. En outre, mes résultats contribuent à mieux comprendre les changements génétiques qui sont impliqués dans l'évolution du système olfactif.

Real-time quantitative PCR assay for measurement of bird telomeres

We present the application of a real-time quantitative PCR assay, previously developed to measure relative telomere length in humans and mice, to two bird species, the zebra finch Taeniopygia guttata and the Alpine swift Apus melba. This technique is based on the PCR amplification of telomeric (TTAGGG)(n) sequences using specific oligonucleotide primers. Relative telomere length is expressed as the ratio (T/S) of telomere repeat copy number (T) to control single gene copy number (S). This method is particularly useful for comparisons of individuals within species, or where the same individuals are followed longitudinally. We used glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as a single control gene. In both species, we validated our PCR measurements of relative telomere length against absolute measurements of telomere length determined by the conventional method of quantifying telomere terminal restriction fragment (TRF) lengths using both the traditional Southern blot analysis (Alpine swifts) and in gel hybridization (zebra finches). As found in humans and mice, telomere lengths in the same sample measured by TRF and PCR were well correlated in both the Alpine swift and the zebra finch.. Hence, this PCR assay for measurement of bird telomeres, which is fast and requires only small amounts of genomic DNA, should open new avenues in the study of environmental factors influencing variation in telomere length, and how this variation translates into variation in cellular and whole organism senescence.

Molecular profiling of non-small cell lung cancer by histologic subtype

Background: A substantial proportion of NSCLC has been shown to harbour specific molecular alterations affecting tumour proliferation and resulting in sensitivity to inhibition of the corresponding activated oncogenic pathway by targeted therapies. Comprehensive tumor profiling can diagnose such alterations and may identify new alterations opening additional treatment options for all distinct NSCLC subtypes. Methods: Over 6,700 non-small cell lung cancer cases referred to Caris Life Sciences between 2009 and 2014 were evaluated; clinical diagnoses and detailed tumor pathology were collected from referring physicians. Specific profiling was performed per physician request and included a combination of sequencing (Sanger, NGS or pyrosequencing), protein expression (IHC), gene amplification/rearrangement (CISH or FISH), and/or RNA fragment analysis within potential cancer-related genes and pathways. Results: Patients were grouped into cohorts according to histological subtype - adenocarcinoma (AD) (n=4,286), squamous cell carcinoma (SCC) (n=1,280), large cell carcinoma (LCC) (n=153) and bronchioalveolar carcinoma (BAC) (n=94). Protein overexpression of cMET (>2+ in >50% cells) was higher in AD (35.9%) compared to other subgroups (12-20%) while RRM1 and TOP2A levels were lower in AD. ALK or ROS1 were rearranged in 5.3% of patients with AD compared to 3.7% of patients with LCC and 1.2% of patients with SCC. EGFR mutations were found at low prevalence in both the LCC (0%) and SCC cohorts (2.8%) compared to 21% in AD. Similar lower rates of BRAF mutations were observed in the LCC and SCC cohorts compared to AD (0%, 1.1% and 5.1%). Pathway analysis showed activating mutations in the ERK pathway in 40% of patients with AD. Only 10-12% of patients with LCC or SCC had activating mutations in the ERK pathway. Conclusions: Despite the limitations of this retrospective series, we report comprehensive profiling of the largest cohort of NSCLC. Tumor profiling reveals that ADs may be more addicted to the ERK pathway than other histological subtypes. Drugs which target cMET may also have most utility in AD. Full analysis by histological subtype and additional correlative data on protein expression, gene copy number and mutations will be presented.

The prognostic and predictive value of selected biomarkers in colorectal cancer

Tissue-based biomarkers are studied to receive information about the pathologic processes and cancer outcome, and to enable development of patient-tailored treatments. The aim of this study was to investigate the potential prognostic and/or predictive value of selected biomarkers in colorectal cancer (CRC). Group IIA secretory phospholipase A2 (IIA PLA2) expression was assessed in 114 samples presenting different phases of human colorectal carcinogenesis. Securin, Ki-67, CD44 variant 6 (CD44v6), aldehyde dehydrogenase 1 (ALDH1) and β-catenin were studied in a material including 227 rectal carcinoma patients treated with short-course preoperative radiotherapy (RT), long-course preoperative (chemo)RT (CRT) or surgery only. Epidermal growth factor receptor (EGFR) gene copy number (GCN), its heterogeneity in CRC tissue, and association with response to EGFR-targeted antibodies cetuximab and panitumumab were analyzed in a cohort of 76 metastatic CRC. IIA PLA2 expression was decreased in invasive carcinomas compared to adenomas, but did not relate to patient survival. High securin expression after long-course (C)RT and high ALDH1 expression in node-negative rectal cancer were independent adverse prognostic factors, ALDH1 specifically in patients treated with adjuvant chemotherapy. The lack of membranous CD44v6 in the rectal cancer invasive front associated with infiltrative growth pattern and the risk of disease recurrence. Heterogeneous EGFR GCN increase predicted benefit from EGFR-targeted antibodies, also in the chemorefractory patient population. In summary, high securin and ALDH1 protein expression independently relate to poor outcome in subgroups of rectal cancer patients, potentially because of resistance to conventional chemotherapeutics. Heterogeneous increase in EGFR GCN was validated to be a promising predictive factor in the treatment of metastatic CRC.

Variability in the Stability and Productivity of Transfected Genes in Chinese Hamster Ovary (CHO) cells

In the field of biologics production, productivity and stability of the transfected gene of interest are two very important attributes that dictate if a production process is viable. To further understand and improve these two traits, we would need to further our understanding of the factors affecting them. These would include integration site of the gene, gene copy number, cell phenotypic variation and cell environment. As these factors play different parts in the development process, they lead to variable productivity and stability of the transfected gene between clones, the well-known phenomenon of “clonal variation”. A study of this phenomenon and how the various factors contribute to it will thus shed light on strategies to improve productivity and stability in the production cell line. Of the four factors, the site of gene integration appears to be one of the most important. Hence, it is proposed that work is done on studying how different integration sites affect the productivity and stability of transfected genes in the development process. For the study to be more industrially relevant, it is proposed that the Chinese Hamster Ovary dhfr-deficient cell line, CHO-DG44, is used as the model system.