Functional variability of rhesus macaque (Macaca mulatta) NAT2 gene for drug-metabolising arylamine N-acetyltransferase 2.

Human NAT2 is a polymorphic pharmacogene encoding for N-acetyltransferase 2, a hepatic enzyme active towards arylamine and arylhydrazine drugs, including the anti-tubercular antibiotic isoniazid. The isoenzyme also modulates susceptibility to chemical carcinogenesis, particularly of the bladder. Human NAT2 represents an ideal model for anthropological investigations into the demographic adaptation of worldwide populations to their xenobiotic environment. Its sequence appears to be subject to positive selection pressures that are population-specific and may be attributed to gene-environment interactions directly associated with exogenous chemical challenges. However, recent evidence suggests that the same evolutionary pattern may not be observed in other primates. Here, we report NAT2 polymorphism in 25 rhesus macaques (Macaca mulatta) and compare the frequencies and functional characteristics of 12 variants. Seven non-synonymous single nucleotide variations (SNVs) were identified, including one nonsense mutation. The missense SNVs were demonstrated to affect enzymatic function in a substrate-dependent manner, albeit more moderately than certain NAT1 SNVs recently characterised in the same cohort. Haplotypic and functional variability of NAT2 was comparable to that previously observed for NAT1 in the same population sample, suggesting that the two paralogues may have evolved under similar selective pressures in the rhesus macaque. This is different to the population variability distribution pattern reported for humans and chimpanzees. Recorded SNVs were also different from those found in other primates. The study contributes to further understanding of NAT2 functional polymorphism in the rhesus macaque, a non-human primate model used in biomedicine and pharmacology, indicating variability in xenobiotic acetylation that could affect drug metabolism.

Evolutionary stasis of the pseudoautosomal boundary in strepsirrhine primates.

Sex chromosomes are typically comprised of a non-recombining region and a recombining pseudoautosomal region. Accurately quantifying the relative size of these regions is critical for sex-chromosome biology both from a functional and evolutionary perspective. The evolution of the pseudoautosomal boundary (PAB) is well documented in haplorrhines (apes and monkeys) but not in strepsirrhines (lemurs and lorises). Here, we studied the PAB of seven species representing the main strepsirrhine lineages by sequencing a male and a female genome in each species and using sex differences in coverage to identify the PAB. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events moved the PAB and shortened the pseudoautosomal region in haplorrhines. Strepsirrhines are well known to have much lower sexual dimorphism than haplorrhines. We suggest that mutations with antagonistic effects between males and females have driven recombination suppression and PAB evolution in haplorrhines.

Population variability of rhesus macaque (Macaca mulatta) NAT1 gene for arylamine N-acetyltransferase 1: Functional effects and comparison with human.

Human NAT1 gene for N-acetyltransferase 1 modulates xenobiotic metabolism of arylamine drugs and mutagens. Beyond pharmacogenetics, NAT1 is also relevant to breast cancer. The population history of human NAT1 suggests evolution through purifying selection, but it is unclear whether this pattern is evident in other primate lineages where population studies are scarce. We report NAT1 polymorphism in 25 rhesus macaques (Macaca mulatta) and describe the haplotypic and functional characteristics of 12 variants. Seven non-synonymous single nucleotide variations (SNVs) were identified and experimentally demonstrated to compromise enzyme function, mainly through destabilization of NAT1 protein and consequent activity loss. One non-synonymous SNV (c.560G > A, p.Arg187Gln) has also been characterized for human NAT1 with similar effects. Population haplotypic and functional variability of rhesus NAT1 was considerably higher than previously reported for its human orthologue, suggesting different environmental pressures in the two lineages. Known functional elements downstream of human NAT1 were also differentiated in rhesus macaque and other primates. Xenobiotic metabolizing enzymes play roles beyond mere protection from exogenous chemicals. Therefore, any link to disease, particularly carcinogenesis, may be via modulation of xenobiotic mutagenicity or more subtle interference with cell physiology. Comparative analyses add the evolutionary dimension to such investigations, assessing functional conservation/diversification among primates.

Humans and Chimpanzees Display Opposite Patterns of Diversity in Arylamine N-Acetyltransferase Genes.

Among the many genes involved in the metabolism of therapeutic drugs, human arylamine N-acetyltransferases (NATs) genes have been extensively studied, due to their medical importance both in pharmacogenetics and disease epidemiology. One member of this small gene family, NAT2, is established as the locus of the classic human acetylation polymorphism in drug metabolism. Current hypotheses hold that selective processes favoring haplotypes conferring lower NAT2 activity have been operating in modern humans' recent history as an adaptation to local chemical and dietary environments. To shed new light on such hypotheses, we investigated the genetic diversity of the three members of the NAT gene family in seven hominid species, including modern humans, Neanderthals and Denisovans. Little polymorphism sharing was found among hominids, yet all species displayed high NAT diversity, but distributed in an opposite fashion in chimpanzees and bonobos (Pan genus) compared to modern humans, with higher diversity in Pan species at NAT1 and lower at NAT2, while the reverse is observed in humans. This pattern was also reflected in the results returned by selective neutrality tests, which suggest, in agreement with the predicted functional impact of mutations detected in non-human primates, stronger directional selection, presumably purifying selection, at NAT1 in modern humans, and at NAT2 in chimpanzees. Overall, the results point to the evolution of divergent functions of these highly homologous genes in the different primate species, possibly related to their specific chemical/dietary environment (exposome) and we hypothesize that this is likely linked to the emergence of controlled fire use in the human lineage.

Comparative analysis of xenobiotic metabolising N-acetyltransferases from ten non-human primates as in vitro models of human homologues.

Xenobiotic metabolising N-acetyltransferases (NATs) perform biotransformation of drugs and carcinogens. Human NAT1 is associated with endogenous metabolic pathways of cells and is a candidate drug target for cancer. Human NAT2 is a well-characterised polymorphic xenobiotic metabolising enzyme, modulating susceptibility to drug-induced toxicity. Human NATs are difficult to express to high purification yields, complicating large-scale production for high-throughput screens or use in sophisticated enzymology assays and crystallography. We undertake comparative functional investigation of the NAT homologues of ten non-human primates, to characterise their properties and evaluate their suitability as models of human NATs. Considering the amount of generated recombinant protein, the enzymatic activity and thermal stability, the NAT homologues of non-human primates are demonstrated to be a much more effective resource for in vitro studies compared with human NATs. Certain NAT homologues are proposed as better models, such as the NAT1 of macaques Macaca mulatta and M. sylvanus, the NAT2 of Erythrocebus patas, and both NAT proteins of the gibbon Nomascus gabriellae which show highest homology to human NATs. This comparative investigation will facilitate in vitro screens towards discovery and optimisation of candidate pharmaceutical compounds for human NAT isoenzymes, while enabling better understanding of NAT function and evolution in primates.

How do rivers, geographic distance, and dispersal behavior influence genetic structure in two sympatric New World monkeys?

Dispersal, one of the major factors affecting the gene flow between populations, shapes the spatial distribution of genetic diversity within species. Alouatta macconnelli and Saguinus midas are two Neotropical monkey species that sympatrically inhabit the Guiana shield in northern Amazonia and are likely to differ in their dispersal behavior and vagility. We took advantage of their sympatry to investigate, over a fine geographical scale (∼50 km long), the relationship between spatial genetic structure, on the one hand, and geographical features and the species' dispersal behavior on the other. A total of 84 A. macconnelli individuals from 25 social units and 76 S. midas individuals from 19 social units were genotyped for nine microsatellite markers. Both species displayed high genetic diversity and allelic richness. However, patterns of genetic structure differed between the two species. In A. macconnelli, no genetic substructuring was observed, while in S. midas we detected significant structuring, but this structuring was not correlated with geographical features, such as the location of individuals relative to the river and/or the distance between them. Instead, the geographical distribution of genetic variation observed for each species is predominantly explained by each species' dispersal pattern. We identified bisexual dispersal for both species, but with significant differences, either in the distance or in the rate of dispersal, between species and sexes. Genetic relatedness within social units was higher in S. midas than in A. macconnelli: gene flow between social units seems limited in S. midas, especially for females, while high dispersal characterizes A. macconnelli, where females seem to disperse at lower rate but at a longer distance than males.

Spatial pattern of genetic diversity and selection in the MHC class II DRB of three Neotropical bat species.

BACKGROUND: Although bats are natural reservoirs of many pathogens, few studies have been conducted on the genetic variation and detection of selection in major histocompatibility complex (MHC) genes. These genes are critical for resistance and susceptibility to diseases, and host-pathogen interactions are major determinants of their extensive polymorphism. Here we examined spatial patterns of diversity of the expressed MHC class II DRB gene of three sympatric Neotropical bats, Carollia perspicillata and Desmodus rotundus (Phyllostomidae), and Molossus molossus (Molossidae), all of which use the same environments (e.g., forests, edge habitats, urban areas). Comparison with neutral marker (mtDNA D-loop) diversity was performed at the same time. RESULTS: Twenty-three DRB alleles were identified in 19 C. perspicillata, 30 alleles in 35 D. rotundus and 20 alleles in 28 M. molossus. The occurrence of multiple DRB loci was found for the two Phyllostomidae species. The DRB polymorphism was high in all sampling sites and different signatures of positive selection were detected depending on the environment. The patterns of DRB diversity were similar to those of neutral markers for C. perspicillata and M. molossus. In contrast, these patterns were different for D. rotundus for which a geographical structure was highlighted. A heterozygote advantage was also identified for this species. No recombination or gene conversion event was found and phylogenetic relationships showed a trans-species mode of evolution in the Phyllostomids. CONCLUSIONS: This study of MHC diversity demonstrated the strength of the environment and contrasting pathogen pressures in shaping DRB diversity. Differences between positively selected sites identified in bat species highlighted the potential role of gut microbiota in shaping immune responses. Furthermore, multiple geographic origins and/or population admixtures observed in C. perspicillata and M. molossus populations acted as an additional force in shaping DRB diversity. In contrast, DRB diversity of D. rotundus was shaped by environment rather than demographic history.

Genetic diversity and structuring across the range of a widely distributed ladybird: focus on rear-edge populations phenotypically divergent.

Population genetics and phenotypic structures are often predicted to vary along the geographic range of a species. This phenomenon would be accentuated for species with large range areas, with discontinuities and marginal populations. We herein compare the genetic patterns of central populations of Coccinella septempunctata L. with those of two phenotypically differentiated populations considered as rear-edge populations and subspecies based on phenotype (Algeria and Japan). According to the central-marginal model and expected characteristics of rear-edge populations, we hypothesize that these rear-edge populations have (1) a reduced genetic diversity, resulting from their relative isolation over long periods of time, (2) a higher population genetic differentiation, explained by low contemporary gene flow levels, and (3) a relationship between genetic diversity characteristics and phenotypes, due to historical isolation and/or local adaptation. Based on genotyping of 28 populations for 18 microsatellite markers, several levels of regional genetic diversity and differentiation are observed between and within populations, according to their localization: low within-population genetic diversity and higher genetic differentiation of rear-edge populations. The genetic structuring clearly dissociates the Algerian and Eastern Asia populations from the others. Geographical patterns of genetic diversity and differentiation support the hypothesis of the central-marginal model. The pattern observed is in agreement with the phenotypic structure across species range. A clear genetic break between populations of Algeria, the Eastern Asia, and the remaining populations is a dominant feature of the data. Differential local adaptations, absence of gene flow between marginal and central populations, and/or incapacity to mate after colonization, have contributed to their distinct genotypic and phenotypic characteristics.

Complete mitochondrial genome of the gray mouse lemur, Microcebus murinus (Primates, Cheirogaleidae).

We report the high-coverage complete mitochondrial genome sequence of the gray mouse lemur Microcebus murinus. The sequencing has been performed on an Illumina Hiseq 2500 platform, with a genome skimming strategy. The total length of this mitogenome is 16 963 bp, containing 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and 1 non-coding region (D-loop region). The genome organization, nucleotide composition and codon usage are similar to those reported from other primate's mitochondrial genomes. The complete mitochondrial genome sequence reported here will be useful for comparative genomics studies in primates.

Rapid birth-and-death evolution of the xenobiotic metabolizing NAT gene family in vertebrates with evidence of adaptive selection.

BACKGROUND: The arylamine N-acetyltransferases (NATs) are a unique family of enzymes widely distributed in nature that play a crucial role in the detoxification of aromatic amine xenobiotics. Considering the temporal changes in the levels and toxicity of environmentally available chemicals, the metabolic function of NATs is likely to be under adaptive evolution to broaden or change substrate specificity over time, making NATs a promising subject for evolutionary analyses. In this study, we trace the molecular evolutionary history of the NAT gene family during the last ~450 million years of vertebrate evolution and define the likely role of gene duplication, gene conversion and positive selection in the evolutionary dynamics of this family. RESULTS: A phylogenetic analysis of 77 NAT sequences from 38 vertebrate species retrieved from public genomic databases shows that NATs are phylogenetically unstable genes, characterized by frequent gene duplications and losses even among closely related species, and that concerted evolution only played a minor role in the patterns of sequence divergence. Local signals of positive selection are detected in several lineages, probably reflecting response to changes in xenobiotic exposure. We then put a special emphasis on the study of the last ~85 million years of primate NAT evolution by determining the NAT homologous sequences in 13 additional primate species. Our phylogenetic analysis supports the view that the three human NAT genes emerged from a first duplication event in the common ancestor of Simiiformes, yielding NAT1 and an ancestral NAT gene which in turn, duplicated in the common ancestor of Catarrhini, giving rise to NAT2 and the NATP pseudogene. Our analysis suggests a main role of purifying selection in NAT1 protein evolution, whereas NAT2 was predicted to mostly evolve under positive selection to change its amino acid sequence over time. These findings are consistent with a differential role of the two human isoenzymes and support the involvement of NAT1 in endogenous metabolic pathways. CONCLUSIONS: This study provides unequivocal evidence that the NAT gene family has evolved under a dynamic process of birth-and-death evolution in vertebrates, consistent with previous observations made in fungi.

A DNA metabarcoding study of a primate dietary diversity and plasticity across its entire fragmented range.

In tropical regions, most primary ecosystems have been replaced by mosaic landscapes in which species must cope with a large shift in the distribution of their habitat and associated food resources. Primates are particularly vulnerable to habitat modifications. Most species persist in small fragments surrounded by complex human-mediated matrices whose structure and connectivity may strongly influence their dispersal and feeding behavior. Behavioral plasticity appears to be a crucial parameter governing the ability of organisms to exploit the resources offered by new matrix habitats and thus to persist in fragmented habitats. In this study, we were interested in the dietary plasticity of the golden-crowned sifaka (Propithecus tattersalli), an endangered species of lemur, found only in the Daraina region in north-eastern Madagascar. We used a DNA-based approach combining the barcoding concept and Illumina next-generation sequencing to (i) describe the species diet across its entire range and (ii) evaluate the influence of landscape heterogeneity on diet diversity and composition. Faeces from 96 individuals were sampled across the entire species range and their contents were analyzed using the trnL metabarcoding approach. In parallel, we built a large DNA reference database based on a checklist of the plant species of the Daraina region. Our results suggest that golden-crowned sifakas exhibit remarkable dietary diversity with at least 130 plant species belonging to 80 genera and 49 different families. We highlighted an influence of both habitat type and openness on diet composition suggesting a high flexibility of foraging strategies. Moreover, we observed the presence of numerous cultivated and naturalized plants in the faeces of groups living in forest edge areas. Overall, our findings support our initial expectation that P. tattersalli is able to cope with the current level of alteration of the landscape and confirm our previous results on the distribution and the dispersal ability of this species.

The evolutionary landscape of cytosolic microbial sensors in humans.

Host-pathogen interactions are generally initiated by host recognition of microbial components or danger signals triggered by microbial invasion. This recognition involves germline-encoded microbial sensors or pattern-recognition receptors (PRRs). By studying the way in which natural selection has driven the evolution of these microbial sensors in humans, we can identify genes playing an essential role and distinguish them from other, more redundant genes. We characterized the sequence diversity of the NOD-like receptor family, including the NALP and NOD/IPAF subfamilies, in various populations worldwide and compared this diversity with that of other PRR families, such as Toll-like receptors (TLRs) and RIG-I-like receptors (RLRs). We found that most NALPs had evolved under strong selective constraints, suggesting that their functions are essential and possibly much broader than previously thought. Conversely, most NOD/IPAF subfamily members were subject to more relaxed selective constraints, suggesting greater redundancy. Furthermore, some NALP genes, including NLRP1, NLRP14, and CIITA, were found to have evolved adaptively. We identified those variants conferring a selective advantage on some human populations as the most likely targets of positive selection. More generally, the strength of selection differed considerably between the major families of microbial sensors. Endosomal TLRs and most NALPs were found to evolve under stronger purifying selection than most NOD/IPAF subfamily members and cell-surface TLRs and RLRs, suggesting some degree of redundancy in the signaling pathways triggered by these molecules. This study provides novel perspectives and experimentally testable hypotheses concerning the relative biological relevance of the various families of microbial sensors in humans.

Genetic data suggest a natural prehuman origin of open habitats in northern Madagascar and question the deforestation narrative in this region.

The impact of climate change and anthropogenic deforestation on biodiversity is of growing concern worldwide. Disentangling how past anthropogenic and natural factors contributed to current biome distribution is thus a crucial issue to understand their complex interactions on wider time scales and to improve predictions and conservation strategies. This is particularly important in biodiversity hotspots, such as Madagascar, dominated by large open habitats whose origins are increasingly debated. Although a dominant narrative argues that Madagascar was originally entirely covered by woodlands, which were destroyed by humans, a number of recent studies have suggested that past climatic fluctuations played a major role in shaping current biome distributions well before humans arrived. Here, we address the question of the origin of open habitats in the Daraina region in northern Madagascar, using a multiproxy approach combining population genetics modeling and remote-sensing analyses. We show that (i) contrary to most regions of Madagascar, the forest cover in Daraina remained remarkably stable over the past 60 y, and (ii) the golden-crowned sifaka (Propithecus tattersalli), a forest-dwelling lemur, underwent a strong population contraction before the arrival of the first humans, hence excluding an anthropogenic cause. Prehuman Holocene droughts may have led to a significant increase of grasslands and a reduction in the species' habitat. This contradicts the prevailing narrative that land cover changes are necessarily anthropogenic in Madagascar but does not preclude the later role played by humans in other regions in which recent lemur bottlenecks have been observed.

The selective footprints of viral pressures at the human RIG-I-like receptor family.

The RIG-I-like receptors (RLRs)--RIG-I, IFIH1 (or MDA5) and LGP2--are thought to be key actors in the innate immune system, as they play a major role in sensing RNA viruses in the cytosol of host cells. Despite the increasingly recognized importance of the RLR family in antiviral immunity, no population genetic studies have yet attempted to compare the evolutionary history of its different members in humans. Here, we characterized the levels of naturally occurring genetic variation in the RLRs in a panel of individuals of different ethnic origins, to assess to what extent natural selection has acted on this family of microbial sensors. Our results show that amino acid-altering variation at RIG-I, particularly in the helicase domain, has been under stronger evolutionary constraint than that at IFIH1 and LGP2, reflecting an important role for RIG-I in sensing numerous RNA viruses and/or functional constraints related to the binding of viral substrates. Such evolutionary constraints have been much more relaxed at IFIH1 and LGP2, which appear to have evolved adaptively in specific human populations. Notably, we identified several mutations showing signatures of positive selection, including two non-synonymous polymorphisms in IFIH1 (R460H and R843H) and one in LGP2 (Q425R), suggesting a selective advantage related to the sensing of RNA viruses by IFIH and to the regulatory functions of LGP2. In light of the fact that some of these mutations have been associated with altered risks of developing autoimmune disorders, our study provides an additional example of the evolutionary conflict between infection and autoimmunity.

Arylamine N-acetyltransferase 2 (NAT2) genetic diversity and traditional subsistence: a worldwide population survey.

Arylamine N-acetyltransferase 2 (NAT2) is involved in human physiological responses to a variety of xenobiotic compounds, including common therapeutic drugs and exogenous chemicals present in the diet and the environment. Many questions remain about the evolutionary mechanisms that have led to the high prevalence of slow acetylators in the human species. Evidence from recent surveys of NAT2 gene variation suggests that NAT2 slow-causing variants might have become targets of positive selection as a consequence of the shift in modes of subsistence and lifestyle in human populations in the last 10,000 years. We aimed to test more extensively the hypothesis that slow acetylation prevalence in humans is related to the subsistence strategy adopted by the past populations. To this end, published frequency data on the most relevant genetic variants of NAT2 were collected from 128 population samples (14,679 individuals) representing different subsistence modes and dietary habits, allowing a thorough analysis at both a worldwide and continent scale. A significantly higher prevalence of the slow acetylation phenotype was observed in populations practicing farming (45.4%) and herding (48.2%) as compared to populations mostly relying on hunting and gathering (22.4%) (P = 0.0007). This was closely mirrored by the frequency of the slow 590A variant that was found to occur at a three-fold higher frequency in food producers (25%) as compared to hunter-gatherers (8%). These findings are consistent with the hypothesis that the Neolithic transition to subsistence economies based on agricultural and pastoral resources modified the selective regime affecting the NAT2 acetylation pathway. Furthermore, the vast amount of data collected enabled us to provide a comprehensive and up-to-date description of NAT2 worldwide genetic diversity, thus building up a useful resource of frequency data for further studies interested in epidemiological or anthropological research questions involving NAT2.

Landscape genetics of an endangered lemur (Propithecus tattersalli) within its entire fragmented range.

Habitat fragmentation may strongly reduce individuals' dispersal among resource patches and hence influence population distribution and persistence. We studied the impact of landscape heterogeneity on the dispersal of the golden-crowned sifaka (Propithecus tattersalli), an endangered social lemur species living in a restricted and highly fragmented landscape. We combined spatial analysis and population genetics methods to describe population units and identify the environmental factors which best predict the rates and patterns of genetic differentiation within and between populations. We used non-invasive methods to genotype 230 individuals at 13 microsatellites in all the main forest fragments of its entire distribution area. Our analyses suggest that the Manankolana River and geographical distance are the primary structuring factors, while a national road crossing the region does not seem to impede gene flow. Altogether, our results are in agreement with a limited influence of forest habitat connectivity on gene flow patterns (except for North of the species' range), suggesting that dispersal is still possible today among most forest patches for this species. Within forest patches, we find that dispersal is mainly among neighbouring social groups, hence confirming previous behavioural observations.

Spatial variation in density and total size estimates in fragmented primate populations: the golden-crowned sifaka (Propithecus tattersalli).

The golden-crowned sifaka (Propithecus tattersalli) is an endangered lemur species found only in the Daraina region, a very restricted area in north-eastern Madagascar. Its forest habitat is highly fragmented and expected to suffer from significant changes in the near future. The species is poorly known and only one census study, carried out in 2000, has ever been published. It is thus crucial to update the conservation status of the golden-crowned sifaka before major anthropogenic environmental changes take place. Using the line-transect approach, we estimated the species density in the main forest fragments located in both the peripheral and central parts of the distribution range, including both protected and unprotected areas. In parallel, we tried to determine whether an edge effect could be detected by comparing densities at different distances from the forest edges. We found important variation of sifaka densities among forest fragments. The total species abundance is thus difficult to determine, but we estimated that it is likely to be over 18,000, two to three times higher than previously thought. However, our data also suggested that most P. tattersalli live in forests located in the central part of the distribution range and that the estimated densities in the central part were high (>80 individuals/km(2)). Two forest fragments, found to host a large part of the total population, are currently outside the managed area and their incorporation to the managed area is strongly recommended. Lastly, as expected for a folivorous and not heavily hunted species, our results are consistent with the hypothesis that this species does not experience a clear edge effect, at least during the first half of the dry season. This could be due to a high resiliency to habitat fragmentation or to the fact that fragmentation has been going on for some time.

A large panel of microsatellite markers for genetic studies in the infra-order catarrhini.

Many genetic studies on catarrhines use microsatellite markers that were isolated from human DNA. A large number of these markers have been characterized in the great apes, macaques and baboons. However, there are few or no markers available for other members of this group. In this study, an extensive literature search was performed to find microsatellite markers that had been successfully amplified across a range of catarrhine species. These conserved loci can provide a valuable starting point for characterizing loci in other catarrhines. Finally, microsatellite markers were tested in a range of species that are not well represented in the literature.

Compound evolutionary history of the rhesus macaque MHC class I B region revealed by microsatellite analysis and localization of retroviral sequences.

In humans, the single polymorphic B locus of the major histocompatibility complex is linked to the microsatellite MIB. In rhesus macaques, however, haplotypes are characterized by the presence of unique combinations of multiple B genes, which may display different levels of polymorphism. The aim of the study was to shed light on the evolutionary history of this highly complex region. First, the robustness of the microsatellite MIB-linked to almost half of the B genes in rhesus macaques (Mamu-B)-for accurate B haplotyping was studied. Based on the physical map of an established haplotype comprising 7 MIB loci, each located next to a certain Mamu-B gene, two MIB loci, MIB1 and MIB6, were investigated in a panel of MHC homozygous monkeys. MIB1 revealed a complex genotyping pattern, whereas MIB6 analysis resulted in the detection of one or no amplicon. Both patterns are specific for a given B haplotype, show Mendelian segregation, and even allow a more precise haplotype definition than do traditional typing methods. Second, a search was performed for retroelements that may have played a role in duplication processes as observed in the macaque B region. This resulted in the description of two types of duplicons. One basic unit comprises an expressed Mamu-B gene, adjacent to an HERV16 copy closely linked to MIB. The second type of duplicon comprises a Mamu-B (pseudo)gene, linked to a truncated HERV16 structure lacking its MIB segment. Such truncation seems to coincide with the loss of B gene transcription. Subsequent to the duplication processes, recombination between MIB and Mamu-B loci appears to have occurred, resulting in a hyperplastic B region. Thus, analysis of MIB in addition to B loci allows deciphering of the compound evolutionary history of the class I B region in Old World monkeys.

Assessing natural introgression in 2 biomedical model species, the rhesus macaque (Macaca mulatta) and the long-tailed macaque (Macaca fascicularis).

Rhesus macaque (Macaca mulatta) and long-tailed macaque (Macaca fascicularis) are the 2 most commonly used primate model species in biomedical sciences. Although morphological studies have revealed a weak hybridization at the interspecific contact zone, in the north of Indochina, a molecular study has suggested an ancient introgression from rhesus to long-tailed macaque into the Indo-Chinese peninsula. However, the gene flow between these 2 taxa has never been quantified using genetic data and theoretical models. In this study, we have examined genetic variation within and between the parapatric Chinese rhesus macaque and Indo-Chinese long-tailed macaque populations, using 13 autosomal, 5 sex-linked microsatellite loci and mitochondrial DNA sequence data. From these data, we assessed genetic structure and estimated gene flow using a Bayesian clustering approach and the "Isolation with Migration" model. Our results reveal a weak interspecific genetic differentiation at both autosomal and sex-linked loci, suggesting large population sizes and/or gene flow between populations. According to the Bayesian clustering, Chinese rhesus macaque is a highly homogeneous gene pool that contributes strongly to the current Indo-Chinese long-tailed macaque genetic makeup, whether or not current admixture is assumed. Coalescent simulations, which integrated the characteristics of the loci, pointed out 1) a higher effective population size in rhesus macaque, 2) no mitochondrial gene flow, and 3) unilateral and male-mediated nuclear gene flow of approximately 10 migrants per generation from rhesus to long-tailed macaque. These patterns of genetic structure and gene flow suggest extensive ancient introgression from Chinese rhesus macaque into the Indo-Chinese long-tailed macaque population.