Evolutionary history of the mtDNA 9-bp deletion in Chinese populations and its relevance to the peopling of east and southeast Asia

In total, 1218 Chinese from twelve ethnic groups and nine Han geographic groups were screened for the mtDNA 9-bp deletion motif. The frequency of the 9-bp deletion in all samples was 14.7% but ranged from 0% to 32% in the various ethnic groups. Three individuals had a triplication of the 9-bp segment. Phylogenetic and demographic analyses of the mtDNA hypervariable segment 1 (HVS 1) sequences suggest that the 9-bp deletion occurred more than once in China. The majority of the Chinese deletion haplotypes (about 90%) have a common origin as a mutational event following an initial expansion of modern humans in eastern Asia. Other deletion haplotypes and the three haplotypes with a 9-bp triplication may have arisen independently in the Chinese, presumably by replication error. HVS1 haplotype analysis suggests two possible migration routes of the 9-bp deletion in east and southeast Asia. Both migrations originated in China with one route leading to the Pacific Islands via Taiwan, the other to southeast Asia and possibly the Nicobar Islands. Along both routes of peopling, a decrease in HVSI diversity of the mtDNA haplotypes is observed. The "Polynesian motif (16217T/C, 16247A/G, and 16261C/T)" and the 16140T/C, 16266C/A, or C/G polymorphisms appear specific to each migration route.     

Different population histories of the Mundari- and Mon-Khmer-speaking Austro-Asiatic tribes inferred from the mtDNA 9-bp deletion/insertion polymorphism in Indian populations

Length variation in the human mtDNA intergenic region between the cytochrome oxidase II (COII) and tRNA lysine (tRNA^lys) genes has been widely studied in world populations. Specifically, Austronesian populations of the Pacific and Austro-Asiatic populations of southeast Asia most frequently carry the 9-bp deletion in that region implying their shared common ancestry in haplogroup B. Furthermore, multiple independent origins of the 9-bp deletion at the background of other mtDNA haplogroups has been shown in populations of Africa, Europe, Australia, and India. We have analyzed 3293 Indian individuals belonging to 58 populations, representing different caste, tribal, and religious groups, for the length variation in the 9-bp motif. The 9-bp deletion (one copy) and insertion (three copies) alleles were observed in 2.51% (2.15% deletion and 0.36% insertion) of the individuals. The maximum frequency of the deletion (45.8%) was observed in the Nicobarese in association with the haplogroup B5a D-loop motif that is common throughout southeast Asia. The low polymorphism in the D-loop sequence of the Nicobarese B5a samples suggests their recent origin and a founder effect, probably involving migration from southeast Asia. Interestingly, none of the 302 (except one Munda sample, which has 9-bp insertion) from Mundari-speaking Austro-Asiatic populations from the Indian mainland showed the length polymorphism of the 9-bp motif, pointing either to their independent origin from the Mon-Khmeric-speaking Nicobarese or to an extensive admixture with neighboring Indo-European-speaking populations. Consistent with previous reports, the Indo-European and Dravidic populations of India showed low frequency of the 9-bp deletion/insertion. More than 18 independent origins of the deletion or insertion mutation could be inferred in the phylogenetic analysis of the D-loop sequences.     

The world-wide distribution of allele frequencies at the human dopamine D4 receptor locus

The dopamine D4 receptor gene (DRD4) has an expressed polymorphism in the third exon that may have functional relevance. The polymorphism exists at two levels. At the higher level there is an imperfect tandem repeat of 48 base pairs (bp) coding for 16 amino acids; alleles have been identified with 2 (32 amino acids) to 10 (160 amino acids) repeats. The imperfect nature of the repeats is responsible for a more subtle level of variation since alleles with the same number of repeats can differ in the exact sequences or in the order of the variants of the 48-bp unit. We have undertaken a global survey of this expressed polymorphism as one approach to understanding the evolutionary significance and origins of the polymorphism as well as understanding what selective forces, if any, may be operating at this locus. As the first step, we have determined the repeat number genotype of the DRD4 repeat polymorphism in 1,327 individuals from 36 different populations. The allele frequencies differ considerably among the different populations. The 4-repeat allele was the most prevalent (global mean allele frequency = 64.3%) and appeared in every population with a frequency ranging from 0.16 to 0.96. The 7-repeat allele was the second most common (global mean = 20.6%), appearing quite frequently in the Americas (mean frequency = 48.3%) but only occasionally in East and South Asia (mean frequency = 1.9%). The 2-repeat allele was the third most common (global mean frequency = 8.2%) and was quite frequent in East and South Asia (mean frequency = 18.1%) while uncommon in the Americas (mean frequency = 2.9%) and Africa (mean frequency = 1.7%). The universality of the polymorphism with only three common repeat-number alleles (4, 7, and 2) indicates that the polymorphism is ancient and arose before the global dispersion of modern humans. The diversity of actual allele frequencies for this expressed polymorphism among different populations emphasizes the importance of population considerations in the design and interpretation of any association studies carried out with this polymorphism. Received: 18 July 1995 / Revised: 18 December 1995     

Mitochondrial sequence variation in the Guahibo Amerindian population from Venezuela

New data were obtained on mitochondrial DNA (mtDNA) from Guahibo from Venezuela, a group so far not studied using molecular data. A population sample (n = 59) was analyzed for mtDNA variation in two control-region hypervariable segments (HV1 and HV2) by sequencing. The presence or absence of a 9-bp polymorphism in the COII/tRNA(Lys) region was studied by direct amplification and electrophoretic identification. Thirty-eight variable sites were detected in regions HV1 and HV2, defining 26 mtDNA lineages; 23.7% of these were present in a single individual. The 9-bp deletion was found in 3.39% of individuals. Nucleotide and haplotype diversities were relatively high compared with other New World populations. The identified sequence haplotypes were classified into four major haplogroups (A-D) according to previous studies, with high frequencies for A (47.46%) and C (49.15%), low frequency for B (3.39%), and an absence of D.     

Bovine spongiform encephalopathy associated insertion/deletion polymorphisms of the prion protein gene in the four beef cattle breeds from North China

Two insertion/deletion (indel) polymorphisms of the prion protein gene (PRNP), a 23-bp indel in the putative promoter region and a 12-bp indel within intron I, are associated with the susceptibility to bovine spongiform encephalopathy (BSE) in cattle. In the present study, the polymorphism frequencies of the two indels in four main beef cattle breeds (Hereford, Simmental, Black Angus, and Mongolian) from North China were studied. The results showed that the frequencies of deletion genotypes and alleles of 23- and 12-bp indels were lower, whereas the frequencies of insertion genotypes and alleles of the two indels were higher in Mongolian cattle than in the other three cattle breeds. In Mongolian cattle, the 23-bp insertion / 12-bp insertion was the major haplotype, whereas in Hereford, Simmental, and Black Angus cattle, the 23-bp deletion / 12-bp deletion was the major haplotype. These results demonstrated that Mongolian cattle could be more resistant to BSE, compared with the other three cattle breeds, because of its relatively low frequencies of deletion genotypes and alleles of 23- and 12-bp indel polymorphisms. Thus, this race could be important for selective breeding to improve resistance against BSE in this area.     

贵州瑶族3支系Y-DNA及线粒体DNA序列多态性分析

采用PCR－RFLP技术,通过观察由12个单核苷酸多态位点(SNPs)组成的Y染色体单倍型及由9个多态位点组成的线粒体DNA单倍型在贵州瑶族中的分布,分析贵州瑶族父系及母系遗传结构,探讨其起源及迁徙。结果显示,97份男性样本分别属于H7、H8、H9、H11 4种Y-DNA单倍型,苗瑶语系特异Y-DNA单倍型H7的平均频率为92.4%;通过对线粒体DNA基因分型,得到8种单倍型,可归入B4、B5、D4、D5和N*单倍型类群中,CoⅡ/tRNA^Lys区域间的9bp缺失平均频率为58.2%。结果提示贵州瑶族父系遗传结构单一,具有典型的苗瑶族群特征,又存在与其他族群的融合。母系遗传结构相对复杂,9 bp缺失是贵州瑶族的母系遗传结构特征。      

Structure and diversity of the mitochondrial gene pool of the aboriginal population of Tuva and Buriatia from restriction polymorphism data

The populations of Tuvinians (N = 36) and Buryats (N = 105) were characterized by using the data on mitochondrial DNA (mtDNA) polymorphism. The gene pools of both ethnic groups possessed the mtDNA types belonging to the four main haplogroups, A, B, C, and D, found only in the indigenous populations of Asia and America. The total frequencies of the A, B, C, and D haplogroups in Tuvinians and Buryats were 72.3% and 52.4%, respectively. These values, along with the frequency for Altai populations (57.2%), were highest in the Asian populations studied, indicating that the populations Southern and Eastern Siberia can be considered as ancestral relatives to the ethnic groups of the New World. Analysis of the mtDNA region V polymorphism showed the presence of 9-bp deletion and 4-bp insertion in both populations with frequencies respectively of 13.9 and 5.56% in Tuvinians and 4.8 and 1.9% in Buryats. The frequency of the +AvaII/8249 variant was 11.1% in Tuvinians and 3.81% in Buryats. Analysis of the association between the region V deletion-insertion polymorphism and certain restriction haplogroups pointed to repeated and independent emergence of the 4-bp insertion in Siberia.     

Role of mtDNA haplogroups in COPD susceptibility in a southwestern Han Chinese population

The interplay of a complex genetic basis with the environmental factors of chronic obstructive pulmonary disease (COPD) may account for the differences in individual susceptibility to COPD. Mitochondrial DNA (mtDNA) contributes to an individual's ability to resist oxidation, an important determinant that affects COPD susceptibility. To investigate whether mtDNA haplogroups play important roles in COPD susceptibility, the frequencies of mtDNA haplogroups and an 822-bp mtDNA deletion in 671 COPD patients and 724 control individuals from southwestern China were compared. Multivariate logistic regression analysis revealed that, whereas mtDNA haplogroups A and M7 might be associated with an increased risk for COPD (OR=1.996, 95% CI=1.149-2.831, p=0.006, and OR=1.754, 95% CI=1.931-2.552, p=0.021, respectively), haplogroups F, D, and M9 might be associated with a decreased risk for COPD in this population (OR=0.554, 95% CI=0.390-0.787, p=0.001; OR=0.758, 95% CI=0.407-0.965, p=0.002; and OR=0.186, 95% CI=0.039-0.881, p=0.034, respectively). Additionally, the increased frequency of the 822-bp mtDNA deletion in male cigarette-smoking subjects among COPD patients and controls of haplogroup D indicated that haplogroup D might increase an individual's susceptibility to DNA damage from external reactive oxygen species derived from heavy cigarette smoking. We conclude that haplogroups A and M7 might be risk factors for COPD, whereas haplogroups D, F, and M9 might decrease the COPD risk in this Han Chinese population.     

Allelic variation of the dopamine receptor D4 gene polymorphic region in gibbons

We examined the tandem repeat sequence of the dopamine receptor D4 (DRD4) gene in 73 individuals derived from 8 species of gibbons (genusHylobates) in an attempt to assess the variability of this gene in gibbon species.H. syndactylus (subgenusSymphalangus) andH. concolor (subgenusNomascus), which were inferred to have diverged at an early time within the family Hylobatidae, shared only long repeat (7–8) alleles. On the other hand, DRD4 was highly polymorphic in gibbons of the subgenusHylobates, with 4-, 5-, 6-, 7-, and 8-repeat alleles being recognized. In this subgenus, 4- and 5-repeat alleles were found in the species distributed mainly in the southern islands such as Sumatra, Java, and Borneo but not in the species inhabiting the Asian continent. Sequence analysis indicated that the repeat structure of the gibbon DRD4 gene was quite complex but most of the 48-bp units could be classified into several groups across the species based on sequence similarities. However, the sequence of the 7-repeat allele ofH. muelleri was unique, since the repeat units had low similarities to other units of gibbons.     

Relative frequencies of polymorphisms of variation in Block 2 repeats and 5' recombinant types of Plasmodium falciparum msp1 alleles

The mechanisms producing the genetic polymorphism at Plasmodium falciparum merozoite surface antigen-1 locus (pfmsp1) include the insertion and deletion of the different type of dimorphic Block 2 9-nucleotide repeat units as well as the intragenic recombination. To study relative occurrence frequencies of these two distinct mechanisms, we have developed a sensitive PCR strategy to identify both 5' recombinant types and the number of Block 2 repeats from the same sample. This method can specifically detect the target 5' recombinant type (Blocks 2-6) at the sensitivity of 1-4 copies of the pfmsp1. Applying the new method to field isolates from the Solomon Islands enabled us to identify six different 5' recombinant types and variation in Block 2 repeat number in three of them, thus distinguishing 10 different alleles. Distribution of these alleles in local three villages in the study area suggests that frequencies of variation in the number of Block 2 9-bp repeats and recombination events within Blocks 2-6 are mutually independent and the frequency of repeat variation is relatively high as compared to that of recombination events at the pfmsp1 locus in P. falciparum populations from the Solomon Islands.     

Maternal footprints of Southeast Asians in North India

We have analyzed 7,137 samples from 125 different caste, tribal and religious groups of India and 99 samples from three populations of Nepal for the length variation in the COII/tRNA(Lys) region of mtDNA. Samples showing length variation were subjected to detailed phylogenetic analysis based on HVS-I and informative coding region sequence variation. The overall frequencies of the 9-bp deletion and insertion variants in South Asia were 1.9 and 0.6%, respectively. We have also defined a novel deep-rooting haplogroup M43 and identified the rare haplogroup H14 in Indian populations carrying the 9-bp deletion by complete mtDNA sequencing. Moreover, we redefined haplogroup M6 and dissected it into two well-defined subclades. The presence of haplogroups F1 and B5a in Uttar Pradesh suggests minor maternal contribution from Southeast Asia to Northern India. The occurrence of haplogroup F1 in the Nepalese sample implies that Nepal might have served as a bridge for the flow of eastern lineages to India. The presence of R6 in the Nepalese, on the other hand, suggests that the gene flow between India and Nepal has been reciprocal.     

The 9-bp deletion in region V of mitochondrial DNA: evidence of mutation recurrence

A deletion of one of the two copies of a 9-bp direct repeat sequence (CCCCCTCTA) in region V of mitochondrial DNA has previously been used as a polymorphic anthropological marker for people of east Asian origin, and to a lesser extent, in Oceanian and African populations. We report the presence of the 9-bp deletion in homoplasmy in skeletal muscle fibers and lymphocytes of a Spanish Caucasian individual. Other mitochondrial DNA polymorphisms associated with the 9-bp deletion characteristic of other populations were not present. Our results suggest that the 9-bp deletion probably originated independently in the maternal lineage of the propositus, and that it can thus be described as a recurrent mutation.     

Comparison of the waxy locus sequence from a non-waxy strain and two waxy mutants of spontaneous and artificial origins in barley

Molecular characterization of 3 alleles of the waxy gene from a non-waxy strain "Shikoku hadaka No. 84" (SH84), an indigenous waxy strain "Mochimugi D" (MMD), and an artificial waxy mutant strain "Shikoku hadaka No. 97" (SH97) of barley (Hordeum vulgare ssp. vulgare) was performed via a PCR direct sequencing strategy. The 3 haplotypes were analyzed in terms of single nucleotide polymorphisms, insertion/deletion mutations, and simple sequence repeat polymorphisms. In comparison with the barley non-waxy gene sequence deposited in the public DNA database, 110 polymorphic sites were found in the 5,190-bp sequenced region of the non-waxy strain SH84. A 418-bp deletion in the 5' non-coding sequence was identified in the indigenous waxy strain MMD. Except for the deletion in the promoter region, the spontaneous mutant wax allele and non-waxy allele were identical. Such highly conserved sequences provide evidence for the recent occurrence of a deletion event in the cultivated barley gene pool. Compared to the original variety SH84, induced waxy mutant SH97 had a base substitution of a C to T in the exon 5, which converting Gln-89 of the wild-type gene into a stop codon, suggesting the involvement of a nonsense-mediated mRNA decay. These results will be helpful for understanding the mechanism of the variable amylose content in waxy cultivars of cereal species.     

A minisatellite and a microsatellite polymorphism within 1.5 kb at the human muscle glycogen phosphorylase (PYGM) locus can be amplified by PCR and have combined informativeness of PIC 0.95.

We sequenced a genomic clone (pMCMP1), previously reported to detect a VNTR polymorphism at the PYGM locus, and found a dinucleotide repeat segment (CA)14(GA)25 and a complex (AT)-repeat-rich segment containing 63 repeats spanning 160 bp. Resolution of PCR-amplified genomic DNA from the (CA)(GA) repeat region on DNA sequencing gels revealed a highly informative polymorphism with alleles differing by 2-bp intervals and ranging in size from 156 to 190 bp. Among three racial groups, a total of 18 alleles were observed. Fourteen alleles were observed in Caucasians (PIC 0.89), 12 alleles in American Blacks (PIC 0.89), and 9 alleles in Pima Indians (PIC 0.73). PCR amplification of the (AT) repeat region and resolution of the products on DNA sequencing gels revealed a complex variable length polymorphism with alleles distributed in size from 367 to 970 bp. Twenty-eight alleles were found in American Blacks (PIC 0.94), 6 alleles in Pima Indians (PIC 0.70), and 11 alleles in Caucasians (PIC 0.71). Comparison of the previously described VNTR RFLP alleles visualized by Southern hybridization to the PCR products described in this report demonstrated that the polymorphism described in both assays was identical. However, a larger number of alleles could be detected from the PCR-amplified products. Combined informativeness, PIC 0.95, for the two polymorphisms was determined from haplotype analysis of 100 Caucasian chromosomes. Therefore, for genotyping purposes, informativeness is maximized from using both polymorphisms.     

Deletions in Plasmid Pbr322: Replication Slippage Involving Leading and Lagging Strands

We test here whether a class of deletions likely to result from errors during DNA replication arise preferentially during synthesis of either the leading or the lagging DNA strand. Deletions were obtained by reversion of particular insertion mutant alleles of the pBR322 amp gene. The alleles contain insertions of palindromic DNAs bracketed by 9-bp direct repeats of amp sequence; in addition, bp 2 to 5 in one arm of the palindrome form a direct repeat with 4 bp of adjoining amp sequence. Prior work had shown that reversion to Ampr results from deletions with endpoints in the 8- or 4-bp repeat, and that the 4-bp repeats are used preferentially because one of them is in the palindrome. To test the role of leading and lagging strand synthesis in deletion formation, we reversed the direction of replication of the amp gene by inverting the pBR322 replication origin, and also constructed new mutant alleles with a 4-bp repeat starting counterclockwise rather than clockwise of the insertion. In both cases the 4-bp repeats were used preferentially as deletion endpoints. A model is presented in which deletions arise during elongation of the strand that copies the palindrome before the adjoining 4-bp repeat, and in which preferential use of the 4-bp repeats independent of the overall direction of replication implies that deletions arise during syntheses of both leading and lagging strands.     

AVPR1A and SLC6A4 polymorphisms in choral singers and non-musicians: a gene association study

Amateur choral singing is a common pastime and worthy of study, possibly conferring benefits to health and social behaviour. Participants might be expected to possess musical ability and share some behavioural characteristics. Polymorphisms in genes concerned with serotonergic neurotransmission are associated with both behaviour and musical aptitude. Those investigated previously include the variable number tandem repeats RS1, RS3 and AVR in the AVPR1A (arginine vasopressin receptor 1a) gene and STin2 in the SLC6A4 (solute carrier family 6 [neurotransmitter transporter, serotonin], member 4) gene, as well as the SLC6A4 promoter region polymorphism, 5-HTTLPR. We conducted a genetic association study on 523 participants to establish whether alleles at these polymorphisms occur more commonly in choral singers than in those not regularly participating in organised musical activity (non-musicians). We also analysed tagging single nucleotide polymorphisms (SNPs) for AVPR1A and SLC6A4 to determine whether other variants in these genes were associated with singer/non-musician status. At the STin2 polymorphism, overall association with singer/non-musician status was evident at P = 0.006. The 9-repeat (P = 0.04) and 12-repeat (P = 0.04) alleles were more common in singers and the 10-repeat allele less so (P = 0.009). Odds ratios were 0.73 (95% CI 0.57-0.94) for the 10-repeat allele and 2.47 (95% CI 0.88-6.94) for the rarer 9-repeat allele. No overall association was detected at P<0.05 between any other polymorphism and singer/non-musician status. Our null findings with respect to RS3, RS1 and AVR, polymorphisms associated with musical ability by other authors, suggest that choir membership may depend partly on factors other than musical ability. In a related musical project involving one participating choir, a new 40-part unaccompanied choral work, "Allele", was composed and broadcast on national radio. In the piece, each singer's part incorporated their personal RS3 genotype.     

Polymorphism analysis of tyrosine hydroxylase (TH) in military working dogs

Canine and human behavior are shaped by similar evolutionary processes, yet the identification of the behavioral phenotype is often difficult. A widely used method relies on breed stereotypes provided by experts such as dog trainers. To reveal a valid association between behavior and genetic factors, an association study of behavioral phenotyping and genotyping is essential. We screened for variable number of tandem repeat (VNTR) polymorphisms in intron 4 of the tyrosine hydroxylase gene in military working dogs (belonging to pass and fail groups based on the results of an in-training examination conducted by the drillmaster), which were scored based on possessiveness, audaciousness, concentration, and motor ability by qualification examination. We first characterized each genotype by sequencing, in which the 1/1 type consists of a single copy of a 36-bp sequence and the 2/2 type is a duplicated form of the 36-bp repeat unit. The 1/2 alleles showed a single nucleotide change as a heteroduplex, which generated a PCR product of similar size as that of the 1/1-182-bp. The military working dogs showed the 2/2 type of VNTR and heteroduplex. For the pass group, two dogs possessed 2/2 type (40 %), whereas three dogs were of the heteroduplex type (60 %). However, all members of fail group showed the 2/2 type (100 %). These data indicate that repeat polymorphisms with behavioral phenotyping can identify military working dogs that would pass or fail the in-training examination.     

Analysis of mtDNA variation in African populations reveals the most ancient of all human continent-specific haplogroups.

mtDNA sequence variation was examined in 140 Africans, including Pygmies from Zaire and Central African Republic (C.A.R.) and Mandenkalu, Wolof, and Pular from Senegal. More than 76% of the African mtDNAs (100% of the Pygmies and 67.3% of the Senegalese) formed one major mtDNA cluster (haplogroup L) defined by an African-specific HpaI site gain at nucleotide pair (np) 3592. Additional mutations subdivided haplogroup L into two subhaplogroups, each encompassing both Pygmy and Senegalese mtDNAs. A novel 12-bp homoplasmic insertion in the intergenic region between tRNA(Tyr) and cytochrome oxidase I (COI) genes was also observed in 17.6% of the Pygmies from C.A.R. This insertion is one of the largest observed in human mtDNAs. Another 25% of the Pygmy mtDNAs harbored a 9-bp deletion between the cytochrome oxidase II (COII) and tRNA(Lys) genes, a length polymorphism previously reported in non-African populations. In addition to haplogroup L, other haplogroups were observed in the Senegalese. These haplogroups were more similar to those observed in Europeans and Asians than to haplogroup L mtDNAs, suggesting that the African mtDNAs without the HpaI np 3592 site could be the ancestral types from which European and Asian mtDNAs were derived. Comparison of the intrapopulation sequence divergence in African and non-African populations confirms that African populations exhibit the largest extent of mtDNA variation, a result that further supports the hypothesis that Africans represent the most ancient human group and that all modern humans have a common and recent African origin. The age of the total African variation was estimated to be 101,000-133,000 years before present (YBP), while the age of haplogroup L was estimated at 98,000-130,000 YBP. These values substantially exceed the ages of all Asian- and European-specific mtDNA haplogroups.     

High-throughput analysis of informative CYP2D6 compound haplotypes

We describe a high-throughput protocol for detecting key polymorphisms in the drug-metabolizing enzyme gene CYP2D6 and a number of linked microsatellites that is both fast and relatively inexpensive to perform. This approach employs GeneScan technology to enable a researcher to determine rapidly the status of seven simple nucleotide polymorphisms in CYP2D6 and also to assay repeat number variation at five closely linked dinucleotide microsatellite loci. The method requires only three PCRs and two GeneScan runs per sample. We anticipate that this will be of value to researchers in three different ways: (1) rapid discrimination of common CYP2D6 alleles, (2) high-resolution haplotyping for association studies involving chromosome 22q13.1 using microsatellite variation, and (3) generation of compound haplotypes for investigating the evolution of CYP2D6 variation. We also report compound haplotype frequencies for an Ashkenazi Jewish and a British sample.     

Y-chromosome polymorphisms of the domestic Bactrian camel in China

Single-nucleotide polymorphisms (SNPs), microsatellites and copy number variation (CNV) were studied on the Y chromosome to understand the paternal origin and phylogenetic relationships for resource protection, rational development and utilization of the domestic Bactrian camel in China. Our sample set consisted of 94 Chinese domestic Bactrian camels from four regions (Inner Mongolia, Gansu, Qinghai and Xinjiang), we screened 29 Y-chromosome-specific loci for SNPs, analysed 40 bovine-derived microsatellite loci and measured CNVs of HSFY and SRY through Sanger sequencing, automated fluorescence-based microsatellite analysis and quantitative real-time PCR, respectively. A multicopy gene, SRY, was first found, and sequence variation was only detected in SRY in a screen of 29 loci in 13 DNA pools of individual camels. In addition, a TG repeat in the USP9Y gene was identified as the first polymorphic microsatellite in the camel Y chromosome, whereas microsatellite based on bovine sequences were not detected. The frequency of each allele varied among different populations. For the Nanjiang, Hexi and Alashan populations, a 243-bp allele was found. For the Sunite population, 241-bp, 243-bp and 247-bp alleles were detected, and the frequencies of these alleles were \(22.2\%\), \(44.5\%\) and \(33.3\%\), respectively; 241-bp and 243-bp alleles were found in other populations. Finally, CNVs in two Y-chromosomal genes were detected; CNV for HSFY and SRY ranged from 1 to 3 and from 1 to 9, respectively.