Frequent copy number variants in a cohort of Mexican-Mestizo individuals
Molecular Cytogenetics volume 16, Article number: 2 (2023)
The human genome presents variation at distinct levels, copy number variants (CNVs) are DNA segments of variable lengths that range from several base pairs to megabases and are present at a variable number of copies in human genomes. Common CNVs have no apparent influence on the phenotype; however, some rare CNVs have been associated with phenotypic traits, depending on their size and gene content. CNVs are detected by microarrays of different densities and are generally visualized, and their frequencies analysed using the HapMap as default reference population. Nevertheless, this default reference is inadequate when the samples analysed are from people from Mexico, since population with a Hispanic genetic background are minimally represented. In this work, we describe the variation in the frequencies of four common CNVs in Mexican-Mestizo individuals.
In a cohort of 147 unrelated Mexican-Mestizo individuals, we found that the common CNVs 2p11.2 (99.6%), 8p11.22 (54.5%), 14q32.33 (100%), and 15q11.2 (71.1%) appeared with unexpectedly high frequencies when contrasted with the HapMap reference (ChAS). Yet, while when comparing to an ethnically related reference population, these differences were significantly reduced or even disappeared.
The findings in this work contribute to (1) a better description of the CNVs characteristics of the Mexican Mestizo population and enhance the knowledge of genome variation in different ethnic groups. (2) emphasize the importance of contrasting CNVs identified in studied individuals against a reference group that—as best as possible—share the same ethnicity while keeping this relevant information in mind when conducting CNV studies at the population or clinical level.
The human genome presents distinct variations: single nucleotide variants, insertion‒deletion of a few nucleotides, repetitive sequences of a variable number of nucleotides, and structural variants. The human genome is 3.1 giga base pairs (Gb) in size, distributed in 23 pairs of chromosomes so that every individual inherits one copy from each parent and has two copies or copy number (CN) CN = 2 of every locus in their cells (except for the X and Y loci in males). In 2004, two independent groups, Iafrate and Sebat [1, 2], described for the first time the genome-wide presence of large-scale copy number variations in the human genome; these variants involve gains or losses of several to hundreds of kilobases (kb) of genomic DNA among phenotypically normal individuals, and these copy number variants (CNVs) are an important source of human genomic variation . Common CNVs have a minor allele frequency (MAF) > 5% and no apparent influence on the phenotype. However, some rare CNVs, with MAF < 1%, have been associated with phenotypic traits , and a few of them have shown clinical relevance depending on their size, gene content, or when overlap with genes that manifest haploinsufficiency (when CNV is CN = 1) or triplosensitivity (when CNV is CN > 2) [5, 6]. In some occasions, their interaction with additional genetic or environmental factors may influence whether CNVs have a detectable phenotypic effect [7, 8]. CNV represents a significant proportion of the total genetic variability in all human populations [9,10,11]. A CNV was initially defined as a DNA segment larger than one kb with a variable copy number compared to a reference genome . Nevertheless, this term has extended to all quantitative variations in the genome, including tandem repeats, deletions, and duplications .
The accurate detection and interpretation of CNVs are essential for research and clinical diagnostics because CNVs may be a general population variant or associated with pathology . Therefore, it has become increasingly important to detect the presence of common and harmless CNVs within different ethnic groups to avoid misinterpretation of pathological variants.
The Database of Genomic Variants (DGV) was initiated to provide a publicly accessible, comprehensive, and curated catalog of genomic variants . The DGV comprises the CNVs and structural variations found in the genomes of control individuals from worldwide populations. However, only three studies in this database included individuals with Mexican ancestry, and they are the same in all three studies [16,17,18]. Furthermore, currently it only comprises individuals born in the US with parents and grandparents of Mexican origin, representing only 3% of the total Mexican population .
In general, this could represent an issue since the contribution of the ethnic group is diluted when the entire group is used as a reference.
Therefore, it is necessary to report information on this type, especially for underrepresented populations.
Several methodologies are commonly used for CNV genome-wide study; one of the most robust methods is microarray analysis, mainly when an array contains a high density of probes. Most CNV studies are performed with single nucleotide polymorphism (SNP) microarrays, which, in addition to providing information on CNV in the genome, also provide information at the allele level, highlighting the presence of regions/runs of homozygosity called long contiguous stretches of homozygosity (LCSH) and providing essential clues regarding parental relatedness (consanguinity), uniparental disomy, chromosomal recombination or rearrangements [20, 21].
This work describes the frequency of four common CNVs found using SNP microarrays in unrelated Mexican-Mestizo individuals, who were healthy or had aneuploidies, and in whom structural alterations in their genome was not suspected. The study highlights the frequency variation of these common CNVs depending on the ethnic origin of the population used as reference for comparison.
It is essential to recognize the ethnic differences in the distribution and frequency of CNVs, not only to recognize their contribution to structural chromosomal variation but also to ease recognition of possible associations between CNVs and phenotypic characteristics that may or may not be pathogenic.
We studied the presence of CNV in a population of 147 Mexican-Mestizo individuals, 50 phenotypically normal (MXM, 30 females and 20 males) and 97 patients with aneuploidies (MXM_A, 41 with trisomy 21, 28 with monosomy X, 11 with trisomy 13 and 17 with trisomy 18), which are part of a more extensive study that aims to find genomic differences among these groups. In this first study, we show four CNVs in common among all these groups, with a very high frequency in our population when contrasted with distant populations, as is the mix of the HapMap reference that includes African, Asian, and Caucasian individuals.
Long contiguous stretches of homozygosity
To recognize the presence of inbreeding and consanguinity in our sample population, we analyzed the regions with homozygosity or long contiguous stretches of homozygosity (LCSH) with Affymetrix Genome-Wide Human SNP Array 6.0 (GW SNP 6.0). Eighteen of the 50 healthy controls (MXM) studied showed at least one autosomal region with LCSH. In the aneuploidy group (MXM_A), 73/97 presented at least one autosomal region with LCSH.
In addition, we analyzed 99 CEL files from 99 blood samples of Mexican women with breast cancer (MXM_WBC) present in public databases, where 75/99 presented at least one autosomal region with LCSH.
Neither of these LCSH regions overlaps with the CNV reported in this study, and more than 42% are common LCSHs, such as the 2q11.1q11.2, 11p11.2p12.1, 16p11.2p11.1, and 20q11.21q11.23 regions previously reported [22,23,24]. The inbreeding coefficient (F) calculated in the three groups supports the absence of inbreeding or consanguinity (Table 1).
We deduced the ancestry in our MXM and MXM_A groups compared with the Mexican ancestry in Los Angeles CA, USA (MXL) included in the HapMap project (64 Mexican individuals born in the US). First, we noticed that our population (MXM and MXM_A) has almost the same genetic content of European descent from Utah, USA (CEU), Native American (NAT), and Yoruba from Nigeria (YRI), then those studied for the HapMap, with NAT being the most significant contribution, even greater than for MXL, 0.488 vs 0.678 and 0.618, then of Caucasians and finally Africans (Fig. 1).
Copy number variants with high frequency in Mexican-Mestizo population
We found CNV in almost all individuals and in all the chromosomes; however, we observed four polymorphisms that occurred in chromosomes 2 (2.9 megabases (Mb)), 8 (0.15 Mb), 14 (1.02 Mb), and 15 (0.76 Mb) with a high frequency in the general population of Mexican-Mestizo (MXM), in Mexican-Mestizo with aneuploidy (MXM_A), and in the MXM_WBC group (Figs. 2 and 3) (Table 2).
In region 2p11.2 for MXM healthy tissues, there was a gain of 3 and 4 copies, CN = 3 in 58% and CN = 4 in 41%, respectively, and only one case (0.67%) was CN = 2 (Fig. 2a); for MXM_A, 21% of cases presented with CN = 3 and 79% with CN = 4, with an overall copy gain across all groups of 99.6% (Fig. 2b, Table 2).
For chromosome 8p11.2, almost half of the cases, 42% of MXM and 51% of MXM_A, have CN = 2, but 49% MXM and 41% MXM_A have losses with CN = 1, and the rest (9% and 8%, respectively) have null CN = 0. The frequency of CN ≠ 2 across all the groups is 54.5% (Fig. 2, Table 2).
The CNV located in 14q32.33 is present as CN = 3 in 55% of MXM healthy samples and 62% in MXM_A samples and as CN = 4 in 45% of MXM and 38% of MXM_A subjects; there is not a single individual who shows CN = 2, 100% of individuals presented with gain in this region (Table 2, Fig. 2a and b).
Finally, in 15q11.2, approximately 70% of the samples had either gains or losses: 40% and 38% of the MXM healthy tissues and MXM_A samples, respectively, showed loss CN = 1, while 31.5% of MXM and 31% of MXM_A had gain with CN = 3, giving a frequency of CN ≠ 2 across all the groups of 71.1%.
The common CNV in MXM healthy subjects and tissues observed utilizing HapMap as a population of reference is shown in Fig. 3a.
Validation of CN by qPCR
We performed qPCR to corroborate the CN obtained with GW SNP Array 6.0 and analyzed the gDNA for a representative sample of the individuals studied, including nine MXM individuals and seven MXM_A individuals. With this method, 15 out of 16 samples (93.75%) in locus 2p11.2 showed CN from 4 to 17 copies (ChAS only detects up to CN = 4). For chromosome 8, locus 8p11.22 A total of 14/16 (87.5%) samples showed < 2 copies. For 14q32.33, we corroborated the gain of the CNV in 13/16 (81.25%) with 3 to 5 copies; for chromosome 15, we found both losses and gains by GW, and we found 4 patients with loss and 6 with copy gains (25% and 37.5%, respectively) (Table 3, Fig. 3).
CRMA v2 analysis
The CNV found with the Chromosome Analysis Suite (ChAS) was validated with an additional CNV calling method CRMA v2.
We declare the raw data set conformed by 50 MXM CEL files (Genome-Wide Human SNP Array 6.0).
This method could compare the copy number estimates in 14q32.33 for each sample. Figure 4 shows the CNV calling for chromosome 14 in subject C12 in the A) C12 vs MXM reference. B) C12 vs Spanish population (IBS) reference. C) C12 vs CEU reference. D) CN in the C12 vs. YRI reference. E) CN in C12 vs. two populations included in the HapMap reference (CEU + YRI). It is important to note that the difference in CN in our sample depends on the reference used.
It is clear that in A and B, the 14q32.33 region does not show changes with respect to the Mexican and Spanish reference population, emphasizing our Mestizo race formed by Spanish and Native American groups. In contrast, comparisons with CEU, YRI and CEU + YRI make evident a gain of probes located in the chr14:106,078,230–107,100,266 region, showing our racial differences with these populations.
In the tested samples, we found four high-frequency CNVs in the following loci: 2p11.2, 8p11.22, 14q32.33, and 15q11.2. No genes were found at 2p11.2 or 15q11.2, while the CNV at 8p11.22 included the genes ADAM5 (HGNC: 212) and ADAM3A (HGNC: 209). The 14q32.33 locus comprises the genes KIAA0125 (HGNC: 9834), ADAM6 (HGNC: 8755), LINC00226 (HGNC: 338,004), IGHD (HGNC:5480), IGHG3 (HGNC:5527), and IGHG1 (HGNC:5525); none of these genes has been associated with Mendelian diseases (OMIM) related to genetic dosage (https://search.clinicalgenome.org/kb/gene-dosage/region/ISCA-37477), and there are no regulatory sequences that could modify gene expression in any of the regions of the four CNV; neither of our CNV reported overlap with LCSH regions.
Additionally, by searching these specific regions in databases such as the sSMC (small supernumerary marker chromosomes) database http://cs-tl.de/DB/CA/sSMC/0-Start.html and the chromosomal heteromorphisms database http://cs-tl.de/DB/CA/HCM/0-Start.html, we found no evidence that a number other than CN = 2 in the 2p11.2, 8p11.22 and 15q11.2 regions or 14q32.33 duplications is associated with pathology, indicating that numerical variants in these regions have no impact on the phenotype, probably due to the absence of dosage-sensitive genes .
The high frequency of some variants could be due to the existence of inbreeding or consanguinity in the population. However, according to our results on the percentage of homozygosity of the LCSH regions and the degree of inbreeding (F value) , there is no indication of inbreeding or consanguinity in the studied groups.
CNVs are genomic variants that confer high variability among individuals and have been recognized for several years [1,2,3]. The first CNV studies were performed in the HapMap collection, which comprises 270 subjects divided into four populations: 90 Yoruba from Ibadan, Nigeria (YRI), 90 European descent from Utah, USA (CEU), 45 Japanese from Tokyo, Japan (JPT) and 45 Han Chinese from Beijing, China (CHB) . In that study, the gDNA was analyzed with 500 K Affymetrix SNP arrays, and the CNV was determined using the software provided by Affymetrix. No individual presented CNVs at 2p11.2, 8p11.22, 14q32.33, or 15q11.2. In 2010, a study of large CNVs (1 kb) performed in 450 samples of HapMap (180 CEU, 180 YRI, 45 JPT and 45 CHB) reported CNV at 14q32.33 with a frequency less than 1% and in 2p11.2, 8p11.22 and 15q11.2 with a frequency near 5% .
Later in 2015, the phase 3 1000 genome study included 2504 healthy individuals, only 64 of them (2.5% of the total sample analyzed) were individuals of Mexican descent, and they all were living in Los Angeles, California (MXL)  This study in which Hispanic genomes were underrepresented showed low frequencies for CNVs at 2p11.2, 8p11.22, 14q32.33 and 15q11.2, coincident with the data previously reported by Conrad et al..
Now we know that the CNV rate of occurrence has wide differences according to the ethnic group that is studied, as it occurs with other polymorphic traits such as SNPs or microsatellites . Furthermore, frequencies may vary depending on the microarray platform used and the CNV calling algorithms [20, 28].
The results presented here are based on the Affymetrix microarray platform and were analyzed using Affymetrix software, with Affymetrix annotations NetAffx 33.1, which is the default software provided with this commercial platform. This workflow is the one being used to find CNV for clinical and research studies in Mexico and in several other countries, however, in this software there is little representation of the Latin American population.
In fact, there are studies reporting association of some of the polymorphisms studied, to pathological conditions [29,30,31]; one of these studies was conducted in Brazil and described the association of CNV gain of the 14q32.33 region with dental tumors , the reported copy gain in this 14q32.33 locus is as compared to the default reference genome without taking into account possible ethnic differences of the studied population, an ethnic matched reference genome would not have shown gain in this region. Indeed in 2020, Godoy et al. identified a CNV obtained from three different microarray platforms from a Brazilian population to conform the Brazilian CNV database. They found that a 14q32.33 gain was present in 97.8% of the samples studied , similar to that found in our study in 100% of the samples. Interestingly, the origins of the Brazilian population and ours have in common the mixture of the Iberian population that conquered us centuries ago.
The foregoing was addressed in this work by CNV calling (CRMA), where a Spanish population was used as reference, with this better-suited ethnical matched reference the gain of the CNV at 14q32.33 was reduced, and even did not appear as copy number gain (Fig. 4a and b).
The findings presented in this study contribute to the description of the frequency of common CNV in the Mexican-Mestizo population. Because microarrays are widely used tools in diagnostic and research contexts, this work is a clear example of why CNVs analysis must be carried out using an ethnically appropriate reference population. The significance of this is being recognized more wildly, as a consequence of this there is a conscient effort to include more diverse ethnic groups in the wildly used reference databases like HapMap . This inclusive policy is fundamental to have a more accurate representation of the human genome.
We found four CNVs, 2p11.2, 8p11.22, 14q32.33, and 15q11.2, with a high frequency in the Mexican-Mestizo population when contrasted with the HapMap population of reference (ChAS), while when using an ethnically related population as a reference, the differences were reduced or disappeared, highlighting the importance of analyzing the CNVs of the studied individuals with a reference group that (as far as possible) shares the same ethnicity. These findings contribute to a better description of the CNVs characteristics of the Mexican-Mestizo population and enhance the knowledge of genome variation in different ethnic groups.
Population, samples, and DNA extraction
We included blood samples from subjects with Mexican-Mestizo origin who consented to participating in this study. Fifty were healthy, unrelated volounteers from the general population. All subjects were over 18 years of age, had no history of genetic diseases and at the time of sample donation were not suffering from any infectious disease nor taking any medications (MXM, 30 females and 20 males). All participants had a normal karyotype 46,XX in females and 46,XY in males (20 corroborated by G banding in 25 metaphases, and 30 by chromosome microarray). We also included 97 aneuploid patients (MXM_A, 41 with trisomy 21, 28 with monosomy X, 11 with trisomy 13 and 17 with trisomy 18), with karyotype showing regular trisomy as well as X-monosomy without mosaicism, which was corroborated by interphase FISH (1000 cells analyzed).
Genomic DNA was extracted from the blood samples obtained from the participants with the saline precipitation method (Gentra Puregene Kit, QIAGEN, Venlo, Limburg, NL).
Genome-wide Human SNP arrays
We analyzed the gDNA with the Affymetrix GeneChip® Genome-Wide Human SNP Array 6.0 (Santa Clara, CA, USA). This array contains 906,600 SNP probes and 946,000 nonpolymorphic oligonucleotides; the median intermarker distance over all 1.8 million SNP and copy number markers combined is less than 700 bases.
The procedures for DNA digestion, ligation, PCR amplification, fragmentation, labeling, denaturing and hybridization into the array were performed in 147 DNA samples (two DNA samples were not included because they did not pass quality controls) according to the protocols provided by the supplier. Arrays were stained and washed in the Affymetrix GeneChip Fluidic Station 450 and scanned using an Affymetrix GeneChip Scanner 3000 7G (Affymetrix, Santa Clara, CA, USA). We analyzed the files obtained with the appropriate bioinformatics tools.
Long contiguous stretches of homozygosity
We visualized the long contiguous stretches of homozygosity (LCSH) in Chromosome Analysis Suite (ChAS) software version 4.1, provided by Affymetrix (Affymetrix, Santa Clara, CA, USA). For the analysis, we used the NetAffx 33 hg19 annotation files (http://www.affymetrix.com). For LCSH > 3 Mb, the analysis configuration was set at LOH with marker count = 50 and size = 3000 kb, and for LCSH > 5 Mb, it was set at marker count = 50 and size = 5000 kb.
Estimation of the coefficient of inbreeding (F)
Individual inbreeding coefficients (F) were estimated using LCSH > 3 Mb data; F was the total length of autosomal LCSH in kb divided by the total autosomal size covered by the Genome-Wide Human SNP Array 6.0 (2,881,033,286 kb for hg19). We report the average of the F value for each group. An F value of 0.25 could reflect a first-degree parental relationship, 0.125 a second-degree relationship, 0.0625 a third-degree relationship and 0.03125 a fourth-degree relationship .
To compare our studied population, we used the following CEL files obtained from Affymetrix GeneChip® Genome-Wide Human SNP Array 6.0 (Santa Clara, CA, USA):
(a) Ninety-nine files from healthy tissues obtained from women with breast cancer (MXM_WBC)  from the GEO open database (GSE87048).
Together with the MXM, they will be the Healthy tissues.
(b) 30 Spanish (IBS) files downloaded from GEO data set GSE67047 
(c) 30 Caucasian (CEU) and 10 Yoruba (YRI), part of the population studied in the HapMap project, CEL files kindly provided by Affymetrix/Thermo Fisher.
The MAP and PED files obtained by Genotyping Console software (Affymetrix) were used for the ancestry analysis; we eliminated two files that did not pass quality control.
The global ancestry of 49/50 healthy Mexican-Mestizo files (MXM) and 96/97 individuals with aneuploidy (MXM_A) was deduced through a supervised maximum likelihood ADMIXTURE approach from K = 2 to K = 3 ancestral components and compared with the global ancestry of the group MXL (Mexican ancestry from Los Angeles California US) (n = 64). For this analysis, genotypic frequencies of Northern European (CEU) (n = 99), Yoruba (YRI) (n = 108), and Native American (NAT) (n = 43), reported in the 1000 genomes project (1KGP), were considered as reference parental populations for demographic and historical reasons.
Copy number variant (CNV) calling
CNVs were visualized with Affymetrix software (Affymetrix, Santa Clara, CA, USA). GeneChip Command Console (AGCC) software was used to generate the CEL files and ARR files from each microarray scanned. The CEL and ARR files were analyzed with Genotyping Console software (Affymetrix, Santa Clara, CA, USA) to obtain the CN data files (CNCHP) MAP and PED files for ancestry. CNCHP data were analyzed using Chromosome Analysis Suite (ChAS) software version 4.1 (Affymetrix, Santa Clara, CA, USA).
The annotation file used in our analysis can be found on the Affymetrix website, listed as NetAffx 33.1 (hg19). The reporting threshold of the copy number was set at "High Resolution" settings, which gave us a total of gains or losses of 100 kb with a marker count ≥ 50. The copy number variants were compared among all samples. The CN analysis in ChAS software (http://www.affymetrix.com) has two ways of showing the results: one as a table that includes both coordinates and CN for each region; or a graphical interphase that displays the results of CN in an image of each of the 23 chromosomes, highlighting the regions with CN > 2 in blue arrowhead, regions with CN = 1 in red arrowhead, or CN regions with copy number (CN = 2), which are only marked in lines (Fig. 5).
To corroborate the changes in the CN obtained by GW SNP 6.0, we analyzed the absolute changes in copy number with qPCR in a representative sample. Oligonucleotides primers were designed at the four loci with CN ≠ 2 and for a distal locus on the same chromosome, where CN = 2 was used as a copy number control (Table 4). qPCR was performed with LightCycler ® Master (Roche Diagnostics) using hydrolysis TaqMan® probes (Universal Probes, Roche) with 100 ng of DNA, 45 amplification cycles and a single quantification. For those chromosomal regions in which it is possible to have a CN = 2 control DNA, the method proposed by Livak of 2−DDCt was used ; for the loci where all subjects showed the change in the CN, the 2−DCt method was used (Table 3) .
For an additional CNV calling method, we used copy number estimation using robust multichip analysis, and this method provides full-resolution raw total copy number estimates by preprocessing and probe summarization. CRMA v2 is available in Bioconductor  and implemented in R (http://www.aroma-project.org/vignettes/CRMAv2/).
Briefly, the CEL files were analyzed for quality after the program calibrated for crosstalk between allele probe pairs, followed by normalization for 25-mer nucleotide-position probe sequence effects. Next, we analyzed the performance of robust probe-summarization, normalized the PCR fragment-length effects on summary signals, and finally calculated the full-resolution total copy numbers. When a standard reference is used, it is often the average of a pool of samples CnChipEffectFile.
The reference control samples for CRMA v2 analysis were a) the same 50 healthy individuals of MXM; b) 30 IBS, 30 CEU and 10 YRI. All groups were analyzed with the same pipeline.
We made the following comparisons: (a) MXM vs. MXM reference; (b) MXM vs. IBS reference; (c) MXM vs. CEU reference; (d) MXM vs. YRI reference; and (e) MXM vs. HapMap reference (CEU + YRI).
We reviewed publicly available databases to classify the CNVs found in our samples and the genes included:
Database of Genomic Variants (DGV, http://projects.tcag.ca/variation/).
Gene ontology (http://geneontology.org/).
The sSMC database (http://cs-tl.de/DB/CA/sSMC/0-Start.html).
The heteromorphims database (http://cs-tl.de/DB/CA/HCM/0-Start.html).
Availability of data and materials
The data that support the findings of this study are not publicly available because one part of them is available from the Instituto Nacional de Ciencias Genómicas (México), but restrictions apply to the availability of these data, which were used under license for the current study. Data are, however, available from the authors upon reasonable request, with permission of Instituto Nacional de Ciencias Genómicas (México) and following the approved ethics committee guidelines.
Utah residents (CEPH) with Northern and Western European ancestry (1KGP population)
Han Chinese from Beijing, China
Copy number variant
The Database of Genomic Variants
- GW SNP 6.0:
Genome-Wide Human SNP Array 6.0, Affymetrix
Homo sapiens (human) genome assembly GRCh37
Iberian population in Spain
Instituto Nacional de Pediatría
Long contiguous stretches of homozygosity
Minor allele frequency
Mexican Ancestry in Los Angeles CA, USA (1 KGP population)
Mexican Mestizo healthy controls
Mexican Mestizo with aneuploidy
Mexican-Woman with breast cancer (blood samples)
Native Americans (1KGP population)
Online Mendelian Inheritance in Man
Quantitative polymerase chain reaction
Small Supernumerary Marker Chromosomes
Single nucleotide polymorphism
Yoruba in Ibadan, Nigeria (1KGP population)
Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C. Detection of large-scale variation in the human genome. Nat Genet. 2004;36(9):949–51. https://doi.org/10.1038/ng1416.
Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Månér S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M. Large-scale copy number polymorphism in the human genome. Science. 2004;305(5683):525–8. https://doi.org/10.1126/science.1098918.
Feuk L, Carson AR, Scherer SW. Structural variation in the human genome. Nat Rev Genet. 2006;7:85–97. https://doi.org/10.1038/nrg1767.
Hujoel MLA, Sherman MA, Barton AR, Mukamel RE, Sankaran VJ, Terao C, Loh PR. Influences of rare copy-number variation on human complex traits. Cell Cell. 2022;185:4233–48. https://doi.org/10.1016/j.cell.2022.09.028.
Hehir-Kwa JY, Pfundt R, Veltman JA, de Leeuw N. Pathogenic or not? Assessing the clinical relevance of copy number variants. Clin Genet. 2013;84(5):415–21. https://doi.org/10.1111/cge.12242.
Rooney ER, Andersen EF, Cherry AM, Kantarci S, Kearney H, Patel A, Raca G, Ritter DI, South ST, Thorland EC, Pineda-Alvarez D, Aradhya S, Martin CL. Technical standards for the interpretation and reporting of constitutional copy-number variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen). Genet Med. 2020;22:245–57. https://doi.org/10.1038/s41436-019-0686-8.
Wellcome Trust Case Control Consortium, Craddock N, Hurles ME, Cardin N, Pearson RD, Plagnol V, Robson S, Vukcevic D, Barnes C, Conrad DF, Giannoulatou E, Holmes C, Marchini JL, Stirrups K, Tobin MD, Wain LV, Yau C, Aerts J, Ahmad T, Andrews TD, Arbury H, Attwood A, Auton A, Ball SG, Balmforth AJ, Barrett JC, Barroso I, Barton A, Bennett AJ, Bhaskar S, Blaszczyk K, Bowes J, Brand OJ, Braund PS, Bredin F, Breen G, Brown MJ, Bruce IN, Bull J, Burren OS, Burton J, Byrnes J, Caesar S, Clee CM, Coffey AJ, Connell JM, Cooper JD, Dominiczak AF, Downes K, Drummond HE, Dudakia D, Dunham A, Ebbs B, Eccles D, Edkins S, Edwards C, Elliot A, Emery P, Evans DM, Evans G, Eyre S, Farmer A, Ferrier IN, Feuk L, Fitzgerald T, Flynn E, Forbes A, Forty L, Franklyn JA, Freathy RM, Gibbs P, Gilbert P, Gokumen O, Gordon-Smith K, Gray E, Green E, Groves CJ, Grozeva D, Gwilliam R, Hall A, Hammond N, Hardy M, Harrison P, Hassanali N, Hebaishi H, Hines S, Hinks A, Hitman GA, Hocking L, Howard E, Howard P, Howson JM, Hughes D, Hunt S, Isaacs JD, Jain M, Jewell DP, Johnson T, Jolley JD, Jones IR, Jones LA, Kirov G, Langford CF, Lango-Allen H, Lathrop GM, Lee J, Lee KL, Lees C, Lewis K, Lindgren CM, Maisuria-Armer M, Maller J, Mansfield J, Martin P, Massey DC, McArdle WL, McGuffin P, McLay KE, Mentzer A, Mimmack ML, Morgan AE, Morris AP, Mowat C, Myers S, Newman W, Nimmo ER, O'Donovan MC, Onipinla A, Onyiah I, Ovington NR, Owen MJ, Palin K, Parnell K, Pernet D, Perry JR, Phillips A, Pinto D, Prescott NJ, Prokopenko I, Quail MA, Rafelt S, Rayner NW, Redon R, Reid DM, Renwick, Ring SM, Robertson N, Russell E, St Clair D, Sambrook JG, Sanderson JD, Schuilenburg H, Scott CE, Scott R, Seal S, Shaw-Hawkins S, Shields BM, Simmonds MJ, Smyth DJ, Somaskantharajah E, Spanova K, Steer S, Stephens J, Stevens HE, Stone MA, Su Z, Symmons DP, Thompson JR, Thomson W, Travers ME, Turnbull C, Valsesia A, Walker M, Walker NM, Wallace C, Warren-Perry M, Watkins NA, Webster J, Weedon MN, Wilson AG, Woodburn M, Wordsworth BP, Young AH, Zeggini E, Carter NP, Frayling TM, Lee C, McVean G, Munroe PB, Palotie A, Sawcer SJ, Scherer SW, Strachan DP, Tyler-Smith C, Brown MA, Burton PR, Caulfield MJ, Compston A, Farrall M, Gough SC, Hall AS, Hattersley AT, Hill AV, Mathew CG, Pembrey M, Satsangi J, Stratton MR, Worthington J, Deloukas P, Duncanson A, Kwiatkowski DP, McCarthy MI, Ouwehand W, Parkes M, Rahman N, Todd JA, Samani NJ, Donnelly P. Genome-wide association study of CNV in 16,000 cases of eight common diseases and 3000 shared controls. Nature. 2010; 464(7289):713–20. https://doi.org/10.1038/nature08979
Mazina V, Gerdts J, Trinh S, Ankenman K, Ward T, Dennis MY, Girirajan S, Eichler EE, Bernier R. Epigenetics of autism-related impairment: copy number variation and maternal infection. J Dev Behav Pediatr. 2015;36(2):61.
Ionita-Laza I, Rogers AJ, Lange C, Raby BA, Lee C. Genetic association analysis of copy number variation (CNV) in human disease pathogenesis. Genomics. 2009;93(1):22–6. https://doi.org/10.1016/j.ygeno.2008.08.012.
Zarrei M, MacDonald JR, Merico D, Scherer SW. A copy number variation map of the human genome. Nat Rev Genet. 2015;16:172–83. https://doi.org/10.1038/nrg3871.
Freeman JL, Perry GH, Feuk L, Redon R, McCarroll SA, Altshuler DM, Aburatani H, Jones KW, Tyler-Smith C, Hurles ME, Carter NP, Scherer SW, Lee C. Copy number variation: new insights in genome diversity. Genome Res. 2006;16(8):949–61. https://doi.org/10.1101/gr.3677206.
Scherer SW, Lee C, Birney E, Altshuler DM, Eichler EE, Carter NP, Hurles ME, Feuk L. Challenges and standards in integrating surveys of structural variation. Nat Genet. 2007;39(7):7–15. https://doi.org/10.1038/ng2093.
Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, Fitzgerald T, Hu M, Ihm CH, Kristiansson K, Macarthur DG, Macdonald JR, Onyiah I, Pang AW, Robson S, Stirrups K, Valsesia A, Walter K, Wei J; Wellcome Trust Case Control Consortium, Tyler-Smith C, Carter NP, Lee C, Scherer SW, Hurles ME. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464(7289):704–12. https://doi.org/10.1038/nature08516
Beckmann JS, Estivill X, Antonarakis SE. Copy number variants and genetic traits: closer to the resolution of phenotypic to genotypic variability. Nat Rev Genet. 2007;8(8):639–46. https://doi.org/10.1038/nrg2149.
MacDonald JR, Ziman R, Yuen RK, Feuk L, Scherer SW. The Database of Genomic Variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 2014;42:D986–92. https://doi.org/10.1093/nar/gkt958.
Sudmant PH, Kitzman JO, Antonacci F, Alkan C, Malig M, Tsalenko A, Sampas N, Bruhn L, Shendure J. 1000 Genomes Project, Eichler EE. Diversity of human copy number variation a646nd multicopy genes. Science. 2010;330(6004):641–6. doi:https://doi.org/10.1126/science.1197005
Campbell CD, Sampas N, Tsalenko A, Sudmant PH, Kidd JM, Malig M, Vu TH, Vives L, Tsang P, Bruhn L, Eichler EE. Population-genetic properties of differentiated human copy-number polymorphisms. Am J Hum Genet. 2011;88(3):317–32. https://doi.org/10.1016/j.ajhg.2011.02.004.
Sudmant PH, Rausch T, Gardner EJ, Handsaker RE, Abyzov A, Huddleston J, Zhang Y, Ye K, Jun G, Fritz MH, Konkel MK, Malhotra A, Stütz AM, Shi X, Casale FP, Chen J, Hormozdiari F, Dayama G, Chen K, Malig M, Chaisson MJP, Walter K, Meiers S, Kashin S, Garrison E, Auton A, Lam HYK, Mu XJ, Alkan C, Antaki D, Bae T, Cerveira E, Chines P, Chong Z, Clarke L, Dal E, Ding L, Emery S, Fan X, Gujral M, Kahveci F, Kidd JM, Kong Y, Lameijer EW, McCarthy S, Flicek P, Gibbs RA, Marth G, Mason CE, Menelaou A, Muzny DM, Nelson BJ, Noor A, Parrish NF, Pendleton M, Quitadamo A, Raeder B, Schadt EE, Romanovitch M, Schlattl A, Sebra R, Shabalin AA, Untergasser A, Walker JA, Wang M, Yu F, Zhang C, Zhang J, Zheng-Bradley X, Zhou W, Zichner T, Sebat J, Batzer MA, McCarroll SA; 1000 Genomes Project Consortium, Mills RE, Gerstein MB, Bashir A, Stegle O, Devine SE, Lee C, Eichler EE, Korbel JO. An integrated map of structural variation in 2,504 human genomes. Nature. 2015;526(7571):75–81. doi:https://doi.org/10.1038/nature15394
IGSR: The International Genome Sample Resource. Supporting open human variation data. https://www.internationalgenome.org/data-portal/population Citation: Susan Fairley, Ernesto Lowy-Gallego, Emily Perry, Paul Flicek, The International Genome Sample Resource (IGSR) collection of open human genomic variation resources. Nucl Acids Res 48(D1)2020:941–947. https://doi.org/10.1093/nar/gkz836
Haraksingh RR, Abyzov A, Urban AE. Comprehensive performance comparison of high-resolution array platforms for genome-wide Copy Number Variation (CNV) analysis in humans. BMC Genom. 2017;18(1):321. https://doi.org/10.1186/s12864-017-3658-x.
Bruce S, Leinonen R, Lindgren CM, Kivinen K, Dahlman-Wright K, Lipsanen-Nyman M, Hannula-Jouppi K, Kere J. Global analysis of uniparental disomy using high density genotyping arrays. J Med Genet. 2005;42(11):847–51. https://doi.org/10.1136/jmg.2005.032367.
Pajusalu S, Žilina O, Yakoreva M, Tammur P, Kuuse K, Mölter-Väär T, Nõukas M, Reimand T, Õunap K. The diagnostic utility of single long contiguous stretches of homozygosity in patients without parental consanguinity. Mol Syndromol. 2015;6:135–40. https://doi.org/10.1159/000438776.
Wang JC, Ross L, Mahon LW, Owen R, Hemmat M, Wang BT, El Naggar M, Kopita KA, Randolph LM, Chase JM, Matas Aguilera MJ, Siles JL, Church JA, Hauser N, Shen JJ, Jones MC, Wierenga KJ, Jiang Z, Haddadin M, Boyar FZ, Anguiano A, Strom CM, Sahoo T. Regions of homozygosity identified by oligonucleotide SNP arrays: evaluating the incidence and clinical utility. Eur J Hum Genet. 2015;23(5):663–71. doi:https://doi.org/10.1038/ejhg.2014.153
Chaves TF, Oliveira LF, Ocampos M, Barbato IT, de Luca GR, Barbato-Filho JH, de Camargo-Pinto LL, Bernardi P, Maris AF. Long contiguous stretches of homozygosity detected by chromosomal microarrays (CMA) in patients with neurodevelopmental disorders in the South of Brazil. BMC Med Genomics. 2019;12:50. https://doi.org/10.1186/s12920-019-0496-5.
Liehr T. Introduction. In: Small supernumerary marker chromosomes (sSMC). Berlin: Springer; 2012. https://doi.org/10.1007/978-3-642-20766-2_1
Kearney HM, Kearney JB, Conlin LK. Diagnostic implications of excessive homozygosity detected by SNP-based microarrays: consanguinity, uniparental disomy, and recessive single-gene mutations. Clin Lab Med. 2011;31(4):595–613. https://doi.org/10.1016/j.cll.2011.08.003.
Li J, Yang T, Wang L, Yan H, Zhang Y, Guo Y, Pan F, Zhang Z, Peng Y, Zhou Q, He L, Zhu X, Deng H, Levy S, Papasian CJ, Drees BM, Hamilton JJ, Recker RR, Cheng J, Deng HW. Whole genome distribution and ethnic differentiation of copy number variation in caucasian and asian populations. PLoS ONE. 2009;4(11):7958. https://doi.org/10.1371/journal.pone.0007958.
Pinto D, Darvishi K, Shi X, Rajan D, Rigler D, Fitzgerald T, Lionel AC, Thiruvahindrapuram B, Macdonald JR, Mills R, Prasad A, Noonan K, Gribble S, Prigmore E, Donahoe PK, Smith RS, Park JH, Hurles ME, Carter NP, Lee C, Scherer SW, Feuk L. Comprehensive assessment of array-based platforms and calling algorithms for detection of copy number variants. Nat Biotechnol. 2011;29(6):512–20. https://doi.org/10.1038/nbt.1852.
Purwar N, Tiwari P, Mathur N, Sharma H, Sahlot R, Garg U, Sharma B, Saxena A, Mathur SK. Higher CNV frequencies in chromosome 14 of girls with turner syndrome phenotype. J Clin Endocrinol Metab. 2021;106(12):e4935–55. https://doi.org/10.1210/clinem/dgab572.
Ambele MA, van Zyl A, Pepper MS, van Heerden MB, van Heerden WFP. Amplification of 3q26.2, 5q14.3, 8q24.3, 8q22.3, and 14q32.33 Are Possible Common Genetic Alterations in Oral Cancer Patients. Front Oncol. 2020;10:683. https://doi.org/10.3389/fonc.2020.00683
Wijesiriwardhana P, Wettasinghe K, Dissanayeke VHW. Copy Number Variants Captured by the Array Comparative Genomic Hybridization in a Cohort of Patients Affected with Hereditary Colorectal Cancer in Sri Lanka: The First CNV Analysis Study of the Hereditary Colorectal Cancer in the Sri Lankan Population. Asian Pac J Cancer Prev. 2021;22(6):1957–1966. https://doi.org/10.31557/APJCP.2021.22.6.1957.
Diniz MG, França JA, Vilas-Boas FAS, de Souza FTA, Calin GA, Gomez RS, Sousa SF, Gomes CC. The long noncoding RNA KIAA0125 is upregulated in ameloblastomas. Pathol Res Pract. 2019;215(3):466–9. https://doi.org/10.1016/j.prp.2018.12.030.
Godoy VCSM, Teixeira BF, Colovati M, Rodrigues OJH, Moysés- Oliveira M, Melaragno MI. Copy number variation (CNV) identification, interpretation, and database from Brazilian patients. Genet Mol Biol. 2020;43(4):e20190218. https://doi.org/10.1590/1678-4685-GMB-2019-0218
GSE87048 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE87048Rebollar-Vega R, Romero-Córdoba S, Rodríguez-Cuevas S, Maffuz-Aziz A, Hidalgo-Miranda A. Analysis of somatic DNA copy number alterations and frequency of breast cancer intrinsic subtypes from Mexican women [copy number]. Contributor(s) Rebollar-Vega R, Romero-Córdoba S, Rodríguez-Cuevas S, Maffuz-Aziz A, Hidalgo-Miranda A Citation(s) Romero-Cordoba SL, Salido-Guadarrama I, Rebollar-Vega R, Bautista-Piña V, Dominguez-Reyes C, Tenorio-Torres A, Villegas-Carlos F, Fernández-López JC, Uribe-Figueroa L, Alfaro-Ruiz L, Hidalgo-Miranda A. Comprehensive omic characterization of breast cancer in Mexican-Hispanic women. Nat Commun. 2021;12(1):2245. https://doi.org/10.1038/s41467-021-22478-5
GSE67047 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE67047. Genome-wide association and epistatic studies in sporadic medullary and juvenile papillary thyroid carcinomas. Bleda M, Garcia-Alonso L, Medina I, Gonzalez CY, Fernandez RM, Nuñez-Torres R, Luzon-Toro B, Torroglosa A, Marba M, del Valle Enguix-Riego M, Montaner D, Antiñolo G, Borrego S, Dopazo J. Citation: Luzón-Toro B, Bleda M, Navarro E, García-Alonso L, Ruiz-Ferrer M, Medina I, Martín-Sánchez M, Gonzalez CY, Fernández RM, Torroglosa A, Antiñolo G, Dopazo J, Borrego S. Identification of epistatic interactions through genome-wide association studies in sporadic medullary and juvenile papillary thyroid carcinomas. BMC Med Genomics. 2015;8:83. https://doi.org/10.1186/s12920-015-0160-7
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
Bengtsson H, Wirapati P, Speed TP. A single-array preprocessing method for estimating full-resolution raw copy numbers from all Affymetrix genotyping arrays including GenomeWideSNP 5 & 6. Bioinformatics. 2009;25(17):2149–56. https://doi.org/10.1093/bioinformatics/btp371.
We would like to thank all the individuals for their consent to participate in this study. We would like to acknowledge Benilde García de Teresa for english revision, Fabiola Morales and Raúl Mojica, Unidad de Genotipificación y Análisis de la Expresión, Instituto Nacional de Medicina Genómica, for help with microarray performance and Jorge Yañez and Paul Gaytan from the Unidad de Síntesis y Secuenciación, Instituto de Biotecnologia, UNAM, for primer synthesis.
CONACyT-FONCICYT 95419, CONACyT-FOSISS 142040 and Federal Resources from the Instituto Nacional de Pediatría 084/2010 and 2020/043. SS had a scholar fellowship CONACyT 100254.
Ethics approval and consent to participate
This study was performed in accordance with the Declaration of Helsinki and was carried out with protocols approved by the Institutional Ethics and Research Committees from the Instituto Nacional de Pediatría, project numbers 084/2010 and 2020/043. Informed consent was obtained from all subjects involved in the study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Sánchez, S., Juárez, U., Domínguez, J. et al. Frequent copy number variants in a cohort of Mexican-Mestizo individuals. Mol Cytogenet 16, 2 (2023). https://doi.org/10.1186/s13039-022-00631-z