Response to Wyckelsma et al.: Loss of α-actinin-3 during human evolution provides superior cold resilience and muscle heat generation.

The common loss-of-function mutation R577X in the structural muscle protein ACTN3 emerged as a potential target of positive selection from early studies and has been the focus of insightful physiological work suggesting a significant impact on muscle metabolism. Adaptation to cold climates has been proposed as a key adaptive mechanism explaining its global allele frequency patterns. Here, we re-examine this hypothesis analyzing modern (n = 3,626) and ancient (n = 1,651) genomic data by using allele-frequency as well as haplotype homozygosity-based methods. The presented results are more consistent with genetic drift rather than selection in cold climates as the main driver of the ACTN3 R577X frequency distribution in human populations across the world. This Matters Arising paper is in response to Wyckelsma et al. (2021),1 published in The American Journal of Human Genetics. See also the response by Wyckelsma et al. (2022),2 published in this issue.

Mitogenomics of the Koryaks and Evens of the northern coast of the Sea of Okhotsk.

Due to the geographical proximity of the northern coast of the Sea of Okhotsk and Kamchatka Peninsula to the Beringia, the indigenous populations of these territories are of great interest for elucidating the human settlement history of northern Asia and America. Meanwhile, there is a clear shortage of genetic studies of the indigenous populations of the northern coast of the Sea of Okhotsk. Here, in order to examine their fine-scale matrilineal genetic structure, ancestry and relationships with neighboring populations, we analyzed 203 complete mitogenomes (174 of which are new) from population samples of the Koryaks and Evens of the northern coast of the Sea of Okhotsk and the Chukchi of the extreme northeast Asia. The patterns observed underscore the reduced level of genetic diversity found in the Koryak, Even, and Chukchi populations, which, along with the high degree of interpopulation differentiation, may be the result of genetic drift. Our phylogeographic analysis reveals common Paleo-Asiatic ancestry for 51.1% of the Koryaks and 17.8% of the Evens. About third of the mitogenomes found in the Koryaks and Evens might be considered as ethno-specific, as these are virtually absent elsewhere in North, Central and East Asia. Coalescence ages of most of these lineages coincide well with the emergence and development of the Tokarev and Old Koryak archaeological cultures associated with the formation of the Koryaks, as well as with the period of separation and split of the North Tungusic groups migrated northwards from the Lake Baikal or the Amur River area.

A Machine-Learning-Based Approach to Prediction of Biogeographic Ancestry within Europe.

Data obtained with the use of massive parallel sequencing (MPS) can be valuable in population genetics studies. In particular, such data harbor the potential for distinguishing samples from different populations, especially from those coming from adjacent populations of common origin. Machine learning (ML) techniques seem to be especially well suited for analyzing large datasets obtained using MPS. The Slavic populations constitute about a third of the population of Europe and inhabit a large area of the continent, while being relatively closely related in population genetics terms. In this proof-of-concept study, various ML techniques were used to classify DNA samples from Slavic and non-Slavic individuals. The primary objective of this study was to empirically evaluate the feasibility of discerning the genetic provenance of individuals of Slavic descent who exhibit genetic similarity, with the overarching goal of categorizing DNA specimens derived from diverse Slavic population representatives. Raw sequencing data were pre-processed, to obtain a 1200 character-long binary vector. A total of three classifiers were used-Random Forest, Support Vector Machine (SVM), and XGBoost. The most-promising results were obtained using SVM with a linear kernel, with 99.9% accuracy and F1-scores of 0.9846-1.000 for all classes.

Mitogenomics of modern Mongolic-speaking populations.

Here, we present a comprehensive data set of 489 complete mitogenomes (211 of which are new) from four Mongolic-speaking populations (Mongols, Barghuts, Khamnigans, and Buryats) to investigate their matrilineal genetic structure, ancestry and relationship with other ethnic groups. We show that along with very high levels of genetic diversity and lack of genetic differentiation, Mongolic-speaking populations exhibit strong genetic resemblance to East Asian populations of Chinese, Japanese, and Uyghurs. Phylogeographic analysis of complete mitogenomes reveals the presence of different components in the gene pools of modern Mongolic-speaking populations-the main East Eurasian component is represented by mtDNA lineages of East Asian, Siberian and autochthonous (the Baikal region/Mongolian) ancestry, whereas the less pronounced West Eurasian component can be ascribed to Europe and West Asia/Caucasus. We also observed that up to one third of the mtDNA subhaplogroups identified in Mongolic-speaking populations can be considered as Mongolic-specific with the coalescence age of most of them not exceeding 1.7 kya. This coincides well with the population size growth which started around 1.1 kya and is detectable only in the Bayesian Skyline Plot constructed based on Mongolic-specific mitogenomes. Our data suggest that the genetic structure established during the Mongol empire is still retained in present-day Mongolic-speaking populations.

Complete mitogenome data for the Serbian population: the contribution to high-quality forensic databases.

Mitochondrial genome (mtDNA) is a valuable resource in resolving various human forensic casework. The usage of variability of complete mtDNA genomes increases their discriminatory power to the maximum and enables ultimate resolution of distinct maternal lineages. However, their wider employment in forensic casework is nowadays limited by the lack of appropriate reference database. In order to fill in the gap in the reference data, which, considering Slavic-speaking populations, currently comprises only mitogenomes of East and West Slavs, we present mitogenome data for 226 Serbians, representatives of South Slavs from the Balkan Peninsula. We found 143 (sub)haplogroups among which West Eurasian ones were dominant. The percentage of unique haplotypes was 85%, and the random match probability was as low as 0.53%. We support previous findings on both high levels of genetic diversity in the Serbian population and patterns of genetic differentiation among this and ten studied European populations. However, our high-resolution data supported more pronounced genetic differentiation among Serbians and two Slavic populations (Russians and Poles) as well as expansion of the Serbian population after the Last Glacial Maximum and during the Migration period (fourth to ninth century A.D.), as inferred from the Bayesian skyline analysis. Phylogenetic analysis of haplotypes found in Serbians contributed towards the improvement of the worldwide mtDNA phylogeny, which is essential for the interpretation of the mtDNA casework.

Insights into matrilineal genetic structure, differentiation and ancestry of Armenians based on complete mitogenome data.

Distinctive peculiarities of Armenians such as their millennia-long genetic isolation and strong national identity attract a keen interest while studying the demographic history of the West Asia. Here, to examine their fine-scale matrilineal genetic structure, ancestry and relationships with neighboring populations, we analyzed 536 complete mitogenomes (141 of which are novel) from 8 geographically different Armenian populations, covering the whole stretch of historical Armenia. The observed patterns highlight a remarkable degree of matrilineal genetic heterogeneity and weak population structuring of Armenians. Moreover, our phylogeographic analysis reveals common ancestries for some mtDNA lineages shared by West Asians, Transcaucasians, Europeans, Central Asians and Armenians. About third of the mtDNA subhaplogroups found in Armenian gene pool might be considered as Armenian-specific, as these are virtually absent elsewhere in Europe, West Asia and Transcaucasia. Coalescence ages of most of these lineages do not exceed 3.1 kya and coincide well with the population size growth started around 1.8-2.8 kya detectable only in the Bayesian Skyline Plots based on the Armenian-specific mtDNA haplotypes.

Whole mitochondrial genome diversity in two Hungarian populations.

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

Mitogenomic diversity and differentiation of the Buryats.

In this paper we present a results of first comprehensive study of the complete mitogenomes in the Buryats with regard to their belonging to the main regional (eastern and western Buryats); tribal (Khori, Ekhirid, Bulagad, and Khongodor), and ethno-territorial (Aginsk, Alar, Balagansk, Barguzin, Ida, Khorinsk, Kuda, Selenga, Verkholensk, Olkhon, Tunka, and Shenehen Buryats) groups. The analysis of molecular variation performed using regional, tribal, and ethno-territorial divisions of the Buryats showed lack of genetic differentiation at all levels. Nonetheless, the complete mitogenome analysis revealed a very high level of genetic diversity in the Buryats which is the highest among Siberian populations and comparable to that in populations of eastern and western Asia. The AMOVA and MDS analyses results imply to a strong genetic similarity between the Buryats and eastern Asian populations of Chinese and Japanese, suggesting their origin on the basis of common maternal ancestry components. Several new Buryat-specific branches of haplogroup G (G2a2a, G2a1i, G2a5a) display signals of dispersals dating to 2.6-6.6 kya with a possible origin in eastern Asia, thus testifying Bronze Age and Neolithic arrival of ancestral eastern Asian component to the South Siberia region.

The genetic legacy of the expansion of Turkic-speaking nomads across Eurasia.

The Turkic peoples represent a diverse collection of ethnic groups defined by the Turkic languages. These groups have dispersed across a vast area, including Siberia, Northwest China, Central Asia, East Europe, the Caucasus, Anatolia, the Middle East, and Afghanistan. The origin and early dispersal history of the Turkic peoples is disputed, with candidates for their ancient homeland ranging from the Transcaspian steppe to Manchuria in Northeast Asia. Previous genetic studies have not identified a clear-cut unifying genetic signal for the Turkic peoples, which lends support for language replacement rather than demic diffusion as the model for the Turkic language's expansion. We addressed the genetic origin of 373 individuals from 22 Turkic-speaking populations, representing their current geographic range, by analyzing genome-wide high-density genotype data. In agreement with the elite dominance model of language expansion most of the Turkic peoples studied genetically resemble their geographic neighbors. However, western Turkic peoples sampled across West Eurasia shared an excess of long chromosomal tracts that are identical by descent (IBD) with populations from present-day South Siberia and Mongolia (SSM), an area where historians center a series of early Turkic and non-Turkic steppe polities. While SSM matching IBD tracts (> 1cM) are also observed in non-Turkic populations, Turkic peoples demonstrate a higher percentage of such tracts (p-values ≤ 0.01) compared to their non-Turkic neighbors. Finally, we used the ALDER method and inferred admixture dates (~9th-17th centuries) that overlap with the Turkic migrations of the 5th-16th centuries. Thus, our results indicate historical admixture among Turkic peoples, and the recent shared ancestry with modern populations in SSM supports one of the hypothesized homelands for their nomadic Turkic and related Mongolic ancestors.

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations.

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6-23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.

Mitochondrial super-haplogroup U diversity in Serbians.

BACKGROUND: Available mitochondrial (mtDNA) data demonstrate genetic differentiation among South Slavs inhabiting the Balkan Peninsula. However, their resolution is insufficient to elucidate the female-specific aspects of the genetic history of South Slavs, including the genetic impact of various migrations which were rather common within the Balkans, a region having a turbulent demographic history. AIM: The aim was to thoroughly study complete mitogenomes of Serbians, a population linking westward and eastward South Slavs. SUBJECTS AND METHODS: Forty-six predominantly Serbian super-haplogroup U complete mitogenomes were analysed phylogenetically against ∼4000 available complete mtDNAs of modern and ancient Western Eurasians. RESULTS: Serbians share a number of U mtDNA lineages with Southern, Eastern-Central and North-Western Europeans. Putative Balkan-specific lineages (e.g. U1a1c2, U4c1b1, U5b3j, K1a4l and K1a13a1) and lineages shared among Serbians (South Slavs) and West and East Slavs were detected (e.g. U2e1b1, U2e2a1d, U4a2a, U4a2c, U4a2g1, U4d2b and U5b1a1). CONCLUSION: The exceptional diversity of maternal lineages found in Serbians may be associated with the genetic impact of both autochthonous pre-Slavic Balkan populations whose mtDNA gene pool was affected by migrations of various populations over time (e.g. Bronze Age pastoralists) and Slavic and Germanic newcomers in the early Middle Ages.

Y chromosome haplotype diversity in Mongolic-speaking populations and gene conversion at the duplicated STR DYS385a,b in haplogroup C3-M407.

Y chromosome microsatellite (Y-STR) diversity has been studied in different Mongolic-speaking populations from South Siberia, Mongolia, North-East China and East Europe. The results obtained indicate that the Mongolic-speaking populations clustered into two groups, with one group including populations from eastern part of South Siberia and Central Asia (the Buryats, Barghuts and Khamnigans) and the other group including populations from western part of Central Asia and East Europe (the Mongols and Kalmyks). High frequency of haplogroup C3-M407 (>50%) is present in the Buryats, Barghuts and Khamnigans, whereas in the Mongols and Kalmyks its frequency is much lower. In addition, two allelic combinations in DYS385a,b loci of C3-M407 haplotypes have been observed: the combination 11,18 (as well as 11,17 and 11,19) is frequent in different Mongolic-speaking populations, but the 11,11 branch is present mainly in the Kalmyks and Mongols. Results of locus-specific sequencing suggest that the action of gene conversion is a more likely explanation for origin of homoallelic 11,11 combination. Moreover, analysis of median networks of Y-STR haplotypes demonstrates that at least two gene conversion events can be revealed-one of them has probably occurred among the Mongols, and the other event occurred in the Barghuts. These two events give an average gene conversion rate range of 0.24-7.1 × 10(-3) per generation.

Mitochondrial DNA perspective of Serbian genetic diversity.

Although south-Slavic populations have been studied to date from various aspects, the population of Serbia, occupying the central part of the Balkan Peninsula, is still genetically understudied at least at the level of mitochondrial DNA (mtDNA) variation. We analyzed polymorphisms of the first and the second mtDNA hypervariable segments (HVS-I and HVS-II) and informative coding-region markers in 139 Serbians to shed more light on their mtDNA variability, and used available data on other Slavic and neighboring non-Slavic populations to assess their interrelations in a broader European context. The contemporary Serbian mtDNA profile is consistent with the general European maternal landscape having a substantial proportion of shared haplotypes with eastern, central, and southern European populations. Serbian population was characterized as an important link between easternmost and westernmost south-Slavic populations due to the observed lack of genetic differentiation with all other south-Slavic populations and its geographical positioning within the Balkan Peninsula. An increased heterogeneity of south Slavs, most likely mirroring turbulent demographic events within the Balkan Peninsula over time (i.e., frequent admixture and differential introgression of various gene pools), and a marked geographical stratification of Slavs to south-, east-, and west-Slavic groups, were also found. A phylogeographic analyses of 20 completely sequenced Serbian mitochondrial genomes revealed not only the presence of mtDNA lineages predominantly found within the Slavic gene pool (U4a2a*, U4a2a1, U4a2c, U4a2g, HV10), supporting a common Slavic origin, but also lineages that may have originated within the southern Europe (H5*, H5e1, H5a1v) and the Balkan Peninsula in particular (H6a2b and L2a1k).

Western Eurasian ancestry in modern Siberians based on mitogenomic data.

BACKGROUND: Although the genetic heritage of aboriginal Siberians is mostly of eastern Asian ancestry, a substantial western Eurasian component is observed in the majority of northern Asian populations. Traces of at least two migrations into southern Siberia, one from eastern Europe and the other from western Asia/the Caucasus have been detected previously in mitochondrial gene pools of modern Siberians. RESULTS: We report here 166 new complete mitochondrial DNA (mtDNA) sequences that allow us to expand and re-analyze the available data sets of western Eurasian lineages found in northern Asian populations, define the phylogenetic status of Siberian-specific subclades and search for links between mtDNA haplotypes/subclades and events of human migrations. From a survey of 158 western Eurasian mtDNA genomes found in Siberia we estimate that nearly 40% of them most likely have western Asian and another 29% European ancestry. It is striking that 65 of northern Asian mitogenomes, i.e. ~41%, fall into 19 branches and subclades which can be considered as Siberian-specific being found so far only in Siberian populations. From the coalescence analysis it is evident that the sequence divergence of Siberian-specific subclades was relatively small, corresponding to only 0.6-9.5 kya (using the complete mtDNA rate) and 1-6 kya (coding region rate). CONCLUSIONS: The phylogeographic analysis implies that the western Eurasian founders, giving rise to Siberian specific subclades, may trace their ancestry only to the early and mid-Holocene, though some of genetic lineages may trace their ancestry back to the end of Last Glacial Maximum (LGM). We have not found the modern northern Asians to have western Eurasian genetic components of sufficient antiquity to indicate traces of pre-LGM expansions.

Y-chromosome diversity in the Kalmyks at the ethnical and tribal levels.

The Mongolic-speaking Kalmyks currently inhabiting the steppes of the Volga region have Central Asian ancestry and are organized into the tribal groups. The genetic relationships among these tribes and their origin have remained obscure. We analyzed 17 short tandem repeat and 44 binary polymorphisms of Y-chromosome in 426 individuals mainly from three major tribes of the Kalmyks (the Torguuds, Dörwöds and Khoshuuds). Among these tribes, the Dörwöds and Torguuds, as well as the Kalmyks collectively as an ethnic group, showed relatively close genetic affinities to each other and to the Mongols and Altaian Kazakhs, whereas the Khoshuuds were clearly separated from all of them, gathering with the Manchu, Tibetans or Evenks (depending on the algorithm used to calculate genetic distances). The genetic results also indicate that paternal gene flow from East Europeans to the Kalmyks is very little, despite their cohabitation in the North Caspian Steppe during the last 380 years. The occurrence of unique cluster of N1c-Tat haplotypes in the Khoshuuds, which dates to about 340 years and is likely to have East European ancestry, is considered as a result of interethnic contacts occurred soon after the appearance of the Kalmyk tribes in the Volga-Ural region.

Phylogeny and genetic history of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870) inferred from complete mitochondrial genomes.

We assessed phylogeny of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870), the most northern ectothermic, terrestrial vertebrate in Eurasia, by sequence analysis of complete mitochondrial genomes in 26 specimens from different localities (China, Khabarovsk region, Sakhalin, Yakutia, Magadan region, Chukotka, Kamchatka, Ural, European part of Russia). In addition, a complete mitochondrial genome of the Schrenck salamander, Salamandrella schrenckii, was determined for the first time. Bayesian phylogenetic analysis of the entire mtDNA genomes of S. keyserlingii demonstrates that two haplotype clades, AB and C, radiated about 1.4 million years ago (Mya). Bayesian skyline plots of population size change through time show an expansion around 250 thousand years ago (kya) and then a decline around the Last Glacial Maximum (25 kya) with subsequent restoration of population size. Climatic changes during the Quaternary period have dramatically affected the population genetic structure of the Siberian salamanders. In addition, complete mtDNA sequence analysis allowed us to recognize that the vast area of Northern Eurasia was colonized only by the Siberian salamander clade C1b during the last 150 kya. Meanwhile, we were unable to find evidence of molecular adaptation in this clade by analyzing the whole mitochondrial genomes of the Siberian salamanders.

Y-chromosome variation in Tajiks and Iranians.

AIM: The purpose of this study was to characterize Y-chromosome diversity in Tajiks from Tajikistan and in Persians and Kurds from Iran. METHOD: Y-chromosome haplotypes were identified in 40 Tajiks, 77 Persians and 25 Kurds, using 12 short tandem repeats (STR) and 18 binary markers. RESULTS: High genetic diversity was observed in the populations studied. Six of 12 haplogroups were common in Persians, Kurds and Tajiks, but only three haplogroups (G-M201, J-12f2 and L-M20) were the most frequent in all populations, comprising together ~60% of the Y-chromosomes in the pooled data set. Analysis of genetic distances between Y-STR haplotypes revealed that the Kurds showed a great distance to the Iranian-speaking populations of Iran, Afghanistan and Tajikistan. The presence of Indian-specific haplogroups L-M20, H1-M52 and R2a-M124 in both Tajik samples from Afghanistan and Tajikistan demonstrates an apparent genetic affinity between Tajiks from these two regions. CONCLUSIONS: Despite the marked similarities between Y-chromosome gene pools of Iranian-speaking populations, there are differences between them, defined by many factors, including geographic and linguistic relationships.

On the Y-chromosome haplogroup C3c classification.

As there are ambiguities in classification of the Y-chromosome haplogroup C3c, relatively frequent in populations of Northern Asia, we analyzed all three haplogroup-defining markers M48, M77 and M86 in C3-M217-individuals from Siberia, Eastern Asia and Eastern Europe. We have found that haplogroup C3c is characterized by the derived state at M48, whereas mutations at both M77 and M86 define subhaplogroup C3c1. The branch defined by M48 alone would belong to subhaplogroup C3c*, characteristic for some populations of Central and Eastern Siberia, such as Koryaks, Evens, Evenks and Yukaghirs. Subhaplogroup C3c* individuals could be considered as remnants of the Neolithic population of Siberia, based on the age of C3c*-short tandem repeat variation amounting to 4.5 ± 2.4 thousand years.

The Y-chromosome C3* star-cluster attributed to Genghis Khan's descendants is present at high frequency in the Kerey clan from Kazakhstan.

To verify the possibility that the Y-chromosome C3* star-cluster attributed to Genghis Khan and his patrilineal descendants is relatively frequent in the Kereys, who are the dominant clan in Kazakhstan and in Central Asia as a whole, polymorphism of the Y-chromosome was studied in Kazakhs, represented mostly by members of the Kerey clan. The Kereys showed the highest frequency (76.5%) of individuals carrying the Y-chromosome variant known as C3* star-cluster ascribed to the descendants of Genghis Khan. C3* star-cluster haplotypes were found in two subclans, Abakh-Kereys and Ashmaily-Kereys, diverged about 20-22 generations ago according to the historical data. Median network of the Kerey star-cluster haplotypes at 17 STR loci displays a bipartite structure, with two subclusters defined by the only difference at the DYS448 locus. Noteworthy is a strong correspondence of these subclusters with the Kerey subclans affiliation. The data obtained suggest that the Kerey clan appears to be the largest known clan in the world descending from a common Y-chromosome ancestor. Possible ways of Genghis Khan's relationship to the Kereys are discussed.

An association between copy number variation of enhancer involved in craniofacial development and biogeographic ancestry.

Human facial morphology is a combination of many complex traits and is determined by a large number of genes and enhancers. Here, we report a Copy Number Variation (CNV) study of enhancer hs1431 in populations of Central European and South Siberian ancestry. Central European samples included 97 Poles, while South Siberian samples included 78 Buryats and 27 Tuvinians. CNVs were detected by real-time PCR, using ViiA™ 7 Real-Time PCR System (Applied Biosystems). We revealed significant differences in CNV of hs1431 enhancer between Polish and Buryat population (p=0.0378), but not between Central European and South Siberian population (p=0.1225). Our results suggest that an increase in copy number variation of hs1431 enhancer is associated with biogeographic ancestry. However, this result needs extending and replicating in larger cohorts. This is the first study revealing the presence of copy number variation of enhancer hs1431 in humans.