Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline MHC-DRB peptide-binding region

Comparative Genomics of the Major Histocompatibility Complex (MHC) of Felids

This review summarizes the current knowledge on the major histocompatibility complex (MHC) of the family Felidae. This family comprises an important domestic species, the cat, as well as a variety of free-living felids, including several endangered species. As such, the Felidae have the potential to be an informative model for studying different aspects of the biological functions of MHC genes, such as their role in disease mechanisms and adaptation to different environments, as well as the importance of genetic diversity for conservation issues in free-ranging or captive populations. Despite this potential, the current knowledge on the MHC in the family as a whole is fragmentary and based mostly on studies of the domestic cat and selected species of big cats. The overall structure of the domestic cat MHC is similar to other mammalian MHCs following the general scheme “centromere-MHC class I-MHC class III-MHC class II” with some differences in the gene contents. An unambiguously defined orthologue of the non-classical class I HLA-E gene has not been identified so far and the class II DQ and DP genes are missing or pseudogenized, respectively. A comparison with available genomes of other felids showed a generally high level of structural and sequence conservation of the MHC region. Very little and fragmentary information on in vitro and/or in vivo biological functions of felid MHC genes is available. So far, no association studies have indicated effects of MHC genetic diversity on a particular disease. No information is available on the role of MHC class I molecules in interactions with Natural Killer (NK) cell receptors or on the putative evolutionary interactions (co-evolution) of the underlying genes. A comparison of complex genomic regions encoding NK cell receptors (the Leukocyte Receptor Complex, LRC and the Natural Killer Cell Complex, NKC) in the available felid genomes showed a higher variability in the NKC compared to the LRC and the MHC regions. Studies of the genetic diversity of domestic cat populations and/or specific breeds have focused mainly on DRB genes. Not surprisingly, higher levels of MHC diversity were observed in stray cats compared to pure breeds, as evaluated by DRB sequencing as well as by MHC-linked microsatellite typing. Immunogenetic analysis in wild felids has only been performed on MHC class I and II loci in tigers, Namibian leopards and cheetahs. This information is important as part of current conservation tasks to assess the adaptive potential of endangered wild species at the human-wildlife interface, which will be essential for preserving biodiversity in a functional ecosystem.

Low genetic variation in the MHC class II DRB gene and MHC-linked microsatellites in endangered island populations of the leopard cat (Prionailurus bengalensis) in Japan

Isolated populations of the leopard cat (Prionailurus bengalensis) on Tsushima and Iriomote islands in Japan are classified as subspecies P. b. euptilurus and P. b. iriomotensis, respectively. Because both populations have decreased to roughly 100, an understanding of their genetic diversity is essential for conservation. We genotyped MHC class II DRB exon 2 and MHC-linked microsatellite loci to evaluate the diversity of MHC genes in the Tsushima and Iriomote cat populations. We detected ten and four DRB alleles in these populations, respectively. A phylogenetic analysis showed DRB alleles from both populations to be closely related to those in other felid DRB lineages, indicating trans-species polymorphism. The MHC-linked microsatellites were more polymorphic in the Tsushima than in the Iriomote population. The MHC diversity of both leopard cat populations is much lower than in the domestic cat populations on these islands, probably due to inbreeding associated with founder effects, geographical isolation, or genetic drift. Our results predict low resistance of the two endangered populations to new pathogens introduced to the islands.

Retention of functional variation despite extreme genomic erosion: MHC allelic repertoires in the Lynx genus

BACKGROUND

Demographic bottlenecks erode genetic diversity and may increase endangered species' extinction risk via decreased fitness and adaptive potential. The genetic status of species is generally assessed using neutral markers, whose dynamic can differ from that of functional variation due to selection. The MHC is a multigene family described as the most important genetic component of the mammalian immune system, with broad implications in ecology and evolution. The genus Lynx includes four species differing immensely in demographic history and population size, which provides a suitable model to study the genetic consequences of demographic declines: the Iberian lynx being an extremely bottlenecked species and the three remaining ones representing common and widely distributed species. We compared variation in the most variable exon of the MHCI and MHCII-DRB loci among the four species of the Lynx genus.

RESULTS

The Iberian lynx was characterised by lower number of MHC alleles than its sister species (the Eurasian lynx). However, it maintained most of the functional genetic variation at MHC loci present in the remaining and genetically healthier lynx species at all nucleotide, amino acid, and supertype levels.

CONCLUSIONS

Species-wide functional genetic diversity can be maintained even in the face of severe population bottlenecks, which caused devastating whole genome genetic erosion. This could be the consequence of divergent alleles being retained across paralogous loci, an outcome that, in the face of frequent gene conversion, may have been favoured by balancing selection.

Genomic structure and expression pattern of MHC IIα and IIβ genes reveal an unusual immune trait in lined seahorse Hippocampus erectus

The major histocompatibility complex (MHC) genes are crucial in the adaptive immune system, and the gene duplication of MHC in animals can generally result in immune flexibility. In this study, we found that the lined seahorse (Hippocampus erectus) has only one gene copy number (GCN) of MHC IIα and IIβ, which is different from that in other teleosts. Together with the lack of spleen and gut-associated lymphatic tissue (GALT), the seahorse may be referred to as having a partial but natural "immunodeficiency". Highly variable amino acid residues were found in the IIα and IIβ domains, especially in the α1 and β1 domains with 9.62% and 8.43% allelic variation, respectively. Site models revealed seven and ten positively selected positions in the α1 and β1 domains, respectively. Real-time PCR experiments showed high expression levels of the MHC II genes in intestine (In), gill (Gi) and trunk kidney (TK) and medium in muscle (Mu) and brood pouch (BP), and the expression levels were significantly up-regulated after bacterial infection. Specially, relative higher expression level of both MHC IIα and IIβ was found in Mu and BP when compared with other fish species, in which MHC II is expressed negligibly in Mu. These results indicate that apart from TK, Gi and In, MU and BP play an important role in the immune response against pathogens in the seahorse. In conclusion, high allelic variation and strong positive selection in PBR and relative higher expression in MU and BP are speculated to partly compensate for the immunodeficiency.

Characterisation of major histocompatibility complex class I genes in Japanese Ranidae frogs

The major histocompatibility complex (MHC) is a key component of adaptive immunity in all jawed vertebrates, and understanding the evolutionary mechanisms that have shaped these genes in amphibians, one of the earliest terrestrial tetrapods, is important. We characterised MHC class I variation in three common Japanese Rana species (Rana japonica, Rana ornativentris and Rana tagoi tagoi) and identified a total of 60 variants from 21 individuals. We also found evolutionary signatures of gene duplication, recombination and balancing selection (including trans-species polymorphism), all of which drive increased MHC diversity. A unique feature of MHC class I from these three Ranidae species includes low synonymous differences per site (d_S) within species, which we attribute to a more recent diversification of these sequences or recent gene duplication. The resulting higher d_N/d_S ratio relative to other anurans studied could be related to stronger selection pressure at peptide binding sites. This is one of the first studies to investigate MHC in Japanese amphibians and permits further exploration of the polygenetic factors associated with resistance to infectious diseases.

Genetic diversity and differentiation of the rhesus macaque (Macaca mulatta) population in western Sichuan, China, based on the second exon of the major histocompatibility complex class II DQB (MhcMamu-DQB1) alleles

BACKGROUND

Rhesus macaques living in western Sichuan, China, have been separated into several isolated populations due to habitat fragmentation. Previous studies based on the neutral or nearly neutral markers (mitochondrial DNA or microsatellites) showed high levels of genetic diversity and moderate genetic differentiation in the Sichuan rhesus macaques. Variation at the major histocompatibility complex (MHC) loci is widely accepted as being maintained by balancing selection, even with a low level of neutral variability in some species. However, in small and isolated or bottlenecked populations, balancing selection may be overwhelmed by genetic drift. To estimate microevolutionary forces acting on the isolated rhesus macaque populations, we examined genetic variation at Mhc-DQB1 loci in 119 wild rhesus macaques from five geographically isolated populations in western Sichuan, China, and compared the levels of MHC variation and differentiation among populations with that previously observed at neutral microsatellite markers.

RESULTS

23 Mamu-DQB1 alleles were identified in 119 rhesus macaques in western Sichuan, China. These macaques exhibited relatively high levels of genetic diversity at Mamu-DQB1. The Hanyuan population presented the highest genetic variation, whereas the Heishui population was the lowest. Analysis of molecular variance (AMOVA) and pairwise FST values showed moderate genetic differentiation occurring among the five populations at the Mhc-DQB1 locus. Non-synonymous substitutions occurred at a higher frequency than synonymous substitutions in the peptide binding region. Levels of MHC variation within rhesus macaque populations are concordant with microsatellite variation. On the phylogenetic tree for the rhesus and crab-eating macaques, extensive allele or allelic lineage sharing is observed between the two species.

CONCLUSIONS

Phylogenetic analyses confirm the apparent trans-species model of evolution of the Mhc-DQB1 genes in these macaques. Balancing selection plays an important role in sharing allelic lineages between species, but genetic drift may share balancing selection dominance to maintain MHC diversity. Great divergence at neutral or adaptive markers showed that moderate genetic differentiation had occurred in rhesus macaque populations in western Sichuan, China, due to the habitat fragmentation caused by long-term geographic barriers and human activity. The Heishui population should be paid more attention for its lowest level of genetic diversity and relatively great divergence from others.

Polymorphisms and tissue expression of the feline leukocyte antigen class I loci FLAI-E, FLAI-H, and FLAI-K

Cytotoxic CD8+ T-cell immunosurveillance for intracellular pathogens, such as viruses, is controlled by classical major histocompatibility complex (MHC) class Ia molecules, and ideally, these antiviral T-cell populations are defined by the specific peptide and restricting MHC allele. Surprisingly, despite the utility of the cat in modeling human viral immunity, little is known about the feline leukocyte antigen class I complex (FLAI). Only a few coding sequences with uncertain locus origin and expression patterns have been reported. Of 19 class I genes, three loci--FLAI-E, FLAI-H, and FLAI-K--are predicted to encode classical molecules, and our objective was to evaluate their status by analyzing polymorphisms and tissue expression. Using locus-specific, PCR-based genotyping, we amplified 33 FLAI-E, FLAI-H, and FLAI-K alleles from 12 cats of various breeds, identifying, for the first time, alleles across three distinct loci in a feline species. Alleles shared the expected polymorphic and invariant sites in the α1/α2 domains, and full-length cDNA clones possessed all characteristic class Ia exons. Alleles could be assigned to a specific locus with reasonable confidence, although there was evidence of potentially confounding interlocus recombination between FLAI-E and FLAI-K. Only FLAI-E, FLAI-H, and FLAI-K origin alleles were amplified from cDNAs of multiple tissue types. We also defined hypervariable regions across these genes, which permitted the assignment of names to both novel and established alleles. As predicted, FLAI-E, FLAI-H, and FLAI-K fulfill the major criteria of class Ia genes. These data represent a necessary prerequisite for studying epitope-specific antiviral CD8+ T-cell responses in cats.

Identification and characterization of the major histocompatibility complex class II DQB (MhcMath-DQB1) alleles in Tibetan macaques (Macaca thibetana)

Tibetan macaque (Macaca thibetana), an endangered primate species endemic to China, have been used as experimental animal model for various human diseases. Major histocompatibility complex (MHC) genes play a crucial role in the susceptibility and/or resistance to many human diseases, but little is known about Tibetan macaques. To gain an insight into the MHC background and to facilitate the experimental use of Tibetan macaques, the second exon of Mhc-DQB1 gene was sequenced in a cohort of wild Tibetan macaques living in the Sichuan province of China. A total of 23 MhcMath-DQB1 alleles were identified for the first time, illustrating a marked allelic polymorphism at the DQB1 locus for these macaques. Most of the sequences (74%) observed in this study belong to DQB1*06 (9 alleles) and DQB1*18 (8 alleles) lineages, and the rest (26%) belong to DQB1*15 (3 alleles) and DQB1*17 (3 alleles) lineages. The most frequent alleles detected among these macaques were MhcMath-DQB1*15:02:02 (17.9%), followed by Math-DQB1*06:06, 17:03 and 18:01, which were detected in 9 (16.1%) of the monkeys, respectively. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the peptide-binding region, suggesting balancing selection for maintaining polymorphisms at the MHC class II DQB1 locus. Phylogenetic analyses confirms the trans-species model of evolution of the Mhc-DQB1 genes in non-human primates, and in particular, the extensive allele sharing is observed between Tibetan and other macaque species.

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.

Diversity and evolutionary patterns of immune genes in free-ranging Namibian leopards (Panthera pardus pardus)

The genes of the major histocompatibility complex (MHC) are a key component of the mammalian immune system and have become important molecular markers for fitness-related genetic variation in wildlife populations. Currently, no information about the MHC sequence variation and constitution in African leopards exists. In this study, we isolated and characterized genetic variation at the adaptively most important region of MHC class I and MHC class II-DRB genes in 25 free-ranging African leopards from Namibia and investigated the mechanisms that generate and maintain MHC polymorphism in the species. Using single-stranded conformation polymorphism analysis and direct sequencing, we detected 6 MHC class I and 6 MHC class II-DRB sequences, which likely correspond to at least 3 MHC class I and 3 MHC class II-DRB loci. Amino acid sequence variation in both MHC classes was higher or similar in comparison to other reported felids. We found signatures of positive selection shaping the diversity of MHC class I and MHC class II-DRB loci during the evolutionary history of the species. A comparison of MHC class I and MHC class II-DRB sequences of the leopard to those of other felids revealed a trans-species mode of evolution. In addition, the evolutionary relationships of MHC class II-DRB sequences between African and Asian leopard subspecies are discussed.

MHC class I and MHC class II DRB gene variability in wild and captive Bengal tigers (Panthera tigris tigris)

Bengal tigers are highly endangered and knowledge on adaptive genetic variation can be essential for efficient conservation and management. Here we present the first assessment of allelic variation in major histocompatibility complex (MHC) class I and MHC class II DRB genes for wild and captive tigers from India. We amplified, cloned, and sequenced alpha-1 and alpha-2 domain of MHC class I and beta-1 domain of MHC class II DRB genes in 16 tiger specimens of different geographic origin. We detected high variability in peptide-binding sites, presumably resulting from positive selection. Tigers exhibit a low number of MHC DRB alleles, similar to other endangered big cats. Our initial assessment-admittedly with limited geographic coverage and sample size-did not reveal significant differences between captive and wild tigers with regard to MHC variability. In addition, we successfully amplified MHC DRB alleles from scat samples. Our characterization of tiger MHC alleles forms a basis for further in-depth analyses of MHC variability in this illustrative threatened mammal.