Evolutionary dynamics of segmental duplications from human Y-chromosomal euchromatin/heterochromatin transition regions

Reconstruction of Segmental Duplication Rates and Associated Genomic Features by Network Analysis

Segmental duplications are long genomic duplications that are fixed in a genome. Segmental duplications play an important evolutionary role because entire genes can be duplicated along with regulatory sequences. The ancestral segmental duplications of the human lineage gave rise to genes that are involved in the development of the human brain and provided sites for further genomic rearrangements. While some duplicated loci have been extensively studied, the universal principles and biological factors underlying the spread of segmental duplications remain unclear. Here, we represent segmental duplications in a network, with edges corresponding to duplication events and nodes corresponding to affected genomic sites. This representation allowed us to estimate how many duplications had occurred at each locus, and thereby enabling the prediction of genomic features associated with increased duplication rates. Our comprehensive study of genomic features associated with duplications and those associated with increased duplication rates allowed us to identify several biological factors affecting a segmental duplication process. In our study, we describe genomic features associated with increased duplication rates, three signatures of the duplication process and associations of segmental duplications with different classes of high-copy repeats. Furthermore our method is readily implemented and can easily be applied to segmental duplications of other genomes to build a network of segmental duplications or to predict real duplication events.

Genetic variation for three Y-STR loci: DYS390, DYS518, and DYS643

Y chromosome short tandem repeats (Y-STRs) are commonly used to analyze male-specific DNA. Although biallelic patterns due to duplication events have been detected at some loci, Y-STRs generally appear as a single peak except for DYS385 because the Y chromosome is haploid. STR loci in regions of segmental duplication by homologous recombination on the Y chromosome exhibit double allelic peaks, rather than single peaks. In this study, we report a bi- and triallelic pattern observed simultaneously in DYS390, DYS518, and DYS643. A bi- and triallelic pattern has not previously been observed simultaneously for these three loci. We also identified the copy number variation in the region including these loci by the microarray-based analysis. Given the peak balance pattern, the copy number variation, and the close position of these three loci on the Y chromosome, we consider that this phenomenon is caused by a segmental duplication in the euchromatin region. By ruling out mixed samples, a common interpretation of multiple peaks, these results have practical implications for the interpretation of Y-STR results in forensics analyses.

Common Bean Subtelomeres Are Hot Spots of Recombination and Favor Resistance Gene Evolution

Subtelomeres of most eukaryotes contain fast-evolving genes usually involved in adaptive processes. In common bean (Phaseolus vulgaris), the Co-2 anthracnose resistance (R) locus corresponds to a cluster of nucleotide-binding-site leucine-rich-repeat (NL) encoding sequences, the prevalent class of plant R genes. To study the recent evolution of this R gene cluster, we used a combination of sequence, genetic and cytogenetic comparative analyses between common bean genotypes from two distinct gene pools (Andean and Mesoamerican) that diverged 0.165 million years ago. Co-2 is a large subtelomeric cluster on chromosome 11 comprising from 32 (Mesoamerican) to 52 (Andean) NL sequences embedded within khipu satellite repeats. Since the recent split between Andean and Mesoamerican gene pools, the Co-2 cluster has experienced numerous gene-pool specific NL losses, leading to distinct NL repertoires. The high proportion of solo-LTR retrotransposons indicates that the Co-2 cluster is located in a hot spot of unequal intra-strand homologous recombination. Furthermore, we observe large segmental duplications involving both Non-Homologous End Joining and Homologous Recombination double-strand break repair pathways. Finally, the identification of a Mesoamerican-specific subtelomeric sequence reveals frequent interchromosomal recombinations between common bean subtelomeres. Altogether, our results highlight that common bean subtelomeres are hot spots of recombination and favor the rapid evolution of R genes. We propose that chromosome ends could act as R gene incubators in many plant genomes.

The Y chromosomes of the great apes

The great apes (orangutans, gorillas, chimpanzees, bonobos and humans) descended from a common ancestor around 13 million years ago, and since then their sex chromosomes have followed very different evolutionary paths. While great-ape X chromosomes are highly conserved, their Y chromosomes, reflecting the general lability and degeneration of this male-specific part of the genome since its early mammalian origin, have evolved rapidly both between and within species. Understanding great-ape Y chromosome structure, gene content and diversity would provide a valuable evolutionary context for the human Y, and would also illuminate sex-biased behaviours, and the effects of the evolutionary pressures exerted by different mating strategies on this male-specific part of the genome. High-quality Y-chromosome sequences are available for human and chimpanzee (and low-quality for gorilla). The chromosomes differ in size, sequence organisation and content, and while retaining a relatively stable set of ancestral single-copy genes, show considerable variation in content and copy number of ampliconic multi-copy genes. Studies of Y-chromosome diversity in other great apes are relatively undeveloped compared to those in humans, but have nevertheless provided insights into speciation, dispersal, and mating patterns. Future studies, including data from larger sample sizes of wild-born and geographically well-defined individuals, and full Y-chromosome sequences from bonobos, gorillas and orangutans, promise to further our understanding of population histories, male-biased behaviours, mutation processes, and the functions of Y-chromosomal genes.

Multicopy gene family evolution on primate Y chromosomes

Background

The primate Y chromosome is distinguished by a lack of inter-chromosomal recombination along most of its length, extensive gene loss, and a prevalence of repetitive elements. A group of genes on the male-specific portion of the Y chromosome known as the “ampliconic genes” are present in multiple copies that are sometimes part of palindromes, and that undergo a form of intra-chromosomal recombination called gene conversion, wherein the nucleotides of one copy are homogenized by those of another. With the aim of further understanding gene family evolution of these genes, we collected nucleotide sequence and gene copy number information for several species of papionin monkey. We then tested for evidence of gene conversion, and developed a novel statistical framework to evaluate alternative models of gene family evolution using our data combined with other information from a human, a chimpanzee, and a rhesus macaque.

Results

Our results (i) recovered evidence for several novel examples of gene conversion in papionin monkeys and indicate that (ii) ampliconic gene families evolve faster than autosomal gene families and than single-copy genes on the Y chromosome and that (iii) Y-linked singleton and autosomal gene families evolved faster in humans and chimps than they do in the other Old World Monkey lineages we studied.

Conclusions

Rapid evolution of ampliconic genes cannot be attributed solely to residence on the Y chromosome, nor to variation between primate lineages in the rate of gene family evolution. Instead other factors, such as natural selection and gene conversion, appear to play a role in driving temporal and genomic evolutionary heterogeneity in primate gene families.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2187-8) contains supplementary material, which is available to authorized users.

FISH and array CGH characterization of de novo derivative Y chromosome (Yq duplication and partial Yp deletion) in an azoospermic male

This study presents a 28-year-old infertile male who was referred to the cytogenetic laboratory for chromosomal analysis after 4 years of regular unprotected intercourse in whom non-obstructive azoospermia was revealed. Standard cytogenetic G-banding was performed on metaphase spreads and a de-novo karyotype 46,X,der(Y)(q11.22;p11.3) was identified. This analysis was followed by flourescence in-situ hybridization(FISH) and array comparative genomic hybridization (aCGH). Finally, the patient's karyotype was identified as 46,X,der(Y)(qter→q11.221::p11.31→qter).ish der(Y) (qter+,pter-,SHOX+,SRY+,Ycen+,DYZ3+;DYZ1+,qter+).arrYq11.221q12(14,448,863-59,288,511) x2, Yp11.32p11.31(104,062-266,388) x0. It is proposed that de-novo derivative monocentric Y chromosome with duplicated region Y qter→q11.221::p11.31→qter with partial deletion of Yp PAR1 region most probably can perturb the conjugation of sex chromosomes during first meiotic division of spermatogenic arrested differentiation (development).

Distribution of segmental duplications in the context of higher order chromatin organisation of human chromosome 7

BACKGROUND

Segmental duplications (SDs) are not evenly distributed along chromosomes. The reasons for this biased susceptibility to SD insertion are poorly understood. Accumulation of SDs is associated with increased genomic instability, which can lead to structural variants and genomic disorders such as the Williams-Beuren syndrome. Despite these adverse effects, SDs have become fixed in the human genome. Focusing on chromosome 7, which is particularly rich in interstitial SDs, we have investigated the distribution of SDs in the context of evolution and the three dimensional organisation of the chromosome in order to gain insights into the mutual relationship of SDs and chromatin topology.

RESULTS

Intrachromosomal SDs preferentially accumulate in those segments of chromosome 7 that are homologous to marmoset chromosome 2. Although this formerly compact segment has been re-distributed to three different sites during primate evolution, we can show by means of public data on long distance chromatin interactions that these three intervals, and consequently the paralogous SDs mapping to them, have retained their spatial proximity in the nucleus. Focusing on SD clusters implicated in the aetiology of the Williams-Beuren syndrome locus we demonstrate by cross-species comparison that these SDs have inserted at the borders of a topological domain and that they flank regions with distinct DNA conformation.

CONCLUSIONS

Our study suggests a link of nuclear architecture and the propagation of SDs across chromosome 7, either by promoting regional SD insertion or by contributing to the establishment of higher order chromatin organisation themselves. The latter could compensate for the high risk of structural rearrangements and thus may have contributed to their evolutionary fixation in the human genome.

Conserved noncoding elements follow power-law-like distributions in several genomes as a result of genome dynamics

Conserved, ultraconserved and other classes of constrained elements (collectively referred as CNEs here), identified by comparative genomics in a wide variety of genomes, are non-randomly distributed across chromosomes. These elements are defined using various degrees of conservation between organisms and several thresholds of minimal length. We here investigate the chromosomal distribution of CNEs by studying the statistical properties of distances between consecutive CNEs. We find widespread power-law-like distributions, i.e. linearity in double logarithmic scale, in the inter-CNE distances, a feature which is connected with fractality and self-similarity. Given that CNEs are often found to be spatially associated with genes, especially with those that regulate developmental processes, we verify by appropriate gene masking that a power-law-like pattern emerges irrespectively of whether elements found close or inside genes are excluded or not. An evolutionary model is put forward for the understanding of these findings that includes segmental or whole genome duplication events and eliminations (loss) of most of the duplicated CNEs. Simulations reproduce the main features of the observed size distributions. Power-law-like patterns in the genomic distributions of CNEs are in accordance with current knowledge about their evolutionary history in several genomes.

Association of microsatellite pairs with segmental duplications in insect genomes

Background

Segmental duplications (SDs), also known as low-copy repeats, are DNA sequences of length greater than 1 kb which are duplicated with a high degree of sequence identity (greater than 90%) causing instability in genomes. SDs are generally found in the genome as mosaic forms of duplicated sequences which are generated by a two-step process: first, multiple duplicated sequences are aggregated at specific genomic regions, and then, these primary duplications undergo multiple secondary duplications. However, the mechanism of how duplicated sequences are aggregated in the first place is not well understood.

Results

By analyzing the distribution of microsatellite sequences among twenty insect species in a genome-wide manner it was found that pairs of microsatellites along with the intervening sequences were duplicated multiple times in each genome. They were found as low copy repeats or segmental duplications when the duplicated loci were greater than 1 kb in length and had greater than 90% sequence similarity. By performing a sliding-window genomic analysis for number of paired microsatellites and number of segmental duplications, it was observed that regions rich in repetitive paired microsatellites tend to get richer in segmental duplication suggesting a “rich-gets-richer” mode of aggregation of the duplicated loci in specific regions of the genome. Results further show that the relationship between number of paired microsatellites and segmental duplications among the species is independent of the known phylogeny suggesting that association of microsatellites with segmental duplications may be a species-specific evolutionary process. It was also observed that the repetitive microsatellite pairs are associated with gene duplications but those sequences are rarely retained in the orthologous genes between species. Although some of the duplicated sequences with microsatellites as termini were found within transposable elements (TEs) of Drosophila, most of the duplications are found in the TE-free and gene-free regions of the genome.

Conclusion

The study clearly suggests that microsatellites are instrumental in extensive sequence duplications that may contribute to species-specific evolution of genome plasticity in insects.

Variable patterns of Y chromosome homology in Akodontini rodents (Sigmodontinae): a phylogenetic signal revealed by chromosome painting

The Akodontini is the second most speciose tribe of sigmodontine rodents, one of the most diverse groups of neotropical mammals. Molecular phylogenetic analyses are discordant regarding the interrelationships of genera, with low support for some clades. However, two clades are concordant, one (clade A) with Akodon sensu strictu (excluding Akodon serrensis), "Akodon" serrensis, Bibimys, Deltamys, Juscelinomys, Necromys, Oxymycterus, Podoxymys, Thalpomys and Thaptomys, and another (clade B) with Blarinomys, Brucepattersonius, Kunsia, Lenoxus and Scapteromys. Here, we present chromosome painting using Akodon paranaensis (APA) Y paint, after suppression of simple repetitive sequences, on ten Akodontini genera. Partial Y chromosome homology, in addition to the homology already reported on the Akodon genus, was detected on the Y chromosomes of "A." serrensis, Thaptomys, Deltamys, Necromys and Thalpomys and on Y and X chromosomes in Oxymycterus. In Blarinomys, Brucepattersonius, Scapteromys and Kunsia, no APA Y signal was observed using different hybridization conditions; APA X paint gave positive signals only on the X chromosome in all genera. The Y chromosome homology was variable in size and positioning among the species studied as follow: (1) whole acrocentric Y chromosome in Akodon and "A." serrensis, (2) Yp and pericentromeric region in submetacentric Y of Necromys and Thaptomys, (3) pericentromeric region in acrocentric Y of Deltamys, (4) distal Yq in the acrocentric Y chromosome of Thalpomys and (5) proximal Yq in the acrocentric Y and Xp in the basal clade A genus Oxymycterus. The results suggest that the homology involves pairing (pseudoautosomal) and additional regions that have undergone rearrangement during divergence. The widespread Y homology represents a phylogenetic signal in Akodontini that provides additional evidence supporting the monophyly of clade A. The findings also raise questions about the evolution of the pseudoautosomal region observed in Oxymycterus. The Y chromosomes of these closely related species seem to have undergone dynamic rearrangements, including restructuring and reduction of homologous segments. Furthermore, the changes observed may indicate progressive attrition of the Y chromosome in more distantly related species.

Y-Chromosome variation in hominids: intraspecific variation is limited to the polygamous chimpanzee

Background

We have previously demonstrated that the Y-specific ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y) varied with respect to copy number and position among chimpanzees (Pan troglodytes). In comparison, seven Y-chromosomal lineages of the bonobo (Pan paniscus), the chimpanzee's closest living relative, showed no variation. We extend our earlier comparative investigation to include an analysis of the intraspecific variation of these genes in gorillas (Gorilla gorilla) and orangutans (Pongo pygmaeus), and examine the resulting patterns in the light of the species' markedly different social and mating behaviors.

Methodology/Principal Findings

Fluorescence in situ hybridization analysis (FISH) of DAZ and CDY in 12 Y-chromosomal lineages of western lowland gorilla (G. gorilla gorilla) and a single lineage of the eastern lowland gorilla (G. beringei graueri) showed no variation among lineages. Similar findings were noted for the 10 Y-chromosomal lineages examined in the Bornean orangutan (Pongo pygmaeus), and 11 Y-chromosomal lineages of the Sumatran orangutan (P. abelii). We validated the contrasting DAZ and CDY patterns using quantitative real-time polymerase chain reaction (qPCR) in chimpanzee and bonobo.

Conclusion/Significance

High intraspecific variation in copy number and position of the DAZ and CDY genes is seen only in the chimpanzee. We hypothesize that this is best explained by sperm competition that results in the variant DAZ and CDY haplotypes detected in this species. In contrast, bonobos, gorillas and orangutans—species that are not subject to sperm competition—showed no intraspecific variation in DAZ and CDY suggesting that monoandry in gorillas, and preferential female mate choice in bonobos and orangutans, probably permitted the fixation of a single Y variant in each taxon. These data support the notion that the evolutionary history of a primate Y chromosome is not simply encrypted in its DNA sequences, but is also shaped by the social and behavioral circumstances under which the specific species has evolved.

Y chromosomal variation tracks the evolution of mating systems in chimpanzee and bonobo

The male-specific regions of the Y chromosome (MSY) of the human and the chimpanzee (Pan troglodytes) are fully sequenced. The most striking difference is the dramatic rearrangement of large parts of their respective MSYs. These non-recombining regions include ampliconic gene families that are known to be important for male reproduction,and are consequently under significant selective pressure. However, whether the published Y-chromosomal pattern of ampliconic fertility genes is invariable within P. troglodytes is an open but fundamental question pertinent to discussions of the evolutionary fate of the Y chromosome in different primate mating systems. To solve this question we applied fluorescence in situ hybridisation (FISH) of testis-specific expressed ampliconic fertility genes to metaphase Y chromosomes of 17 chimpanzees derived from 11 wild-born males and 16 bonobos representing seven wild-born males. We show that of eleven P. troglodytes Y-chromosomal lines, ten Y-chromosomal variants were detected based on the number and arrangement of the ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y)-a so-far never-described variation of a species' Y chromosome. In marked contrast, no variation was evident among seven Y-chromosomal lines of the bonobo, P. paniscus, the chimpanzee's closest living relative. Although, loss of variation of the Y chromosome in the bonobo by a founder effect or genetic drift cannot be excluded, these contrasting patterns might be explained in the context of the species' markedly different social and mating behaviour. In chimpanzees, multiple males copulate with a receptive female during a short period of visible anogenital swelling, and this may place significant selection on fertility genes. In bonobos, however, female mate choice may make sperm competition redundant (leading to monomorphism of fertility genes), since ovulation in this species is concealed by the prolonged anogenital swelling, and because female bonobos can occupy high-ranking positions in the group and are thus able to determine mate choice more freely.

Y chromosome STR allelic and haplotype diversity in five ethnic Tamil populations from Tamil Nadu, India

We have analyzed 17 Y chromosomal STR loci in a population sample of 154 unrelated male individuals of the Tamil ethnic group residing in the state of Tamil Nadu, Southern India using AmpFlSTR(R) Yfiler PCR amplification kit. The population samples consist of the following castes: Kongu Gounder (KOG), Nadar Hindu (NAH), Agamudayar (AGA), Parayar (PAR) and other Tamil individuals (MCT) of mixed castes. A total of 152 unique haplotypes were identified among the 154 individuals studied. The haplotype diversity was found to be 0.9935 or higher for all the five groups. The results of population pairwise Fst p values indicate no statistically significant differentiation between the five populations in this study, but the results were highly significant when compared with 12 other global populations (p<0.05). Comparison of populations in this study with other national and global populations using Principal co-ordinate analysis (PCA) using Rst distance matrix indicates a delineation of all the Indian populations from other unrelated populations.

The quest for Y-chromosomal markers - methodological strategies for mammalian non-model organisms

Tracing maternal and paternal lineages independently to explore breeding systems and dispersal strategies in natural populations has been high on the wish-list of evolutionary biologists. As males are the heterogametic sex in mammals, such sex-specific patterns can be indirectly observed when Y chromosome polymorphism is combined with mitochondrial sequence information. Over the past decade, Y-chromosomal markers applied to human populations have revealed remarkable differences in the demographic history and behaviour between the sexes. However, with a few exceptions, genetic data tracing the paternal line are lacking in most other mammalian species. This deficit can be attributed to the difficulty of developing Y-specific genetic markers in non-model organisms and the general low levels of polymorphisms observed on the Y chromosome. Here, we present an overview of the currently employed strategies for developing paternal markers in mammals. Moreover, we review the practical feasibility and requirements of various methodological strategies and highlight their future prospects when combined with new molecular techniques such as next generation sequencing.

Genomic analyses of sex chromosome evolution

Sex chromosomes have evolved multiple times in many taxa. The recent explosion in the availability of whole genome sequences from a variety of organisms makes it possible to investigate sex chromosome evolution within and across genomes. Comparative genomic studies have shown that quite distant species may share fundamental properties of sex chromosome evolution, while very similar species can evolve unique sex chromosome systems. Furthermore, within-species genomic analyses can illuminate chromosome-wide sequence and expression polymorphisms. Here, we explore recent advances in the study of vertebrate sex chromosomes achieved using genomic analyses.

Power-laws in the genomic distribution of coding segments in several organisms: an evolutionary trace of segmental duplications, possible paleopolyploidy and gene loss

Large-scale features of the spatial arrangement of protein-coding segments (PCS) are investigated by means of the inter-PCS spacers' size distributions, which have been found to follow power-laws. Linearity in double-logarithmic scale extends to several orders of magnitude in the genomes of organisms as disparate as mammals, insects and plants. This feature is also present in the most compact eukaryotic genomes and in half of the examined bacteria, despite their very limited non-coding space. We have tried to determine the sequence of events in the course of genomes' evolution which may account for the formation of the observed size distributions. The proposed mechanism essentially includes two types of events: (i) segmental duplications (and possibly paleopolyploidy), and (ii) the subsequent loss of most of the duplicated genes. It is shown by computer simulations that the formulated scenario generates power-law-like inter-PCS spacers' size distributions, which remain robust for a variety of parameter choices, even if insertion of external sequences, such as viruses or proliferating retroelements is included. Moreover, power-laws are preserved after most of the non-coding DNA has been removed, thus explaining the finding of this pattern in genomes as compact as that of Takifugu rubripes.

Complex evolution of a Y-chromosomal double homeobox 4 (DUX4)-related gene family in hominoids

The human Y chromosome carries four human Y-chromosomal euchromatin/heterochromatin transition regions, all of which are characterized by the presence of interchromosomal segmental duplications. The Yq11.1/Yq11.21 transition region harbours a peculiar segment composed of an imperfectly organized tandem-repeat structure encoding four members of the double homeobox (DUX) gene family. By comparative fluorescence in situ hybridization (FISH) analysis we have documented the primary appearance of Y-chromosomal DUX genes (DUXY) on the gibbon Y chromosome. The major amplification and dispersal of DUXY paralogs occurred after the gibbon and hominid lineages had diverged. Orthologous DUXY loci of human and chimpanzee show a highly similar structural organization. Sequence alignment survey, phylogenetic reconstruction and recombination detection analyses of human and chimpanzee DUXY genes revealed the existence of all copies in a common ancestor. Comparative analysis of the circumjacent beta-satellites indicated that DUXY genes and beta-satellites evolved in concert. However, evolutionary forces acting on DUXY genes may have induced amino acid sequence differences in the orthologous chimpanzee and human DUXY open reading frames (ORFs). The acquisition of complete ORFs in human copies might relate to evolutionary advantageous functions indicating neo-functionalization. We propose an evolutionary scenario in which an ancestral tandem array DUX gene cassette transposed to the hominoid Y chromosome followed by lineage-specific chromosomal rearrangements paved the way for a species-specific evolution of the Y-chromosomal members of a large highly diverged homeobox gene family.

Evolutionary analysis of the highly dynamic CHEK2 duplicon in anthropoids

BACKGROUND

Segmental duplications (SDs) are euchromatic portions of genomic DNA (> or = 1 kb) that occur at more than one site within the genome, and typically share a high level of sequence identity (>90%). Approximately 5% of the human genome is composed of such duplicated sequences. Here we report the detailed investigation of CHEK2 duplications. CHEK2 is a multiorgan cancer susceptibility gene encoding a cell cycle checkpoint kinase acting in the DNA-damage response signalling pathway. The continuous presence of the CHEK2 gene in all eukaryotes and its important role in maintaining genome stability prompted us to investigate the duplicative evolution and phylogeny of CHEK2 and its paralogs during anthropoid evolution.

RESULTS

To study CHEK2 duplicon evolution in anthropoids we applied a combination of comparative FISH and in silico analyses. Our comparative FISH results with a CHEK2 fosmid probe revealed the single-copy status of CHEK2 in New World monkeys, Old World monkeys and gibbons. Whereas a single CHEK2 duplication was detected in orangutan, a multi-site signal pattern indicated a burst of duplication in African great apes and human. Phylogenetic analysis of paralogous and ancestral CHEK2 sequences in human, chimpanzee and rhesus macaque confirmed this burst of duplication, which occurred after the radiation of orangutan and African great apes. In addition, we used inter-species quantitative PCR to determine CHEK2 copy numbers. An amplification of CHEK2 was detected in African great apes and the highest CHEK2 copy number of all analysed species was observed in the human genome. Furthermore, we detected variation in CHEK2 copy numbers within the analysed set of human samples.

CONCLUSION

Our detailed analysis revealed the highly dynamic nature of CHEK2 duplication during anthropoid evolution. We determined a burst of CHEK2 duplication after the radiation of orangutan and African great apes and identified the highest CHEK2 copy number in human. In conclusion, our analysis of CHEK2 duplicon evolution revealed that SDs contribute to inter-species variation. Furthermore, our qPCR analysis led us to presume CHEK2 copy number variation in human, and molecular diagnostics of the cancer susceptibility gene CHEK2 inside the duplicated region might be hampered by the individual-specific set of duplicons.