Extensive DRB region diversity in cynomolgus macaques: recombination as a driving force

The DR region of primate species is generally complex and displays diversity concerning the number and combination of distinct types of DRB genes present per region configuration. A highly variable short tandem repeat (STR) present in intron 2 of nearly all primate DRB genes can be utilized as a quick and accurate high through-put typing procedure. This approach resulted previously in the description of unique and haplotype-specific DRB-STR length patterns in humans, chimpanzees, and rhesus macaques. For the present study, a cohort of 230 cynomolgus monkeys, including self-sustaining breeding groups, has been examined. MtDNA analysis showed that most animals originated from the Indonesian islands, but some are derived from the mainland, south and north of the Isthmus of Kra. Haplotyping and subsequent sequencing resulted in the detection of 118 alleles, including 28 unreported ones. A total of 49 Mafa-DRB region configurations were detected, of which 28 have not yet been described. Humans and chimpanzees possess a low number of different DRB region configurations in concert with a high degree of allelic variation. In contrast, however, allelic heterogeneity within a given Mafa-DRB configuration is even less frequently observed than in rhesus macaques. Several of these region configurations appear to have been generated by recombination-like events, most probably propagated by a retroviral element mapping within DRB6 pseudogenes, which are present on the majority of haplotypes. This undocumented high level of DRB region configuration-associated diversity most likely represents a species-specific strategy to cope with various pathogens

Definition of Mafa-A and -B haplotypes in pedigreed cynomolgus macaques (Macaca fascicularis)

The major histocompatibility complex (MHC) class I B gene/allelic repertoire was investigated in a pedigreed population of cynomolgus macaques of mixed Indonesian/Malaysian origin. The Mafa-B alleles detected in this cohort are mostly specific for a given geographic area, and only a small number of alleles appears to be shared with other populations. This suggests the fast evolution of Mafa-B alleles due to adaptation to new environments. In contrast to humans, the B locus in Old World monkeys displays extensive copy number variation. The Mafa-B and previously defined -A gene combinations segregate in families and thus allowed the definition of extended haplotypes. In many cases it was possible to assign a particular Mafa-I allele to one of these Mafa-A/B haplotypes as well. The presence of a large number of stable haplotypes in this cohort of animals, which was pedigreed for up to eight generations, looks promising for developing discriminative MHC typing tools that are less cumbersome. Furthermore, the discovery of 53 unreported Mafa-B sequences expands the lexicon of alleles significantly, and may help in understanding the complex organisation of the macaque B region

DR haplotype diversity of the cynomolgus macaque as defined by its transcriptome

The DR region of particular primate species may display allelic polymorphism and gene copy number variation (region configuration polymorphism). The sum of these distinct types of polymorphism is defined as complexity. To date, however, the DR region of cynomolgus macaques (Macaca fascicularis) has been poorly defined. Transcriptome analysis of a pedigreed colony, comprising animals from Indonesia and Indochina, revealed a total of 15 Mafa-DRA and 57 DRB alleles, specifying 28 different region configurations. The DRA alleles can be divided into two distinct lineages. One lineage is polymorphic, but the majority of the amino acid replacements map to the leader peptide. The second lineage is at best oligomorphic, and segregates with one specific Mafa-DRB allele. The number of Mafa-DRB genes ranges from two to five per haplotype. Due to the presence of pseudogenes, however, each haplotype encodes only one to three bona fide DRB transcripts. Depending on the region configuration in which the Mafa-DRB gene is embedded, identical alleles may display differential transcription levels. Region configurations appear to have been generated by recombination-like events. When genes or gene segments are relocated, it seems plausible that they may be placed in the context of distinct transcription control elements. As such, DRB region-related transcription level differences may add an extra layer of polymorphism to this section of the adaptive immune system

Improved Xenobiotic Metabolism and Reduced Susceptibility to Cancer in Gluten-Sensitive Macaques upon Introduction of a Gluten-Free Diet

A non-human primate (NHP) model of gluten sensitivity was employed to study the gene perturbations associated with dietary gluten changes in small intestinal tissues from gluten-sensitive rhesus macaques (Macaca mulatta).Stages of remission and relapse were accomplished in gluten-sensitive animals by administration of gluten-free (GFD) and gluten-containing (GD) diets, as described previously. Pin-head-sized biopsies, obtained non-invasively by pediatric endoscope from duodenum while on GFD or GD, were used for preparation of total RNA and gene profiling, using the commercial Rhesus Macaque Microarray (Agilent Technologies),targeting expression of over 20,000 genes.When compared with normal healthy control, gluten-sensitive macaques showed differential gene expressions induced by GD. While observed gene perturbations were classified into one of 12 overlapping categories--cancer, metabolism, digestive tract function, immune response, cell growth, signal transduction, autoimmunity, detoxification of xenobiotics, apoptosis, actin-collagen deposition, neuronal and unknown function--this study focused on cancer-related gene networks such as cytochrome P450 family (detoxification function) and actin-collagen-matrix metalloproteinases (MMP) genes.A loss of detoxification function paralleled with necessity to metabolize carcinogens was revealed in gluten-sensitive animals while on GD. An increase in cancer-promoting factors and a simultaneous decrease in cancer-preventing factors associated with altered expression of actin-collagen-MMP gene network were noted. In addition, gluten-sensitive macaques showed reduced number of differentially expressed genes including the cancer-associated ones upon withdrawal of dietary gluten. Taken together, these findings indicate potentially expanded utility of gluten-sensitive rhesus macaques in cancer research

Genomic plasticity of the MHC class I A region in rhesus macaques: extensive haplotype diversity at the population level as revealed by microsatellites

The Mamu-A genes of the rhesus macaque show different degrees of polymorphism, transcription level variation, and differential haplotype distribution. Per haplotype, usually one “major” transcribed gene is present, A1 (A7), in various combinations with “minor” genes, A2 to A6. In silico analysis of the physical map of a heterozygous animal revealed the presence of similar Mamu-A regions consisting of four duplication units, but with dissimilar positions of the A1 genes on both haplotypes, and in combination with different minor genes. Two microsatellites, D6S2854 and D6S2859, have been selected as potential tools to characterize this complex region. Subsequent analysis of a large breeding colony resulted in the description of highly discriminative patterns, displaying copy number variation in concert with microsatellite repeat length differences. Sequencing and segregation analyses revealed that these patterns are unique for each Mamu-A haplotype. In animals of Indian, Burmese, and Chinese origin, 19, 15, or 9 haplotypes, respectively, could be defined, illustrating the occurrence of differential block duplications and subsequent rearrangements by recombination. The haplotypes can be assigned to 12 unique combinations of genes (region configurations). Although most configurations harbor two transcribed A genes, one or three genes per haplotype are also present. Additionally, haplotypes lacking an A1 gene or with an A1 duplication appear to exist. The presence of different transcribed A genes/alleles in monkeys from various origins may have an impact on differential disease susceptibilities. The high-throughput microsatellite technique will be a valuable tool in animal selection for diverse biomedical research projects

Blood ties: ABO is a trans-species polymorphism in primates

The ABO histo-blood group, the critical determinant of transfusion incompatibility, was the first genetic polymorphism discovered in humans. Remarkably, ABO antigens are also polymorphic in many other primates, with the same two amino acid changes responsible for A and B specificity in all species sequenced to date. Whether this recurrence of A and B antigens is the result of an ancient polymorphism maintained across species or due to numerous, more recent instances of convergent evolution has been debated for decades, with a current consensus in support of convergent evolution. We show instead that genetic variation data in humans and gibbons as well as in Old World Monkeys are inconsistent with a model of convergent evolution and support the hypothesis of an ancient, multi-allelic polymorphism of which some alleles are shared by descent among species. These results demonstrate that the ABO polymorphism is a trans-species polymorphism among distantly related species and has remained under balancing selection for tens of millions of years, to date, the only such example in Hominoids and Old World Monkeys outside of the Major Histocompatibility Complex.Comment: 45 pages, 4 Figures, 4 Supplementary Figures, 5 Supplementary Table

Hidden Patterns of Anti-HLA Class I Alloreactivity Revealed Through Machine Learning

Detection of alloreactive anti-HLA antibodies is a frequent and mandatory test before and after organ transplantation to determine the antigenic targets of the antibodies. Nowadays, this test involves the measurement of fluorescent signals generated through antibody–antigen reactions on multi-beads flow cytometers. In this study, in a cohort of 1,066 patients from one country, anti-HLA class I responses were analyzed on a panel of 98 different antigens. Knowing that the immune system responds typically to “shared” antigenic targets, we studied the clustering patterns of antibody responses against HLA class I antigens without any a priori hypothesis, applying two unsupervised machine learning approaches. At first, the principal component analysis (PCA) projections of intralocus specific responses showed that anti-HLA-A and anti-HLA-C were the most distantly projected responses in the population with the anti-HLA-B responses to be projected between them. When PCA was applied on the responses against antigens belonging to a single locus, some already known groupings were confirmed while several new crossreactive patterns of alloreactivity were detected. Anti-HLA-A responses projected through PCA suggested that three cross-reactive groups accounted for about 70% of the variance observed in the population, while anti-HLA-B responses were mainly characterized by a distinction between previously described Bw4 and Bw6 cross-reactive groups followed by several yet undocumented or poorly described ones. Furthermore, anti-HLA-C responses could be explained by two major cross-reactive groups completely overlapping with previously described C1 and C2 allelic groups. A second featurebased analysis of all antigenic specificities, projected as a dendrogram, generated a robust measure of allelic antigenic distances depicting bead-array defined cross reactive groups. Finally, amino acid combinations explaining major population specific crossreactive groups were described. The interpretation of the results was based on the current knowledge of the antigenic targets of the antibodies as they have been characterized either experimentally or computationally and appear at the HLA epitope registry

Patterns of 1,748 Unique Human Alloimmune Responses Seen by Simple Machine Learning Algorithms

Allele specific antibody response against the polymorphic system of HLA is the allogeneic response marker determining the immunological risk for graft acceptance before and after organ transplantation and therefore routinely studied during the patient’s workup. Experimentally, bead bound antigen- antibody reactions are detected using a special multicolor flow cytometer (Luminex). Routinely for each sample, antibody responses against 96 different HLA antigen groups are measured simultaneously and a 96-dimensional immune response vector is created. Under a common experimental protocol, using unsupervised clustering algorithms, we analyzed these immune intensity vectors of anti HLA class II responses from a dataset of 1,748 patients before or after renal transplantation residing in a single country. Each patient contributes only one serum sample in the analysis. A population view of linear correlations of hierarchically ordered fluorescence intensities reveals patterns in human immune responses with striking similarities with the previously described CREGs but also brings new information on the antigenic properties of class II HLA molecules. The same analysis affirms that “public” anti-DP antigenic responses are not correlated to anti DR and anti DQ responses which tend to cluster together. Principal Component Analysis (PCA) projections also demonstrate ordering patterns clearly differentiating anti DP responses from anti DR and DQ on several orthogonal planes. We conclude that a computer vision of human alloresponse by use of several dimensionality reduction algorithms rediscovers proven patterns of immune reactivity without any a priori assumption and might prove helpful for a more accurate definition of public immunogenic antigenic structures of HLA molecules. Furthermore, the use of Eigen decomposition on the Immune Response generates new hypotheses that may guide the design of more effective patient monitoring tests

A snapshot of the Mamu-B genes and their allelic repertoire in rhesus macaques of Chinese origin

The major histocompatibility complex class I gene repertoire was investigated in a large panel of rhesus macaques of Chinese origin. As observed in Indian animals, subjects of Chinese derivation display Mamu-B gene copy number variation, and the sum of expressed genes varies among haplotypes. In addition, these genes display differential transcription levels. The majority of the Mamu-B alleles discovered during this investigation appear to be unique for the population studied. Only one particular Mamu-B haplotype is shared between Indian and Chinese animals, and it must have been present in the progenitor stock. Hence, the data highlight the fact that most allelic polymorphism, and most of the Mamu-B haplotypes themselves, are of relatively recent origin and were most likely generated after the separation of the Indian and Chinese rhesus macaque populations

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