Characterization of MHC class II genes from an ancient reptile lineage, Sphenodon (tuatara)

Genome-wide diversity and MHC characterisation in a critically endangered freshwater turtle susceptible to disease

Small, isolated populations are often vulnerable to increased inbreeding and genetic drift, both of which elevate the risk of extinction. The Bellinger River turtle (Myuchelys georgesi) is a critically endangered species endemic to a single river catchment in New South Wales, Australia. The only extant wild population, along with the breeding program, face significant threats from viral outbreaks, most notably a nidovirus outbreak in 2015 that led to a 90% population decline. To enhance our understanding of genomic characteristics in the species, including genome-wide and functional gene diversity, we re-sequenced, assembled, and analysed 31 re-sequenced genomes for pure M. georgesi (N = 31). We manually annotated the major histocompatibility complex (MHC), identifying five MHC class I and ten MHC class II genes and investigated genetic diversity across both classes in M. georgesi. Our results showed that genome-wide diversity is critically low in pure M. georgesi, contexualised through comparison with opportunistically sampled backcross animals-offspring of F1 hybrids (M. georgesi × Emydura macquarii) backcrossed to pure M. georgesi (N = 4). However, the variation observed within the core MHC region of pure M. georgesi, extending across scaffold 10, exceeded that of all other macrochromosomes. Additionally, no significant short-term changes in either genome-wide or immunogenetic diversity were detected following the 2015 nidovirus outbreak (before; N = 19, after; N = 12). Demographic history reconstructions indicated a sustained, long-term decline in effective population size since the last interglacial period, accompanied by more recent steep declines. These patterns suggested that prolonged isolation and reduced population size have significantly influenced the dynamics of genome-wide diversity. It is likely that contemporary stressors, including the recent nidovirus outbreak, are acting on an already genetically depleted population. This study offers new insights into genome-wide and immune gene diversity, including immune gene annotation data with broader implications for testudines. These findings provide crucial information to support future management strategies for the species.

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.

Avian MHC Evolution in the Era of Genomics: Phase 1.0

Birds are a wonderfully diverse and accessible clade with an exceptional range of ecologies and behaviors, making the study of the avian major histocompatibility complex (MHC) of great interest. In the last 20 years, particularly with the advent of high-throughput sequencing, the avian MHC has been explored in great depth in several dimensions: its ability to explain ecological patterns in nature, such as mating preferences; its correlation with parasite resistance; and its structural evolution across the avian tree of life. Here, we review the latest pulse of avian MHC studies spurred by high-throughput sequencing. Despite high-throughput approaches to MHC studies, substantial areas remain in need of improvement with regard to our understanding of MHC structure, diversity, and evolution. Recent studies of the avian MHC have nonetheless revealed intriguing connections between MHC structure and life history traits, and highlight the advantages of long-term ecological studies for understanding the patterns of MHC variation in the wild. Given the exceptional diversity of birds, their accessibility, and the ease of sequencing their genomes, studies of avian MHC promise to improve our understanding of the many dimensions and consequences of MHC variation in nature. However, significant improvements in assembling complete MHC regions with long-read sequencing will be required for truly transformative studies.

The effects of historical fragmentation on major histocompatibility complex class II β and microsatellite variation in the Aegean island reptile, Podarcis erhardii

The major histocompatibility complex (MHC) plays a key role in disease resistance and is the most polymorphic gene region in vertebrates. Although habitat fragmentation is predicted to lead to a loss in MHC variation through drift, the impact of other evolutionary forces may counter this effect. Here we assess the impact of selection, drift, migration, and recombination on MHC class II and microsatellite variability in 14 island populations of the Aegean wall lizard Podarcis erhardii. Lizards were sampled from islands within the Cyclades (Greece) formed by rising sea levels as the last glacial maximum approximately 20,000 before present. Bathymetric data were used to determine the area and age of each island, allowing us to infer the corresponding magnitude and timing of genetic bottlenecks associated with island formation. Both MHC and microsatellite variation were positively associated with island area, supporting the hypothesis that drift governs neutral and adaptive variation in this system. However, MHC but not microsatellite variability declined significantly with island age. This discrepancy is likely due to the fact that microsatellites attain mutation‐drift equilibrium more rapidly than MHC. Although we detected signals of balancing selection, recombination and migration, the effects of these evolutionary processes appeared negligible relative to drift. This study demonstrates how land bridge islands can provide novel insights into the impact of historical fragmentation on genetic diversity as well as help disentangle the effects of different evolutionary forces on neutral and adaptive diversity.

Major Histocompatibility Complex Genes Map to Two Chromosomes in an Evolutionarily Ancient Reptile, the Tuatara Sphenodon punctatus

Major histocompatibility complex (MHC) genes are a central component of the vertebrate immune system and usually exist in a single genomic region. However, considerable differences in MHC organization and size exist between different vertebrate lineages. Reptiles occupy a key evolutionary position for understanding how variation in MHC structure evolved in vertebrates, but information on the structure of the MHC region in reptiles is limited. In this study, we investigate the organization and cytogenetic location of MHC genes in the tuatara (Sphenodon punctatus), the sole extant representative of the early-diverging reptilian order Rhynchocephalia. Sequencing and mapping of 12 clones containing class I and II MHC genes from a bacterial artificial chromosome library indicated that the core MHC region is located on chromosome 13q. However, duplication and translocation of MHC genes outside of the core region was evident, because additional class I MHC genes were located on chromosome 4p. We found a total of seven class I sequences and 11 class II β sequences, with evidence for duplication and pseudogenization of genes within the tuatara lineage. The tuatara MHC is characterized by high repeat content and low gene density compared with other species and we found no antigen processing or MHC framework genes on the MHC gene-containing clones. Our findings indicate substantial differences in MHC organization in tuatara compared with mammalian and avian MHCs and highlight the dynamic nature of the MHC. Further sequencing and annotation of tuatara and other reptile MHCs will determine if the tuatara MHC is representative of nonavian reptiles in general.

De novo sequence assembly and characterisation of a partial transcriptome for an evolutionarily distinct reptile, the tuatara (Sphenodon punctatus)

Background

The tuatara (Sphenodon punctatus) is a species of extraordinary zoological interest, being the only surviving member of an entire order of reptiles which diverged early in amniote evolution. In addition to their unique phylogenetic placement, many aspects of tuatara biology, including temperature-dependent sex determination, cold adaptation and extreme longevity have the potential to inform studies of genome evolution and development. Despite increasing interest in the tuatara genome, genomic resources for the species are still very limited. We aimed to address this by assembling a transcriptome for tuatara from an early-stage embryo, which will provide a resource for genome annotation, molecular marker development and studies of development and adaptation in tuatara.

Results

We obtained 30 million paired-end 50 bp reads from an Illumina Genome Analyzer and assembled them with Velvet and Oases using a range of kmers. After removing redundancy and filtering out low quality transcripts, our transcriptome dataset contained 32911 transcripts, with an N50 of 675 and a mean length of 451 bp. Almost 50% (15965) of these transcripts could be annotated by comparison with the NCBI non-redundant (NR) protein database or the chicken, green anole and zebrafish UniGene sets. A scan of candidate genes and repetitive elements revealed genes involved in immune function, sex differentiation and temperature-sensitivity, as well as over 200 microsatellite markers.

Conclusions

This dataset represents a major increase in genomic resources for the tuatara, increasing the number of annotated gene sequences from just 60 to almost 16,000. This will facilitate future research in sex determination, genome evolution, local adaptation and population genetics of tuatara, as well as inform studies on amniote evolution.

Polymorphism and Balancing Selection of MHC Class II DAB Gene in 7 Selective Flounder (Paralichthys olivaceus) Families

In order to determine the genetic variation of the MHC class IIB exon2 allele in the offspring, 700 fry from seven families of Japanese flounder challenged with V. anguillarum were studied, and different mortality rates were found in those families. Five to ten surviving and dead fry from each of the seven families were selected to study the MHC class II B exon2 gene with PCR and a direct sequencing method. One hundred and sixteen different exon2 sequences were found and 116 different alleles were identified, while a minimum of four loci were revealed in the MHC class II B exon2 gene. The ratio (d_N/d_S) of nonsynonymous substitution (d_N) to synonymous substitutions (d_S) in the peptide-binding region (PBR) of the MHC class IIB gene was 6.234, which indicated that balancing selection is acting on the MHC class IIB genes. The MHC IIB alleles were thus being passed on to their progeny. Some alleles were significantly more frequent in surviving than dead individuals. All together our data suggested that the alleles Paol-DAB*4301, Paol-DAB*4601, Paol-DAB*4302, Paol-DAB*3803, and Paol-DAB*4101 were associated with resistance to V. anguillarum in flounder.

Genome evolution in Reptilia, the sister group of mammals

The genomes of birds and nonavian reptiles (Reptilia) are critical for understanding genome evolution in mammals and amniotes generally. Despite decades of study at the chromosomal and single-gene levels, and the evidence for great diversity in genome size, karyotype, and sex chromosome diversity, reptile genomes are virtually unknown in the comparative genomics era. The recent sequencing of the chicken and zebra finch genomes, in conjunction with genome scans and the online publication of the Anolis lizard genome, has begun to clarify the events leading from an ancestral amniote genome--predicted to be large and to possess a diverse repeat landscape on par with mammals and a birdlike sex chromosome system--to the small and highly streamlined genomes of birds. Reptilia exhibit a wide range of evolutionary rates of different subgenomes and, from isochores to mitochondrial DNA, provide a critical contrast to the genomic paradigms established in mammals.

Characterization and evolution of MHC class II B genes in Galápagos marine iguanas (Amblyrhynchus cristatus)

Major histocompatibility complex (MHC) class II molecules play a key role in the adaptive immune system of vertebrates. Class II B genes appear to evolve in a very different manner in mammals and birds. Orthology is commonly observed among mammal loci, while genes tend to cluster phylogenetically within bird species. Here we present class II B data from a representative of another major group of amniotes, the squamates (i.e. lizards, snakes, amphisbaenians), with the ultimate goal of placing mammalian and avian MHC evolution into a broader context. In this study, eight class II B cDNA sequences were obtained from the Galápagos marine iguana (Amblyrhynchus cristatus) which were divided into five locus groups, Amcr-DAB1 through -DAB5, based on similarities along most of the coding and noncoding portions of the transcribed gene. All marine iguana sequences were monophyletic with respect to class II genes from other vertebrates indicating that they originated from a common ancestral locus after squamates split from other reptiles. The beta-1 domain, which is involved in antigen binding, exhibited signatures of positive selection as well as interlocus gene conversion in both long and short tracts-a pattern also observed in birds and fish, but not in mammals. On the other hand, the beta-2 domain was divergent between gene groups, which is characteristic of mammals. Based on these results, we preliminarily show that squamate class II B genes have been shaped by a unique blend of evolutionary forces that have been observed in differing degrees in other vertebrates.

Species-specific evolution of class I MHC genes in iguanas (order: Squamata; subfamily: Iguaninae)

Over the last few decades, the major histocompatibility complex (MHC) has emerged as a model for understanding the influence of natural selection on genetic diversity in populations as well as for investigating the genetic basis of host resistance to pathogens. However, many vertebrate taxa remain underrepresented in the field of MHC research, preventing its application to studies of disease, evolution, and conservation genetics in these groups. This is particularly true for squamates, which are by far the most diversified order of non-avian reptiles but have not been the subject of any recent MHC studies. In this paper, we present MHC class I complementary DNA data from three squamate species in the subfamily Iguaninae (iguanas): the Galápagos marine iguana (Amblyrhynchus cristatus), the Galápagos land iguana (Conolophus subcristatus), and the green iguana (Iguana iguana). All sequences obtained are related to the few published class I genes from other squamates. There is evidence for multiple loci in each species, and the conserved alpha-3 domain appears to be evolving in a species-specific manner. Conversely, there is some indication of shared polymorphism between species in the peptide-binding alpha-1 and alpha-2 domains, suggesting that these two regions have different phylogenetic histories. The great similarity between alpha-3 sequences in marine iguanas in particular suggests that concerted evolution is acting to homogenize class I loci within species. However, while less likely, the data are also compatible with a birth and death model of evolution.

Two patterns of variation among MHC class I loci in Tuatara (Sphenodon punctatus)

The genes of the major histocompatibility complex (MHC) are a central component of the immune system in vertebrates and have become important markers of functional, fitness-related genetic variation. We have investigated the evolutionary processes that generate diversity at MHC class I genes in a large population of an archaic reptile species, the tuatara (Sphenodon punctatus), found on Stephens Island, Cook Strait, New Zealand. We identified at least 2 highly polymorphic (UA type) loci and one locus (UZ) exhibiting low polymorphism. The UZ locus is characterized by low nucleotide diversity and weak balancing selection and may be either a nonclassical class I gene or a pseudogene. In contrast, the UA-type alleles have high nucleotide diversity and show evidence of balancing selection at putative peptide-binding sites. Twenty-one different UA-type genotypes were identified among 26 individuals, suggesting that the Stephens Island population has high levels of MHC class I variation. UA-type allelic diversity is generated by a mixture of point mutation and gene conversion. As has been found in birds and fish, gene conversion obscures the genealogical relationships among alleles and prevents the assignment of alleles to loci. Our results suggest that the molecular mechanisms that underpin MHC evolution in nonmammals make locus-specific amplification impossible in some species.

Characterization, polymorphism, and evolution of MHC class II B genes in birds of prey

During the last decade, the major histocompatibility complex (MHC) has received much attention in the fields of evolutionary and conservation biology because of its potential implications in many biological processes. New insights into the gene structure and evolution of MHC genes can be gained through study of additional lineages of birds not yet investigated at the genomic level. In this study, we characterized MHC class II B genes in five families of birds of prey (Accipitridae, Pandionidae, Strigidae, Tytonidae, and Falconidae). Using PCR approaches, we isolated genomic MHC sequences up to 1300 bp spanning exons 1 to 3 in 26 representatives of each raptor lineage, finding no stop codons or frameshift mutations in any coding region. A survey of diversity across the entirety of exon 2 in the lesser kestrel Falco naumanni reported 26 alleles in 21 individuals. Bayesian analysis revealed 21 positively selected amino acid sites, which suggests that the MHC genes described here are functional and probably expressed. Finally, through interlocus comparisons and phylogenetic analysis, we also discuss genetic evidence for concerted and transspecies evolution in the raptor MHC.

Molecular characterization of MHC class II in a nonmodel anuran species, the fire-bellied toad Bombina bombina

While the anuran Xenopus comprises one of the best characterized nonmammalian taxa regarding the major histocompatibility complex (MHC), the organization of this gene complex has never been studied in other anurans, and information on amphibian MHC (other than Xenopus) is generally very scarce. Here, we describe the characterization of the first MHC class II B cDNA sequences from a nonmodel anuran species, the European fire-bellied toad (Bombina bombina). We isolated two transcript sequences differing substantially in amino acid composition and length within the beta2 domain. To investigate the variability of the peptide binding region in this species, we sequenced a 158-bp large fragment from wild B. bombina (n = 20) and identified eight distinct alleles. All substitutions but one were nonsynonymous, and many of the highly polymorphic sites corresponded with amino acid positions known to be involved in antigen binding. The level of variation we found in B. bombina was similar compared to that previously found in a comparable sample of a wild urodelan species, Ambystoma tigrinum, and to that found in Xenopus laevis. Based on the cDNA data and the individual's allelic diversity, we conclude that Bombina possesses at least two class II B loci. With our new beta1 primers, we were able to generate sequences in other species of anurans. We provide here a first phylogenetic analysis of this gene in amphibians.

Combining immunological and molecular data to assess phylogenetic relations of some Greek Podarcis species

Most recent molecular studies revealed the phylogeny of Greek Podarcis species, which for years remained elusive, due to discordant data produced from various chromosomal, complement fixation and protein studies. In this report, we analyzed cellular immune responses of spleen-derived lymphocytes from six allopatric Podarcis species encountered in Greece, by assessing two-way mixed lymphocyte reaction (MLR)-induced proliferation. On the basis of stimulation indices (S.I.) as determined from cultures set up from xenogeneic splenocytes coincubated in pairs, we generated a phylogenetic tree, fully consistent with the phylogenetic relationships of Podarcis as determined by parallel analyses based on partial mitochondrial (mt) DNA sequences. Although the exact mechanisms triggering lymphocyte responses in lizard two-way xenogeneic MLR are not fully understood, our results show the potential use of cell-mediated immune responses as an additional approach to mtDNA analysis, for species delimitation within specific lizard taxa.