Natural selection at major histocompatibility complex loci of vertebrates

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.

Developing a crop- wild-reservoir pathogen system to understand pathogen evolution and emergence

Crop pathogens reduce yield and contribute to global malnourishment. Surveillance not only detects presence/absence but also reveals genetic diversity, which can inform our understanding of rapid adaptation and control measures. An often neglected aspect is that pathogens may also use crop wild relatives as alternative hosts. This study develops the beet (Beta vulgaris) rust (Uromyces beticola) system to explore how crop pathogens evolve to evade resistance using a wild reservoir. We test predictions that crop selection will drive virulence gene differentiation and affect rates of sex between crop- and wild-host rust populations. We sequenced, assembled, and annotated the 588 Mb beet rust genome, developed a novel leaf peel pathogen DNA extraction protocol, and analysed genetic diversity in 42 wild and crop isolates. We found evidence for two populations: one containing exclusively wild-host isolates; the other containing all crop-host isolates, plus five wild isolates. Effectors showed greater diversity in the exclusively wild population and greater differentiation between populations. Preliminary evidence suggests the rates of sexual reproduction may differ between populations. This study highlights how differences in pathogen populations might be used to identify genes important for survival on crops and how reproduction might impact adaptation. These findings are relevant to all crop-reservoir systems and will remain unnoticed without comparison to wild reservoirs.

Characterization, expression, and polymorphism of MHC II α and MHC II β in Sichuan taimen (Hucho bleekeri)

The major histocompatibility complex (MHC) is involved in antigen presentation and plays an essential role in regulating immune function. In the present study, we identified two MHC class II genes and investigated their potential roles in Hucho bleekeri. The MHC II α and MHC II β of H. bleekeri had typical leading peptides, extracellular domains, connecting peptides, transmembrane region, and cytoplasmic region. Amino acid sequence comparison revealed that MHC II of H. bleekeri has high homology with other vertebrates, among which homology with salmonid fish was the highest. Phylogenetic analysis showed that H. bleekeri MHC II clustered with salmonid fish; moreover they clustered with orthologous genes of other fish, whereas mammalian MHC II clustered into a separate branch. Tissue distribution analysis revealed MHC II was widely expressed in all tested tissues, with both MHC II α and MHC II β highly expressed in the spleen, gill, kidney, and hindgut. After lipopolysaccharide (LPS) and polyinosinic:polycytidylic acid (poly(I:C)) stimulation, the expression of MHC II in the head kidney and spleen of H. bleekeri was significantly upregulated. Compared with MHC II α, MHC II β acted faster in response to the stimulation. Polymorphism analysis of MHC II revealed that all the different alleles belonged to the same major type, and very limited polymorphisms were found in H. bleekeri MHC II α and II β. Selection pressure analysis showed signs of weak and non-significant positive selection in the MHC II α and MHC II β extracellular region. Our study reveals the potential role of MHC II in the immune response of H. bleekeri and provides a reference for studying the evolutionary model of teleost MHC II.

Expansion of MHC-IIB Has Constrained the Evolution of MHC-IIA in Passerines

The major histocompatibility complex (MHC) is central in adaptive immunity, with the highly polymorphic MHC genes encoding antigen-presenting molecules. Two MHC class II (MHC-II) loci, DA1 and DA2, predate the radiation of extant birds and persist throughout much of the avian phylogeny. Within each locus, the MHC-II molecules are encoded by A-genes (DAA) and B-genes (DAB), which are arranged in A-B dyads. However, in passerines (order Passeriformes), the DA2 locus has been lost, and the ancestral A-B dyad at the DA1 locus has been replaced by a putatively single A-gene (DAA1) and an array of highly polymorphic B-genes (DAB1). In this study, we genotyped the DAA1 gene of 15 passerine species and confirmed that passerines possess just one copy of DAA1. We then compared selection patterns in DAA1 between passerines and nonpasserines and found that exon 2, which encodes the antigen-presenting domain, has been subject to weaker positive selection and stronger negative selection in passerines compared with nonpasserines. Additional comparisons showed that the patterns of selection in the passerine DAA1 gene are unlikely to be related to the loss of the DA2 locus. Instead, our findings suggest that the expansion of DAB1 (MHC-IIB) has imposed an evolutionary constraint on the passerine DAA1 (MHC-IIA) gene. We speculate that this constraint may be the result of each DAA1 chain forming heterodimers with many different DAB1 chains.

Moderate Genetic Diversity of MHC Genes in an Isolated Small Population of Black-and-White Snub-Nosed Monkeys (Rhinopithecus bieti)

Genetic diversity is an essential indicator that echoes the natural selection and environmental adaptation of a species. Isolated small populations are vulnerable to genetic drift, inbreeding, and limited gene flow; thus, assessing their genetic diversity is critical in conservation. In this study, we studied the genetic diversity of black-and-white snub-nosed monkeys (Rhinopithecus bieti) using neutral microsatellites and five adaptive major histocompatibility complex (MHC) genes. Two DQA1 alleles, two DQB1 alleles, two DRB1 alleles, two DRB5 alleles, and three DPB1 alleles were isolated from a population. The results indicate that neutral microsatellites demonstrate a high degree of heterozygosity and polymorphism, while adaptive MHC genes display a high degree of heterozygosity and moderate polymorphism. The results also show that balancing selection has prominently influenced the MHC diversity of the species during evolution: (1) significant positive selection is identified at several amino acid sites (primarily at and near antigen-binding sites) of the DRB1, DRB5, and DQB1 genes; (2) phylogenetic analyses display the patterns of trans-species evolution for all MHC loci. This study provides valuable genetic diversity insights into black-and-white snub-nosed monkeys, which dwell at the highest altitude and have experienced the harshest environmental selection of all primates globally since the Pleistocene. Such results provide valuable scientific evidence and a reference for making or amending conservation strategies for this endangered primate species.

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle)

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

Evolution of T cells in the cancer-resistant naked mole-rat

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αβT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

BACKGROUND

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

RESULTS

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages.

CONCLUSIONS

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Assessment of human leukocyte antigen-based neoantigen presentation to determine pan-cancer response to immunotherapy

Despite the central role of human leukocyte antigen class I (HLA-I) in tumor neoantigen presentation, quantitative determination of presentation capacity remains elusive. Based on a pooled pan-cancer genomic dataset of 885 patients treated with immune checkpoint inhibitors (ICIs), we developed a score integrating the binding affinity of neoantigens to HLA-I, as well as HLA-I allele divergence, termed the HLA tumor-Antigen Presentation Score (HAPS). Patients with a high HAPS were more likely to experience survival benefit following ICI treatment. Analysis of the tumor microenvironment indicated that the antigen presentation pathway was enriched in patients with a high HAPS. Finally, we built a neural network incorporating factors associated with neoantigen production, presentation, and recognition, which exhibited potential for differentiating cancer patients likely to benefit from ICIs. Our findings highlight the clinical utility of evaluating HLA-I tumor antigen presentation capacity and describe how ICI response may depend on HLA-mediated immunity.

Limited associations between MHC diversity and reproductive success in a bird species with biparental care

The selective pressure from pathogens on individuals can have direct consequences on reproduction. Genes from the major histocompatibility complex (MHC) are central to the vertebrate adaptive immune system and pathogen resistance. In species with biparental care, each sex has distinct reproductive roles and levels of investment, and due to a trade-off with immunity, one can expect different selective regimes acting upon the MHC of each parent. Here, we addressed whether couples combine each other's variation at MHC loci to increase their breeding success. Specifically, we used a 23-year dataset from a barn owl population (Tyto alba) to understand how MHC class Iα and IIβ functional divergence and supertypes of each parent were associated with clutch size and fledging success. We did not detect associations between MHC diversity and supertypes with the clutch size or with the fledging success. In addition, to understand the relative contribution from the MHC of the genetic parents and the social parents, we analyzed the fledging success using only a cross-fostered dataset. We found several associations of weak-to-moderate effect sizes between the father's MHC and fledging success: (i) lower MHC-Iα divergence in the genetic father increases fledging success, which might improve paternal care during incubation, and (ii) one and two MHC-IIβ DAB2 supertypes in the social father decrease and increase, respectively, fledging success, which may affect the paternal care after hatching. Furthermore, fledging success increased when both parents did not carry MHC-IIβ DAB1 supertype 2, which could suggest conditional effects of this supertype. Although our study relied on a substantial dataset, we showed that the associations between MHC diversity and reproductive success remain scarce and of complex interpretation in the barn owl. Moreover, our results highlighted the need to incorporate more than one proxy of reproductive success and several MHC classes to capture more complex associations.

Molecular Characteristics, Functional Definitions, and Regulatory Mechanisms for Cross-Presentation Mediated by the Major Histocompatibility Complex: A Comprehensive Review

The major histocompatibility complexes of vertebrates play a key role in the immune response. Antigen-presenting cells are loaded on MHC I molecules, which mainly present endogenous antigens; when MHC I presents exogenous antigens, this is called cross-presentation. The discovery of cross-presentation provides an important theoretical basis for the study of exogenous antigens. Cross-presentation is a complex process in which MHC I molecules present antigens to the cell surface to activate CD8+ T lymphocytes. The process of cross-representation includes many components, and this article briefly outlines the origins and development of MHC molecules, gene structures, functions, and their classical presentation pathways. The cross-presentation pathways of MHC I molecules, the cell lines that support cross-presentation, and the mechanisms of MHC I molecular transporting are all reviewed. After more than 40 years of research, the specific mechanism of cross-presentation is still unclear. In this paper, we summarize cross-presentation and anticipate the research and development prospects for cross-presentation.

Bw4 ligand and direct T-cell receptor binding induced selection on HLA A and B alleles

Introduction

The HLA region is the hallmark of balancing selection, argued to be driven by the pressure to present a wide variety of viral epitopes. As such selection on the peptide-binding positions has been proposed to drive HLA population genetics. MHC molecules also directly binds to the T-Cell Receptor and killer cell immunoglobulin-like receptors (KIR).

Methods

We here combine the HLA allele frequencies in over six-million Hematopoietic Stem Cells (HSC) donors with a novel machine-learning-based method to predict allele frequency.

Results

We show for the first time that allele frequency can be predicted from their sequences. This prediction yields a natural measure for selection. The strongest selection is affecting KIR binding regions, followed by the peptide-binding cleft. The selection from the direct interaction with the KIR and TCR is centered on positively charged residues (mainly Arginine), and some positions in the peptide-binding cleft are not associated with the allele frequency, especially Tyrosine residues.

Discussion

These results suggest that the balancing selection for peptide presentation is combined with a positive selection for KIR and TCR binding.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

Background

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

Results

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages.

Conclusions

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Identification of long non-coding RNA MSTRG.5748.1 and MSTRG.7894.1 from Megalobrama amblycephala and their potential roles in innate immunity

Megalobrama amblycephala is one of the most economically important freshwater fish in China, and the bacterial septicemia caused by Aeromonas hydrophila is a serious threat to the breeding industry of M. amblycephala. Unfortunately, the characterization of long noncoding RNA (lncRNA) in response to A. hydrophila infection has not been performed in M. amblycephala. To better understand the biological significance of lncRNA in the immune system, we identified two lncRNA, named MSTRG.5748.1 and MSTRG.7894.1, as playing critical roles in the antibacterial response of M. amblycephala. After separating the nucleus and cytoplasm of the hepatocytes from M. amblycephala, cellular localization of MSTRG.5748.1 and MSTRG.7894.1 was performed to predict their functions. The results showed that MSTRG.5748.1 was mainly expressed in the nucleus, suggesting that its functions are mostly to regulate the expression of downstream genes through epistasis and transcription. MSTRG.7894.1 existed in both the nucleus and cytoplasm, which indicated that it has many regulatory modes. qPCR analysis showed that MSTRG.5748.1 and MSTRG.7894.1 were expressed in the immune-related organs of M. amblycephala, and significantly changed in the liver after A. hydrophila infection. RNA-seq analysis revealed that differentially expressed genes (DEGs) were mainly enriched in antigen processing and presentation via MHC class I, RIG-I-like receptor (RLR) signaling pathway, and IFN-related pathway, and a large number of pathway-related genes were significantly regulated after lncRNA overexpression in muscle cell of M. amblycephala. Overexpression of MSTRG.5748.1 and MSTRG.7894.1 significantly inhibited the expression of STING and IFN, significantly upregulated muscle cell viability, and promoted cell proliferation by targeting STING and IFN.

Paternity data reveal high MHC diversity among sires in a polygynandrous, egalitarian primate

Evidence from human and nonhuman primates suggests that females avoid breeding with close kin and may choose mates based on MHC diversity, which can improve offspring survival. In despotic societies, female mate choice may be hindered by male sexual coercion, but in egalitarian societies, females may be less constrained. Among northern muriquis-an egalitarian, polygynandrous primate with male philopatry-analyses of new data on paternity and variation at microsatellite and MHC loci, combined with behavioural and life-history data, revealed that sires showed higher MHC diversity than expected by chance and were never close kin of dams, consistent with predictions of female mate choice and close inbreeding avoidance. However, females did not differentially reproduce with males who were more distantly related to them or more dissimilar at the MHC than expected by chance, nor with those who had more MHC alleles distinct from their own. The lack of male dominance may permit females to identify and reproduce preferentially with non-offspring males and with males who are more diverse at the MHC. Nonetheless, the absence of disassortative mating at the MHC and neutral loci suggests that female mate choice may be limited by other factors impacting male fertilization success.

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

Neutral Forces and Balancing Selection Interplay to Shape the Major Histocompatibility Complex Spatial Patterns in the Striped Hamster in Inner Mongolia: Suggestive of Broad-Scale Local Adaptation

BACKGROUND

The major histocompatibility complex (MHC) plays a key role in the adaptive immune response to pathogens due to its extraordinary polymorphism. However, the spatial patterns of MHC variation in the striped hamster remain unclear, particularly regarding the relative contribution of the balancing selection in shaping MHC spatial variation and diversity compared to neutral forces.

METHODS

In this study, we investigated the immunogenic variation of the striped hamster in four wild populations in Inner Mongolia which experience a heterogeneous parasitic burden. Our goal was to identify local adaptation by comparing the genetic structure at the MHC with that at seven microsatellite loci, taking into account neutral processes.

RESULTS

We observed significant variation in parasite pressure among sites, with parasite burden showing a correlation with temperature and precipitation. Molecular analysis revealed a similar co-structure between MHC and microsatellite loci. We observed lower genetic differentiation at MHC loci compared to microsatellite loci, and no correlation was found between the two.

CONCLUSIONS

Overall, these results suggest a complex interplay between neutral evolutionary forces and balancing selection in shaping the spatial patterns of MHC variation. Local adaptation was not detected on a small scale but may be applicable on a larger scale.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Cell-free protein synthesis of CD1E and B2M protein and in vitro interaction

CD1E, one of the most important glycolipid antigens on T cell membranes, is required for glycolipid antigen presentation on the cell surface. Cell-based recombinant expression systems have many limitations for synthesizing transmembrane proteins such as CD1E, including low protein yields and miss folding. To overcome these challenges, here we successfully synthesized high-quality soluble CD1E using an E.coli cell-free protein synthesis system (CFPS) with the aid of detergent. Following purification by Ni²⁺ affinity chromatography, we were able to obtain CD1E with ≥90% purity. Furthermore, we used the string website to predict the protein interaction network of CD1E and identified a potential binding partner━B2M. Similarly, we synthesized soluble B2M in the E.coli CFPS. Finally, we verified the interaction between CD1E and B2M by using Surface Plasmon Resonance (SPR). Taken together, the methods described here provide an alternative way to obtain active transmembrane protein and may facilitate future structural and functional studies on CD1E.

Divergent allele advantage in the MHC and amphibian emerging infectious disease

Genetic variation in the major histocompatibility complex (MHC) may be associated with resistance to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). The pathogen originated in Asia, then spread worldwide, causing amphibian population declines and species extinctions. We compared the expressed MHC IIβ1 alleles of a Bd-resistant toad species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian frog species, Litoria caerulea. We found at least six expressed MHC IIβ1 loci in each of the two species. Amino acid diversity encoded by these MHC alleles was similar between species, but the genetic divergence of those alleles known for broader pathogen-derived peptide binding was greater in the Bd-resistant species. In addition, we found a potentially rare allele in one resistant individual from the Bd-susceptible species. Deep next-generation sequencing recovered approximately triple the genetic resolution accessible from traditional cloning-based genotyping. Targeting more than one MHC IIβ1 expressed locus enables us to better understand how host MHC may adapt to emerging infectious diseases.

Malaria-driven adaptation of MHC class I in wild bonobo populations

The malaria parasite Plasmodium falciparum causes substantial human mortality, primarily in equatorial Africa. Enriched in affected African populations, the B*53 variant of HLA-B, a cell surface protein that presents peptide antigens to cytotoxic lymphocytes, confers protection against severe malaria. Gorilla, chimpanzee, and bonobo are humans' closest living relatives. These African apes have HLA-B orthologs and are infected by parasites in the same subgenus (Laverania) as P. falciparum, but the consequences of these infections are unclear. Laverania parasites infect bonobos (Pan paniscus) at only one (TL2) of many sites sampled across their range. TL2 spans the Lomami River and has genetically divergent subpopulations of bonobos on each side. Papa-B, the bonobo ortholog of HLA-B, includes variants having a B*53-like (B07) peptide-binding supertype profile. Here we show that B07 Papa-B occur at high frequency in TL2 bonobos and that malaria appears to have independently selected for different B07 alleles in the two subpopulations.

The HLA-DRB1*07 Allele Is Associated with Interstitial Lung Abnormalities (ILA) and Subpleural Location in a Mexican Mestizo Population

Interstitial lung abnormalities (ILA) are defined as the presence of different patterns of increased lung density, including ground glass attenuation and reticular opacities on chest high-resolution computed tomography (HRCT). In this study, we included 90 subjects with ILA and 189 healthy controls (HC) from our Aging Lung Program. We found that subjects with ILA are older, have a significant smoking history, and have worse pulmonary function than HC (p < 0.05). When we evaluated the allele frequencies of the human leukocyte antigen (HLA) system, we found that HLA-DRB1*07 was associated with a higher risk for ILA (p < 0.05, OR = 1.95, 95% CI = 1.06-3.57). When we compared subjects with subpleural ILA vs. HC, the association with HLA-DRB1*07 became stronger than the whole ILA group (p < 0.05, OR = 2.29, 95% CI = 1.24-4.25). Furthermore, subjects with subpleural ILA and central ILA display differences in allele frequencies with HLA-DRB1*14 (3.33% vs. 13.33%, p < 0.05) and *15 (3.33% vs. 20%, p < 0.05). Our findings indicate that the HLA-DRB1*07 allele contributes to the risk of ILA, especially those of subpleural locations.

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.

Next-generation sequencing: unraveling genetic mechanisms that shape cancer immunotherapy efficacy

Immunity is governed by fundamental genetic processes. These processes shape the nature of immune cells and set the rules that dictate the myriad complex cellular interactions that power immune systems. Everything from the generation of T cell receptors and antibodies, control of epitope presentation, and recognition of pathogens by the immunoediting of cancer cells is, in large part, made possible by core genetic mechanisms and the cellular machinery that they encode. In the last decade, next-generation sequencing has been used to dissect the complexities of cancer immunity with potent effect. Sequencing of exomes and genomes has begun to reveal how the immune system recognizes “foreign” entities and distinguishes self from non-self, especially in the setting of cancer. High-throughput analyses of transcriptomes have revealed deep insights into how the tumor microenvironment affects immunotherapy efficacy. In this Review, we discuss how high-throughput sequencing has added to our understanding of how immune systems interact with cancer cells and how cancer immunotherapies work.

Molecular Mechanisms Underlying Vibrio Tolerance in Ruditapes philippinarum Revealed by Comparative Transcriptome Profiling

The clam Ruditapes philippinarum is an important species in the marine aquaculture industry in China. However, in recent years, the aquaculture of R. philippinarum has been negatively impacted by various bacterial pathogens. In this study, the transcriptome libraries of R. philippinarum showing different levels of resistance to challenge with Vibrio anguillarum were constructed and RNA-seq was performed using the Illumina sequencing platform. Host immune factors were identified that responded to V. anguillarum infection, including C-type lectin domain, glutathione S-transferase 9, lysozyme, methyltransferase FkbM domain, heat shock 70 kDa protein, Ras-like GTP-binding protein RHO, C1q, F-box and BTB/POZ domain protein zf-C2H2. Ten genes were selected and verified by RT-qPCR, and nine of the gene expression results were consistent with those of RNA-seq. The lectin gene in the phagosome pathway was expressed at a significantly higher level after V. anguillarum infection, which might indicate the role of lectin in the immune response to V. anguillarum. Comparing the results from R. philippinarum resistant and nonresistant to V. anguillarum increases our understanding of the resistant genes and key pathways related to Vibrio challenge in this species. The results obtained here provide a reference for future immunological research focusing on the response of R. philippinarum to V. anguillarum infection.

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.

Selection balancing at innate immune genes: adaptive polymorphism maintenance in Toll-like receptors

Balancing selection is a classic mechanism for maintaining variability in immune genes involved in host–pathogen interactions. However, it remains unclear how widespread the mechanism is across immune genes other than the major histocompatibility complex (MHC). Although occasional reports suggest that balancing selection (heterozygote advantage, negative frequency-dependent selection, and fluctuating selection) may act on other immune genes, the current understanding of the phenomenon in non-MHC immune genes is far from solid. In this review, we focus on Toll-like receptors (TLRs), innate immune genes directly involved in pathogen recognition and immune response activation, as there is a growing body of research testing the assumptions of balancing selection in these genes. After reviewing infection- and fitness-based evidence, along with evidence based on population allelic frequencies and heterozygosity levels, we conclude that balancing selection maintains variation in TLRs, though it tends to occur under specific conditions in certain evolutionary lineages rather than being universal and ubiquitous. Our review also identifies key gaps in current knowledge and proposes promising areas for future research. Improving our understanding of host–pathogen interactions and balancing selection in innate immune genes are increasingly important, particularly regarding threats from emerging zoonotic diseases.

Using de novo genome assembly and high-throughput sequencing to characterize the MHC region in a non-model bird, the Eurasian coot

Genes of the Major Histocompatibility Complex (MHC) form a key component of vertebrate adaptive immunity, as they code for molecules which bind antigens of intra- and extracellular pathogens (MHC class I and II, respectively) and present them to T cell receptors. In general, MHC genes are hyper-polymorphic and high MHC diversity is often maintained within natural populations (via balancing selection) and within individuals (via gene duplications). Because of its complex architecture with tandems of duplicated genes, characterization of MHC region in non-model vertebrate species still poses a major challenge. Here, we combined de novo genome assembly and high-throughput sequencing to characterize MHC polymorphism in a rallid bird species, the Eurasian coot Fulica atra. An analysis of genome assembly indicated high duplication rate at MHC-I, which was also supported by targeted sequencing of peptide-binding exons (at least five MHC-I loci genotyped). We found high allelic richness at both MHC-I and MHC-II, although signature of diversifying selection and recombination (gene conversion) was much stronger at MHC-II. Our results indicate that Eurasian coot retains extraordinary polymorphism at both MHC classes (when compared to other non-passerine bird species), although they may be subject to different evolutionary mechanism.

Immunogenetic variation shapes the gut microbiome in a natural vertebrate population

BACKGROUND

The gut microbiome (GM) can influence many biological processes in the host, impacting its health and survival, but the GM can also be influenced by the host's traits. In vertebrates, Major Histocompatibility Complex (MHC) genes play a pivotal role in combatting pathogens and are thought to shape the host's GM. Despite this-and the documented importance of both GM and MHC variation to individual fitness-few studies have investigated the association between the GM and MHC in the wild.

RESULTS

We characterised MHC class I (MHC-I), MHC class II (MHC-II) and GM variation in individuals within a natural population of the Seychelles warbler (Acrocephalus sechellensis). We determined how the diversity and composition of the GM varied with MHC characteristics, in addition to environmental factors and other host traits. Our results show that the presence of specific MHC alleles, but not MHC diversity, influences both the diversity and composition of the GM in this population. MHC-I alleles, rather than MHC-II alleles, had the greatest impact on the GM. GM diversity was negatively associated with the presence of three MHC-I alleles (Ase-ua3, Ase-ua4, Ase-ua5), and one MHC-II allele (Ase-dab4), while changes in GM composition were associated with the presence of four different MHC-I alleles (Ase-ua1, Ase-ua7, Ase-ua10, Ase-ua11). There were no associations between GM diversity and TLR3 genotype, but GM diversity was positively correlated with genome-wide heterozygosity and varied with host age and field period.

CONCLUSIONS

These results suggest that components of the host's immune system play a role in shaping the GM of wild animals. Host genotype-specifically MHC-I and to a lesser degree MHC-II variation-can modulate the GM, although whether this occurs directly, or indirectly through effects on host health, is unclear. Importantly, if immune genes can regulate host health through modulation of the microbiome, then it is plausible that the microbiome could also influence selection on immune genes. As such, host-microbiome coevolution may play a role in maintaining functional immunogenetic variation within natural vertebrate populations. Video abstract.

Sequence Diversity, Locus Structure, and Evolutionary History of the SpTransformer Genes in the Sea Urchin Genome

The generation of large immune gene families is often driven by evolutionary pressure exerted on host genomes by their pathogens, which has been described as the immunological arms race. The SpTransformer (SpTrf) gene family from the California purple sea urchin, Strongylocentrotus purpuratus, is upregulated upon immune challenge and encodes the SpTrf proteins that interact with pathogens during an immune response. Native SpTrf proteins bind both bacteria and yeast, and augment phagocytosis of a marine Vibrio, while a recombinant SpTrf protein (rSpTrf-E1) binds a subset of pathogens and a range of pathogen associated molecular patterns. In the sequenced sea urchin genome, there are four SpTrf gene clusters for a total of 17 genes. Here, we report an in-depth analysis of these genes to understand the sequence complexities of this family, its genomic structure, and to derive a putative evolutionary history for the formation of the gene clusters. We report a detailed characterization of gene structure including the intron type and UTRs with conserved transcriptional start sites, the start codon and multiple stop codons, and locations of polyadenylation signals. Phylogenetic and percent mismatch analyses of the genes and the intergenic regions allowed us to predict the last common ancestral SpTrf gene and a theoretical evolutionary history of the gene family. The appearance of the gene clusters from the theoretical ancestral gene may have been driven by multiple duplication and deletion events of regions containing SpTrf genes. Duplications and ectopic insertion events, indels, and point mutations in the exons likely resulted in the extant genes and family structure. This theoretical evolutionary history is consistent with the involvement of these genes in the arms race in responses to pathogens and suggests that the diversification of these genes and their encoded proteins have been selected for based on the survival benefits of pathogen binding and host protection.

Urbanization processes drive divergence at the major histocompatibility complex in a common waterbird

Urban sprawl is one of the most common landscape alterations occurring worldwide, and there is a growing list of species that are recognised to have adapted to urban life. To be successful, processes of urban colonization by wildlife require a broad spectrum of phenotypic (e.g., behavioural or physiological) adjustments, but evidence for genetic adaptations is much scarcer. One hypothesis proposes that different pathogen-driven selective pressures between urban and non-urban landscapes leads to adaptations in host immune genes. Here, we examined urbanization-related differentiation at the key pathogen-recognition genes of vertebrate adaptive immunity-the major histocompatibility complex (MHC)-in a common waterbird, the Eurasian coot (Fulica atra). Samples were collected from an old urban population (established before the 1950s), a new urban population (established in the 2000s), and two rural populations from central Poland. We found strong significant divergence (as measured with Jost's D) at the MHC class II between the old urban population and the remaining (new urban and rural) populations. Also, there was a moderate, but significant divergence at the MHC between the new urban population and two rural populations, while no divergence was found between the two rural populations. The total number of MHC alleles and the number of private (population-specific) MHC alleles was lower in old urban populations, as compared to the rural ones. These patterns of differentiation at the MHC were not consistent with patterns found for neutral genetic markers (microsatellites), which showed few differences between the populations. Our results indicate that MHC allele composition depended on the level of anthropogenic disturbance and the time which passed since urban colonization, possibly due to the processes of genotype sorting and local adaptation. As such, our study contributes to the understanding of genetic mechanisms associated with urbanization processes in wildlife.

Integrated analysis reveals prognostic value of HLA-I LOH in triple-negative breast cancer

Background

Triple-negative breast cancers (TNBCs), especially those non-immune-inflamed tumors, have a poor prognosis and limited therapies. Human leukocyte antigen (HLA)-I not only contributes to antitumor immune response and the phenotype of the tumor microenvironment, but also is a negative predictor of outcomes after immunotherapy. However, the importance of HLA functional status in TNBCs remains poorly understood.

Methods

Using the largest original multiomics datasets on TNBCs, we systematically characterized the HLA-Ⅰ status of TNBCs from the perspective of HLA-Ⅰ homogeneity and loss of heterozygosity (LOH). The prognostic significance of HLA-I status was measured. To explain the potential mechanism of prognostic value in HLA-Ⅰ status, the mutational signature, copy number alteration, neoantigen and intratumoral heterogeneity were measured. Furthermore, the correlation between HLA-Ⅰ functional status and the tumor immune microenvironment was analyzed.

Results

LOH and homogeneity in HLA-I accounted for 18% and 21% of TNBCs, respectively. HLA-I LOH instead of HLA-I homogeneity was an independent prognostic biomarker in TNBCs. In particular, for patients with non-immune-inflamed tumors, HLA-I LOH indicated a worse prognosis than HLA-I non-LOH. Furthermore, integrated genomic and transcriptomic analysis showed that HLA-I LOH was accompanied by upregulated scores of mutational signature 3 and homologous recombination deficiency scores, which implied the failure of DNA double-strand break repair. Moreover, HLA-I LOH had higher mutation and neoantigen loads and more subclones than HLA-I non-LOH. These results indicated that although HLA-I LOH tumors with failure of DNA double-strand break repair were prone to produce neoantigens, their limited capacity for antigen presentation finally contributed to poor immune selection pressure.

Conclusion

Our study illustrates the genomic landscape of HLA-I functional status and stresses the prognostic significance of HLA-I LOH in TNBCs. For “cold” tumors in TNBCs, HLA-I LOH indicated a worse prognosis than HLA-I non-LOH.

BoLA-DRB3 gene haplotypes show divergence in native Sudanese cattle from taurine and indicine breeds

Autochthonous Sudanese cattle breeds, namely Baggara for beef and Butana and Kenana for dairy, are characterized by their adaptive characteristics and high performance in hot and dry agro-ecosystems. They are thus used largely by nomadic and semi-nomadic pastoralists. We analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus linked to the immune response, for the indigenous cattle of Sudan and in the context of the global cattle repository. Blood samples (n = 225) were taken from three indigenous breeds (Baggara; n = 113, Butana; n = 60 and Kenana; n = 52) distributed across six regions of Sudan. Nucleotide sequences were genotyped using the sequence-based typing method. We describe 53 alleles, including seven novel alleles. Principal component analysis (PCA) of the protein pockets implicated in the antigen-binding function of the MHC complex revealed that pockets 4 and 9 (respectively) differentiate Kenana-Baggara and Kenana-Butana breeds from other breeds. Venn analysis of Sudanese, Southeast Asian, European and American cattle breeds with 115 alleles showed 14 were unique to Sudanese breeds. Gene frequency distributions of Baggara cattle showed an even distribution suggesting balancing selection, while the selection index (ω) revealed the presence of diversifying selection in several amino acid sites along the BoLA-DRB3 exon 2 of these native breeds. The results of several PCA were in agreement with clustering patterns observed on the neighbor joining (NJ) trees. These results provide insight into their high survival rate for different tropical diseases and their reproductive capacity in Sudan's harsh environment.

Nucleotide sequence variation of the major histocompatibility complex class II DQA1 gene in different cattle breeds from Nigeria and New Zealand

The major histocompatibility complex (MHC) plays a role in immune response. Among other activities, the bovine MHC genes (BoLA) trigger immune responses, including the activation of antibody-producing B-cells. In this study, White Fulani (n = 24), Red Bororo (n = 5) and Holstein-White × Fulani-cross (n = 11) cattle from Nigeria, and New Zealand Holstein-Friesian × Jersey-cross (n = 40) cattle were used to investigate variability in exon 2 of BoLA-DQA1. Ten alleles were identified using a PCR-Single Strand Conformation Polymorphism (SSCP) approach and their nucleotide sequences confirmed by DNA sequencing. A total of 12.60 % of all nucleotide positions analysed were revealed to be variable and two novel BoLA-DQA1 alleles are reported here for the first time.

Genetic assessment of BoLA-DRB3 polymorphisms by comparing Bangladesh, Ethiopian, and Korean cattle

Attributable to their major function in pathogen recognition, the use of bovine leukocyte antigens (BoLA) as disease markers in immunological traits in cattle is well established. However, limited report exists on polymorphism of the BoLA gene in zebu cattle breeds by high resolution typing methods. Thus, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 100 animals (Boran, n = 13; Sheko, n = 20; Fogera, n = 16; Horro, n = 19), Hanwoo cattle (n = 18) and Bangladesh Red Chittagong zebu (n = 14). Out of the 59 detected alleles, 43 were already deposited under the Immuno Polymorphism Database for major histocompatibility complex (IPD-MHC) while 16 were unique to this study. Assessment of the level of genetic variability at the population and sequence levels with genetic distance in the breeds considered in this study showed that Zebu breeds had a gene diversity score greater than 0.752, nucleotide diversity score greater than 0.152, and mean number of pairwise differences higher than 14, being very comparable to those investigated for other cattle breeds. Regarding neutrality tests analyzed, we investigated that all the breeds except Hanwoo had an excess number of alleles and could be expected from a recent population expansion or genetic hitchhiking. Howbeit, the observed heterozygosity was not significantly (p < 0.05) higher than the expected heterozygosity. The Hardy Weinberg equilibrium (HWE) analysis revealed non-significant excess of heterozygote animals, indicative of plausible over-dominant selection. The pairwise FST values suggested a low genetic variation among all the breeds (FST = 0.056; p < 0.05), besides the rooting from the evolutionary or domestication history of the cattle. No detached clade was observed in the evolutionary divergence study of the BoLA-DRB3 gene, inferred from the phylogenetic tree based on the maximum likelihood model. The investigation herein indicated the clear differences in BoLA-DRB3 gene variability between African and Asian cattle breeds.

Genetic diversity of major histocompatibility complex class I genes in Zootoca vivipara

The Major Histocompatibility Complex (MHC), as a family of highly polymorphic genes associated with immunity in the genome of the vertebrate, has become an important indicator for assessing the evolutionary potential of wildlife. In order to better protect Zootoca vivipara in the Greater Khingan Range and Lesser Khingan Range, to understand the genetic structure of Z. vivipara, and to explore the mechanism and phylogenetic relationship of the gene polymorphisms, the MHC molecular marker method was used to analyze Z. vivipara population. Forty-seven alleles were obtained from four populations. The four populations were highly polymorphic, rich in genetic information, and had significant genetic diversity. There were certain inbreeding phenomena. There was a high degree of genetic differentiation among populations, which was caused by genetic drift and natural selection. The sequence undergoes genetic duplication and recombination. The existence of trans-species polymorphism was found in the constructed phylogenetic tree. The present study provides a theoretical basis for species protection of Z. vivipara.

Evolution of MHC class I genes in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Major histocompatibility complex (MHC) genes have been widely studied to assess the immunological fitness and evolutionary adaptation of animal populations. Among the Canidae, the raccoon dog's adventurous nature, omnivorous behavior, and high variability of intracellular pathogens make it ideal to study selection on MHC class I in a non-model canid species. Here, we examined allelic diversity and evolutionary patterns of MHC class I genes in the raccoon dog (Nyctereutes procyonoides). We identified 48 novel MHC class I alleles from 31 raccoon dogs from Japan and Russia. Some alleles were geographically restricted, whereas others were widely distributed across the species' range. The rate of non-synonymous substitutions was greater than that of synonymous substitutions for both exon 2 and exon 3 encoding α1 and α2 domains, respectively, in the α chain of the MHC class I protein. Positively selected sites at the amino acid level were evident in both the α1 and α2 domains, and a recombination breakpoint was found in exon 3. Bayesian phylogenetic trees showed no evidence of trans-species polymorphism (TSP) with alleles from carnivoran species in other families but did detect TSP between raccoon dogs and the domestic dog, Canis familiaris, indicative of long-term balancing selection in canids. Our results indicate that the extensive allelic diversity of MHC class I in Japanese and Russian raccoon dogs has been influenced and maintained by pathogen-driven positive selection, recombination, and long-term balancing selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13364-021-00561-y.

Response of the ileum transcriptome to fructo-oligosaccharides in Taiping chickens

The aim of this study was to investigate the effects of fructo-oligosaccharide (FOS) supplementation intake of Taiping chickens (Gallus gallus domesticus) and its stimulating effects on ileum. 120 healthy chickens were randomly divided into two groups; control group (CT) and fructo-oligosaccharides group (FOS). At the 60th day of age, ileum mucosa of three chickens per group were collected and performed transcriptome profiling of Taiping chicken ileum mucosa using the Hiseq™ 2500 sequencing platform. Compared with CT group, 50 genes were differentially expressed in the FOS group. Ten of the differently expressed genes were further validated by RT-qPCR. In addition, gene ontology and Kyoto encyclopedia of genes and genomes analyses revealed that these differentially expressed genes were mainly enriched to drug metabolism-cytochrome P450, metabolism of xenobiotics by cytochrome P450, retinol metabolism, fat digestion and absorption, herpes simplex infection and valine, leucine and isoleucine biosynthesis. The results of this study provided the help to our understanding application of fructo-oligosaccharides in indigenous chicken production and provide a theoretical basis for the genetic development of indigenous chickens.

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

Chromosome level assembly reveals a unique immune gene organization and signatures of evolution in the common pheasant

The common pheasant Phasianus colchicus, belonging to the order Galliformes and family Phasianidae, is the most widespread species. Despite a long history of captivity, the domestication of this bird is still at a preliminary stage. Recently, the demand for accelerating its transformation to poultry for meat and egg production has been increasing. In this study, we assembled high quality, chromosome scale genome of the common pheasant by using PacBio long reads, next-generation short reads, and Hi-C technology. The primary assembly has contig N50 size of 1.33 Mb and scaffold N50 size of 59.46 Mb, with a total size of 0.99 Gb, resolving most macrochromosomes into single scaffolds. A total of 23,058 genes and 10.71 Mb interspersed repeats were identified, constituting 30.31% and 10.71% of the common pheasant genome, respectively. Our phylogenetic analysis revealed that the common pheasant shared common ancestors with turkey about 24.7-34.5 million years ago (Ma). Rapidly evolved gene families, as well as branch-specific positively selected genes, indicate that calcium-related genes are potentially related to the adaptive and evolutionary change of the common pheasant. Interestingly, we found that the common pheasant has a unique major histocompatibility complex B locus (MHC-B) structure: three major inversions occurred in the sequence compared with chicken MHC-B. Furthermore, we detected signals of selection in five breeds of domestic common pheasant, several of which are production-oriented.

Probabilistic Estimation of Identity by Descent Segment Endpoints and Detection of Recent Selection

Most methods for fast detection of identity by descent (IBD) segments report identity by state segments without any quantification of the uncertainty in the endpoints and lengths of the IBD segments. We present a method for determining the posterior probability distribution of IBD segment endpoints. Our approach accounts for genotype errors, recent mutations, and gene conversions which disrupt DNA sequence identity within IBD segments, and it can be applied to large cohorts with whole-genome sequence or SNP array data. We find that our method's estimates of uncertainty are well calibrated for homogeneous samples. We quantify endpoint uncertainty for 77.7 billion IBD segments from 408,883 individuals of white British ancestry in the UK Biobank, and we use these IBD segments to find regions showing evidence of recent natural selection. We show that many spurious selection signals are eliminated by the use of unbiased estimates of IBD segment endpoints and a pedigree-based genetic map. Eleven of the twelve regions with the greatest evidence for recent selection in our scan have been identified as selected in previous analyses using different approaches. Our computationally efficient method for quantifying IBD segment endpoint uncertainty is implemented in the open source ibd-ends software package.

Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae

Background

The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family.

Results

Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3.

Conclusions

The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

MHC Genotyping by SSCP and Amplicon-Based NGS Approach in Chamois

Genes of the major histocompatibility complex (MHC) code for cell surface proteins essential for adaptive immunity. They show the most outstanding genetic diversity in vertebrates, which has been connected with various fitness traits and thus with the long-term persistence of populations. In this study, polymorphism of the MHC class II DRB locus was investigated in chamois with Single-Strand Conformation Polymorphism (SSCP)/Sanger genotyping and Ion Torrent S5 next-generation sequencing (NGS). From eight identified DRB variants in 28 individuals, five had already been described, and three were new, undescribed alleles. With conventional SSCP/Sanger sequencing, we were able to detect seven alleles, all of which were also detected with NGS. We found inconsistencies in the individual genotypes between the two methods, which were mainly caused by allelic dropout in the SSCP/Sanger method. Six out of 28 individuals were falsely classified as homozygous with SSCP/Sanger analysis. Overall, 25% of the individuals were identified as genotyping discrepancies between the two methods. Our results show that NGS technologies are better performing in sequencing highly variable regions such as the MHC, and they also have a higher detection capacity, thus allowing a more accurate description of the genetic composition, which is crucial for evolutionary and population genetic studies.

Adult survival in migratory caribou is negatively associated with MHC functional diversity

Genes of the major histocompatibility complex (MHC) are involved in acquired immunity in vertebrates. Only a few studies have investigated the fitness consequences of MHC gene diversity in wild populations. Here, we looked at the association between annual survival and body mass and MHC-DRB exon 2 (MHC-DRB) genetic diversity, obtained from high-throughput sequencing, in two declining migratory caribou (Rangifer tarandus) herds. To disentangle the potential direct and general effects of MHC-DRB genetic diversity, we compared different indices of diversity that were either based on DNA-sequence variation or on physicochemical divergence of the translated peptides, thereby covering a gradient of allelic-to-functional diversity. We found that (1) body mass was not related to MHC-DRB diversity or genotype, and (2) adult survival probability was negatively associated with point accepted mutation distance, a corrected distance that considers the likelihood of each amino acid substitution to be accepted by natural selection. In addition, we found no evidence of fluctuating selection over time on MHC-DRB diversity. We concluded that direct effects were involved in the negative relationship between MHC functional diversity and survival, although the mechanism underlying this result remains unclear. A possible explanation could be that individuals with higher MHC diversity suffer higher costs of immunity (immunopathology). Our results suggest that genetic diversity is not always beneficial even in genes that are likely to be strongly shaped by balancing selection.

Population-Specific Genetic and Expression Differentiation in Europeans

Much of the enormous phenotypic variation observed across human populations is thought to have arisen from events experienced as our ancestors peopled different regions of the world. However, little is known about the genes involved in these population-specific adaptations. Here, we explore this problem by simultaneously examining population-specific genetic and expression differentiation in four human populations. In particular, we derive a branch-based estimator of population-specific differentiation in four populations, and apply this statistic to single-nucleotide polymorphism and RNA-seq data from Italian, British, Finish, and Yoruban populations. As expected, genome-wide estimates of genetic and expression differentiation each independently recapitulate the known relationships among these four human populations, highlighting the utility of our statistic for identifying putative targets of population-specific adaptations. Moreover, genes with large copy number variations display elevated levels of population-specific genetic and expression differentiation, consistent with the hypothesis that gene duplication and deletion events are key reservoirs of adaptive variation. Further, many top-scoring genes are well-known targets of adaptation in Europeans, including those involved in lactase persistence and vitamin D absorption, and a handful of novel candidates represent promising avenues for future research. Together, these analyses reveal that our statistic can aid in uncovering genes involved in population-specific genetic and expression differentiation, and that such genes often play important roles in a diversity of adaptive and disease-related phenotypes in humans.

HLA Heterozygote Advantage against HIV-1 Is Driven by Quantitative and Qualitative Differences in HLA Allele-Specific Peptide Presentation

Pathogen-mediated balancing selection is regarded as a key driver of host immunogenetic diversity. A hallmark for balancing selection in humans is the heterozygote advantage at genes of the human leukocyte antigen (HLA), resulting in improved HIV-1 control. However, the actual mechanism of the observed heterozygote advantage is still elusive. HLA heterozygotes may present a broader array of antigenic viral peptides to immune cells, possibly resulting in a more efficient cytotoxic T-cell response. Alternatively, heterozygosity may simply increase the chance to carry the most protective HLA alleles, as individual HLA alleles are known to differ substantially in their association with HIV-1 control. Here, we used data from 6,311 HIV-1-infected individuals to explore the relative contribution of quantitative and qualitative aspects of peptide presentation in HLA heterozygote advantage against HIV. Screening the entire HIV-1 proteome, we observed that heterozygous individuals exhibited a broader array of HIV-1 peptides presented by their HLA class I alleles. In addition, viral load was negatively correlated with the breadth of the HIV-1 peptide repertoire bound by an individual's HLA variants, particularly at HLA-B. This suggests that heterozygote advantage at HLA-B is at least in part mediated by quantitative peptide presentation. We also observed higher HIV-1 sequence diversity among HLA-B heterozygous individuals, suggesting stronger evolutionary pressure from HLA heterozygosity. However, HLA heterozygotes were also more likely to carry certain HLA alleles, including the highly protective HLA-B*57:01 variant, indicating that HLA heterozygote advantage ultimately results from a combination of quantitative and qualitative effects in antigen presentation.

Divergent Allele Advantage at Human MHC Genes: Signatures of Past and Ongoing Selection

The highly polymorphic genes of the major histocompatibility complex (MHC) play a key role in adaptive immunity. Divergent allele advantage, a mechanism of balancing selection, is proposed to contribute to their exceptional polymorphism. It assumes that MHC genotypes with more divergent alleles allow for broader antigen-presentation to immune effector cells, by that increasing immunocompetence. However, the direct correlation between pairwise sequence divergence and the corresponding repertoire of bound peptides has not been studied systematically across different MHC genes. Here, we investigated this relationship for five key classical human MHC genes (human leukocyte antigen; HLA-A, -B, -C, -DRB1, and -DQB1), using allele-specific computational binding prediction to 118,097 peptides derived from a broad range of human pathogens. For all five human MHC genes, the genetic distance between two alleles of a heterozygous genotype was positively correlated with the total number of peptides bound by these two alleles. In accordance with the major antigen-presentation pathway of MHC class I molecules, HLA-B and HLA-C alleles showed particularly strong correlations for peptides derived from intracellular pathogens. Intriguingly, this bias coincides with distinct protein compositions between intra- and extracellular pathogens, possibly suggesting adaptation of MHC I molecules to present specifically intracellular peptides. Eventually, we observed significant positive correlations between an allele’s average divergence and its population frequency. Overall, our results support the divergent allele advantage as a meaningful quantitative mechanism through which pathogen-mediated selection leads to the evolution of MHC diversity.