新一代测序的拷贝数变异检测算法研究与设计

基于不同的测序技术,基因拷贝数变异的检测方法有多种,但时间复杂度较高,而新一代测序技术的发展为基因拷贝数变异检测的研究开辟了新领域。通过仿真实验、置换检验设计出一种新的基于新一代测序的拷贝数变异检测算法。不同于其它算法,本算法无需参考样本,通过直接研究比对后的序列以及reads与拷贝数的关系,来研究检测拷贝数变异,实验结果表明在时间复杂度上能提高50%以上的运算速度,这对今后拷贝数与疾病的研究具有重要意义。     

拷贝数变异: 基因组多样性的新形式

基因拷贝数变异是指DNA片段大小范围从kb到Mb的亚微观突变, 是一可能具有致病性、良性或未知临床意义的基因组改变。Fosmid末端配对序列比较策略、比较基因组杂交芯片是当前较多使用的检测手段。染色体非等位的同源重排、非同源突变和非b DNA结构是造成基因组拷贝数变异的重要原因。拷贝数变异可导致不同程度的基因表达差异, 对正常表型的构成及疾病的发生发展具有一定作用。文章在总结基因拷贝数变异的认识过程和研究策略的基础上, 分析了拷贝数变异的形成和作用机制, 介绍了第一代人类基因组拷贝数变异图谱, 阐述了拷贝数变异研究的临床意义, 提示在探索疾病相关的遗传变异时不能错失拷贝数变异这一基因组多样性的新形式。     

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

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

高通量测序数据的基因组拷贝数变异检测方法综述

拷贝数变异是指基因组中发生大片段的DNA序列的拷贝数增加或者减少。根据现有的研究可知,拷贝数变异是多种人类疾病的成因,与其发生与发展机制密切相关。高通量测序技术的出现为拷贝数变异检测提供了技术支持,在人类疾病研究、临床诊疗等领域,高通量测序技术已经成为主流的拷贝数变异检测技术。虽然不断有新的基于高通量测序技术的算法和软件被人们开发出来,但是准确率仍然不理想。本文全面地综述基于高通量测序数据的拷贝数变异检测方法,包括基于reads深度的方法、基于双末端映射的方法、基于拆分read的方法、基于从头拼接的方法以及基于上述4种方法的组合方法,深入探讨了每类不同方法的原理,代表性的软件工具以及每类方法适用的数据以及优缺点等,并展望未来的发展方向。     

CNV及其与癌症间的关系

拷贝数变异(copy number variation,CNV)是指基因组发生1 Kb 以上的DNA片段的增添、缺失或重排。癌症的早期诊断与治疗一直是本世纪亟待解决的难题。CNV的相关研究为人类健康和疾病的治疗提供了宝贵的见解。目前,CNV的研究引发了人们对疾病的新探索,尤其体现在与遗传物质息息相关的疾病(例如,癌症)的病因研究、临床诊断、新药研发和治疗。该文主要综述了CNV的研究方法、形成机制以及其与癌症间的联系,以期推动癌症相关研究的发展。     

基因拷贝数变异与人类疾病

拷贝数变异(copy number variation,CNV)是人类遗传多样性的一类重要形式。在前期的研究中,人们通过寡核苷酸分型、比较基因组杂交以及测序等技术手段,在人类基因组中鉴定出了大量拷贝数变异位点。这些变异可能是由于基因组重组或复制过程中的差错而产生。CNV在人群中的覆盖率远远高于寡核苷酸多态性(single nucleotide polymorphism,SNP),它们可以通过多种机制改变基因的表达水平,如基因剂量效应、基因断裂-融合效应,以及远距调控效应,进而引起多种人类复杂疾病。认识基因组中的拷贝数变异对于我们更好地认识基因与疾病的关系、遗传-环境因素的相互作用,以及基因组变异与物种进化的关系具有重要的意义。     

复杂疾病全基因组关联研究进展—— 研究设计和遗传标记

实现全基因组关联研究(Genome-wide association study, GWA)在数年前还是遗传学家们的梦想, 如今它已经变成了现实。自2005年Science杂志报道了第一项有关年龄相关性(视网膜)黄斑变性全基因组关联研究研究以来, 有关与复杂疾病的全基因组关联研究如雨后春笋般层出不穷。文中介绍了近两年来全基因组关联研究在复杂疾病研究领域内的主要发现、全基因组关联研究设计原理、遗传标记的选择、比较及相关商品信息。最后介绍了人类基因组拷贝数变异的研究进展, 总结了人类全基因组关联研究所取得成就和存在的问题, 并对全基因组关联研究未来的研究重点和要解决的问题进行了展望。     

拷贝数变异与疾病的关系及其在动物抗病育种中的应用前景

摘要: 拷贝数变异(Copy number variations, CNVs)主要指大于1 kb以上的DNA片段的缺失、插入、重复等。CNVs广泛存在于人类和其他哺乳动物的基因组中。文章主要介绍了CNVs对人类疾病的影响及其检测技术, 并对CNVs在动物抗病育种中的应用前景进行了展望。由于拷贝数变异对抗病性和易感性的影响至关重要, 因此采用生物技术手段有望将其运用于家畜标记辅助选择、QTL精细定位以及动物优良抗病品种培育当中。     

基因组拷贝数变异及其突变机理与人类疾病

拷贝数变异(Copy number variation,CNV)是由基因组发生重排而导致的,一般指长度为1 kb以上的基因组大片段的拷贝数增加或者减少,主要表现为亚显微水平的缺失和重复。CNV是基因组结构变异(Structural variation,SV)的重要组成部分。CNV位点的突变率远高于SNP(Single nucleotide polymorphism),是人类疾病的重要致病因素之一。目前,用来进行全基因组范围的CNV研究的方法有:基于芯片的比较基因组杂交技术(array-based comparative genomic hybridization,aCGH)、SNP分型芯片技术和新一代测序技术。CNV的形成机制有多种,并可分为DNA重组和DNA错误复制两大类。CNV可以导致呈孟德尔遗传的单基因病与罕见疾病,同时与复杂疾病也相关。其致病的可能机制有基因剂量效应、基因断裂、基因融合和位置效应等。对CNV的深入研究,可以使我们对人类基因组的构成、个体间的遗传差异、以及遗传致病因素有新的认识。     

一种基于二代测序拷贝数变异检测的新方法

区域捕获测序是针对基因组特定区段如对MHC(Major histocompatibility complex)区域、外显子区域等测序的有效手段,但是由于捕获测序中探针设计不均匀而造成区域内测序深度变异很大,因此,与基于全基因组的测序数据相比,其拷贝数变异的检测难度更大.目前已经出现了捕获测序下拷贝数变异(copy number variations,CNV)的检测方法,但对CNV的检测准确性仍然很低,特别是对于低频率CNV来说效果极差.因此,本研究开发了一个新的拷贝数变异检测方法,其特点是:(1)以区域内划分的区间为单位检测区间内的CNV,而不是直接对每个个体检测CNV;(2)全面利用群体内所有个体信息,通过区间内read深度在群体的分布规律来检测CNV的分离规律,假设区间内只有1个CNV,那么区间内的read深度将服从三峰的混合正态分布.将该方法应用于21 327个银屑病个体区域捕获测序的CNV检测中,结果表明,XHMM,ExomeDepth和本方法跟金标准重叠的窗口总数与金标准总窗口数的百分比(即重叠率)分别是7％、18％和62％.与XHMM和ExomeDepth相比,新方法在区间内CNV检测覆盖度可以分别提高55个百分点和44个百分点.本研究完善拷贝数变异检测方法,为疾病的诊断治疗提供一定的理论依据.     

Variation of CNV distribution in five different ethnic populations

Recent studies have revealed a new type of variation in the human genome encompassing relatively large genomic segments ( approximately 100 kb-2.5 Mb), commonly referred to as copy number variation (CNV). The full nature and extent of CNV and its frequency in different ethnic populations is still largely unknown. In this study we surveyed a set of 12 CNVs previously detected by array-CGH. More than 300 individuals from five different ethnic populations, including three distinct European, one Asian and one African population, were tested for the occurrence of CNV using multiplex ligation-dependent probe amplification (MLPA). Seven of these loci indeed showed CNV, i.e., showed copy numbers that deviated from the population median. More precise estimations of the actual genomic copy numbers for (part of) the NSF gene locus, revealed copy numbers ranging from two to at least seven. Additionally, significant inter-population differences in the distribution of these copy numbers were observed. These data suggest that insight into absolute DNA copy numbers for loci exhibiting CNV is required to determine their potential contribution to normal phenotypic variation and, in addition, disease susceptibility.     

Genetic Copy Number Variation and General Cognitive Ability

Differences in genomic structure between individuals are ubiquitous features of human genetic variation. Specific copy number variants (CNVs) have been associated with susceptibility to numerous complex psychiatric disorders, including attention-deficit-hyperactivity disorder, autism-spectrum disorders and schizophrenia. These disorders often display co-morbidity with low intelligence. Rare chromosomal deletions and duplications are associated with these disorders, so it has been suggested that these deletions or duplications may be associated with differences in intelligence. Here we investigate associations between large (≥500kb), rare (<1% population frequency) CNVs and both fluid and crystallized intelligence in community-dwelling older people. We observe no significant associations between intelligence and total CNV load. Examining individual CNV regions previously implicated in neuropsychological disorders, we find suggestive evidence that CNV regions around SHANK3 are associated with fluid intelligence as derived from a battery of cognitive tests. This is the first study to examine the effects of rare CNVs as called by multiple algorithms on cognition in a large non-clinical sample, and finds no effects of such variants on general cognitive ability.     

Maize Inbreds Exhibit High Levels of Copy Number Variation (CNV) and Presence/Absence Variation (PAV) in Genome Content

Following the domestication of maize over the past ∼10,000 years, breeders have exploited the extensive genetic diversity of this species to mold its phenotype to meet human needs. The extent of structural variation, including copy number variation (CNV) and presence/absence variation (PAV), which are thought to contribute to the extraordinary phenotypic diversity and plasticity of this important crop, have not been elucidated. Whole-genome, array-based, comparative genomic hybridization (CGH) revealed a level of structural diversity between the inbred lines B73 and Mo17 that is unprecedented among higher eukaryotes. A detailed analysis of altered segments of DNA conservatively estimates that there are several hundred CNV sequences among the two genotypes, as well as several thousand PAV sequences that are present in B73 but not Mo17. Haplotype-specific PAVs contain hundreds of single-copy, expressed genes that may contribute to heterosis and to the extraordinary phenotypic diversity of this important crop.     

Copy Number Variation in Thai Population

Copy number variation (CNV) is a major genetic polymorphism contributing to genetic diversity and human evolution. Clinical application of CNVs for diagnostic purposes largely depends on sufficient population CNV data for accurate interpretation. CNVs from general population in currently available databases help classify CNVs of uncertain clinical significance, and benign CNVs. Earlier studies of CNV distribution in several populations worldwide showed that a significant fraction of CNVs are population specific. In this study, we characterized and analyzed CNVs in 3,017 unrelated Thai individuals genotyped with the Illumina Human610, Illumina HumanOmniexpress, or Illumina HapMap550v3 platform. We employed hidden Markov model and circular binary segmentation methods to identify CNVs, extracted 23,458 CNVs consistently identified by both algorithms, and cataloged these high confident CNVs into our publicly available Thai CNV database. Analysis of CNVs in the Thai population identified a median of eight autosomal CNVs per individual. Most CNVs (96.73%) did not overlap with any known chromosomal imbalance syndromes documented in the DECIPHER database. When compared with CNVs in the 11 HapMap3 populations, CNVs found in the Thai population shared several characteristics with CNVs characterized in HapMap3. Common CNVs in Thais had similar frequencies to those in the HapMap3 populations, and all high frequency CNVs (>20%) found in Thai individuals could also be identified in HapMap3. The majorities of CNVs discovered in the Thai population, however, were of low frequency, or uniquely identified in Thais. When performing hierarchical clustering using CNV frequencies, the CNV data were clustered into Africans, Europeans, and Asians, in line with the clustering performed with single nucleotide polymorphism (SNP) data. As CNV data are specific to origin of population, our population-specific reference database will serve as a valuable addition to the existing resources for the investigation of clinical significance of CNVs in Thais and related ethnicities.     

Comparative Studies of Copy Number Variation Detection Methods for Next-Generation Sequencing Technologies

Copy number variation (CNV) has played an important role in studies of susceptibility or resistance to complex diseases. Traditional methods such as fluorescence in situ hybridization (FISH) and array comparative genomic hybridization (aCGH) suffer from low resolution of genomic regions. Following the emergence of next generation sequencing (NGS) technologies, CNV detection methods based on the short read data have recently been developed. However, due to the relatively young age of the procedures, their performance is not fully understood. To help investigators choose suitable methods to detect CNVs, comparative studies are needed. We compared six publicly available CNV detection methods: CNV-seq, FREEC, readDepth, CNVnator, SegSeq and event-wise testing (EWT). They are evaluated both on simulated and real data with different experiment settings. The receiver operating characteristic (ROC) curve is employed to demonstrate the detection performance in terms of sensitivity and specificity, box plot is employed to compare their performances in terms of breakpoint and copy number estimation, Venn diagram is employed to show the consistency among these methods, and F-score is employed to show the overlapping quality of detected CNVs. The computational demands are also studied. The results of our work provide a comprehensive evaluation on the performances of the selected CNV detection methods, which will help biological investigators choose the best possible method.     

Reliability and Stability of Contingent Negative Variation

Reliability parameters of a test indicate the stability (and quality) of the test itself. Reliability coefficients greater than 0.70 suggest an attribute as being sufficiently stable over time to be characterized as a trait. Reliability parameters of contingent negative variation (CNV) amplitudes in 27 healthy individuals were determined using a test-retest design. CNV was recorded at Cz, with an interstimulus interval of 3 s, on 2 separate occasions: initial session and 10 days later. Correlation coefficients between the 2 recording sessions were 0.675 for the total-CNV (tCNV), 0.855 for the early component (iCNV), 0.631 for the late component (lCNV), and 0.420 for the post-imperative negative variation (PINV). Statistical retest parameters for Spearman Brown were 0.806 for tCNV, 0.922 for iCNV, 0.774 for lCNV, and 0.655 for PINV. The iCNV, more than the other parameters, remained stable over the period of 10 days. It is suggested that the described standardized CNV recording procedure ensures reproducible and stable results in healthy subjects.     

A Definitive Haplotype Map as Determined by Genotyping Duplicated Haploid Genomes Finds a Predominant Haplotype Preference at Copy-Number Variation Events

The majority of complete hydatidiform moles (CHMs) harbor duplicated haploid genomes that originate from sperm. This makes CHMs more advantageous than conventional diploid cells for determining haplotypes of SNPs and copy-number variations (CNVs), because all of the genetic variants in a CHM genome are homozygous. Here we report SNP and CNV haplotype structures determined by analysis of 100 CHMs from Japanese subjects via high-density DNA arrays. The obtained haplotype map should be useful as a reference for the haplotype structure of Asian populations. We resolved common CNV regions (merged CNV segments across the examined samples) into CNV events (clusters of CNV segments) on the basis of mutual overlap and found that the haplotype backgrounds of different CNV events within the same CNV region were predominantly similar, perhaps because of inherent structural instability.     

Extensive Variation in Gene Copy Number at the Killer Immunoglobulin-Like Receptor Locus in Humans

Killer immunoglobulin-like receptors (KIRs) are involved in the regulation of natural killer cell cytotoxicity. Within the human genome seventeen KIR genes are present, which all contain a large number of allelic variants. The high level of homology among KIR genes has hampered KIR genotyping in larger cohorts, and determination of gene copy number variation (CNV) has been difficult. We have designed a multiplex ligation-dependent probe amplification (MLPA) technique for genotyping and CNV determination in one single assay and validated the results by next-generation sequencing and with a KIR gene-specific short tandem repeat assay. In this way, we demonstrate in a cohort of 120 individuals a high level of CNV for all KIR genes except for the framework genes KIR3DL3 and KIR3DL2. Application of our MLPA assay in segregation analyses of families from the Centre d’Etude du Polymorphisme Humaine, previously KIR-genotyped by classical techniques, confirmed an earlier reported duplication and resulted in the identification of a novel duplication event in one of these families. In summary, our KIR MLPA assay allows rapid and accurate KIR genotyping and CNV detection, thus rendering improved transplantation programs and oncology treatment feasible, and enables more detailed studies on the role of KIRs in human (auto)immunity and infectious disease.     

Cortical Habituation Deficit in Tinnitus Sufferers: Contingent Negative Variation as an Indicator of Duration of the Disease

Cortical attention and habituation parameters are altered in patients suffering from tinnitus. The aim of the study was to quantify cortical attention and habituation parameters in tinnitus patients by recording the contingent negative variation (CNV) response and to correlate amplitudes of different CNV parameters with duration of disease. Twenty patients suffering from tinnitus (median: 44?years) and twenty age- and sex-matched healthy controls (median: 41?years) were tested by a CNV paradigm. We recorded overall CNV, initial CNV, and terminal CNV and calculated habituation slopes. All CNV parameters were Spearman-correlated with individual duration of disease. Highly significant between groups differences emerged in total (tinnitus: -8.4?uV vs. controls: -3.8?uV), initial (-11.2 vs. -6.0?uV), and terminal CNV (-11.9 vs. -6.5?uV) demonstrating higher negative amplitudes in tinnitus patients. Habituation differed in total and terminal CNV, indicating missing habituation in tinnitus patients. Overall CNV (??=?-.365) and initial CNV (??=?-.529) showed a medium Spearman correlation with duration of disease. We conclude that the correlation between duration of tinnitus and the initial CNV amplitudes indicates an altered state of cortical excitability that can also be observed in more negative CNV-amplitudes in tinnitus patients. We assume that this state indicates a chronicity process in tinnitus disease.     

Genetic Variation of the Fc Gamma Receptor 3B Gene and Association with Rheumatoid Arthritis

Background

Fc gamma receptors (FcγRs) play a crucial role in immunity by linking IgG antibody-mediated responses with cellular effector and regulatory functions. Genetic variants in these receptors have been previously identified as risk factors for several chronic inflammatory conditions. The present study aimed to investigate the presence of copy number variations (CNVs) in the FCGR3B gene and its potential association with the autoimmune disease rheumatoid arthritis (RA).

Methodology/Principal Findings

CNV of the FCGR3B gene was studied using Multiplex Ligation Dependent Probe Amplification (MLPA) in 518 Dutch RA patients and 304 healthy controls. Surprisingly, three independent MLPA probes targeting the FCGR3B promoter measured different CNV frequencies, with probe#1 and #2 measuring 0 to 5 gene copies and probe#3 showing little evidence of CNV. Quantitative-PCR correlated with the copy number results from MLPA probe#2, which detected low copy number (1 copy) in 6.7% and high copy number (≥3 copies) in 9.4% of the control population. No significant difference was observed between RA patients and the healthy controls, neither in the low copy nor the high copy number groups (p-values = 0.36 and 0.71, respectively). Sequencing of the FCGR3B promoter region revealed an insertion/deletion (indel) that explained the disparate CNV results of MLPA probe#1. Finally, a non-significant trend was found between the novel -256A>TG indel and RA (40.7% in healthy controls versus 35.9% in RA patients; P = 0.08).

Conclusions/Significance

The current study highlights the complexity and poor characterization of the FCGR3B gene sequence, indicating that the design and interpretation of genotyping assays based on specific probe sequences must be performed with caution. Nonetheless, we confirmed the presence of CNV and identified novel polymorphisms in the FCGR3B gene in the Dutch population. Although no association was found between RA and FCGR3B CNV, the possible protective effect of the -256A>TG indel polymorphism must be addressed in larger studies.