Effects of geographic origin on captive Macaca mulatta mitochondrial DNA variation

Partial sequences from mitochondrial (mt) 12S and 16S rRNA genes were analyzed to characterize diversity among captive rhesus macaques (Macaca mulatta) originating from various geographic regions. Several nested clades, defined by closely related haplotypes, were identified, suggesting considerable genetic subdivision, probably relics from heterogeneous origins, founder effects, and genetic drift, followed by breeding isolation. The rhesus matrilineages from India differed discretely and markedly from Chinese matrilineages; approximately 90% of the genetic heterogeneity among the combined samples of Indian and Chinese rhesus macaques studied here was due to country of origin. In addition, mtDNA sequences from macaques of China were more diverse than those from rhesus macaques of India, an outcome consistent with China's greater subspecies diversity and with nuclear genotype distributions. Otherwise, the distribution of mtDNA variation within rhesus macaques of China, and especially within those of India, exhibited far less structure and did not conform to a simple isolation-by-distance model. As the demand for genetically heterogeneous and well-characterized rhesus macaques for biomedical-based research increases, mtDNA haplotypes can be useful for genetically defining, preserving maximal levels of genetic diversity within, and confirming the geographic origin of captive breeding groups of rhesus macaques.     

Mitochondrial DNA variation within and among regional populations of longtail macaques (Macaca fascicularis) in relation to other species of the fascicularis group of macaques

An 835 base pair (bp) fragment of mitochondrial DNA (mtDNA) was sequenced to characterize genetic variation within and among 1,053 samples comprising five regional populations each of longtail macaques (Macaca fascicularis) and rhesus macaques (Macaca mulatta), and one sample each of Japanese (M. fuscata) and Taiwanese (M. cyclopis) macaques. The mtDNA haplotypes of longtail macaques clustered in two large highly structured clades (Fas1 and Fas2) of a neighbor-joining tree that were reciprocally monophyletic with respect to those representing rhesus macaques, Japanese macaques, and Taiwanese macaques. Both clades exhibited haplotypes of Indonesian and Malaysian longtail macaques widely dispersed throughout them; however, longtail macaques from Indochina, Philippines, and Mauritius each clustered in a separate well-defined clade together with one or a few Malaysian and/or Indonesian longtail macaques, suggesting origins on the Sunda shelf. Longtail macaques from Malaysia and Indonesia were far more genetically diverse, and those from Mauritius were far less diverse than any other population studied. Nucleotide diversity between mtDNA sequences of longtail macaques from different geographic regions is, in some cases, greater than that between Indian and Chinese rhesus macaques. Approximately equal amounts of genetic diversity are due to differences among animals in the same regional population, different regional populations, and different species. A greater proportion of genetic variance was explained by interspecies differences when Japanese and Taiwanese macaques were regarded as regional populations of rhesus macaques than when they were treated as separate species. Rhesus macaques from China were more closely related to both Taiwanese and Japanese macaques than to their own conspecifics from India.     

Phylogenetic analysis of chinese rhesus macaques (Macaca mulatta) based on mitochondrial control region sequences

Between one and six subspecies of Chinese rhesus macaques (Macaca mulatta) have been proposed based on morphological differences and/or their geographic distribution. In this study, a 489 base pair fragment of the mitochondrial control region was amplified from 230 DNA samples collected from rhesus macaques in the Sichuan province in Western China. The fragment was then sequenced and aligned with 208 sequences from wild rhesus macaques, sampled throughout the species' geographic range in China downloaded from GenBank. Phylogenetic analysis of the 182 unique sequences identified among these samples divided Chinese rhesus macaques into two western haplogroups (haplogroups A and B) and three older eastern haplogroups (haplogroups C, D, and E), whose differentiation probably occurred during the penultimate glacial event. During the warming after the penultimate glacial event, haplogroups A, B, and E rapidly expanded and a relatively young subhaplogroup of haplogroup E, E', limited to Southern China but shared with Vietnamese rhesus macaques, was reintroduced from Indochina during the last glacial event. One haplotype most closely related to subhaplogroup E' probably represents the isolation of Hainan Island, to where it is restricted, from the mainland by the formation of the Qiongzhou Strait approximately 8,500 years ago. The distribution of haplogroups both informs the phylogeographic history of dispersal of Chinese rhesus macaques and has implications for their suitability as animal models in biomedical research.     

Genetic characterization of wild and captive rhesus macaques in China

The genetic structures of wild and captive rhesus macaque populations within China were compared by analyzing the mtDNA sequences of 203 captive-bred Chinese rhesus macaques with 77 GenBank sequences from wild-caught animals trapped throughout China. The genotypes of 22 microsatellites of captive Chinese rhesus macaques were also compared with those of captive Indian animals. The Chinese population is significantly differentiated from the Indian population and is more heterogeneous. Thus, compared with Indian rhesus macaques the phenotypic variance of traits with high heritability will be inflated in Chinese animals. Our data suggest that the western Chinese provinces have more subdivided populations than the eastern and southern Chinese provinces. The southern Chinese populations are the least structured and might have been more recently established. Human-mediated interbreeding among captive Chinese populations has occurred, implying that Chinese breeding strategies can influence the interpretation of biomedical research in the USA.     

Geographic Patterns of Genetic Variation in a Broadly Distributed Marine Vertebrate: New Insights into Loggerhead Turtle Stock Structure from Expanded Mitochondrial DNA Sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (∼800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.     

High frequency of triplicatedα-globin genes in tropical primates,crab-eating macaques (Macaca fascicularis), chimpanzees (Pan troglodytes), and orang-utans (Pongo pygmaeus)

Chimpanzees and orang-utans had triplicatedα-globin gene haplotypes in the frequency of 0.80 and 0.20, respectively. Homozygous duplicated haplotypes could not be found in any of the 44 chimpanzees examined. Chimpanzees having homozygous triplicated haplotypes have greater numbers of red blood cells than those chimpanzees heterozygous for the duplicate and triplicate haplotype. Crab-eating macaques in Malayan peninsula of Thailand had triplicated haplotypes occurring in frequencies ranging from 0.13 to 0.50. On the other hand, triplicated haplotypes occurred in very low frequencies (0–0.07) in crab-eating macaques in the northern and eastern part of Thailand as well as in rhesus macaques from India and China, and in Japanese macaques.     

Identification of Country of Origin and Admixture Between Indian and Chinese Rhesus Macaques

We describe a restriction analysis that distinguishes between rhesus macaques of unmixed Indian and Chinese ancestry and between western and eastern Chinese ancestry. We amplified a 254-bp fragment of mitochondrial DNA (mtDNA) that contains restriction sites hypothesized to be diagnostic of country of origin for samples from 534 and 567 individuals alleged to be of solely Indian or solely Chinese origin, respectively. After digestion with the MaeIII, SmlI, and BccI restriction enzymes, the polymerase chain reaction (PCR) products of only 3 of the 1101 samples exhibited restriction patterns uncharacteristic of their alleged country of origin. A sample comprising 392 of these rhesus macaques was genotyped for 24 nuclear microsatellite (STR) loci. Principal coordinates analysis confirmed marked genetic similarity of regional populations within each country but a substantial difference between Indian and Chinese rhesus macaques. Using STRUCTURE (Pritchard and Wen, 2003),we assigned probabilities of Chinese and Indian ancestry to each sample based on its STR genotypes. We assigned all the unmixed rhesus macaques to their correct countries of origin with probabilities >0.95. We constructed an artificial sample of 1st-generation hybrid Indian/Chinese rhesus macaques by randomly sampling from the genotypes of Indian and Chinese individuals. STRUCTURE assigned robabilities of Chinese and Indian ancestry to hybrids that closely corresponded with the proportions of alleles in that sample drawn from unmixed Chinese and Indian rhesus macaques.     

A Large‐Scale SNP‐Based Genomic Admixture Analysis of the Captive Rhesus Macaque Colony at the California National Primate Research Center

Some breeding facilities in the United States have crossbred Chinese and Indian rhesus macaque (Macaca mulatta) founders either purposefully or inadvertently. Genetic variation that reflects geographic origins among research subjects has the potential to influence experimental outcomes. The use of animals from different geographic regions, their hybrids, and animals of varying degrees of kinship in an experiment can obscure treatment effects under study because high interanimal genetic variance can increase phenotypic variance among the research subjects. The intent of this study, based on a broad genomic analysis of 2,808 single nucleotide polymorphisms (SNPs), is to ensure that only animals estimated to be of pure Indian or Chinese ancestry, based on both demographic and genetic information, are used as sources of infants for derivation and expansion of the California National Primate Research Center's (CNPRC) super‐Specific Pathogen Free (SSPF) rhesus macaque colony. Studies of short tandem repeats (STRs) in Indian and Chinese rhesus macaques have reported that heterozygosity of STRs is higher in Chinese rhesus macaques than in Indian rhesus macaques. The present study shows that heterozygosity of SNPs is actually higher in Indian than in Chinese rhesus macaques and that the Chinese SSPF rhesus macaque colony is far less differentiated from their founders compared to the Indian‐origin animals. The results also reveal no evidence of recent gene flow from long‐tailed and pig‐tailed macaques into the source populations of the SSPF rhesus macaques. This study indicates that many of the long‐tailed macaques held in the CNPRC are closely related individuals. Most polymorphisms shared among the captive rhesus, long‐tailed, and pig‐tailed macaques likely predate the divergence among these groups. Am. J. Primatol. 74:747‐757, 2012. © 2012 Wiley Periodicals, Inc.     

Regional differences in the distribution of the sub-Saharan, West Eurasian, and South Asian mtDNA lineages in Yemen

Despite its key location for population movements out of and back into Africa, Yemen has not yet been sampled on a regional level for an investigation of sub-Saharan, West Eurasian, and South Asian genetic contributions. In this study, we present mitochondrial DNA (mtDNA) data for regionally distinct Yemeni populations that reveal different distributions of mtDNA lineages. An extensive database of mtDNA sequences from North and East African, Middle Eastern and Indian populations was analyzed to provide a context for the regional Yemeni mtDNA datasets. The groups of western Yemen appear to be most closely related to Middle Eastern and North African populations, while the eastern Yemeni population from Hadramawt is most closely related to East Africa. Furthermore, haplotype matches with Africa are almost exclusively confined to West Eurasian R0a haplogroup in southwestern Yemen, although more sub-Saharan L-type matches appear in more northern Yemeni populations. In fact, Yemeni populations have the highest frequency of R0a haplotypes detected to date, thus Yemen or southern Arabia may be the site of the initial expansion of this haplogroup. Whereas two variants of the sub-Saharan haplogroup M1 were detected only in southwestern Yemen close to the Bab el-Mandeb Strait, different non-African M haplotypes were detected at low frequencies (approximately 2%) in western parts of the country and at a higher frequency (7.5%) in the Hadramawt. We conclude that the Yemeni gene pool is highly stratified both regionally and temporally and that it has received West Eurasian, Northeast African, and South Asian gene flow.     

4040 SNPs for genomic analysis in the rhesus macaque (Macaca mulatta)

Although the rhesus macaque (Macaca mulatta) is commonly used for biomedical research and becoming a preferred model for translational medicine, quantification of genome-wide variation has been slow to follow the publication of the genome in 2007. Here we report the properties of 4040 single nucleotide polymorphisms discovered and validated in Chinese and Indian rhesus macaques from captive breeding colonies in the United States. Frequency-matched measures of linkage disequilibrium were much greater in the Indian sample. Although the majority of polymorphisms were shared between the two populations, rare alleles were over twice as common in the Chinese sample. Indian rhesus had higher rates of heterozygosity, as well as previously undetected substructure, potentially due to admixture from Burma in wild populations and demographic events post-captivity.     

Characterization of 47 MHC class I sequences in Filipino cynomolgus macaques

Cynomolgus macaques (Macaca fascicularis) provide increasingly common models for infectious disease research. Several geographically distinct populations of these macaques from Southeast Asia and the Indian Ocean island of Mauritius are available for pathogenesis studies. Though host genetics may profoundly impact results of such studies, similarities and differences between populations are often overlooked. In this study we identified 47 full-length MHC class I nucleotide sequences in 16 cynomolgus macaques of Filipino origin. The majority of MHC class I sequences characterized (39 of 47) were unique to this regional population. However, we discovered eight sequences with perfect identity and six sequences with close similarity to previously defined MHC class I sequences from other macaque populations. We identified two ancestral MHC haplotypes that appear to be shared between Filipino and Mauritian cynomolgus macaques, notably a Mafa-B haplotype that has previously been shown to protect Mauritian cynomolgus macaques against challenge with a simian/human immunodeficiency virus, SHIV_89.6P. We also identified a Filipino cynomolgus macaque MHC class I sequence for which the predicted protein sequence differs from Mamu-B*17 by a single amino acid. This is important because Mamu-B*17 is strongly associated with protection against simian immunodeficiency virus (SIV) challenge in Indian rhesus macaques. These findings have implications for the evolutionary history of Filipino cynomolgus macaques as well as for the use of this model in SIV/SHIV research protocols. Kevin J. Campbell and Ann M. Detmer contributed equally to this work.     

Comparison of the frequencies of major histocompatibility (MHC) class-II DQA1 and DQB1 alleles in Indian and Chinese rhesus macaques (Macaca mulatta)

The major histocompatibility complex (MHC) comprises related gene families, some of which are highly polymorphic, whose protein products mediate immune response. Rhesus macaques (Macaca mulatta) are a vital animal model for research in human diseases and are native to regions extending from Afghanistan in the west to the Eastern Plains of China and from Peking to the north, southward through islands of Southeast Asia. The distributions of MHC class-II Mamu DQA1 and Mamu DQB1 alleles in two groups of domestically bred rhesus macaques of Indian and Chinese origin and the Mamu DQA1 genotypes of a small number of Burmese rhesus macaques were compared. Major allelic differences were observed between the Indian and Chinese rhesus macaques, and gene diversity decreased from east to west. These and other intra-specific genetic differences among regional populations of rhesus macaques might influence the outcome of biomedical research in which they are used as subjects, and illustrate the importance of completely genetically characterizing subjects used as animal models in biomedical research.     

Genetic characterization of Indian-origin and Chinese-origin rhesus macaques (Macaca mulatta)

Genetic differences between Indian-origin and Chinese-origin rhesus macaques are as great as those between some primate species and can influence the results of experiments in which both are used as animal models for the study of the same human diseases. Unfortunately, many breeding facilities do not know with certainty the origin of the founders of their rhesus breeding colonies. Here I summarize the most definitive of the genetic traits among the microsatellite (STR) loci and mitochondrial DNA sequences that my laboratory previously reported to characterize Indian-origin and Chinese-origin rhesus macaques and then estimate the frequencies of these traits and their reliability as indicators of country of origin. The expression of diagnostic traits at two or more of four different unlinked loci provides a nearly 100% reliability in distinguishing rhesus macaques of Indian and Chinese origin.     

Genetic relationship of the melon fly,Bactrocera cucurbitae (Diptera: Tephritidae) inferred from mitochondrial DNA

Abstract The melon fruit fly, Bactrocera cucurbitae (Coquillett) (Diptera: Tephritidae), has been the subject of worldwide quarantine and management efforts due to its widespread agricultural impact and potential for rapid range expansion. From its presumed native distribution in India, this species has spread throughout the hot‐humid regions of the world. We provide information that reveals population structure, invasion history and population connectivity from 23 locations covering nine countries based on DNA sequences of the mitochondrial cytochrome oxidase I (COI) gene. Forty‐two polymorphic sites were described among 38 haplotypes. The most common haplotype, H1, was observed in 73% of the samples distributed among all populations. Highest genetic diversity was seen within populations, and no isolation‐by‐distance was detected. The western regions (Nepal, Bangladesh, Thailand, Burma and China‐west) showed higher haplotype diversity than eastern regions (China‐east). China‐Yunnan showed highest levels of genetic diversity in China. Haplotype diversity decreased with longitude from west to east. Together, these analyses suggest that B. cucurbitae has expanded from west to east within a limited geographic scale and recently invaded China through Yunnan Province.     

Mitochondrial DNA variation of domestic sheep (Ovis aries) in Kenya

The history of domestic sheep (Ovis aries) in Africa remains largely unknown. After being first introduced from the Near East, sheep gradually spread through the African continent with pastoral societies. The eastern part of Africa was important either for the first diffusion of sheep southward or for putative secondary introductions from the Arabian Peninsula or southern Asia. We analysed mitochondrial DNA control region sequences of 91 domestic sheep from Kenya and found a high diversity of matrilines from the widespread haplogroup B, whereas only a single individual from haplogroup A was detected. Our phylogeography analyses of more than 500 available mitochondrial DNA sequences also identified ancestral haplotypes that were probably first introduced in Africa and are now widely distributed. Moreover, we found no evidence of an admixture between East and West African sheep. The presence of shared haplotypes in eastern and ancient southern African sheep suggests the possible southward movement of sheep along the eastern part of Africa. Finally, we found no evidence of an extensive introduction of sheep from southern Asia into Africa via the Indian Ocean trade. The overall findings on the phylogeography of East African domestic sheep set the grounds for understanding the origin and subsequent movements of sheep in Africa. The richness of maternal lineages in Kenyan breeds is of prime importance for future conservation and breeding programmes.     

Simple sequence repeat (SSR) polymorphisms for colony management and population genetics in rhesus macaques (Macaca mulatta)

Cross-species amplification of 72 SSR (predominantly tetranucleotide) loci from the DNA of six rhesus macaques of diverse regional origins was conducted using human primers for the polymerase chain reaction (PCR). Thirteen of these primer pairs, which consistently and unambiguously amplified polymorphic fragments from these six samples and which exhibited Mendelian properties, were also used to amplify SSR loci for 176 male rhesus macaques that are founders of six different captive breeding colonies. These include four groups of macaques originating in India and one group of macaques each that originated in China and Thailand. Gene diversity based on the SSR loci provided a reliable estimate of average heterozygosity but was between two and four times higher than that for 9 protein coding loci. Based on the SSR loci, Chinese rhesus were more genetically diverse and unique than were rhesus from India or Thailand, a conclusion not consistent with data based on protein coding loci. The six most informative SSR loci are unlinked and provide probabilities of single parent exclusion and genetic identity exceeding 0.99 and one in one million, respectively, both of which are reasonable standards for colony management purposes. Am. J. Primatol. 42:199–213, 1997. © 1997 Wiley-Liss, Inc.     

Genomic structures and population histories of linguistically distinct tribal groups of India

There are various conflicting hypotheses regarding the origins of the tribal groups of India, who belong to three major language groups--Austro-Asiatic, Dravidian and Tibeto-Burman. To test some of the major hypotheses we designed a genetic study in which we sampled tribal populations belonging to all the three language groups. We used a set of autosomal DNA markers, mtDNA restriction-site polymorphisms (RSPs) and mtDNA hypervariable segment-1 (HVS-1) sequence polymorphisms in this study. Using the unlinked autosomal markers we found that there is a fair correspondence between linguistic and genomic affinities among the Indian tribal groups. We reconstructed mtDNA RSP haplotypes and found that there is extensive haplotype sharing among all tribal populations. However, there is very little sharing of mtDNA HVS-1 sequences across populations, and none across language groups. Haplogroup M is ubiquitous, and the subcluster U2i of haplogroup U occurs in a high frequency. Our analyses of haplogroup and HVS-1 sequence data provides evidence in support of the hypothesis that the Austro-Asiatic speakers are the most ancient inhabitants of India. Our data also support the earlier finding that some of the western Eurasian haplogroups found in India may have been present in India prior to the entry of Aryan speakers. However, we do not find compelling evidence to support the theory that haplogroup M was brought into India on an "out of Africa" wave of migration through a southern exit route from Ethiopia. On the contrary, our data raise the possibility that this haplogroup arose in India and was later carried to East Africa from India.     

MHC class I characterization of Indonesian cynomolgus macaques

Cynomolgus macaques (Macaca fascicularis) are quickly becoming a useful model for infectious disease and transplantation research. Even though cynomolgus macaques from different geographic regions are used for these studies, there has been limited characterization of full-length major histocompatibility complex (MHC) class I immunogenetics of distinct geographic populations. Here, we identified 48 MHC class I cDNA nucleotide sequences in eleven Indonesian cynomolgus macaques, including 41 novel Mafa-A and Mafa-B sequences. We found seven MHC class I sequences in Indonesian macaques that were identical to MHC class I sequences identified in Malaysian or Mauritian macaques. Sharing of nucleotide sequences between these geographically distinct populations is also consistent with the hypothesis that Indonesia was a source of the Mauritian macaque population. In addition, we found that the Indonesian cDNA sequence Mafa-B7601 is identical throughout its peptide binding domain to Mamu-B03, an allele that has been associated with control of Simian immunodeficiency virus (SIV) viremia in Indian rhesus macaques. Overall, a better understanding of the MHC class I alleles present in Indonesian cynomolgus macaques improves their value as a model for disease research, and it better defines the biogeography of cynomolgus macaques throughout Southeast Asia.     

Mitochondrial phylogeography and subspecific variation in the red panda (Ailurus fulgens): implications for conservation

The red panda (Ailurus fulgens) is an endangered species and its present distribution is restricted to isolated mountain ranges in western China (Sichuan, Yunnan, and Tibet provinces) and the Himalayan Mountains chain of Nepal, India, Bhutan, and Burma. To examine the evolutionary history across its current range, and to assess the genetic divergence among current subspecies and population structure among different geographic locations, we sequenced mitochondrial DNA from the control region (CR) and cytochrome (cyt) b gene for 41 individuals in Sichuan, Yunnan, Tibet of China, and Burma. 25 CR haplotypes (10 for cyt b) were identified from 11 geographic locations. Only three haplotypes were shared among sample localities, including one among current subspecies. Nine haplotypes were shared with the study of Su et al. [Mol. Biol. Evol. 18 (2001) 1070]. CR haplotype diversity was high (0.95+/-0.02) and nucleotide diversity among all haplotypes was relatively low (0.018+/-0.009). Phylogenetic confirmed trees show a shallow pattern with very little structure or statistical robustness. The application of two coalescent-based tests for population growth allowed us to interpret this phylogeny as the result of a recent population expansion. Analysis of molecular variance and nested clade analysis failed to detect significant geographic structure in both data sets. The lack of significant differentiation between subspecies does not indicate the presence of evolutionary significant units. We suggest that the present population structure has resulted from habitat fragmentation and expansion from glacial refugia. Due to its habitat requirements it is likely that the red panda has undergone bottlenecks and population expansions several times in the recent past. The present population may exhibit a pattern reminiscent of a relatively recent population expansion.