The Acquisition of an Inheritable 50-bp Deletion in the Human mtDNA Control Region Does Not Affect the mtDNA Copy Number in Peripheral Blood Cells

The mitochondrial DNA (mtDNA) control region is believed to play an important biological role in mtDNA replication. Large deletions in this region are rarely found, but when they do occur they might be expected to interfere with the replication of the molecule, thus leading to a reduction of mtDNA copy number. During a survey for mtDNA sequence variations in 5,559 individuals from the general Chinese population and 2,538 individuals with medical disorders, we identified a 50-bp deletion (m.298_347del50) in the mtDNA control region in a member of a healthy Han Chinese family belonging to haplogroup B4c1b2, as suggested by complete mtDNA genome sequencing. This deletion removes the conserved sequence block II (CSBII; region 299-315) and the replication primer location (region 317-321). However, quantification of the mtDNA copy number in this subject showed a value within a range that was observed in 20 healthy subjects without the deletion. The deletion was detected in the hair samples of the maternal relatives of the subject and exhibited variable heteroplasmy. Our current observation, together with a recent report for a benign 154-bp deletion in the mtDNA control region, suggests that the control of mtDNA replication may be more complex than we had thought. Hum Mutat 31:538-543, 2010. (C) 2010 Wiley-Liss, Inc.

Plasmi copy number quantification by spectrofluorometry

No data (2012)

Quantitative Analysis of Single Nucleotide Polymorphisms within Copy Number Variation

Background Single nucleotide polymorphisms (SNPs) have been used extensively in genetics and epidemiology studies. Traditionally, SNPs that did not pass the Hardy-Weinberg equilibrium (HWE) test were excluded from these analyses. Many investigators have addressed possible causes for departure from HWE, including genotyping errors, population admixture and segmental duplication. Recent large-scale surveys have revealed abundant structural variations in the human genome, including copy number variations (CNVs). This suggests that a significant number of SNPs must be within these regions, which may cause deviation from HWE. Results We performed a Bayesian analysis on the potential effect of copy number variation, segmental duplication and genotyping errors on the behavior of SNPs. Our results suggest that copy number variation is a major factor of HWE violation for SNPs with a small minor allele frequency, when the sample size is large and the genotyping error rate is 0~1%. Conclusions Our study provides the posterior probability that a SNP falls in a CNV or a segmental duplication, given the observed allele frequency of the SNP, sample size and the significance level of HWE testing.

The DUB/USP17 deubiquitinating enzymes: A gene family within a tandemly repeated sequence, is also embedded within the copy number variable Beta-defensin cluste

Background: The DUB/USP17 subfamily of deubiquitinating enzymes were originally identified as immediate early genes induced in response to cytokine stimulation in mice (DUB-1, DUB-1A, DUB-2, DUB-2A). Subsequently we have identified a number of human family members and shown that one of these (DUB-3) is also cytokine inducible. We originally showed that constitutive expression of DUB-3 can block cell proliferation and more recently we have demonstrated that this is due to its regulation of the ubiquitination and activity of the 'CAAX' box protease RCE1.

Results: Here we demonstrate that the human DUB/USP17 family members are found on both chromosome 4p16.1, within a block of tandem repeats, and on chromosome 8p23.1, embedded within the copy number variable betadefensin cluster. In addition, we show that the multiple genes observed in humans and other distantly related mammals have arisen due to the independent expansion of an ancestral sequence within each species. However, it is also apparent when sequences from humans and the more closely related chimpanzee are compared, that duplication events have taken place prior to these species separating.

Conclusions: The observation that the DUB/USP17 genes, which can influence cell growth and survival, have evolved from an unstable ancestral sequence which has undergone multiple and varied duplications in the species examined marks this as a unique family. In addition, their presence within the beta-defensin repeat raises the question whether they may contribute to the influence of this repeat on immune related conditions.

Targets of genome copy number reduction in primary breast cancers identified by integrative genomics

The identification of specific oncogenes and tumor suppressor genes in regions of recurrent aneuploidy is a major challenge of molecular cancer research. Using both oligonucleotide single-nucleotide polymorphism and mRNA expression arrays, we integrated genomic and transcriptional information to identify and prioritize candidate cancer genes in regions of increased and decreased chromosomal copy number in a cohort of primary breast cancers. Confirming the validity of this approach, several regions of previously-known copy number (CN) alterations in breast cancer could be successfully reidentified. Focusing on regions of decreased CN, we defined a prioritized list of eighteen candidate genes, which included ARPIN, FBNI, and LZTSI, previously shown to be associated with cancers in breast or other tissue types, and novel genes such as P29, MORF4LI, and TBCID5. One such gene, the RUNX3 transcription factor, was selected for further study. We show that RUNX3 is present at reduced CNs in proportion to the rest of the tumor genome and that RUNX3 CN reductions can also be observed in a breast cancer series from a different center. Using tissue microarrays, we demonstrate in an independent cohort of over 120 breast tissues that RUNX3 protein is expressed in normal breast epithelium but not fat and stromal tissue, and widely down-regulated in the majority of breast cancers (> 85%). In vitro, RUNX3 overexpression suppressed the invasive potential of MDA-MB-231 breast cancer cells in a matrigel assay. Our results demonstrate the utility of integrative genomic approaches to identify novel potential cancer-related genes in primary tumors. This article contains Supplementary Material available at http:// www.interscience.wiley.com/jpages/1045-2257/suppmat. (c) 2006 Wiley-Liss, Inc.

Genome-wide association of early-onset myocardial infarction with single nucleotide polymorphisms and copy number variants

We conducted a genome-wide association study testing single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) for association with early-onset myocardial infarction in 2,967 cases and 3,075 controls. We carried out replication in an independent sample with an effective sample size of up to 19,492. SNPs at nine loci reached genome-wide significance: three are newly identified (21q22 near MRPS6-SLC5A3-KCNE2, 6p24 in PHACTR1 and 2q33 in WDR12) and six replicated prior observations1-4 (9p21, 1p13 near CELSR2-PSRC1-SORT1, 10q11 near CXCL12, 1q41 in MIA3, 19p13 near LDLR and 1p32 near PCSK9). We tested 554 common copy number polymorphisms (>1% allele frequency) and none met the pre-specified threshold for replication (P < 10-3). We identified 8,065 rare CNVs but did not detect a greater CNV burden in cases compared to controls, in genes compared to the genome as a whole, or at any individual locus. SNPs at nine loci were reproducibly associated with myocardial infarction, but tests of common and rare CNVs failed to identify additional associations with myocardial infarction risk.

Discovering subgroups of patients from DNA copy number data using NMF on compacted matrices

In the study of complex genetic diseases, the identification of subgroups of patients sharing similar genetic characteristics represents a challenging task, for example, to improve treatment decision. One type of genetic lesion, frequently investigated in such disorders, is the change of the DNA copy number (CN) at specific genomic traits. Non-negative Matrix Factorization (NMF) is a standard technique to reduce the dimensionality of a data set and to cluster data samples, while keeping its most relevant information in meaningful components. Thus, it can be used to discover subgroups of patients from CN profiles. It is however computationally impractical for very high dimensional data, such as CN microarray data. Deciding the most suitable number of subgroups is also a challenging problem. The aim of this work is to derive a procedure to compact high dimensional data, in order to improve NMF applicability without compromising the quality of the clustering. This is particularly important for analyzing high-resolution microarray data. Many commonly used quality measures, as well as our own measures, are employed to decide the number of subgroups and to assess the quality of the results. Our measures are based on the idea of identifying robust subgroups, inspired by biologically/clinically relevance instead of simply aiming at well-separated clusters. We evaluate our procedure using four real independent data sets. In these data sets, our method was able to find accurate subgroups with individual molecular and clinical features and outperformed the standard NMF in terms of accuracy in the factorization fitness function. Hence, it can be useful for the discovery of subgroups of patients with similar CN profiles in the study of heterogeneous diseases.

Integration of copy number and transcriptomics provides risk stratification in prostate cancer: A discovery and validation cohort study

BACKGROUND: Understanding the heterogeneous genotypes and phenotypes of prostate cancer is fundamental to improving the way we treat this disease. As yet, there are no validated descriptions of prostate cancer subgroups derived from integrated genomics linked with clinical outcome.

METHODS: In a study of 482 tumour, benign and germline samples from 259 men with primary prostate cancer, we used integrative analysis of copy number alterations (CNA) and array transcriptomics to identify genomic loci that affect expression levels of mRNA in an expression quantitative trait loci (eQTL) approach, to stratify patients into subgroups that we then associated with future clinical behaviour, and compared with either CNA or transcriptomics alone.

FINDINGS: We identified five separate patient subgroups with distinct genomic alterations and expression profiles based on 100 discriminating genes in our separate discovery and validation sets of 125 and 103 men. These subgroups were able to consistently predict biochemical relapse (p = 0.0017 and p = 0.016 respectively) and were further validated in a third cohort with long-term follow-up (p = 0.027). We show the relative contributions of gene expression and copy number data on phenotype, and demonstrate the improved power gained from integrative analyses. We confirm alterations in six genes previously associated with prostate cancer (MAP3K7, MELK, RCBTB2, ELAC2, TPD52, ZBTB4), and also identify 94 genes not previously linked to prostate cancer progression that would not have been detected using either transcript or copy number data alone. We confirm a number of previously published molecular changes associated with high risk disease, including MYC amplification, and NKX3-1, RB1 and PTEN deletions, as well as over-expression of PCA3 and AMACR, and loss of MSMB in tumour tissue. A subset of the 100 genes outperforms established clinical predictors of poor prognosis (PSA, Gleason score), as well as previously published gene signatures (p = 0.0001). We further show how our molecular profiles can be used for the early detection of aggressive cases in a clinical setting, and inform treatment decisions.

INTERPRETATION: For the first time in prostate cancer this study demonstrates the importance of integrated genomic analyses incorporating both benign and tumour tissue data in identifying molecular alterations leading to the generation of robust gene sets that are predictive of clinical outcome in independent patient cohorts.

Copy number variation and Huntington’s disease

Tese de doutoramento, Ciências Biomédicas (Neurociências), Universidade de Lisboa, Faculdade de Medicina, 2014

Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

The Tnt1 family member Retrosol copy number and structure disclose retrotransposon diversification in different Solanum species

Eukaryotic genome expansion/retraction caused by LTR-retrotransposon activity is dependent on the expression of full length copies to trigger efficient transposition and recombination-driven events. The Tnt1 family of retrotransposons has served as a model to evaluate the diversity among closely related elements within Solanaceae species and found that members of the family vary mainly in their U3 region of the long terminal repeats (LTRs). Recovery of a full length genomic copy of Retrosol was performed through a PCR-based approach from wild potato, Solanum oplocense. Further characterization focusing on both LTR sequences of the amplified copy allowed estimating an approximate insertion time at 2 million years ago thus supporting the occurrence of transposition cycles after genus divergence. Copy number of Tnt1-like elements in Solanum species were determined through genomic quantitative PCR whereby results sustain that Retrosol in Solanum species is a low copy number retrotransposon (1-4 copies) while Retrolyc1 has an intermediate copy number (38 copies) in S. peruvianum. Comparative analysis of retrotransposon content revealed no correlation between genome size or ploidy level and Retrosol copy number. The tetraploid cultivated potato with a cellular genome size of 1,715 Mbp harbours similar copy number per monoploid genome than other diploid Solanum species (613-884 Mbp). Conversely, S. peruvianum genome (1,125 Mbp) has a higher copy number. These results point towards a lineage specific dynamic flux regarding the history of amplification/activity of Tnt1-like elements in the genome of Solanum species.

Copy number variants on the X chromosome in women with primary ovarian insufficiency

Objective: To investigate whether submicroscopic copy number variants (CNVs) on the X chromosome can be identified in women with primary ovarian insufficiency (POI), defined as spontaneous secondary amenorrhea before 40 years of age accompanied by follicle-stimulating hormone levels above 40 IU/L on at least two occasions. Design: Analysis of intensity data of single nucleotide polymorphism (SNP) probes generated by genomewide Illumina 370k CNV BeadChips, followed by the validation of identified loci using a custom designed ultra-high-density comparative genomic hybridization array containing 48,325 probes evenly distributed over the X chromosome. Setting: Multicenter genetic cohort study in the Netherlands. Patient(s): 108 Dutch Caucasian women with POI, 97 of whom passed quality control, who had a normal karyogram and absent fragile X premutation, and 235 healthy Dutch Caucasian women as controls. Intervention(s): None. Main Outcome Measure(s): Amount and locus of X chromosomal microdeletions or duplications. Result(s): Intensity differences between SNP probes identify microdeletions and duplications. The initial analysis identified an overrepresentation of deletions in POI patients. Moreover, CNVs in two genes on the Xq21.3 locus (i.e., PCDH11X and TGIF2LX) were statistically significantly associated with the POI phenotype. Mean size of identified CNVs was 262 kb. However, in the validation study the identified putative Xq21.3 deletions samples did not show deviations in intensities in consecutive probes. Conclusion(s): X chromosomal submicroscopic CNVs do not play a major role in Caucasian POI patients. We provide guidelines on how submicroscopic cytogenetic POI research should be conducted. (Fertil Steril (R) 2011;95:1584-8. (C) 2011 by American Society for Reproductive Medicine.)

Copy number variation and transposable elements feature in recent, ongoing adaptation at the Cyp6g1 locus

The increased transcription of the Cyp6g1 gene of Drosophila melanogaster, and consequent resistance to insecticides such as DDT, is a widely cited example of adaptation mediated by cis-regulatory change. A fragment of an Accord transposable element inserted upstream of the Cyp6g1 gene is causally associated with resistance and has spread to high frequencies in populations around the world since the 1940s. Here we report the existence of a natural allelic series at this locus of D. melanogaster, involving copy number variation of Cyp6g1, and two additional transposable element insertions (a P and an HMS-Beagle). We provide evidence that this genetic variation underpins phenotypic variation, as the more derived the allele, the greater the level of DDT resistance. Tracking the spatial and temporal patterns of allele frequency changes indicates that the multiple steps of the allelic series are adaptive. Further, a DDT association study shows that the most resistant allele, Cyp6g1-[BP], is greatly enriched in the top 5% of the phenotypic distribution and accounts for ~16% of the underlying phenotypic variation in resistance to DDT. In contrast, copy number variation for another candidate resistance gene, Cyp12d1, is not associated with resistance. Thus the Cyp6g1 locus is a major contributor to DDT resistance in field populations, and evolution at this locus features multiple adaptive steps occurring in rapid succession.

Recurrent copy number gains of ACVR1 and corresponding transcript overexpression are associated with survival in head and neck squamous cell carcinomas

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)