MHC class I of saltwater crocodiles (Crocodylus porosus): polymorphism and balancing selection

Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles

Characterizing polymorphism at the major histocompatibility complex (MHC) genes is key to understanding the vertebrate immune response to disease. Despite being globally afflicted by the infectious tumour disease fibropapillomatosis (FP), immunogenetic variation in sea turtles is minimally explored. We sequenced the α 1 peptide-binding region of MHC class I genes (162 bp) from 268 juvenile green (Chelonia mydas) and 88 loggerhead (Caretta caretta) sea turtles in Florida, USA. We recovered extensive variation (116 alleles) and trans-species polymorphism. Supertyping analysis uncovered three functional MHC supertypes corresponding to the three well-supported clades in the phylogeny. We found significant evidence of positive selection at seven amino acid sites in the class I exon. Random forest modelling and risk ratio analysis of Ch. mydas alleles uncovered one allele weakly associated with smooth FP tumour texture, which may be associated with disease outcome. Our study represents the first characterization of MHC class I diversity in Ch. mydas and the largest sample of sea turtles used to date in any study of adaptive genetic variation, revealing tremendous genetic variation and high adaptive potential to viral pathogen threats. The novel associations we identified between MHC diversity and FP outcomes in sea turtles further highlight the importance of evaluating genetic predictors of disease, including MHC and other functional markers.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

Effects of Population Bottleneck and Balancing Selection on the Chinese Alligator Are Revealed by Locus-Specific Characterization of MHC Genes

Chinese alligator (Alligator sinensis) is an endangered freshwater crocodilian endemic to China, which experienced a severe bottleneck about 30 years ago. In this study, we developed locus-specific primers to investigate the polymorphism of 3 major histocompatibility complex (MHC) loci in 3 Chinese alligator populations, in combination with 6 neutral microsatellite markers as a contrast. We found the genetic trace for the bottleneck effect on the endangered Chinese alligator: the low allelic diversity (2 alleles at each locus), the low nucleotide substitution rate (no more than 0.009) at all sites, the deviation from Hardy-Weinberg Equilibrium/heterozygote deficiency, and the significant Tajima’s D values, indicating the MHC class I and class II loci being at different stages of bottleneck. We also obtained 3 pieces of evidence for balancing selection on this severely bottlenecked reptile: an obvious excess of nonsynonymous substitutions over synonymous at the antigen-binding positions, the mean synonymous substitution rate of MHC exons significantly higher than mean nucleotide substitution rate of introns, and the differentiation coefficient F_ST of MHC loci significantly lower than that of microsatellite loci. Consequently, we emphasize that the Chinese alligator holds a pretty low adaptive ability and requires scientific conservation strategies to ensure the long-term population development.

Balancing selection and recombination as evolutionary forces caused population genetic variations in golden pheasant MHC class I genes

BACKGROUND

The major histocompatibility complex (MHC) genes are vital partners in the acquired immune processes of vertebrates. MHC diversity may be directly associated with population resistance to infectious pathogens. Here, we screened for polymorphisms in exons 2 and 3 of the IA1 and IA2 genes in 12 golden pheasant populations across the Chinese mainland to characterize their genetic variation levels, to understand the effects of historical positive selection and recombination in shaping class I diversity, and to investigate the genetic structure of wild golden pheasant populations.

RESULTS

Among 339 individual pheasants, we identified 14 IA1 alleles in exon 2 (IA1-E2), 11 IA1-E3 alleles, 27 IA2-E2 alleles, and 28 IA2-E3 alleles. The non-synonymous substitution rate was significantly greater than the synonymous substitution rate at sequences in the IA2 gene encoding putative peptide-binding sites but not in the IA1 gene; we also found more positively selected sites in IA2 than in IA1. Frequent recombination events resulted in at least 9 recombinant IA2 alleles, in accordance with the intermingling pattern of the phylogenetic tree. Although some IA alleles are widely shared among studied populations, large variation occurs in the number of IA alleles across these populations. Allele frequency analysis across 2 IA loci showed low levels of genetic differentiation among populations on small geographic scales; however, significant genetic differentiation was observed between pheasants from the northern and southern regions of the Yangtze River. Both STRUCTURE analysis and F-statistic (F ST ) value comparison classified those populations into 2 major groups: the northern region of the Yangtze River (NYR) and the southern region of the Yangtze River (SYR).

CONCLUSIONS

More extensive polymorphisms in IA2 than IA1 indicate that IA2 has undergone much stronger positive-selection pressure during evolution. Moreover, the recombination events detected between the genes and the intermingled phylogenetic pattern indicate that interlocus recombination accounts for much of the allelic variation in IA2. Analysis of the population differentiation implied that homogenous balancing selection plays an important part in maintaining an even distribution of MHC variations. The natural barrier of the Yangtze River and heterogeneous balancing selection might help shape the NYR-SYR genetic structure in golden pheasants.

High genetic diversity in the endangered and narrowly distributed amphibian species Leptobrachium leishanense

Threatened species typically have a small or declining population size, which make them highly susceptible to loss of genetic diversity through genetic drift and inbreeding. Genetic diversity determines the evolutionary potential of a species; therefore, maintaining the genetic diversity of threatened species is essential for their conservation. In this study, we assessed the genetic diversity of the adaptive major histocompatibility complex (MHC) genes in an endangered and narrowly distributed amphibian species, Leptobrachium leishanense in Southwest China. We compared the genetic variation of MHC class I genes with that observed in neutral markers (5 microsatellite loci and cytochrome b gene) to elucidate the relative roles of genetic drift and natural selection in shaping the current MHC polymorphism in this species. We found a high level of genetic diversity in this population at both MHC and neutral markers compared with other threatened amphibian species. Historical positive selection was evident in the MHC class I genes. The higher allelic richness in MHC markers compared with that of microsatellite loci suggests that selection rather than genetic drift plays a prominent role in shaping the MHC variation pattern, as drift can affect all the genome in a similar way but selection directly targets MHC genes. Although demographic analysis revealed no recent bottleneck events in L. leishanense, additional population decline will accelerate the dangerous status for this species. We suggest that the conservation management of L. leishanense should concentrate on maximizing the retention of genetic diversity through preventing their continuous population decline. Protecting their living habitats and forbidding illegal hunting are the most important measures for conservation of L. leishanense.

Comparative genome analyses reveal distinct structure in the saltwater crocodile MHC

The major histocompatibility complex (MHC) is a dynamic genome region with an essential role in the adaptive immunity of vertebrates, especially antigen presentation. The MHC is generally divided into subregions (classes I, II and III) containing genes of similar function across species, but with different gene number and organisation. Crocodylia (crocodilians) are widely distributed and represent an evolutionary distinct group among higher vertebrates, but the genomic organisation of MHC within this lineage has been largely unexplored. Here, we studied the MHC region of the saltwater crocodile (Crocodylus porosus) and compared it with that of other taxa. We characterised genomic clusters encompassing MHC class I and class II genes in the saltwater crocodile based on sequencing of bacterial artificial chromosomes. Six gene clusters spanning ∼452 kb were identified to contain nine MHC class I genes, six MHC class II genes, three TAP genes, and a TRIM gene. These MHC class I and class II genes were in separate scaffold regions and were greater in length (2-6 times longer) than their counterparts in well-studied fowl B loci, suggesting that the compaction of avian MHC occurred after the crocodilian-avian split. Comparative analyses between the saltwater crocodile MHC and that from the alligator and gharial showed large syntenic areas (>80% identity) with similar gene order. Comparisons with other vertebrates showed that the saltwater crocodile had MHC class I genes located along with TAP, consistent with birds studied. Linkage between MHC class I and TRIM39 observed in the saltwater crocodile resembled MHC in eutherians compared, but absent in avian MHC, suggesting that the saltwater crocodile MHC appears to have gene organisation intermediate between these two lineages. These observations suggest that the structure of the saltwater crocodile MHC, and other crocodilians, can help determine the MHC that was present in the ancestors of archosaurs.

MHC class II diversity of koala (Phascolarctos cinereus) populations across their range

Major histocompatibility complex class II (MHCII) genes code for proteins that bind and present antigenic peptides and trigger the adaptive immune response. We present a broad geographical study of MHCII DA β1 (DAB) and DB β1 (DBB) variants of the koala (Phascolarctos cinereus; n=191) from 12 populations across eastern Australia, with a total of 13 DAB and 7 DBB variants found. We identified greater MHCII variation and, possibly, additional gene copies in koala populations in the north (Queensland and New South Wales) relative to the south (Victoria), confirmed by STRUCTURE analyses and genetic differentiation using analysis of molecular variance. The higher MHCII diversity in the north relative to south could potentially be attributed to (i) significant founder effect in Victorian populations linked to historical translocation of bottlenecked koala populations and (ii) increased pathogen-driven balancing selection and/or local genetic drift in the north. Low MHCII genetic diversity in koalas from the south could reduce their potential response to disease, although the three DAB variants found in the south had substantial sequence divergence between variants. This study assessing MHCII diversity in the koala with historical translocations in some populations contributes to understanding the effects of population translocations on functional genetic diversity.

Selection and trans-species polymorphism of major histocompatibility complex class II genes in the order Crocodylia

Major Histocompatibility Complex (MHC) class II genes encode for molecules that aid in the presentation of antigens to helper T cells. MHC characterisation within and between major vertebrate taxa has shed light on the evolutionary mechanisms shaping the diversity within this genomic region, though little characterisation has been performed within the Order Crocodylia. Here we investigate the extent and effect of selective pressures and trans-species polymorphism on MHC class II α and β evolution among 20 extant species of Crocodylia. Selection detection analyses showed that diversifying selection influenced MHC class II β diversity, whilst diversity within MHC class II α is the result of strong purifying selection. Comparison of translated sequences between species revealed the presence of twelve trans-species polymorphisms, some of which appear to be specific to the genera Crocodylus and Caiman. Phylogenetic reconstruction clustered MHC class II α sequences into two major clades representing the families Crocodilidae and Alligatoridae. However, no further subdivision within these clades was evident and, based on the observation that most MHC class II α sequences shared the same trans-species polymorphisms, it is possible that they correspond to the same gene lineage across species. In contrast, phylogenetic analyses of MHC class II β sequences showed a mixture of subclades containing sequences from Crocodilidae and/or Alligatoridae, illustrating orthologous relationships among those genes. Interestingly, two of the subclades containing sequences from both Crocodilidae and Alligatoridae shared specific trans-species polymorphisms, suggesting that they may belong to ancient lineages pre-dating the divergence of these two families from the common ancestor 85-90 million years ago. The results presented herein provide an immunogenetic resource that may be used to further assess MHC diversity and functionality in Crocodylia.

Evolution of MHC class I in the order Crocodylia

The major histocompatibility complex (MHC) is a dynamic genomic region with an essential role in the adaptive immunity of jawed vertebrates. The evolution of the MHC has been dominated by gene duplication and gene loss, commonly known as the birth-and-death process. Evolutionary studies of the MHC have mostly focused on model species. However, the investigation of this region in non-avian reptiles is still in its infancy. To provide insights into the evolutionary mechanisms that have shaped the diversity of this region in the Order Crocodylia, we investigated MHC class I exon 3, intron 3, and exon 4 across 20 species of the families Alligatoridae and Crocodilidae. We generated 124 DNA sequences and identified 31 putative functional variants as well as 14 null variants. Phylogenetic analyses revealed three gene groups, all of which were present in Crocodilidae but only one in Alligatoridae. Within these groups, variants generally appear to cluster at the genus or family level rather than in species-specific groups. In addition, we found variation in gene copy number and some indication of interlocus recombination. These results suggest that MHC class I in Crocodylia underwent independent events of gene duplication, particularly in Crocodilidae. These findings enhance our understanding of MHC class I evolution and provide a preliminary framework for comparative studies of other non-avian reptiles as well as diversity assessment within Crocodylia.

Strong purifying selection in endogenous retroviruses in the saltwater crocodile (Crocodylus porosus) in the Northern Territory of Australia

Background

Endogenous retroviruses (ERVs) are remnants of exogenous retroviruses that have integrated into the nuclear DNA of a germ-line cell. Here we present the results of a survey into the ERV complement of Crocodylus porosus, the saltwater crocodile, representing 45 individuals from 17 sampling locations in the Northern Territory of Australia. These retroelements were compared with published ERVs from other species of Crocodylia (Crocodilians; alligators, caimans, gharials and crocodiles) as well as representatives from other vertebrates. This study represents one of the first in-depth studies of ERVs within a single reptilian species shedding light on the diversity of ERVs and proliferation mechanisms in crocodilians.

Results

Analyses of the retroviral pro-pol gene region have corroborated the presence of two major clades of ERVs in C. porosus and revealed 18 potentially functional fragments out of the 227 recovered that encode intact pro-pol ORFs. Interestingly, we have identified some patterns of diversification among those ERVs as well as a novel sequence that suggests the presence of an additional retroviral genus in C. porosus. In addition, considerable diversity but low genetic divergence within one of the C. porosus ERV lineages was identified.

Conclusions

We propose that the ERV complement of C. porosus has come about through a combination of recent infections and replication of ancestral ERVs. Strong purifying selection acting on these clades suggests that this activity is recent or still occurring in the genome of this species. The discovery of potentially functional elements is an interesting development that warrants further investigation.