Copy number variation (CNVs) is a relatively new field in genomics and it is defined as a phenomenon in which sections of the genome are repeated and the number of repeats in the genome varies between individuals in the human population. Copy number variation is a type of structural variation, specifically, it is a type of duplication or deletion event that affects a considerable number of base pairs. However, note that although modern genomics research is mostly focused on human genomes, copy number variations also occur in a variety of other organisms including E. coli. Recent research indicates that approximately two thirds of the entire human genome is composed of repeats and 4.8-9.5% of the human genome can be classified as copy number variations. In mammals, copy number variations play an important role in generating necessary variation in the population as well as disease phenotype.
Copy number variations can be generally categorized into two main groups: short repeats and long repeats. However, there are no clear boundaries between the two groups and the classification depends on the nature of the loci of interest. Short repeats include mainly bi-nucleotide repeats (two repeating nucleotides e.g. A-B-A-B-A-B…) and tri-nucleotide repeats. Long repeats include repeats of entire genes. This classification based on size of the repeat is the most obvious type of classification as size is an important factor in examining the types of mechanisms that most likely gave rise to the repeats, hence the likely effects of these repeats on phenotype.
Types of Copy Number Variation and Chromosomal Rearrangements
One of the most well known examples of a short copy number variation is the tri-nucleotide repeat of the CAG base pairs in the Huntingtin gene, the gene that is responsible for the neurological disorder Huntington’s disease. For this particular case, once the CAG tri-nucleotide repeats more than 44-46 times, Huntington’s disease will likely develop in the individual and it will likely be inherited by his or her offspring. Interestingly, the number of repeats of the CAG tri-nucleotide is correlated with the age of onset of Huntington’s disease. These types of short repeats are often thought to be due to errors in polymerase activity during replication including polymerase slippage, template switching, and fork switching which will be discussed in detail later. The short repeat size of these copy number variations lends itself to errors in the polymerase as these repeated regions are prone to misrecognition by the polymerase and replicated regions may be replicated again, leading to extra copies of the repeat. In addition, if these tri-nucleotide repeats are in the same reading frame in the coding portion of a gene, it may lead to a long chain of the same amino acid, possibly creating protein aggregates in the cell, and if these short repeats fall into the non-coding portion of the gene, it may affect gene expression and regulation. On the other hand, a variable number of repeats of entire genes is less commonly identified in the genome. One example of a whole gene repeat is the alpha-amylase 1 gene (AMY1) that encodes alpha-amylase which has a significant copy number variation between different populations with different diets. Although the specific mechanism that allows the AMY1 gene to increase or decrease its copy number is still a topic of debate, some hypotheses suggests that the non-homologous end joining or the microhomology-mediated end joining is likely responsible for these whole gene repeats. Repeats of entire genes has immediate effects on expression of that particular gene, and the fact that the copy number variation of the AMY1 gene has been related to diet is a remarkable example of recent human evolutionary adaptation. Although these are the general groups that copy number variations are grouped into, the exact number of base pairs copy number variations affect depends on the specific loci of interest. Currently, using data from all reported copy number variations, the mean size of copy number variant is around 118kb, and the median is around 18kb.
In terms of the structural architecture of copy number variations, research has suggested and defined hotspot regions in the genome where copy number variations are four times more enriched. These hotspot regions were defined to be regions containing long repeats that are 90-100% similar known as segmental duplications either tandem or interspersed and most importantly, these hotspot regions have an increased rate of chromosomal rearrangement. It was thought that these large-scale chromosomal rearrangements give rise to normal variation and genetic diseases, including copy number variations. Moreover, these copy number variation hotspots are consistent throughout many populations from different continents, implying that these hotspots were either independently acquired by all the populations and passed on through generations, or they were acquired in early human evolution before the populations split, the latter seems more likely. Lastly, spatial biases of the location at which copy number variations are most densely distributed does not seem to occur in the genome. Although it was originally detected by fluorescent in situ hybridization and microsatellite analysis that copy number repeats are localized to regions that are highly repetitive such as telomeres, centromeres, and heterochromatin, recent genome-wide studies have concluded otherwise. Namely, the subtelomeric regions and pericentromeric regions are where most chromosomal rearrangement hotspots are found, and there is no considerable increase in copy number variations in that region. Furthermore, these regions of chromosomal rearrangement hotspots do not have decreased gene numbers, again, implying that there is minimal spatial bias of the genomic location of copy number variations.
Detection and Identification
Copy number variation was initially thought to occupy an extremely small and negligible portion of the genome through cytogenetic observations. Copy number variations were generally associated only with small tandem repeats or specific genetic disorders, therefore, copy number variations were initially only examined in terms of a specific loci. However, breakthroughs in the past decade or so led to an increasing number of highly accurate ways of identifying and studying copy number variations, one of which is the genome-wide association study that allow copy number variations in general to be located and identified in the genome. Copy number variations were originally studied by cytogenetic techniques, which are techniques that allow one to observe the physical structure of the chromosome. One of these techniques is fluorescent in situ hybridization (FISH) which involves inserting fluorescent probes that require a high degree of complementarity in the genome for binding. Comparative genomic hybridization was also commonly used to detect copy number variations by fluorophore visualization and then comparing the length of the chromosomes. One major drawback of these early techniques is that the genomic resolution is relatively low and only large repeats such as whole gene repeats can be detected.
Recent advances in biotechnology gave rise to many important techniques that are of extremely high genomic resolution and as a result, an increasing number of copy number variations in the genome have been reported. One of these advances involve using bacterial artificial chromosome (BAC) array with around 1 megabase of intervals throughout the entire gene, BACs can also detect copy number variations in rearrangement hotspots allowing for the detection of 119 novel copy number variations. Throughout the past decade or so, high throughput genomic sequencing has revolutionized the field of human genomics and in silico studies have been performed to detect copy number variations in the genome. Reference sequences have been compared to other sequences of interest using fosmids by strictly controlling the fosmid clones to be 40kb. Sequencing end reads would provide adequate information to align the reference sequence to the sequence of interest, and any misalignments are easily noticeable thus concluded to be copy number variations within that region of the clone. This type of detection technique offers a high genomic resolution and precise location of the repeat in the genome, and it can also detect other types of structural variation such as inversions. In addition, another way of detecting copy number variation that can ensure high genomic resolution is using single nucleotide polymorphisms (SNP). Since the International HapMap project had begun, common SNPs that occur between four different populations from different continents have been sequenced and located. Due to the abundance of the human SNP data, the direction of detecting copy number variation has changed to utilize these SNPs. Relying on the fact that human recombination is relatively rare and that many recombination events occur in specific regions of the genome known as recombination hotspots, linkage disequilibrium can be used to identify copy number variations. Efforts have been made in associating copy number variations with specific haplotype SNPs by analyzing the linkage disequilibrium, using these associations, one is able to recognize copy number variations in the genome using SNPs as markers. One drawback of this method is that because SNPs in the International HapMap are not optimized for detecting copy number variations, therefore the data is biased towards large copy number variations. Lastly, next-generation sequencing has also been used recently to detect copy number variations with high genomic resolutions. Using whole-genome shot-gun sequencing data, assays have been developed to accurately detect and identify regions of duplications.
Accurately detecting, identifying, and categorizing copy number variations is extremely important because of the complications they bring to DNA sequencing. Traditionally, DNA sequencing relied heavily on sequencing short reads from a large genome and using any overlapping regions of the reads to combine the short reads to form longer reads. This will eventually be mapped together to give the sequence of the entire genome. However, the issues related to copy number variations arise in linking the overlapping regions together. By definition, copy number variation is a region of the genome duplicated a variable number of times in the population and due to the large variation between the number of times portions of the genome are duplicated, when mapping overlapping sequences, it becomes unclear whether or not a region is an overlap or a duplicated region.
The more fascinating question to ask is perhaps not how copy number variations are detected in the genome, but rather how did they arise in the first place. An important idea to bear in mind while considering these mechanisms is that although many theories have been put forward, most of these theories are speculations and conjecture. More importantly, there is no conclusive evidence that correlates a specific copy number variation to a specific mechanism. There are two main types of molecular mechanism for the formation of copy number variations: homologous based and non-homologous based.
One of the most well recognized theory that lead to copy number variations as well as deletions and inversion is non-allelic homologous recombinations. This mechanism suggests that during meiotic recombination, homologous chromosomes pair up and form two ended double-stranded breaks leading to Holliday junctions as usual. However, during the formation of Holliday junctions, the double-stranded breaks are misaligned and the crossover lands in non-allelic positions on the same chromosome. When the Holliday junction is resolved, the unequal crossing over event allows transfer of genetic material between the two homologous chromosomes, and as a result, a portion of the DNA on both the homologies is repeated. Since the repeated regions are no longer segregating independently, the duplicated region of the chromosome is inherited. Another type of homologous recombination based mechanism that can lead to copy number variation is known as break induced replication. When a double stranded break occurs in the genome unexpectedly the cell activates pathways that mediate the repair of the break. Errors in repairing the break, similar to non-allelic homologous recombination, can lead to an increase in copy number of a particular region of the genome. During the repair of a double stranded break, the broken end can invade its homologous chromosome instead of rejoin the original strand. As in the non-allelic homologous recombination mechanism, an extra copy of a particular region is transferred to another chromosome, leading to a duplication event. Furthermore, cohesin proteins are found to aid in the repair system of double stranded breaks through clamping the two ends in close proximity which prevents interchromosomal invasion of the ends. If for any reason, such as activation of ribosomal RNA, cohesin activity is affected then there may be local increase in double stranded break repair errors.
The other class of possible mechanisms that are hypothesized to lead to copy number variations is non-homologous based. To distinguish between this and homologous based mechanisms, one must understand the concept of homology. Homologous pairing of chromosomes involved using DNA strands that are highly similar to each other (~97%) and these strands must be longer than a certain length to avoid short but highly similar pairings. Non-homologous pairings, on the other hand, rely on only few base pairs of similarity between two strands, therefore it is possible for genetic materials to be exchanged or duplicated in the process of non-homologous based double stranded repairs.
One type of non-homologous based mechanism is the non-homologous end joining or micro-homology end joining mechanism. These mechanisms are also involved in repairing double stranded breaks but require no homology or limited micro-homology. When these strands are repaired, oftentimes there are small deletions or insertions added into the repaired strand. It is possible that retrotransposons are inserted into the genome through this repair system. If retrotransposons are inserted into a non-allelic position on the chromosome, meiotic recombination can drive the insertion to be recombined into the same strand as an already existing copy of the same region. Another mechanism is the break-fusion-bridge cycle which involves sister chromatids that have both lost its telomeric region due to double stranded breaks. It is proposed that these sister chromatids will fuse together to form one dicentric chromosome, and then segregate into two different nuclei. Because pulling the dicentric chromosome apart causes a double stranded break, the end regions can fuse to other double stranded breaks and repeat the cycle. The fusion of two sister chromatids can cause inverted duplication and when these events are repeated throughout the cycle, the inverted region will be repeated leading to an increase in copy number. The last mechanism that can lead to copy number variations is polymerase slippage, which is also known as template switching. During normal DNA replication, the polymerase on the lagging strand is required to unclamp and re-clamp the replication region continuously. When small scale repeats in the DNA sequence exist already, the polymerase can be ‘confused’ when it re-clamps to continue replication and instead of clamping to the correct base pairs, it may shift a few base pairs and replicate a portion of the repeated region again. Note that although this has been experimentally observed and is a widely accepted mechanism, the molecular interactions that led to this error remains unknown. In addition, because this type of mechanism requires the polymerase to jump around the DNA strand and it is unlikely that the polymerase can re-clamp at another locus some kilobases apart, therefore this is more applicable to short repeats such as bi-nucleotide or tri-nucleotide repeats.
Amylase is an enzyme in saliva that is responsible for the breakdown of starch into monosaccharides, and one type of amylase is encoded by the alpha-amylase gene (AMY1). The AMY1 locus, as well as the amylase enzyme, is one of the most extensively studied and sequenced gene in the human genome. Its homologs are also found in other primates and therefore it is likely that the primate AMY1 gene is ancestral to the human AMY1 gene and was adapted early in primate evolution. AMY1 is a particularly interesting gene to examine as it is one of the most well studied genes which has wide range of variable numbers of copies throughout different human populations. The AMY1 gene is also one of the few genes that had been studied that displayed convincing evidence which correlates its protein function to its copy number. Copy number is known to alter transcription as well as translation levels of a particular gene, however research has shown that the relationship between protein levels and copy number is variable. In the AMY1 genes of European Americans it is found that the concentration of salivary amylase is closely correlated to the copy number of the AMY1 gene. As a result, it was hypothesized that the copy number of the AMY1 gene is closely correlated with its protein function, which is to digest starch.
The AMY1 gene copy number has been found to be correlated to different levels of starch in diets of different populations. 8 Populations from different continents were categorized into high starch diets and low starch diets and their AMY1 gene copy number was visualized using high resolution FISH and qPCR. It was found that the high starch diet populations which consists of the Japanese, Hadza, and European American populations had a significantly higher (2 times higher) average AMY1 copy number than the low starch diet populations including Biaka, Mbuti, Datog, Yakut populations. It was hypothesized that the levels of starch in one’s regular diet, the substrate for AMY1, can directly affect the copy number of the AMY1 gene. Since it was concluded that the copy number of AMY1 is directly correlated with salivary amylase, the more starch present in the population’s daily diet, the more evolutionarily favorable it is to have multiple copies of the AMY1 gene. The AMY1 gene was the first gene to provide strong evidence for evolution on a molecular genetics level. Moreover, using comparative genomic hybridization, copy number variations of the entire genomes of the Japanese population was compared to that of the Yakut population. It was found that the copy number variation of the AMY1 gene was significantly different from the copy number variation in other genes or regions of the genome, suggesting that the AMY1 gene was under a strong selective pressure that had little or no influence on the other copy number variations. Finally, the variability of length of 783 microsatellites between the two populations were compared to copy number variability of the AMY1 gene. It was found that the AMY1 gene copy number range was larger than that of over 97% of the microsatellites examined. This implies that natural selection played a considerable role in shaping the average number of AMY1 genes in these two populations. However, as only 6 populations were studied, it is important to consider the possibility that there may be other factors in their diet or culture that influenced the AMY1 copy number other than starch.
Although it is unclear when the AMY1 gene copy number began to increase, it is known and confirmed that the AMY1 gene existed in early primates. Chimpanzees, the closest evolutionary relatives to humans, were found to have 2 diploid copies of the AMY1 gene that is identical in length to the human AMY1 gene, which is significantly less than that of humans. On the other hand, bonobos, also a close relative of modern humans, was found to have more than 2 diploid copies of the AMY1 gene. Nonetheless, the bonobo AMY1 genes were sequenced and analyzed, and it was found that the coding sequences of the AMY1 genes were disrupting, which may lead to the production of dysfunctional salivary amylase. It can be inferred from the results that the increase in bonobo AMY1 copy number is likely not correlated to the amount of starch in their diet. It was further hypothesized that the increase in copy number began recently during early hominin evolution as none of the great apes had more than two copies of the AMY1 gene that produced functional protein. In addition, it was speculated that the increase in the AMY1 copy number began around 20,000 years ago when humans shifted from a hunter-gatherer lifestyle to agricultural societies, which was also when humans relied heavily on root vegetables high in starch. This hypothesis, although logical, lacks experimental evidence due to the difficulties in gathering information on the shift of human diets, especially on root vegetables that are high in starch as they cannot be directly observed or tested. Recent breakthroughs in DNA sequencing has allowed researchers to sequence older DNA such as that of Neanderthals to a certain degree of accuracy. Perhaps sequencing Neanderthal DNA can provide a time marker as to when the AMY1 gene copy number increased and offer insight into human diet and gene evolution.
Currently it is unknown which mechanism gave rise to the initial duplication of the amylase gene, and it can imply that the insertion of the retroviral sequences was due to non-homologous end joining, which caused the duplication of the AMY1 gene. However, there is currently no evidence to support this theory and therefore this hypothesis remains conjecture. The recent origin of the multi-copy AMY1 gene implies that depending on the environment, the AMY1 gene copy number can increase and decrease very rapidly relative to genes that do not interact as directly with the environment. The AMY1 gene is an excellent example of how gene dosage affects the survival of an organism in a given environment. The multiple copies of the AMY1 gene gives those who rely more heavily on high starch diets an evolutionary advantage, therefore the high gene copy number persists in the population.
Copy Number Variations, Gene Families, and Natural Selection
Recently, there had been discussion connecting copy number variations to gene families. Gene families are defined as a set of related genes that serve similar functions but have minor temporal or spatial differences and these genes likely derived from one ancestral gene. The main reason copy number variations are connected to gene families is that there is a possibility that genes in a family may have derived from one ancestral gene which got duplicated into different copies. Mutations accumulate through time in the genes and with natural selection acting on the genes, some mutations lead to environmental advantages allowing those genes to be inherited and eventually clear gene families are separated out. An example of a gene family that may have been created due to copy number variations is the globin gene family. The globin gene family is an elaborate network of genes consisting of alpha and beta globin genes including genes that are expressed in both embryos and adults as well as pseudogenes. These globin genes in the globin family are all well conserved and only differ by a small portion of the gene, indicating that they were derived from a common ancestral gene, perhaps due to duplication of the initial globin gene.
Research have shown that copy number variations are significantly more common in genes that encode proteins that directly interact with the environment than proteins that are involved in basic cellular activities. It was suggested that the gene dosage effect accompanying copy number variation may lead to detrimental effects if essential cellular functions are disrupted, therefore proteins involved in cellular pathways are subjected to strong purifying selection. In addition, proteins function together and interact with proteins of other pathways, therefore it is important to view the effects of natural selection on bio-molecular pathways rather than on individual proteins. With that being said, it was found that proteins in the periphery of the pathway are enriched in copy number variations whereas proteins in the center of the pathways are depleted in copy number variations. It was explained that proteins in the periphery of the pathway interact with fewer proteins and so a change in protein dosage affected by a change in copy number may have a smaller effect on the overall outcome of the cellular pathway.
In the past few years, researchers seem to have shifted their focus from detecting, locating, and sequencing copy number variations to in depth analyses of the role of these copy number variations in the human genome and in nature in general. Evidence is needed to further validate the relationship between copy number variations and gene families as well as the role that natural selection plays in shaping these relationships and changes. Furthermore, researchers are also aiming to elucidating the molecular mechanisms involved in copy number variations as it may reveal essential information regarding structural variations in general. Taking a step back, the area of structural variation in the human genome seems to be a rapidly growing research topic. Not only can these research data provide additional evidence for evolution and natural selection, it can also be used to develop treatments for a wide range of genetic diseases.
- Comparative genomics
- Copy number analysis
- Human genome
- Molecular evolution
- Segmental duplication
- Tandem exon duplication
- Virtual Karyotype
- Mccarroll, S. A.; Altshuler, D. M. (2007). "Copy-number variation and association studies of human diseases". Nature Genetics 39: 37–42. doi:10.1038/ng2080. PMID 17597780.
- Sharp, A. J.; Locke, D. P.; Mcgrath, S. D.; Cheng, Z; Bailey, J. A.; Vallente, R. U.; Pertz, L. M.; Clark, R. A.; Schwartz, S.; Segraves, R. (2005). "Segmental Duplications and Copy-Number Variation in the Human Genome". The American Journal of Human Genetics 77 (1): 78–88. doi:10.1086/431652. PMC 1226196. PMID 15918152.
- Taniguchi, Y.; Choi, P. J.; Li, J. W.; Chen, H.; Babu, M.; Hearn, J.; Emili, A.; Xie, X. S. (2010). "Quantifying E. coli Proteome and Transcriptome with Single-Molecule Sensitivity in Single Cells". Science 329 (5991): 533–538. doi:10.1126/science.1188308.
- Koning, A. P. J. D.; Gu, W.; Castoe, T. A.; Bazter, M. A.; Pollock, D. D. (2011). "Repetitive Elements May Comprise Over Two-Thirds of the Human Genome". PLOS Genetics 7 (12): e1002384. doi:10.1371/journal.pgen.1002384. PMC 3228813. PMID 22144907.
- Zarrei, M.; Macdonald, J. R.; Merico, D.; Scherer, S. W. (2015). "A copy number variation map of the human genome". Nature Reviews Genetics 16: 172–183. doi:10.1038/nrg3871. PMID 25645873.
- Hastings, P. J.; Lupski, J. R.; Roseberg, S. M.; Ira, G. (2009). "Mechanisms of change in gene copy number". Nature Reviews Genetics 10: 551–564. doi:10.1038/nrg2593. PMC 2864001. PMID 19597530.
- Macdonald, M.; Ambrose, C. M.; Duyao, M. P.; Myers, R. H.; Lin, C.; Srinidhi, L.; Barnes, G.; Taylor, S. A.; James, M.; Groot, N. (1993). "A novel gene containing a trinucleotide repeat that is expanded and unstable on Huntington's disease chromosomes". Cell 72 (6): 971–983. doi:10.1016/0092-8674(93)90585-e. PMID 8458085.
- Myer, R. H. (2004). "Huntington’s Disease Genetics". NeuroRx 1 (2): 255–262. doi:10.1602/neurorx.1.2.255. PMC 534940. PMID 15717026.
- Albertini, A. M.; Hofer, M.; Calos, M. P.; Miller, J. H. (1982). "On the formation of spontaneous deletions: The importance of short sequence homologies in the generation of large deletions". Cell 29 (2): 319–328. doi:10.1016/0092-8674(82)90148-9.
- Perry, G. H.; Dominy, N. J.; Claw, K. G.; Lee, A. S.; Fiegler, H.; Redon, R.; Werner, J.; Villanae, F. A.; Mountain, J. L.; Misra, R (2007). "Diet and the evolution of human amylase gene copy number variation". Nature Genetics 39 (10): 1256–1260. doi:10.1038/ng2123. PMC 2377015. PMID 17828263.
- Freeman, J. L. (2006). "Copy number variation: New insights in genome diversity". Genome Research 16: 949–961. doi:10.1101/gr.3677206.
- Bailey, J. A.; Gu, Z.; Clark, R. A.; Reinert, K.; Samonte, R. V.; Schwartz, S.; Adams, M. D.; Myers, E. W.; Li, P. W.; Eichler, E. E. (2002). "Recent Segmental Duplications in the Human Genome". Science 297 (5583): 1003–1007. doi:10.1126/science.1072047.
- Jacobs, P. A.; Browne, C.; Gregson, N.; Joyce, C.; White, H. (1992). "Estimates of the frequency of chromosome abnormalities detectable in unselected newborns using moderate levels of banding". Journal of Medical Genetics 29 (2): 103–108. doi:10.1136/jmg.29.2.103. PMC 1015848. PMID 1613759.
- Inuoe, K; Lupski, J. R. (2002). "Molecular mechanism for genetic disorders". Annual Review of Genomics and Human Genetics 3: 199–242. doi:10.1146/annurev.genom.3.032802.120023.
- Iafrate, A. J.; Feuk, L.; Rivera, M. N.; Listewnik, M. L.; Donahoe, P. K.; Qi, Y.; Scherer, S. W.; Lee, C. (2004). "Detection of large-scale variation in the human genome". Nature Genetics 36 (9): 949–951. doi:10.1038/ng1416. PMID 15286789.
- Tuzun, E.; Sharp, A. J.; Bailey, J. A.; Kaul, R.; Morrison, V. A.; Pertz, L. M.; Haugen, E.; Hayden, H.; Albertson, D.; Pinkel, D. (2005). "Fine-scale structural variation of the human genome". Nature Genetic 37 (7): 727–732. doi:10.1038/ng1562. PMID 15895083.
- Conrad, D. F.; Antonarakis, S. E. (2007). "Gene Duplication: A Drive for Phenotypic Diversity and Cause of Human Disease. Annual Review of Genomics and Human Genetics". Annual Review of Genomics and Human Genetics 8: 17–35. doi:10.1146/annurev.genom.8.021307.110233. PMID 17386002.
- Xie, C; Tammi, M. T. (2009). "CNV-seq, a new method to detect copy number variation using high-throughput sequencing". BMC Bioinformatics 10 (8): 80. doi:10.1186/1471-2105-10-80.
- Pagues, F; Haber, J. E. (1999). "Multiple pathways of recombination induced by double-strand breaks in Saccharomyces Cerevisiae". Microbiology Review 63 (2): 349–404. PMC 98970. PMID 10357855.
- Bauters, M.; Esch, H. V.; Friez, M. J.; Boespflug-Tanguy, O.; Zenker, M.; Vianna-Morgante, A. M.; Rosenberg, C.; Ignatius, J.; Raynaud, M.; Hollanders, K. (2008). "Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair". Genome Research 18 (6): 847–858. doi:10.1101/gr.075903.107.
- Kobayashi, T.; Ganley, A. R. (2005). "Recombination Regulation by Transcription-Induced Cohesin Dissociation in rDNA Repeats". Science 309 (5740): 1581–1584. doi:10.1126/science.1116102. PMID 16141077.
- Leiber, M.R. (2007). "The Mechanism of Human Nonhomologous DNA End Joining". Journal of Biological Chemistry 283 (1): 1–5. doi:10.1074/jbc.R700039200. PMID 17999957.
- McClintock, B. (1951). "Chromosome Organization And Genic Expression". Cold Spring Harbor Symposia on Quantitative Biology 16: 13–47. doi:10.1101/sqb.1951.016.01.004. PMID 14942727.
- Smith, C. E.; Llorente, B.; Symington, L. S. (2007). "Template switching during break-induced replication". Nature 447: 102–105. doi:10.1038/nature05723.
- Bi, S.; Liu, L. F. (1994). "recA-independent and recA-dependent Intramolecular Plasmid Recombination". Journal of Molecular Biology 235 (2): 414–423. doi:10.1006/jmbi.1994.1002. PMID 8289271.
- Korbel, J. O.; Kim, P. M.; Chen, X.; Urban, A. E.; Weissman, S.; Snyder, M.; Gerstein, M. B. (2008). "The current excitement about copy-number variation: how it relates to gene duplications and protein families". Current Opinion in Structural Biology 18 (3): 366–374. doi:10.1016/j.sbi.2008.02.005.
- Samuelson, L. C.; Wiebauer, L.; Snow, C. M.; Meisler, M. H. (1990). "Retroviral and pseudogene insertion sites reveal the lineage of human salivary and pancreatic amylase genes from a single gene during primate evolution". Molecular and Cell Biology 10 (6): 2513–2520. doi:10.1128/mcb.10.6.2513. PMC 360608. PMID 1692956.
- Goodman, M.; Koop, B. F.; Czelnusniak, J.; Weiss, M. L. (1984). "The n-globin gene: Its long evolutionary history in the B-globin gene family of mammals". Journal of Molecular Biology 180 (4): 803–823. PMID 6527390.
- Redon, R.; Ishikawa, S.; Fitch, K. R.; Feuk, L.; Perry, G. H.; Andrew, T. D.; Fiegler, H.; Shapero, M. H.; Carson, A. R.; Chen, W. (2006). "Global variation in copy number in the human genome". Nature 444: 444–454. doi:10.1038/nature05329. PMC 2669898. PMID 17122850.
- Kim, P. M.; Korbel, J. O.; Gerstein, M. B. (2007). "Positive selection at the protein network periphery: Evaluation in terms of structural constraints and cellular context". Proceedings of the National Academy of Sciences 104 (51): 20274–20279. doi:10.1073/pnas.0710183104.
- Genome-wide analysis of DNA copy-number changes using cDNA microarrays "Gene amplifications and deletions frequently contribute to tumorigenesis. Characterization of these DNA copy-number changes is important for both the basic understanding of cancer and its diagnosis".
- "Huge genetic variation in healthy people". New Scientist. 7 August 2004.
- "As normal as normal can be". Nature Genetics. 1 September 2004.
- "Human Genome: Patchwork people". Nature. 20 October 2005.
- "Gene duplications may define who you are". New Scientist. 22 November 2006.
- "DNA varies more widely from person to person, Genetic maps reveal". National Geographic. 22 November 2006.
- "Finding the right lenses" (PDF). Nature Genetics. 1 July 2007.
- "Nucleotide-resolution analysis of structural variants using BreakSeq and a breakpoint library" (PDF). Nature Biotechnoloty. 1 January 2010.
- "New Research Sheds Light on Autism's Genetic Causes". Singularity Hub. 15 June 2010.
- Copy Number Variation Project, Sanger Institute
- Integrative annotation platform for copy number variations in humans
- A bibliography on copy number variation
- Database of Genomic Variants, a database of structural variants in the human genome
- Copy Number Variation Detection via High-Density SNP Genotyping
- Oxford Gene Technology
- BioDiscovery Nexus Copy Number
- High-resolution mapping of copy number variations in 2,026 healthy individuals
- The 1000 Genomes Project
- cn.FARMS: a latent variable model to detect copy number variations in microarray data with a low false discovery rate, an R package —software