The importance of immune gene variability (MHC) in evolutionary ecology and conservation

Immune targeting of filarial glutaredoxin through a multi-epitope peptide-based vaccine: A reverse vaccinology approach

Despite the efforts of global programme to eliminate lymphatic filariasis (GPELF), the threat of lymphatic filariasis (LF) still looms over humanity in terms of long-term disabilities, and morbidities across the globe. In light of this situation, investigators have chosen to focus on the development of immunotherapeutics targeting the physiologically important filarial-specific proteins. Glutaredoxin (16.43 kDa) plays a pivotal role in filarial redox biology, serving as a vital contributor. In the context of the intra-host survival of filarial parasites, this antioxidant helps in mitigating the oxidative stress imposed by the host immune system. Given its significant contribution, the development of a vaccine targeting glutaredoxin holds promise as a new avenue for achieving a filaria-free world. Herein, multi-epitope-based vaccine was designed using advanced immunoinformatics approach. Initially, 4B-cell epitopes and 6 T-cell epitopes (4 MHC I and 2 MHC II) were identified from the 146 amino acid long sequence of glutaredoxin of the human filarid, Wuchereria bancrofti. Subsequent clustering of these epitopes with linker peptides finalized the vaccine structure. To boost TLR-mediated innate immunity, TLR-specific adjuvants were incorporated into the designed vaccine. After that, experimental analyses confirm the designed vaccine, Vac4 as anefficient ligand of human TLR5 to elicit protective innate immunity against filarial glutaredoxin. Immune simulation further demonstrated abundant levels of IgG and IgM as crucial contributors in triggering vaccine-induced adaptive responses in the recipients. Hence, to facilitate the validation of immunogenicity of the designed vaccine, Vac4 was cloned in silico in pET28a(+) expression vector for recombinant production. Taken together, our findings suggest that vaccine-mediated targeting of filarial glutaredoxin could be a future option for intervening LF on a global scale.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Organisation and evolution of the major histocompatibility complex class I genes in cetaceans

A quarter of marine mammals are at risk of extinction, with disease and poor habitat quality contributing to population decline. Investigation of the Major Histocompatibility Complex (MHC) provides insight into species' capacity to respond to immune and environmental challenges. The eighteen available cetacean chromosome level genomes were used to annotate MHC Class I loci, and to reconstruct the phylogenetic relationship of the described loci. The highest number of loci was observed in the striped dolphin (Stenella coeruleoalba), while the least was observed in the pygmy sperm whale (Kogia breviceps) and rough toothed dolphin (Steno bredanensis). Of the species studied, Mysticetes had the most pseudogenes. Evolutionarily, MHC Class I diverged before the speciation of cetaceans. Yet, locus one was genomically and phylogenetically similar in many species, persisting over evolutionary time. This characterisation of MHC Class I in cetaceans lays the groundwork for future population genetics and MHC expression studies.

Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

No evidence for a role of MHC class II genotype in the chemical encoding of heterozygosity and relatedness in Antarctic fur seals

Despite decades of research, surprisingly little is known about the mechanism(s) by which an individual's genotype is encoded in odour. Many studies have focused on the role of the major histocompatibility complex (MHC) owing to its importance for survival and mate choice. However, the salience of MHC-mediated odours compared to chemicals influenced by the rest of the genome remains unclear, especially in wild populations where it is challenging to quantify and control for the effects of the genomic background. We addressed this issue in Antarctic fur seals by analysing skin swabs together with full-length MHC DQB II exon 2 sequences and data from 41 genome-wide distributed microsatellites. We did not find any effects of MHC relatedness on chemical similarity and there was also no relationship between MHC heterozygosity and chemical diversity. However, multilocus heterozygosity showed a significant positive association with chemical diversity, even after controlling for MHC heterozygosity. Our results appear to rule out a dominant role of the MHC in the chemical encoding of genetic information in a wild vertebrate population and highlight the need for genome-wide approaches to elucidate the mechanism(s) and specific genes underlying genotype-odour associations.

MHC-I alleles mediate clearance and antibody response to the zoonotic Lassa virus in Mastomys rodent reservoirs

West African Mastomys rodents are the primary reservoir of the zoonotic Lassa virus (LASV). The virus causes haemorrhagic Lassa fever and considerable mortality in humans. To date, the role of Mastomys immunogenetics in resistance to, and persistence of, LASV infections is largely unknown. Here, we investigated the role of Major Histocompatibility Complex class I (MHC-I) on LASV infection status (i.e., active vs. cleared infection, determined via PCR and an immunofluorescence assay on IgG antibodies, respectively) in Mastomys natalensis and M. erythroleucus sampled within southwestern Nigeria. We identified more than 190 and 90 MHC-I alleles by Illumina high throughput-sequencing in M. natalensis and M. erythroleucus, respectively, with different MHC allele compositions and frequencies between LASV endemic and non-endemic sites. In M. natalensis, the MHC allele ManaMHC-I*006 was negatively associated with active infections (PCR-positive) and positively associated with cleared infections (IgG-positive) simultaneously, suggesting efficient immune responses that facilitate LASV clearance in animals carrying this allele. Contrarily, alleles ManaMHC-I*008 and ManaMHC-I*021 in M. natalensis, and MaerMHC-I*008 in M. erythroleucus, were positively associated with active infection, implying susceptibility. Alleles associated with susceptibility shared a glutamic acid at the positively selected codon 57, while ManaMHC-I*006 featured an arginine. There was no link between number of MHC alleles per Mastomys individual and LASV prevalence. Thus, specific alleles, but not MHC diversity per se, seem to mediate antibody responses to viremia. We conclude that co-evolution with LASV likely shaped the MHC-I diversity of the main LASV reservoirs in southwestern Nigeria, and that information on reservoir immunogenetics may hold insights into transmission dynamics and zoonotic spillover risks.

Recent selection created distinctive variability patterns on MHC class II loci in three dolphin species from the Mediterranean Sea

The major histocompatibility complex (MHC) includes highly polymorphic genes involved in antigen presentation, which is crucial for adaptive immune response. They represent fitness related genetic markers particularly informative for populations exposed to environmental challenges. Here we analyse the diversity and evolutionary traits of MHC class II DQA and DQB genes in the dolphins Stenella coeruleoalba and Grampus griseus from the Mediterranean Sea. We found substantial nucleotide and functional diversity, as well as strong evidence of balancing selection indicated by allele and supertype frequencies, Tajima's D statistics and dN/dS tests. The Risso's dolphin, considered the least abundant in the region, showed the effect of divergent allele advantage at the nucleotide and functional-peptide levels. An outstanding polymorphism was found in the striped dolphin, particularly intriguing in the DQA gene where the Ewens-Watterson test detected a selection sweep that occurred in recent history. We hypothesize that morbillivirus, which has recurrently invaded Mediterranean populations over the last decades, exerted the detected selective pressure.

AmpliSAS and AmpliHLA: Web Server and Local Tools for MHC Typing of Non-model Species and Human Using NGS Data

AmpliSAS and AmpliHLA are tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human specific version; it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Both tools are available in AmpliSAT web-server as well as scripts for local/server installation. Here we describe the installation and deployment of AmpliSAS and AmpliHLA Perl scripts and dependencies on a local or a server computer. We will show how to run them in the command line using as examples four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

Genetic Diversity and Sequence Conservation of Peptide-Binding Regions of MHC Class I Genes in Pig, Cattle, Chimpanzee, and Human

Comparative analyses of MHC gene diversity and evolution across different species could offer valuable insights into the evolution of MHC genes. Intra- and inter-species sequence diversity and conservation of 12 classical major histocompatibility complex (MHC) class I genes from cattle, chimpanzees, pigs, and humans was analyzed using 20 representative allelic groups for each gene. The combined analysis of paralogous loci for each species revealed that intra-locus amino-acid sequence variations in the peptide-binding region (PBR) of MHC I genes did not differ significantly between species, ranging from 8.44% for SLA to 10.75% for BoLA class I genes. In contrast, intraspecies differences in the non-PBRs of these paralogous genes were more pronounced, varying from 4.59% for SLA to 16.89% for HLA. Interestingly, the Shannon diversity index and rate of nonsynonymous substitutions for PBR were significantly higher in SLA and BoLA than those in Patr and HLA. Analysis of peptide-binding pockets across all analyzed MHC class I genes of the four species indicated that pockets A and E showed the lowest and highest diversity, respectively. The estimated divergence times suggest that primate and artiodactyl MHC class I genes diverged 60.41 Mya, and BoLA and SLA genes diverged 35.34 Mya. These results offer new insights into the conservation and diversity of MHC class I genes in various mammalian species.

Immunogenetics, sylvatic plague and its vectors: insights from the pathogen reservoir Mastomys natalensis in Tanzania

Yersinia pestis is a historically important vector-borne pathogen causing plague in humans and other mammals. Contemporary zoonotic infections with Y. pestis still occur in sub-Saharan Africa, including Tanzania and Madagascar, but receive relatively little attention. Thus, the role of wildlife reservoirs in maintaining sylvatic plague and spillover risks to humans is largely unknown. The multimammate rodent Mastomys natalensis is the most abundant and widespread rodent in peri-domestic areas in Tanzania, where it plays a major role as a Y. pestis reservoir in endemic foci. Yet, how M. natalensis' immunogenetics contributes to the maintenance of plague has not been investigated to date. Here, we surveyed wild M. natalensis for Y. pestis vectors, i.e., fleas, and tested for the presence of antibodies against Y. pestis using enzyme-linked immunosorbent assays (ELISA) in areas known to be endemic or without previous records of Y. pestis in Tanzania. We characterized the allelic and functional (i.e., supertype) diversity of the major histocompatibility complex (MHC class II) of M. natalensis and investigated links to Y. pestis vectors and infections. We detected antibodies against Y. pestis in rodents inhabiting both endemic areas and areas considered non-endemic. Of the 111 nucleotide MHC alleles, only DRB*016 was associated with an increased infestation with the flea Xenopsylla. Surprisingly, we found no link between MHC alleles or supertypes and antibodies of Y. pestis. Our findings hint, however, at local adaptations towards Y. pestis vectors, an observation that more exhaustive sampling could unwind in the future.

Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus

Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.

Gut Microbiome Profiling of the Endangered Southern Greater Glider (Petauroides volans) after the 2019-2020 Australian Megafire

Studying the gut microbiome can provide valuable insights into animal health and inform the conservation management of threatened wildlife. Gut microbiota play important roles in regulating mammalian host physiology, including digestion, energy metabolism and immunity. Dysbiosis can impair such physiological processes and compromise host health, so it is essential that the gut microbiome be considered in conservation planning. The southern greater glider (Petauroides volans) is an endangered arboreal marsupial that faced widespread habitat fragmentation and population declines following the 2019-2020 Australian bushfire season. This study details baseline data on the gut microbiome of this species. The V3-V4 region of the 16S rRNA gene was amplified from scats collected from individuals inhabiting burnt and unburnt sites across southeastern Australia and sequenced to determine bacterial community composition. Southern greater glider gut microbiomes were characterised by high relative abundances of Firmicutes and Bacteroidota, which is consistent with that reported for other marsupial herbivores. Significant differences in gut microbial diversity and community structure were detected among individuals from different geographic locations. Certain microbiota and functional orthologues were also found to be significantly differentially abundant between locations. The role of wildfire in shaping southern greater glider gut microbiomes was shown, with some significant differences in the diversity and abundance of microbiota detected between burnt and unburnt sites. Overall, this study details the first data on greater glider (Petauroides) gut microbiomes, laying the foundation for future studies to further explore relationships between microbial community structure, environmental stressors and host health.

Diversity and evolution of the MHC class II DRB gene in the Capra sibirica experienced a demographic fluctuation in China

The major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates, and their proteins play a critical role in adaptive immunity for defense against a variety of pathogens. MHC diversity was lost in many species after experiencing a decline in size. To understand the variation and evolution of MHC genes in the Siberian ibex, Capra sibirica, which has undergone a population decline, we analyzed the variation of the second exon of MHC class II DRB genes in samples collected from five geographic localities in Xinjiang, China, that belong to three diverged mitochondrial clades. Consequently, we identified a total of 26 putative functional alleles (PFAs) with 260 bp in length from 43 individuals, and found one (for 27 individuals) to three (for 5 individuals) PFAs per individual, indicating the presence of one or two DRB loci per haploid genome. The Casi-DRB1*16 was the most frequently occurring PFA, Casi-DRB1*22 was found in only seven individuals, 14 PFAs occurred once, 7 PFAs twice, implying high frequency of rare PFAs. Interestingly, more than half (15) of the PFAs were specific to clade I, only two and three PFAs were specific to clades II and III, respectively. So, we assume that the polygamy and sexual segregation nature of this species likely contributed to the allelic diversity of DRB genes. Genetic diversity indices showed that PFAs of clade II were lower in nucleotide, amino acid, and supertype diversity compared to those of the other two clades. The pattern of allele sharing and FST values between the three clades was to some extent in agreement with the pattern observed in mitochondrial DNA divergence. In addition, recombination analyses revealed no evidence for significant signatures of recombination events. Alleles shared by clades III and the other two clades diverged 6 million years ago, and systematic neighbor grids showed Trans-species polymorphism. Together with the PAML and MEME analyses, the results indicated that the DRB gene in C. sibirica evolved under balancing and positive selection. However, by comparison, it can be clearly seen that different populations were under different selective pressures. Our results are valuable in understanding the diversity and evolution of the DRB gene in a mountain living C. sibirica and in making decisions on future long-term protection strategies.

The Impact of Increased Homozygosity on Human Fertility: A Comprehensive Review

This comprehensive review explores the multifaceted relationship between increased homozygosity and human fertility, delving into the genetic, ethical, cultural, and public health dimensions of this complex phenomenon. Homozygosity, characterized by identical alleles at specific gene loci, can result from consanguineous marriages, genetic drift, and population isolation. The review highlights key findings, including the heightened risk of recessive genetic disorders, the implications for immune system diversity, and the influence on complex traits and diseases. It underscores the critical role of genetic counseling in addressing these consequences, considering the ethical implications, and respecting cultural practices. The delicate balance between genetic diversity and cultural norms is emphasized, calling for increased awareness and community engagement. Looking ahead, the review suggests emerging technologies, longitudinal studies, and interdisciplinary research as crucial avenues for further exploration, with the ultimate goal of informing effective public health policies and interventions that safeguard genetic diversity and cultural traditions for future generations.

Assessing micro-macroparasite selective pressures and anthropogenic disturbance as drivers of immune gene diversity in a Neotropical wild cat

Anthropogenic environmental change is reducing available habitat for wild species, providing novel selection pressures such as infectious diseases and causing species to interact in new ways. The potential for emerging infectious diseases and zoonoses at the interface between humans, domestic animals, and wild species is a key global concern. In vertebrates, diversity at the major histocompatibility complex MHC is critical to disease resilience, and its study in wild populations provides insights into eco-evolutionary dynamics that human activities alter. In natural populations, variation at MHC loci is partly maintained by balancing selection, driven by pathogenic selective pressures. We hypothesize that MHC genetic diversity differs between guigna populations inhabiting human-dominated landscapes (higher pathogen pressures) versus more natural habitats (lower pathogen pressures). We predict that MHC diversity in guignas would be highest in human-dominated landscapes compared with continuous forest habitats. We also expected to find higher MHC diversity in guignas infected with micro and macro parasites (higher parasite load) versus non infected guignas. We characterized for the first time the genetic diversity at three MHC class I and II exons in 128 wild guignas (Leopardus guigna) across their distribution range in Chile (32-46° S) and Argentina, representing landscapes with varying levels of human disturbance. We integrated MHC sequence diversity with multiple measures of anthropogenic disturbance and both micro and macro parasite infection data. We also assessed signatures of positive selection acting on MHC genes. We found significantly higher MHC class I diversity in guignas inhabiting landscapes where houses were present, and with lower percentage of vegetation cover, and also in animals with more severe cardiorespiratory helminth infection (richness and intensity) and micro-macroparasite co-infection. This comprehensive, landscape-level assessment further enhances our knowledge on the evolutionary dynamics and adaptive potential of vertebrates in the face of emerging infectious disease threats and increasing anthropogenic impacts.

HLA-G*14 bp indel variant in autism spectrum disorder in a population from southern Brazil

Altered immune response during pregnancy has been associated with ASD susceptibility. HLA-G is expressed by the trophoblast at the maternal/fetal interface and induces allogenic tolerance toward the fetus. A 14-bp insertion in the HLA-G 3'UTR (rs371194629) was associated with reduced levels of HLA-G. We aimed to assess the influence of the HLA-G*14 bp indel variant in ASD susceptibility and symptomatology in a Brazilian admixed sample. The insertion genotype (14 bp+/14 bp+) was firstly associated with hetero aggression, but statistical significance was lost after correction (p = 0.035, pcorrected = 0.35). No association between the HLA-G variant and susceptibility to ASD or differential clinical manifestations were observed.

A suite of selective pressures supports the maintenance of alleles of a Drosophila immune peptide

The innate immune system provides hosts with a crucial first line of defense against pathogens. While immune genes are often among the fastest evolving genes in the genome, in Drosophila, antimicrobial peptides (AMPs) are notable exceptions. Instead, AMPs may be under balancing selection, such that over evolutionary timescales multiple alleles are maintained in populations. In this study, we focus on the Drosophila antimicrobial peptide Diptericin A, which has a segregating amino acid polymorphism associated with differential survival after infection with the Gram-negative bacteria Providencia rettgeri. Diptericin A also helps control opportunistic gut infections by common Drosophila gut microbes, especially those of Lactobacillus plantarum. In addition to genotypic effects on gut immunity, we also see strong sex-specific effects that are most prominent in flies without functional diptericin A. To further characterize differences in microbiomes between different diptericin genotypes, we used 16S metagenomics to look at the microbiome composition. We used both lab reared and wild caught flies for our sequencing and looked at overall composition as well as the differential abundance of individual bacterial families. Overall, we find flies that are homozygous serine for diptericin A are better equipped to survive a systemic infection from P. rettgeri, but in general homozygous arginine flies have a longer lifespan after being fed common gut commensals. Our results suggest a possible mechanism for the maintenance of genetic variation of diptericin A through the complex interactions of sex, systemic immunity, and the maintenance of the gut microbiome.

Local translation in microglial processes is required for efficient phagocytosis

Neurons, astrocytes and oligodendrocytes locally regulate protein translation within distal processes. Here, we tested whether there is regulated local translation within peripheral microglial processes (PeMPs) from mouse brain. We show that PeMPs contain ribosomes that engage in de novo protein synthesis, and these are associated with transcripts involved in pathogen defense, motility and phagocytosis. Using a live slice preparation, we further show that acute translation blockade impairs the formation of PeMP phagocytic cups, the localization of lysosomal proteins within them, and phagocytosis of apoptotic cells and pathogen-like particles. Finally, PeMPs severed from their somata exhibit and require de novo local protein synthesis to effectively surround pathogen-like particles. Collectively, these data argue for regulated local translation in PeMPs and indicate a need for new translation to support dynamic microglial functions.

Conserved Evolution of MHC Supertypes among Japanese Frogs Suggests Selection for Bd Resistance

The chytrid fungus Batrachochytrium dendrobatidis (Bd) is a major threat to amphibians, yet there are no reports of major disease impacts in East Asian frogs. Genetic variation of the major histocompatibility complex (MHC) has been associated with resistance to Bd in frogs from East Asia and worldwide. Using transcriptomic data collated from 11 Japanese frog species (one individual per species), we isolated MHC class I and IIb sequences and validated using molecular cloning. We then compared MHC from Japanese frogs and other species worldwide, with varying Bd susceptibility. Supertyping analysis, which groups MHC alleles based on physicochemical properties of peptide binding sites, identified that all examined East Asian frogs contained at least one MHC-IIb allele belonging to supertype ST-1. This indicates that, despite the large divergence times between some Japanese frogs (up to 145 million years), particular functional properties in the peptide binding sites of MHC-II are conserved among East Asian frogs. Furthermore, preliminary analysis using NetMHCIIpan-4.0, which predicts potential Bd-peptide binding ability, suggests that MHC-IIb ST-1 and ST-2 have higher overall peptide binding ability than other supertypes, irrespective of whether the peptides are derived from Bd, other fungi, or bacteria. Our findings suggest that MHC-IIb among East Asian frogs may have co-evolved under the same selective pressure. Given that Bd originated in this region, it may be a major driver of MHC evolution in East Asian frogs.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Influence of HLA-A, -B, -DR Polymorphisms on the Severity of COVID-19: A Case-Control Study in the Iranian Population

BACKGROUND

As an emerging pandemic disease, COVID-19 encompasses a spectrum of clinical diagnoses, from the common cold to severe respiratory syndrome. Considering the shreds of evidence demonstrating the relationship between human leukocyte antigen (HLA) allele diversity and infectious disease susceptibility, this study was conducted to determine the association of HLA alleles with COVID-19 severity in Iranian subjects.

METHODS

In this case-control study, a total of 200 unrelated individuals (consisting of 100 people with severe COVID-19 and an average age of 55.54 as the case group, and 100 patients with mild COVID-19 with an average age of 48.97 as the control group) were recruited, and HLA typing (Locus A, B, and DR) was performed using the Olerup sequence-specific oligonucleotide (SSO) HLA-typing kit.

RESULTS

Our results showed that HLA-A*11 and HLA-DRB1*14 alleles were more frequently observed in severe COVID-19 cases, while HLA-B*52 was more common in mild cases, which was in agreement with some previous studies.

CONCLUSION

Our results confirmed the evidence for the association of HLA alleles with COVID-19 outcomes. We found that HLA-A*11 and HLA-DRB1*14 alleles may be susceptibility factors for severe COVID-19, while HLA-B*52 may be a protective factor. These findings provide new insight into the pathogenesis of COVID-19 and help patient management.

Composition and function of the Galapagos penguin gut microbiome vary with age, location, and a putative bacterial pathogen

Microbial colonization plays a direct role in host health. Understanding the ecology of the resident microbial community for a given host species is thus an important step for detecting population vulnerabilities like disease. However, the idea of integrating microbiome research into conservation is still relatively new, and wild birds have received less attention in this field than mammals or domesticated animals. Here we examine the composition and function of the gut microbiome of the endangered Galapagos penguin (Spheniscus mendiculus) with the goals of characterizing the normal microbial community and resistome, identifying likely pathogens, and testing hypotheses of structuring forces for this community based on demographics, location, and infection status. We collected fecal samples from wild penguins in 2018 and performed 16S rRNA gene sequencing and whole genome sequencing (WGS) on extracted DNA. 16S sequencing revealed that the bacterial phyla Fusobacteria, Epsilonbacteraeota, Firmicutes, and Proteobacteria dominate the community. Functional pathways were computed from WGS data, showing genetic functional potential primarily focused on metabolism-amino acid metabolism, carbohydrate metabolism, and energy metabolism are the most well-represented functional groups. WGS samples were each screened for antimicrobial resistance, characterizing a resistome made up of nine antibiotic resistance genes. Samples were screened for potential enteric pathogens using virulence factors as indicators; Clostridium perfringens was revealed as a likely pathogen. Overall, three factors appear to be shaping the alpha and beta diversity of the microbial community: penguin developmental stage, sampling location, and C. perfringens. We found that juvenile penguins have significantly lower alpha diversity than adults based on three metrics, as well as significantly different beta diversity. Location effects are minimal, but one site has significantly lower Shannon diversity than the other primary sites. Finally, when samples were grouped by C. perfringens virulence factors, we found dramatic changes in beta diversity based on operational taxonomic units, protein families, and functional pathways. This study provides a baseline microbiome for an endangered species, implicates both penguin age and the presence of a potential bacterial pathogen as primary factors associated with microbial community variance, and reveals widespread antibiotic resistance genes across the population.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

The Association of HLA Alleles and Haplotypes with Age of Disease Onset in Pakistani Psoriatic Patients

INTRODUCTION

The human leukocyte antigen (HLA) region on chromosome 6p21 is well known to carry the most important genetic factors in susceptibility to psoriasis. Different HLA alleles and haplotypes have been reported to be associated with psoriasis in different populations. Psoriasis has a variable age of onset and, based on this, it can be classified into two types; type I with age of onset before 40 years of age and type II with age of onset after 40 years of age. The objective of this study was to determine the association of HLA class I and class II alleles and haplotypes with disease and stratification using age of onset in Pakistani psoriatic patients.

METHODS

A group of 603 individuals (326 cases and 277 controls) were analyzed for HLA class I and II alleles and haplotype association by sequence specific PCR. The association was further analyzed according to the age of onset of the patients.

RESULTS

We found that HLA alleles B*57 and Cw*06:02, DQB1*03:03:02 are strongly associated with early onset psoriasis, while alleles B*15, DRB1*13:02 and DQB1*03:03:02 are associated with late-onset psoriasis. Cw*06:02 allele was not associated with late-onset psoriasis patients. Allele DQB1*03:03:02 had the highest odds ratio in all patients. We found a novel association specifically with late-onset psoriasis samples with the haplotype HLA-A*11; B*15; Cw*04; DRB1*15; DQB1*05 (Pc = 3.60 × 10-7). We also found strong association with previously reported extended haplotype EH-57.1: HLA-B*57; Cw*06:02; DRB1*07:01; DQB1*03:03:02 in all our patients (Pc = 8.34 × 10-07).

CONCLUSION

Our results show that different HLA class I and II alleles and haplotypes are associated with psoriasis at different age of onset. In this study, we have reported novel alleles and haplotype association with late-onset psoriasis. Our data confirm the previous strong associations with HLA alleles and haplotypes and also reports novel alleles and haplotype association in Pakistani psoriasis patients.

Genetic slippage after sex maintains diversity for parasite resistance in a natural host population

Although parasite-mediated selection is a major driver of host evolution, its influence on genetic variation for parasite resistance is not yet well understood. We monitored resistance in a large population of the planktonic crustacean Daphnia magna over 8 years, as it underwent yearly epidemics of the bacterial pathogen Pasteuria ramosa. We observed cyclic dynamics of resistance: Resistance increased throughout the epidemics, but susceptibility was restored each spring when hosts hatched from sexual resting stages. Host resting stages collected across the year showed that largely resistant host populations can produce susceptible sexual offspring. A genetic model of resistance developed for this host-parasite system, based on multiple loci and strong epistasis, is in partial agreement with our findings. Our results reveal that, despite strong selection for resistance in a natural host population, genetic slippage after sexual reproduction can be a strong factor for the maintenance of genetic diversity of host resistance.

Genome-wide scan for potential CD4+ T-cell vaccine candidates in Candida auris by exploiting reverse vaccinology and evolutionary information

Candida auris is a globally emerging fungal pathogen responsible for causing nosocomial outbreaks in healthcare associated settings. It is known to cause infection in all age groups and exhibits multi-drug resistance with high potential for horizontal transmission. Because of this reason combined with limited therapeutic choices available, C. auris infection has been acknowledged as a potential risk for causing a future pandemic, and thus seeking a promising strategy for its treatment is imperative. Here, we combined evolutionary information with reverse vaccinology approach to identify novel epitopes for vaccine design that could elicit CD4+ T-cell responses against C. auris. To this end, we extensively scanned the family of proteins encoded by C. auris genome. In addition, a pathogen may acquire substitutions in epitopes over a period of time which could cause its escape from the immune response thus rendering the vaccine ineffective. To lower this possibility in our design, we eliminated all rapidly evolving genes of C. auris with positive selection. We further employed highly conserved regions of multiple C. auris strains and identified two immunogenic and antigenic T-cell epitopes that could generate the most effective immune response against C. auris. The antigenicity scores of our predicted vaccine candidates were calculated as 0.85 and 1.88 where 0.5 is the threshold for prediction of fungal antigenic sequences. Based on our results, we conclude that our vaccine candidates have the potential to be successfully employed for the treatment of C. auris infection. However, in vivo experiments are imperative to further demonstrate the efficacy of our design.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Microsatellite markers of the major histocompatibility complex genomic region of domestic camels

We identified and characterized 11 polymorphic microsatellite markers suitable for routine testing (three in the MHC class I sub-region, four in MHC class II and four in the MHC class III sub-region) of dromedaries and Bactrian camels. In total, 38 dromedaries and 33 Bactrian camels were genotyped, and interspecific differences were observed in the numbers of alleles and in allelic frequencies, as well as in the observed heterozygosity. These loci may be used as markers to study the adaptive genetic diversity of the MHC region in Old World camels.

Performance Comparison of Different Approaches in Genotyping MHC-DRB: The Contrast between Single-Locus and Multi-Locus Species

Major histocompatibility complex (MHC) genes are widely recognised as valuable markers for wildlife genetic studies given their extreme polymorphism and functional importance in fitness-related traits. Newly developed genotyping methods, which rely on the use of next-generation sequencing (NGS), are gradually replacing traditional cloning and Sanger sequencing methods in MHC genotyping studies. Allele calling in NGS methods remains challenging due to extreme polymorphism and locus multiplication in the MHC coupled with allele amplification bias and the generation of artificial sequences. In this study, we compared the performance of molecular cloning with Illumina and Ion Torrent NGS sequencing in MHC-DRB genotyping of single-locus species (roe deer) and species with multiple DRB loci (red deer) in an attempt to adopt a reliable and straightforward method that does not require complex bioinformatic analyses. Our results show that all methods work similarly well in roe deer, but we demonstrate non-consistency in results across methods in red deer. With Illumina sequencing, we detected a maximum number of alleles in 10 red deer individuals (42), while other methods were somewhat less accurate as they scored 69–81% of alleles detected with Illumina sequencing.

Molecular Characterization of MHC Class I Alpha 1 and 2 Domains in Asian Seabass (Lates calcarifer)

The Asian seabass is of importance both as a farmed and wild animal. With the emergence of infectious diseases, there is a need to understand and characterize the immune system. In humans, the highly polymorphic MHC class I (MHC-I) molecules play an important role in antigen presentation for the adaptive immune system. In the present study, we characterized a single MHC-I gene in Asian seabass (Lates calcarifer) by amplifying and sequencing the MHC-I alpha 1 and alpha 2 domains, followed by multi-sequence alignment analyses. The results indicated that the Asian seabass MHC-I α1 and α2 domain sequences showed an overall similarity within Asian seabass and retained the majority of the conserved binding residues of human leukocyte antigen-A2 (HLA-A2). Phylogenetic tree analysis revealed that the sequences belonged to the U lineage. Mapping the conserved binding residue positions on human HLA-A2 and grass carp crystal structure showed a high degree of similarity. In conclusion, the availability of MHC-I α1 and α2 sequences enhances the quality of MHC class I genetic information in Asian seabass, providing new tools to analyze fish immune responses to pathogen infections, and will be applicable in the study of the phylogeny and the evolution of antigen-specific receptors.

Individual HLA heterogeneity and its implications for cellular immune evasion in cancer and beyond

Structural and functional variability of human leukocyte antigen (HLA) is the foundation for competent adaptive immune responses against pathogen and tumor antigens as it assures the breadth of the presented immune-peptidome, theoretically sustaining an efficient and diverse T cell response. This variability is presumably the result of the continuous selection by pathogens, which over the course of evolution shaped the adaptive immune system favoring the assortment of a hyper-polymorphic HLA system able to elaborate efficient immune responses. Any genetic alteration affecting this diversity may lead to pathological processes, perturbing antigen presentation capabilities, T-cell reactivity and, to some extent, natural killer cell functionality. A highly variable germline HLA genotype can convey immunogenetic protection against infections, be associated with tumor surveillance or influence response to anti-neoplastic treatments. In contrast, somatic aberrations of HLA loci, rearranging the original germline configuration, theoretically decreasing its variability, can facilitate mechanisms of immune escape that promote tumor growth and immune resistance. The purpose of the present review is to provide a unified and up-to-date overview of the pathophysiological consequences related to the perturbations of the genomic heterogeneity of HLA complexes and their impact on human diseases, with a special focus on cancer.

First Steps towards the Development of Epigenetic Biomarkers in Female Cheetahs (Acinonyx jubatus)

Free-ranging cheetahs (Acinonyx jubatus) are generally healthy, whereas cheetahs under human care, such as those in zoological gardens, suffer from ill-defined infectious and degenerative pathologies. These differences are only partially explained by husbandry management programs because both groups share low genetic diversity. However, mounting evidence suggests that physiological differences between populations in different environments can be tracked down to differences in epigenetic signatures. Here, we identified differentially methylated regions (DMRs) between free-ranging cheetahs and conspecifics in zoological gardens and prospect putative links to pathways relevant to immunity, energy balance and homeostasis. Comparing epigenomic DNA methylation profiles obtained from peripheral blood mononuclear cells (PBMCs) from eight free-ranging female cheetahs from Namibia and seven female cheetahs living in zoological gardens within Europe, we identified DMRs of which 22 were hypermethylated and 23 hypomethylated. Hypermethylated regions in cheetahs under human care were located in the promoter region of a gene involved in host-pathogen interactions (KLC1) and in an intron of a transcription factor relevant for the development of pancreatic β-cells, liver, and kidney (GLIS3). The most canonical mechanism of DNA methylation in promoter regions is assumed to repress gene transcription. Taken together, this could indicate that hypermethylation at the promoter region of KLC1 is involved in the reduced immunity in cheetahs under human care. This approach can be generalized to characterize DNA methylation profiles in larger cheetah populations under human care with a more granular longitudinal data collection, which, in the future, could be used to monitor the early onset of pathologies, and ultimately translate into the development of biomarkers with prophylactic and/or therapeutic potential.

Transcriptomic analysis elucidates evolution of the major histocompatibility complex class I in neotropical bats

The Order Chiroptera comprises more than 1,400 species, each with its evolutionary history and under unique selective pressures, among which are the host–pathogen interactions. Bats have coped with complex interactions with a broad spectrum of microbes throughout their evolutionary history, prompting the development of unique adaptations that allow them to co-exist with microbes with pathogenic potential more efficiently than other nonadapted species. In this sense, an extraordinary immune system with unique adaptations has been hypothesized in bats. To explore this, we focused on the major histocompatibility complex (MHC), which plays a crucial role in pathogen recognition and presentation to T cells to trigger the adaptive immune response. We analyzed MHC class I transcripts in five species, each from different families of New World bats. From RNA-seq data, we assembled a partial region of the MHC-I comprising the α1 and α2 domains, which are responsible for peptide binding and recognition. We described five putative functional variants, two of which have two independent insertions at the α2 domain. Our results suggest that this insertion appeared after the divergence of the order Chiroptera and may have an adaptive function in the defense against intracellular pathogens, providing evidence of positive selection and trans-specific polymorphism on the peptide-binding sites.

Oxidative stress markers, immune-regulating cytokines, and the pathological evaluation of sheep co-infected with Oestrus ovis and Coenuruscerebralis

This study aimed to evaluate the immune response using the immune-regulating cytokines as (IL-1β; MHC-I and MHC-II) associated with co-infected sheep Oestrusovis and Coenurus cerebralisas well as oxidative stress markers (malondialdehyde "MDA" and nitric oxide "NO"). So; sheep samples from different regions in Egypt showed different neurological signs, were examined for detection of the cause of the nervous manifestations. Moreover, the O. ovis and C. cerebralis cysts were collected and identified using scanning electron microscopy. The brain tissues were evaluated for different immunological genes such as MHC-I, MHC-II, and Interleukin-1β activity using quantitative real-time PCR (qRT-PCR) techniques, where the infected sheep showed higher MHC-I gene expression (10-fold), higher MHC-II gene expression (peaked at 25-fold), and higher IL-1β gene expression (14-fold) than the control group. The MDA level was significantly increased. Also, stress marker (nitric oxide) levels were significantly higher in infectedsheep than in negative control one. During gross pathology, migrating larvae of O. ovis andC. cerebralis were noticed In such areas, hemorrhages and patches of clotted blood were noticed. cysts with prominent protoscolices were also observed and were attached to the caudal region near the cerebellum.

Human Leukocyte Antigen (HLA) System: Genetics and Association with Bacterial and Viral Infections

The human leukocyte antigen (HLA) system is one of the most crucial host factors influencing disease progression in bacterial and viral infections. This review provides the basic concepts of the structure and function of HLA molecules in humans. Here, we highlight the main findings on the associations between HLA class I and class II alleles and susceptibility to important infectious diseases such as tuberculosis, leprosy, melioidosis, Staphylococcus aureus infection, human immunodeficiency virus infection, coronavirus disease 2019, hepatitis B, and hepatitis C in populations worldwide. Finally, we discuss challenges in HLA typing to predict disease outcomes in clinical implementation. Evaluation of the impact of HLA variants on the outcome of bacterial and viral infections would improve the understanding of pathogenesis and identify those at risk from infectious diseases in distinct populations and may improve the individual treatment.

Characterization, Selection, and Trans-Species Polymorphism in the MHC Class II of Heermann’s Gull (Charadriiformes)

The major histocompatibility complex (MHC) enables vertebrates to cope with pathogens and maintain healthy populations, thus making it a unique set of loci for addressing ecology and evolutionary biology questions. The aim of our study was to examine the variability of Heermann’s Gull MHC class II (MHCIIB) and compare these loci with other Charadriiformes. Fifty-nine MHCIIB haplotypes were recovered from sixty-eight Heermann’s Gulls by cloning, of them, twelve were identified as putative true alleles, forty-five as unique alleles, and two as pseudogenes. Intra and interspecific relationships indicated at least two loci in Heermann’s Gull MHCIIB and trans-species polymorphism among Charadriiformes (coinciding with the documented evidence of two ancient avian MHCIIB lineages, except in the Charadriidae family). Additionally, sites under diversifying selection revealed a better match with peptide-binding sites inferred in birds than those described in humans. Despite the negative anthropogenic activity reported on Isla Rasa, Heermann’s Gull showed MHCIIB variability consistent with population expansion, possibly due to a sudden growth following conservation efforts. Duplication must play an essential role in shaping Charadriiformes MHCIIB variability, buffering selective pressures through balancing selection. These findings suggest that MHC copy number and protected islands can contribute to seabird conservation.

Transcriptome annotation reveals minimal immunogenetic diversity among Wyoming toads, Anaxyrus baxteri

Briefly considered extinct in the wild, the future of the Wyoming toad (Anaxyrus baxteri) continues to rely on captive breeding to supplement the wild population. Given its small natural geographic range and history of rapid population decline at least partly due to fungal disease, investigation of the diversity of key receptor families involved in the host immune response represents an important conservation need. Population decline may have reduced immunogenetic diversity sufficiently to increase the vulnerability of the species to infectious diseases. Here we use comparative transcriptomics to examine the diversity of toll-like receptors and major histocompatibility complex (MHC) sequences across three individual Wyoming toads. We find reduced diversity at MHC genes compared to bufonid species with a similar history of bottleneck events. Our data provide a foundation for future studies that seek to evaluate the genetic diversity of Wyoming toads, identify biomarkers for infectious disease outcomes, and guide breeding strategies to increase genomic variability and wild release successes.

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.

Development of an Immortalized Porcine Fibroblast Cell Panel With Different Swine Leukocyte Antigen Genotypes

Immortalized cell lines are valuable resources to expand the molecular characterization of major histocompatibility complex genes and their presented antigens. We generated a panel of immortalized cell lines by transfecting human telomerase reverse transcriptase (hTERT) into primary fibroblast cells prepared from ear, fetal, and lung tissues of 10 pigs from five breeds and successfully cultured them for 30-45 passages. The cell growth characteristic of the immortalized fibroblasts was similar to that of primary fibroblast, which was unable to form colonies on soft agar. The genotypes of major swine leukocyte antigen (SLA) genes, including three classical class I (SLA-1, -2, and -3) and three class II genes (DQB1, DRB1, and DQA), were determined using high-resolution typing. A total of 58 alleles, including a novel allele for SLA-2, were identified. Each cell line was unique. A cell line derived from a National Institutes of Health miniature pig was homozygous across the six major SLA genes. The expression levels of SLA classical class I genes varied among the cell lines and were slightly upregulated in the immortalized compared to the primary cells based on semiquantitative reverse transcription polymerase chain reaction. The immortalized porcine fibroblast cell lines with diverse SLA haplotypes that were developed in this study have potential to be applied in studies regarding the molecular characteristics and genetic structure of SLA genes and epitope-major histocompatibility complex interactions in pigs.

Transmembrane Helices Are an Over-Presented and Evolutionarily Conserved Source of Major Histocompatibility Complex Class I and II Epitopes

Cytolytic T cell responses are predicted to be biased towards membrane proteins. The peptide-binding grooves of most alleles of histocompatibility complex class I (MHC-I) are relatively hydrophobic, therefore peptide fragments derived from human transmembrane helices (TMHs) are predicted to be presented more often as would be expected based on their abundance in the proteome. However, the physiological reason of why membrane proteins might be over-presented is unclear. In this study, we show that the predicted over-presentation of TMH-derived peptides is general, as it is predicted for bacteria and viruses and for both MHC-I and MHC-II, and confirmed by re-analysis of epitope databases. Moreover, we show that TMHs are evolutionarily more conserved, because single nucleotide polymorphisms (SNPs) are present relatively less frequently in TMH-coding chromosomal regions compared to regions coding for extracellular and cytoplasmic protein regions. Thus, our findings suggest that both cytolytic and helper T cells are more tuned to respond to membrane proteins, because these are evolutionary more conserved. We speculate that TMHs are less prone to mutations that enable pathogens to evade T cell responses.

Whole-genome resequencing of temporally stratified samples reveals substantial loss of haplotype diversity in the highly inbred Scandinavian wolf population

Genetic drift can dramatically change allele frequencies in small populations and lead to reduced levels of genetic diversity, including loss of segregating variants. However, there is a shortage of quantitative studies of how genetic diversity changes over time in natural populations, especially on genome-wide scales. Here, we analyzed whole-genome sequences from 76 wolves of a highly inbred Scandinavian population, founded by only one female and two males, sampled over a period of 30 yr. We obtained chromosome-level haplotypes of all three founders and found that 10%–24% of their diploid genomes had become lost after about 20 yr of inbreeding (which approximately corresponds to five generations). Lost haplotypes spanned large genomic regions, as expected from the amount of recombination during this limited time period. Altogether, 160,000 SNP alleles became lost from the population, which may include adaptive variants as well as wild-type alleles masking recessively deleterious alleles. Although not sampled, we could indirectly infer that the two male founders had megabase-sized runs of homozygosity and that all three founders showed significant haplotype sharing, meaning that there were on average only 4.2 unique haplotypes in the six copies of each autosome that the founders brought into the population. This violates the assumption of unrelated founder haplotypes often made in conservation and management of endangered species. Our study provides a novel view of how whole-genome resequencing of temporally stratified samples can be used to visualize and directly quantify the consequences of genetic drift in a small inbred population.

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.

Molecular and phylogenetic analysis of MHC class I exons 7-8 in a variety of cattle and buffalo breeds

The present study was conducted on the MHC class I (BoLA-A/BuLA-A) gene in Sahiwal, Jersey, Hariana, and Tharparkar breeds of cattle and Murrah, Mehsana, and Bhadawari breeds of buffalo to study the polymorphism. Exons 7-8 of the MHC class I gene was first characterized for polymorphism study in buffalo and the results reveal that this gene has a higher level of nucleotide changes than the cattle. Genes were investigated for polymorphisms in 285 animals of cattle and buffalo breeds. Molecular characterization of the MHC class I (BoLa-A/Bula-A) gene reveals a higher degree of polymorphism at the nucleotide level in cattle and buffalo. Results revealed this region has a higher level of polymorphisms in buffalo as campared to the cattle. Alul restriction patterns were monomorphic except for three different patterns but it was able to illustrate the differences in buffalo and cattle. SSCP analysis of exons 7-8 showed remarkable differences in cattle and buffalo. Sequence analysis revealed more closeness of Murrah breed with crossbred and indigenous cattle than Holstein Friesian. Exon 8 had more deletion and stop codon as compared to exon 7. The investigation confirmed that MHC class I BoLa-A/Bula-A exons 7-8 is highly polymorphic in buffalo as compared to cattle.

Genetic effects of long-term captive breeding on the endangered pygmy hog

Long-term captive populations often accumulate genetic changes that are detrimental to their survival in the wild. Periodic genetic evaluation of captive populations is thus necessary to identify deleterious changes and minimize their impact through planned breeding. Pygmy hog (Porcula salvania) is an endangered species with a small population inhabiting the tall sub-Himalayan grasslands of Assam, India. A conservation breeding program of pygmy hog from six founders has produced a multi-generational captive population destined for reintroduction into the wild. However, the impact of conservation breeding on its genetic diversity remained undocumented. Here, we evaluate temporal genetic changes in 39 pygmy hogs from eight consecutive generations of a captive population using genome-wide SNPs, mitochondrial genomes, and MHC sequences, and explore the relationship between genetic diversity and reproductive success. We find that pygmy hog harbors a very low genome-wide heterozygosity (H) compared to other members of the Suidae family. However, within the captive population we find excess heterozygosity and a significant increase in H from the wild-caught founders to the individuals in subsequent generations due to the selective pairing strategy. The MHC and mitochondrial nucleotide diversities were lower in captive generations compared to the founders with a high prevalence of low-frequency MHC haplotypes and more unique mitochondrial genomes. Further, even though no signs of genetic inbreeding were observed from the estimates of individual inbreeding coefficient F and between individuals (F_IS) in each generation, the kinship coefficient showed a slightly increasing trend in the recent generations, due to a relatively smaller non-random sample size compared to the entire captive population. Surprisingly, male pygmy hogs that had higher heterozygosity also showed lower breeding success. We briefly discuss the implications of our findings in the context of breeding management and recommend steps to minimize the genetic effects of long-term captive breeding.

Improved high-throughput MHC typing for non-model species using long-read sequencing

The major histocompatibility complex (MHC) plays a critical role in the vertebrate immune system. Accurate MHC typing is critical to understanding not only host fitness and disease susceptibility, but also the mechanisms underlying host‐pathogen co‐evolution. However, due to the high degree of gene duplication and diversification of MHC genes, it is often technically challenging to accurately characterise MHC genetic diversity in non‐model species. Here we conducted a systematic review to identify common issues associated with current widely used MHC typing approaches. Then to overcome these challenges, we developed a long‐read based MHC typing method along with a new analysis pipeline. Our approach enables the sequencing of fully phased MHC alleles spanning all key functional domains and the separation of highly similar alleles as well as the removal of technical artefacts such as PCR heteroduplexes and chimeras. Using this approach, we performed population‐scale MHC typing in the Tasmanian devil (Sarcophilus harrisii), revealing previously undiscovered MHC functional diversity in this endangered species. Our new method provides a better solution for addressing research questions that require high MHC typing accuracy. Since the method is not limited by species or the number of genes analysed, it will be applicable for studying not only the MHC but also other complex gene families.