Presequence-Independent Mitochondrial Import of DNA Ligase Facilitates Establishment of Cell Lines with Reduced mtDNA Copy Number

Due to the essential role played by mitochondrial DNA (mtDNA) in cellular physiology and bioenergetics, methods for establishing cell lines with altered mtDNA content are of considerable interest. Here, we report evidence for the existence in mammalian cells of a novel, low- efficiency, presequence-independent pathway for mitochondrial protein import, which facilitates mitochondrial uptake of such proteins as Chlorella virus ligase (ChVlig) and Escherichia coli LigA. Mouse cells engineered to depend on this pathway for mitochondrial import of the LigA protein for mtDNA maintenance had severely (up to >90%) reduced mtDNA content. These observations were used to establish a method for the generation of mouse cell lines with reduced mtDNA copy number by, first, transducing them with a retrovirus encoding LigA, and then inactivating in these transductants endogenous Lig3 with CRISPR-Cas9. Interestingly, mtDNA depletion to an average level of one copy per cell proceeds faster in cells engineered to maintain mtDNA at low copy number. This makes a low-mtDNA copy number phenotype resulting from dependence on mitochondrial import of DNA ligase through presequence-independent pathway potentially useful for rapidly shifting mtDNA heteroplasmy through partial mtDNA depletion.     

Age-Related and Heteroplasmy-Related Variation in Human mtDNA Copy Number

The mitochondrial (mt) genome is present in many copies in human cells, and intra-individual variation in mtDNA sequences is known as heteroplasmy. Recent studies found that heteroplasmies are highly tissue-specific, site-specific, and allele-specific, however the functional implications have not been explored. This study investigates variation in mtDNA copy numbers (mtCN) in 12 different tissues obtained at autopsy from 152 individuals (ranging in age from 3 days to 96 years). Three different methods to estimate mtCN were compared: shotgun sequencing (in 4 tissues), capture-enriched sequencing (in 12 tissues) and droplet digital PCR (ddPCR, in 2 tissues). The highest precision in mtCN estimation was achieved using shotgun sequencing data. However, capture-enrichment data provide reliable estimates of relative (albeit not absolute) mtCNs. Comparisons of mtCN from different tissues of the same individual revealed that mtCNs in different tissues are, with few exceptions, uncorrelated. Hence, each tissue of an individual seems to regulate mtCN in a tissue-related rather than an individual-dependent manner. Skeletal muscle (SM) samples showed an age-related decrease in mtCN that was especially pronounced in males, while there was an age-related increase in mtCN for liver (LIV) samples. MtCN in SM samples was significantly negatively correlated with both the total number of heteroplasmic sites and with minor allele frequency (MAF) at two heteroplasmic sites, 408 and 16327. Heteroplasmies at both sites are highly specific for SM, accumulate with aging and are part of functional elements that regulate mtDNA replication. These data support the hypothesis that selection acting on these heteroplasmic sites is reducing mtCN in SM of older individuals.     

Host Genetic Variants and Gene Expression Patterns Associated with Epstein-Barr Virus Copy Number in Lymphoblastoid Cell Lines

Lymphoblastoid cell lines (LCLs) are commonly used in molecular genetics, supplying DNA for the HapMap and 1000 Genomes Projects, used to test chemotherapeutic agents, and informing the basis of a number of population genetics studies of gene expression. The process of transforming human B cells into LCLs requires the presence of Epstein-Barr virus (EBV), a double-stranded DNA virus which through B-cell immortalisation maintains an episomal virus genome in every cell of an LCL at variable copy numbers. Previous studies have reported that EBV alters host-gene expression and EBV copy number may be under host genetic control. We performed a genome-wide association study of EBV genome copy number in LCLs and found the phenotype to be highly heritable, although no individual SNPs achieved a significant association with EBV copy number. The expression of two host genes (CXCL16 and AGL) was positively correlated and expression of ADARB2 was negatively correlated with EBV copy number in a genotype-independent manner. This study shows an association between EBV copy number and the gene expression profile of LCLs, and suggests that EBV copy number should be considered as a covariate in future studies of host gene expression in LCLs.     

Proteomic Changes Resulting from Gene Copy Number Variations in Cancer Cells

Along the transformation process, cells accumulate DNA aberrations, including mutations, translocations, amplifications, and deletions. Despite numerous studies, the overall effects of amplifications and deletions on the end point of gene expression—the level of proteins—is generally unknown. Here we use large-scale and high-resolution proteomics combined with gene copy number analysis to investigate in a global manner to what extent these genomic changes have a proteomic output and therefore the ability to affect cellular transformation. We accurately measure expression levels of 6,735 proteins and directly compare them to the gene copy number. We find that the average effect of these alterations on the protein expression is only a few percent. Nevertheless, by using a novel algorithm, we find the combined impact that many of these regional chromosomal aberrations have at the protein level. We show that proteins encoded by amplified oncogenes are often overexpressed, while adjacent amplified genes, which presumably do not promote growth and survival, are attenuated. Furthermore, regulation of biological processes and molecular complexes is independent of general copy number changes. By connecting the primary genome alteration to their proteomic consequences, this approach helps to interpret the data from large-scale cancer genomics efforts.     

The Role of Cell Growth-Related Gene Copy Number Variation in Autoimmune Thyroid Disease

Biological Trace Element Research - Autoimmune thyroid disease (AITD) is a recurrent and refractory clinical endocrine disease. Some studies have shown that the incidence of AITD is not only...     

Mitochondrial DNA Copy Number in Peripheral Blood and Melanoma Risk

Mitochondrial DNA (mtDNA) copy number in peripheral blood has been suggested as risk modifier in various types of cancer. However, its influence on melanoma risk is unclear. We evaluated the association between mtDNA copy number in peripheral blood and melanoma risk in 500 melanoma cases and 500 healthy controls from an ongoing melanoma study. The mtDNA copy number was measured using real-time polymerase chain reaction. Overall, mean mtDNA copy number was significantly higher in cases than in controls (1.15 vs 0.99, P<0.001). Increased mtDNA copy number was associated with a 1.45-fold increased risk of melanoma (95% confidence interval: 1.12-1.97). Significant joint effects between mtDNA copy number and variables related to pigmentation and history of sunlight exposure were observed. This study supports an association between increased mtDNA copy number and melanoma risk that is independent on the known melanoma risk factors (pigmentation and history of sunlight exposure).     

mtDNA拷贝数的生物学意义及其调控

线粒体是细胞能量和自由基代谢中心,并在细胞凋亡、钙调控、细胞周期和信号转导中发挥重要作用,维持线粒体功能正常对于细胞正常行使职能意义重大。线粒体的功能与线粒体DNA（mitochondrial DNA,mtDNA）的数量和质量紧密相关,mtDNA的数量即mtDNA拷贝数又受到mtDNA质量的影响,因此mtDNA拷贝数可作为线粒体功能的重要表征。mtDNA拷贝数变异引起线粒体功能紊乱,进而导致疾病发生。本文综述了mtDNA拷贝数变异与神经退行性疾病、心血管疾病、肿瘤等疾病的发生发展和个体衰老之间的关系,以及mtDNA复制转录相关因子、氧化应激、细胞自噬等因素介导mtDNA拷贝数变异的调控机制。以期为进一步深入探究mtDNA拷贝数调控的分子机制,以及未来治疗神经退行性疾病、肿瘤及延缓衰老等提供一定的理论基础。     

Haplotype Phasing and Inheritance of Copy Number Variants in Nuclear Families

DNA copy number variants (CNVs) that alter the copy number of a particular DNA segment in the genome play an important role in human phenotypic variability and disease susceptibility. A number of CNVs overlapping with genes have been shown to confer risk to a variety of human diseases thus highlighting the relevance of addressing the variability of CNVs at a higher resolution. So far, it has not been possible to deterministically infer the allelic composition of different haplotypes present within the CNV regions. We have developed a novel computational method, called PiCNV, which enables to resolve the haplotype sequence composition within CNV regions in nuclear families based on SNP genotyping microarray data. The algorithm allows to i) phase normal and CNV-carrying haplotypes in the copy number variable regions, ii) resolve the allelic copies of rearranged DNA sequence within the haplotypes and iii) infer the heritability of identified haplotypes in trios or larger nuclear families. To our knowledge this is the first program available that can deterministically phase null, mono-, di-, tri- and tetraploid genotypes in CNV loci. We applied our method to study the composition and inheritance of haplotypes in CNV regions of 30 HapMap Yoruban trios and 34 Estonian families. For 93.6% of the CNV loci, PiCNV enabled to unambiguously phase normal and CNV-carrying haplotypes and follow their transmission in the corresponding families. Furthermore, allelic composition analysis identified the co-occurrence of alternative allelic copies within 66.7% of haplotypes carrying copy number gains. We also observed less frequent transmission of CNV-carrying haplotypes from parents to children compared to normal haplotypes and identified an emergence of several de novo deletions and duplications in the offspring.     

Copy Number Variants in the Kallikrein Gene Cluster

The kallikrein gene family (KLK1-KLK15) is the largest contiguous group of protease genes within the human genome and is associated with both risk and outcome of cancer and other diseases. We searched for copy number variants in all KLK genes using quantitative PCR analysis and analysis of inheritance patterns of single nucleotide polymorphisms. Two deletions were identified: one 2235-bp deletion in KLK9 present in 1.2% of alleles, and one 3394-bp deletion in KLK15 present in 4.0% of alleles. Each deletion eliminated one complete exon and created out-of-frame coding that eliminated the catalytic triad of the resulting truncated gene product, which therefore likely is a non-functional protein. Deletion breakpoints identified by DNA sequencing located the KLK9 deletion breakpoint to a long interspersed element (LINE) repeated sequence, while the deletion in KLK15 is located in a single copy sequence. To search for an association between each deletion and risk of prostate cancer (PC), we analyzed a cohort of 667 biopsied men (266 PC cases and 401 men with no evidence of PC at biopsy) using short deletion-specific PCR assays. There was no association between evidence of PC in this cohort and the presence of either gene deletion. Haplotyping revealed a single origin of each deletion, with most recent common ancestor estimates of 3000-8000 and 6000-14 000 years for the deletions in KLK9 and KLK15, respectively. The presence of the deletions on the same haplotypes in 1000 Genomes data of both European and African populations indicate an early origin of both deletions. The old age in combination with homozygous presence of loss-of-function variants suggests that some kallikrein-related peptidases have non-essential functions.     

增加植酸酶基因appA-m的拷贝提高其在巴斯德毕赤酵母的表达量

为了进一步提高植酸酶的发酵效价,降低植酸酶生产成本,对毕赤酵母表达载体pGAPZα-A进行了改造。将表达载体pPIC9的AOX1启动子序列引入pGAPZα-A,使之成为甲醇可诱导型表达载体pAOXZα,插入植酸酶基因appA-m后得到重组载体pAOXZα-appA-m。以染色体上带有一个拷贝的appA-m基因、发酵效价可达到7.5×106IU/mL发酵液的重组酵母菌株74#为受体菌进行转化,在该重组菌株的染色体上的另一位点整合含有植酸酶基因的表达盒,经筛选到高表达植酸酶的重组子。通过PCR进行验证,植酸酶基因被整合到重组酵母的染色体上,且受体菌中原有的植酸酶基因结构未改变。重组菌在5L发酵罐经甲醇诱导120h植酸酶蛋白表达量达到4mg/mL发酵液,酶活性(发酵效价)达到1.2×107IU/mL发酵液以上,较含单拷贝植酸酶基因的受体菌株表达量有较大程度提高。PCR检测及表达量分析证明改良的菌株具有很好的遗传稳定性和表达稳定性。     

An Atlas of the Speed of Copy Number Changes in Animal Gene Families and Its Implications

The notion that gene duplications generating new genes and functions is commonly accepted in evolutionary biology. However, this assumption is more speculative from theory rather than well proven in genome-wide studies. Here, we generated an atlas of the rate of copy number changes (CNCs) in all the gene families of ten animal genomes. We grouped the gene families with similar CNC dynamics into rate pattern groups (RPGs) and annotated their function using a novel bottom-up approach. By comparing CNC rate patterns, we showed that most of the species-specific CNC rates groups are formed by gene duplication rather than gene loss, and most of the changes in rates of CNCs may be the result of adaptive evolution. We also found that the functions of many RPGs match their biological significance well. Our work confirmed the role of gene duplication in generating novel phenotypes, and the results can serve as a guide for researchers to connect the phenotypic features to certain gene duplications.     

Spectrum of EGFR Gene Copy Number Changes and KRAS Gene Mutation Status in Korean Triple Negative Breast Cancer Patients

Anti-epidermal growth factor receptor (EGFR) therapy has been tried in triple negative breast cancer (TNBC) patients without evaluation of molecular and clinical predictors in several randomized clinical studies. Only fewer than 20% of metastatic TNBCs showed response to anti-EGFR therapy. In order to increase the overall response rate, first step would be to classify TNBC into good or poor responders according to oncogenic mutation profiles. This study provides the molecular characteristics of TNBCs including EGFR gene copy number changes and mutation status of EGFR and KRAS gene in Korean TNBC patients. Mutation analysis for EGFR, KRAS, BRAF and TP53 from a total of 105 TNBC tissue samples was performed by direct sequencing, peptide nucleic acid-mediated PCR clamping method and real-time PCR. Copy number changes of EGFR gene were evaluated using multiplex ligation-dependent probe amplification. Out of all 105 TNBCs, 15.2% (16/105) showed EGFR copy number changes. Among them, increased or decreased EGFR copy number was detected in 13 (5 single copy gain, 2 amplification and 4 high-copy number amplification) and 3 cases (3 hemizygous deletion), respectively. The mutation frequencies of KRAS, EGFR and TP53 gene were 1.9% (G12V and G12D), 1.0% (exon 19 del) and 31.4%, respectively. There was no BRAF V600E mutation found. Future studies are needed to evaluate the clinical outcomes of TNBC patients who undergo anti-EGFR therapy according to the genetic status of EGFR.     

整合分析基因表达与拷贝数变异识别癌症的驱动基因及调控子miRNAs

目的:通过整合分析基因的表达与拷贝数变异(CNV)识别癌症的驱动基因及调控子mi RNAs。方法:通过整合基因表达与CNV数据,分别计算了乳腺癌、结肠癌、肺癌、肾癌、膀胱癌、头颈癌六种癌症中mi RNAs的调控得分,提出了一个识别驱动基因和显著调控子mi RNAs的方法。结果:本文研究发现,CNV区域上编码的基因相比于非CNV区域上编码的基因更倾向于受mi RNAs调控。但是,癌相关CNV区域上的基因相比正常CNV区域上的基因更少受mi RNAs调控。本研究识别出了EXOSC4、ZNF7、BOP1等原癌基因,以及mi R-488、mi R-27a、mi R-454等在多种癌症中都起调控作用的调控子mi RNAs。结论:本文的方法为癌症研究带来了新的启发,这些具有调控扩增基因过表达作用的mi RNAs的发现,有助于我们更进一步了解癌基因表达的复杂调控机制,进而推动癌症的诊断、治疗和预后。     

基因拷贝数与染色体位置对酵母表达乙肝病毒融合表面抗原SA-28基因的影响

应用ＦＬＰ重组酶介导的染色体定点整合技术,将带有不同拷贝数的乙肝病毒融合表面抗原ＳＡ－２８基因表达单元的质粒整合在酵母不同的染色体位点,并测定了ＳＡ－２８基因的表达情况,初步研究了基因拷贝数与染色体位置对酵母表达外源基因的影响。结果表明ＳＡ－２８基因在ＨＩＳ３位点整 合时的表达水平随基因拷贝数的增加而提高,遵循基因剂量效应;在某些染色体位点整２合时,插入方向对其表达有不同程度的影响,呈现出明显的染     

Whole Cell Formaldehyde Cross-Linking Simplifies Purification of Mitochondrial Nucleoids and Associated Proteins Involved in Mitochondrial Gene Expression

Mitochondrial DNA/protein complexes (nucleoids) appear as discrete entities inside the mitochondrial network when observed by live-cell imaging and immunofluorescence. This somewhat trivial observation in recent years has spurred research towards isolation of these complexes and the identification of nucleoid-associated proteins. Here we show that whole cell formaldehyde crosslinking combined with affinity purification and tandem mass-spectrometry provides a simple and reproducible method to identify potential nucleoid associated proteins. The method avoids spurious mitochondrial isolation and subsequent multifarious nucleoid enrichment protocols and can be implemented to allow for label-free quantification (LFQ) by mass-spectrometry. Using expression of a Flag-tagged Twinkle helicase and appropriate controls we show that this method identifies many previously identified nucleoid associated proteins. Using LFQ to compare HEK293 cells with and without mtDNA, but both expressing Twinkle-FLAG, identifies many proteins that are reduced or absent in the absence of mtDNA. This set not only includes established mtDNA maintenance proteins but also many proteins involved in mitochondrial RNA metabolism and translation and therefore represents what can be considered an mtDNA gene expression proteome. Our data provides a very valuable resource for both basic mitochondrial researchers as well as clinical geneticists working to identify novel disease genes on the basis of exome sequence data.