Sequence polymorphism and geographical variation at a positively selected MHC-DRB gene in the finless porpoise (Neophocaena phocaenoides): implication for recent differentiation of the Yangtze finless porpoise?

Novel insights on aquatic mammal MHC evolution: Evidence from manatee DQB diversity

The low diversity in marine mammal major histocompatibility complex (MHC) appears to support the hypothesis of reduced pathogen selective pressure in aquatic systems compared to terrestrial environments. However, the lack of characterization of the aquatic and evolutionarily distant Sirenia precludes drawing more generalized conclusions. Therefore, we aimed to characterize the MHC DQB diversity of two manatee species and compare it with those reported for marine mammals. Our results identified 12 and 6 alleles in T. inunguis and T. manatus, respectively. Alleles show high rates of nonsynonymous substitutions, suggesting loci are evolving under positive selection. Among aquatic mammals, Pinnipeda DQB had smaller numbers of alleles, higher synonymous substitution rate, and a dN/dS ratio closer to 1, suggesting it may be evolving under more relaxed selection compared to fully aquatic mammals. This contradicts one of the predictions of the hypothesis that aquatic environments impose reduced pathogen pressure to mammalian immune system. These results suggest that the unique evolutionary trajectories of mammalian MHC may impose challenges in drawing ecoevolutionary conclusions from comparisons across distant vertebrate lineages.

Integrated analysis of blood mRNAs and microRNAs reveals immune changes with age in the Yangtze finless porpoise (Neophocaena asiaeorientalis)

Populations of Yangtze finless porpoises (YFPs) have rapidly declined in recent decades, raising the specter of extinction. In order to protect YFPs, a greater understanding of their biology is needed, including studying how their immune functioning changes with age. Here, we systematically studied the hematologic and biochemical parameters, as well as mRNAs and miRNAs profiles of old, adult, and young YFPs. The lymphocyte (LYMPH), neutrophils (NEUT) and eosinophils (EOS) counts in old YFPs were lower than those in young or adult YFPs. When comparing old to adult YFPs, the latter showed higher expression of genes associated with the innate and adaptive immune systems, including complement components, major histocompatibility complex, interleukins, TNF receptors, and chemokines/cytokines. When comparing old to young YFPs, the most striking difference was in higher toll-like receptor signaling in the latter. When comparing adult to young YFPs, the former exhibited higher expression of genes related to adaptive immunity and the FoxO signaling pathway, but lower expression of genes associated with the PI3K-Akt signaling pathway. Negative miRNA-mRNA interactions were predicted in comparisons of the old and adult (326), old and young (316), adult and young (211) groups. Overall, these results delineate a progression from early innate immune function dominance to adaptive immune function enhancement (young to adult) and deterioration (adult to old), and the changes in miRNAs profile correlate with the effects of age on immune functions. This study is the first to observe the changes of immune function of Yangtze finless porpoise with age using transcriptome method, and the study's findings are of great significance for protecting this endangered species.

The MHC Class Ia Genes in Chenfu's Treefrog (Zhangixalus chenfui) Evolved via Gene Duplication, Recombination, and Selection

Simple Summary

Amphibians, the first terrestrial vertebrates, provide materials for adaptive evolutionary studies, such as the evolution of the major histocompatibility complex (MHC). To date, various MHC evolutionary mechanisms have been identified in frogs, but more research is needed to determine the evolutionary mechanisms of the frog MHC. The main purpose of this study was to evaluate polymorphisms in the MHC class Ia genes of the Chenfu’s Treefrog. The MHC class Ia genes of the Chenfu’s Treefrog have high polymorphism. The mechanisms responsible for the formation of the polymorphisms include gene duplication, recombination, and selection.

Abstract

The molecular mechanisms underlying the evolution of adaptive immunity-related proteins can be deduced by a thorough examination of the major histocompatibility complex (MHC). Currently, in vertebrates, there is a relatively large amount of research on MHCs in mammals and birds. However, research related to amphibian MHC genes and knowledge about the evolutionary patterns is limited. This study aimed to isolate the MHC class I genes from Chenfu’s Treefrog (Zhangixalus chenfui) and reveal the underlying evolutionary processes. A total of 23 alleles spanning the coding region of MHC class Ia genes were identified in 13 individual samples. Multiple approaches were used to test and identify recombination from the 23 alleles. Amphibian MHC class Ia alleles, from NCBI, were used to construct the phylogenetic relationships in MEGA. Additionally, the partition strategy was adopted to construct phylogenetic relationships using MrBayes and MEGA. The sites of positive selection were identified by FEL, PAML, and MEME. In Chenfu’s Treefrog, we found that: (1) recombination usually takes place between whole exons of MHC class Ia genes; (2) there are at least 3 loci for MHC class Ia, and (3) the diversity of genes in MHC class Ia can be attributed to recombination, gene duplication, and positive selection. We characterized the evolutionary mechanisms underlying MHC class Ia genes in Chenfu’s Treefrog, and in so doing, broadened the knowledge of amphibian MHC systems.

Distinct evolution of toll-like receptor signaling pathway genes in cetaceans

BACKGROUND

The relatively rapid spread and diversity of marine pathogens posed an initial and ongoing challenge for cetaceans (whales, dolphins, and porpoises), descendants of terrestrial mammals that transitioned from land to sea approximately 56 million years ago. Toll-like receptors (TLRs) play important roles in regulating immunity against pathogen infections by detecting specific molecular patterns and activating a wide range of downstream signaling pathways. The ever-increasing catalogue of mammalian genomes offers unprecedented opportunities to reveal genetic changes associated with evolutionary and ecological processes.

OBJECTIVE

This study aimed to explore the molecular evolution of TLR signaling pathway genes in cetaceans.

METHODS

Genes involved in the TLR signaling pathway were retrieved by BLAST searches using human coding sequences as queries. We tested each gene for positive selection along the cetacean branches using PAML and Hyphy. Physicochemical property changes of amino acids at all positively selected residues were assessed by TreeSAAP and visualized with WebLogo. Bovine and dolphin TLR4 was assessed using human embryonic kidney cell line HEK293, which lacks TLR4 and its co-receptor MD-2.

RESULTS

We demonstrate that eight TLR signaling pathway genes are under positive selection in cetaceans. These include key genes in the response to Gram-negative bacteria: TLR4, CD14, and LY96 (MD-2). Moreover, 41 out of 65 positively selected sites were inferred to harbor substitution that dramatically changes the physicochemical properties of amino acids, with most of them situated in or adjacent to functional regions. We also found strong evidence that positive selection occurred in the lineage of the Yangtze finless porpoise, likely reflecting relatively recent adaptions to a freshwater milieu. Species-specific differences in TLR4 response were observed between cetacean and terrestrial species. Cetacean TLR4 was significantly less responsive to lipopolysaccharides from a terrestrial E. coli strain, possibly a reflection of the arms race of host-pathogen co-evolution faced by cetaceans in an aquatic environment.

CONCLUSION

This study provides further impetus for studies on the evolution and function of the cetacean immune system.

Divergent Evolution of TRC Genes in Mammalian Niche Adaptation

Mammals inhabit a wide variety of ecological niches, which in turn can be affected by various ecological factors, especially in relation to immunity. The canonical TRC repertoire (TRAC, TRBC, TRGC, and TRDC) codes C regions of T cell receptor chains that form the primary antigen receptors involved in the activation of cellular immunity. At present, little is known about the correlation between the evolution of mammalian TRC genes and ecological factors. In this study, four types canonical of TRC genes were identified from 37 mammalian species. Phylogenetic comparative methods (phyANOVA and PGLS) and selective pressure analyses among different groups of ecological factors (habitat, diet, and sociality) were carried out. The results showed that habitat was the major ecological factor shaping mammalian TRC repertoires. Specifically, trade-off between TRGC numbers and positive selection of TRAC and the balanced evolutionary rates between TRAC and TRDC genes were speculated as two main mechanisms in adaption to habitat and sociality. Overall, our study suggested divergent mechanisms for the evolution of TRCs, prompting mammalian immunity adaptions within diverse niches.

Widespread positive selection on cetacean TLR extracellular domain

Toll like receptors (TLRs), key members of innate immune system, can recognize a wide diversity of pathogens and initiate both innate and adaptive immune responses in vertebrate. Cetaceans must have faced new challenges of pathogens when their terrestrial relatives transitioned from the terrestrial to aquatic environment. Here, we sequenced the extracellular domain (ECD) of 10 TLRs in cetacean lineages because this region involved in the recognition of pathogens. A total of 148 sites ranging between 5-26 codons (0.01%-4.83%) were identified to be robust candidates of positive selection at the ECD of 10 TLRs. In addition, the majority (90.54%) of these positively selected codons were found to have radical amino acid changes, which strengthen the evidence of positive selection. Importantly, more radical amino acid changes in selected sites were enriched in the period of early evolutionary transition from land to semi-aquatic and from semi-aquatic to full-aquatic habitat, which might endow cetaceans with a faster adaptation to new pathogens as they transitioned into novel habitat. Interestingly, similar selective intensity was detected in both viral and non-viral TLRs in cetaceans, which was not in line with previous studies on primates and birds that reported stronger positive selection in non-viral TLRs than in viral TLRs. This result may be explained by the fact that cetaceans might have faced diversity of bacteria and viruses during its transitions from terrestrial to aquatic environment whereas both primates and birds probably being affected by only a restricted number of related viruses due to their homogeneous habitat.

Bridging immunogenetics and immunoproteomics: Model positional scanning library analysis for Major Histocompatibility Complex class II DQ in Tursiops truncatus

The Major Histocompatibility Complex (MHC) is a critical element in mounting an effective immune response in vertebrates against invading pathogens. Studies of MHC in wildlife populations have typically focused on assessing diversity within the peptide binding regions (PBR) of the MHC class II (MHC II) family, especially the DQ receptor genes. Such metrics of diversity, however, are of limited use to health risk assessment since functional analyses (where changes in the PBR are correlated to recognition/pathologies of known pathogen proteins), are difficult to conduct in wildlife species. Here we describe a means to predict the binding preferences of MHC proteins: We have developed a model positional scanning library analysis (MPSLA) by harnessing the power of mixture based combinatorial libraries to probe the peptide landscapes of distinct MHC II DQ proteins. The algorithm provided by NNAlign was employed to predict the binding affinities of sets of peptides generated for DQ proteins. These binding affinities were then used to retroactively construct a model Positional Scanning Library screen. To test the utility of the approach, a model screen was compared to physical combinatorial screens for human MHC II DP. Model library screens were generated for DQ proteins derived from sequence data from bottlenose dolphins from the Indian River Lagoon (IRL) and the Atlantic coast of Florida, and compared to screens of DQ proteins from Genbank for dolphin and three other cetaceans. To explore the peptide binding landscape for DQ proteins from the IRL, combinations of the amino acids identified as active were compiled into peptide sequence lists that were used to mine databases for representation in known proteins. The frequency of which peptide sequences predicted to bind the MHC protein are found in proteins from pathogens associated with marine mammals was found to be significant (p values <0.0001). Through this analysis, genetic variation in MHC (classes I and II) can now be associated with the binding repertoires of the expressed MHC proteins and subsequently used to identify target pathogens. This approach may be eventually applied to evaluate individual population and species risk for outbreaks of emerging diseases.

MHC class I BFIV gene polymorphisms in four Chinese native chicken breeds

The major histocompatibility complex (MHC) includes the most polymorphic genes in vertebrates, and balancing selection has been proposed as a main evolutionary force. Here we present one of the first data sets examining the genetic characteristics of chicken MHC I BFIV molecules in four Chinese native breeds, sourced from different regions in China. In all, 89 BFIV alleles were isolated from 102 individuals sampled, and 13 repeated alleles were observed. No significant correlation was found between genetic differentiation and geographical distance in the phylogenetic tree. BFIV genes exhibited a high level of nucleotide polymorphisms, and most of the polymorphic sites were located in the peptide-binding region (PBR) encoded in exons 2 and 3. A comparison of the three-dimensional structures of PBRs in chicken BFIV and human HLA-A molecules revealed evident structural and functional similarities. The results suggested that MHC I molecules had similar structural features in different species.

Organization and characteristics of the major histocompatibility complex class II region in the Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis)

Little is known about the major histocompatibility complex (MHC) in the genome of Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis) (YFP) or other cetaceans. In this study, a high-quality YFP bacterial artificial chromosome (BAC) library was constructed. We then determined the organization and characterization of YFP MHC class II region by screening the BAC library, followed by sequencing and assembly of positive BAC clones. The YFP MHC class II region consists of two segregated contigs (218,725 bp and 328,435 bp respectively) that include only eight expressed MHC class II genes, three pseudo MHC genes and twelve non-MHC genes. The YFP has fewer MHC class II genes than ruminants, showing locus reduction in DRB, DQA, DQB, and loss of DY. In addition, phylogenic and evolutionary analyses indicated that the DRB, DQA and DQB genes might have undergone birth-and-death evolution, whereas the DQB gene might have evolved under positive selection in cetaceans. These findings provide an essential foundation for future work, such as estimating MHC genetic variation in the YFP or other cetaceans. This work is the first report on the MHC class II region in cetaceans and offers valuable information for understanding the evolution of MHC genome in cetaceans.

Genetic Variation at Exon 2 of the MHC Class II DQB Locus in Blue Whale (Balaenoptera musculus) from the Gulf of California

The genes of the Major Histocompatibility Complex (MHC) play an important role in the vertebrate immune response and are among the most polymorphic genes known in vertebrates. In some marine mammals, MHC genes have been shown to be characterized by low levels of polymorphism compared to terrestrial taxa; this reduction in variation is often explained as a result of lower pathogen pressures in marine habitats. To determine if this same reduction in variation applies to the migratory population of blue whales (Balaenoptera musculus) that occurs in the Gulf of California, we genotyped a 172 bp fragment of exon 2 of the MHC Class II DQB locus for 80 members of this population. Twenty-two putatively functional DQB allotypes were identified, all of which were homologous with DQB sequences from other cetacean species. Up to 5 putative alleles per individual were identified, suggesting that gene duplication has occurred at this locus. Rates of non-synonymous to synonymous substitutions (ω) and maximum likelihood analyses of models of nucleotide variation provided potential evidence of ongoing positive selection at this exon. Phylogenetic analyses of DQB alleles from B. musculus and 16 other species of cetaceans revealed trans-specific conservation of MHC variants, suggesting that selection has acted on this locus over prolonged periods of time. Collectively our findings reveal that immunogenic variation in blue whales is comparable to that in terrestrial mammals, thereby providing no evidence that marine taxa are subject to reduced pathogen-induced selective pressures.

Evolution by selection, recombination, and gene duplication in MHC class I genes of two Rhacophoridae species

BACKGROUND

Comparison of major histocompatibility complex (MHC) genes across vertebrate species can reveal molecular mechanisms underlying the evolution of adaptive immunity-related proteins. As the first terrestrial tetrapods, amphibians deserve special attention because of their exposure to probably increased spectrum of microorganisms compared with ancestral aquatic fishes. Knowledge regarding the evolutionary patterns and mechanisms associated with amphibian MHC genes remains limited. The goal of the present study was to isolate MHC class I genes from two Rhacophoridae species (Rhacophorus omeimontis and Polypedates megacephalus) and examine their evolution.

RESULTS

We identified 27 MHC class I alleles spanning the region from exon 2 to 4 in 38 tree frogs. The available evidence suggests that these 27 sequences all belong to classical MHC class I (MHC Ia) genes. Although several anuran species only display one MHC class Ia locus, at least two or three loci were observed in P. megacephalus and R. omeimontis, indicating that the number of MHC class Ia loci varies among anuran species. Recombination events, which mainly involve the entire exons, played an important role in shaping the genetic diversity of the 27 MHC class Ia alleles. In addition, signals of positive selection were found in Rhacophoridae MHC class Ia genes. Amino acid sites strongly suggested by program to be under positive selection basically accorded with the putative antigen binding sites deduced from crystal structure of human HLA. Phylogenetic relationships among MHC class I alleles revealed the presence of trans-species polymorphisms.

CONCLUSIONS

In the two Rhacophoridae species (1) there are two or three MHC class Ia loci; (2) recombination mainly occurs between the entire exons of MHC class Ia genes; (3) balancing selection, gene duplication and recombination all contribute to the diversity of MHC class Ia genes. These findings broaden our knowledge on the evolution of amphibian MHC systems.

Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (Rhinopithecus roxellana)

BACKGROUND

Small, isolated populations often experience loss of genetic variation due to random genetic drift. Unlike neutral or nearly neutral markers (such as mitochondrial genes or microsatellites), major histocompatibility complex (MHC) genes in these populations may retain high levels of polymorphism due to balancing selection. The relative roles of balancing selection and genetic drift in either small isolated or bottlenecked populations remain controversial. In this study, we examined the mechanisms maintaining polymorphisms of MHC genes in small isolated populations of the endangered golden snub-nosed monkey (Rhinopithecus roxellana) by comparing genetic variation found in MHC and microsatellite loci. There are few studies of this kind conducted on highly endangered primate species.

RESULTS

Two MHC genes were sequenced and sixteen microsatellite loci were genotyped from samples representing three isolated populations. We isolated nine DQA1 alleles and sixteen DQB1 alleles and validated expression of the alleles. Lowest genetic variation for both MHC and microsatellites was found in the Shennongjia (SNJ) population. Historical balancing selection was revealed at both the DQA1 and DQB1 loci, as revealed by excess non-synonymous substitutions at antigen binding sites (ABS) and maximum-likelihood-based random-site models. Patterns of microsatellite variation revealed population structure. FST outlier analysis showed that population differentiation at the two MHC loci was similar to the microsatellite loci.

CONCLUSIONS

MHC genes and microsatellite loci showed the same allelic richness pattern with the lowest genetic variation occurring in SNJ, suggesting that genetic drift played a prominent role in these isolated populations. As MHC genes are subject to selective pressures, the maintenance of genetic variation is of particular interest in small, long-isolated populations. The results of this study may contribute to captive breeding and translocation programs for endangered species.

Character of chicken polymorphic major histocompatibility complex class II alleles of 3 Chinese local breeds

To better understand the major histocompatibility complex (MHC) genetic character of domestic birds, we sequenced and analyzed chicken MHC II (B-L) genes of 3 local chicken breeds, derived from 3 separate areas in China. We amplified cDNA sequences from 105 individuals, accounting for 35 alleles. Some of the same B-LB alleles with a high frequency were found in all samples. The putative B-L α-chain had few polymorphic sites, whereas the B-L β-chain had several polymorphic sites. Most of the mutation positions were located in the B-LB β1 domain encoded by exon 2, especially in the peptide-binding region. This indicated that the highly polymorphic peptide-binding region could potentiate binding diverse antigen epitopes. The comparison of 3-D molecule structures of chicken B-L and human HLA-DR1 revealed a distinctly structural similarity, but the chicken B-L molecule had more polymorphic sites than the human HLA-DR1 molecule, which presumably might be a mechanism to compensate for responding to a wider array of pathogens due to fewer loci for chicken. Moreover, some conserved sites in human and chicken MHC class II molecules reflected their common ancestry and similar functions. These results suggest that the chicken B-L gene showed more polymorphic sites and distinctly dominant trans-breed alleles, potentially to adapt to pathogens.

Considerable MHC diversity suggests that the functional extinction of baiji is not related to population genetic collapse

To further extend our understanding of the mechanism causing the current nearly extinct status of the baiji (Lipotes vexillifer), one of the most critically endangered species in the world, genetic diversity at the major histocompatibility complex (MHC) class II DRB locus was investigated in the baiji. Nine highly divergent DRB alleles were identified in 17 samples, with an average of 28.4 (13.2%) nucleotide difference and 16.7 (23.5%) amino acid difference between alleles. The unexpectedly high levels of DRB allelic diversity in the baiji may partly be attributable to its evolutionary adaptations to the freshwater environment which is regarded to have a higher parasite diversity compared to the marine environment. In addition, balancing selection was found to be the main mechanisms in generating sequence diversity at baiji DRB gene. Considerable sequence variation at the adaptive MHC genes despite of significant loss of neutral genetic variation in baiji genome might suggest that intense selection has overpowered random genetic drift as the main evolutionary forces, which further suggested that the critically endangered or nearly extinct status of the baiji is not an outcome of genetic collapse.

Relationship between Angiotensin-Converting Enzyme Insertion/Deletion Gene Polymorphism and Susceptibility of Minimal Change Nephrotic Syndrome: A Meta-Analysis

Aim. This meta-analysis was performed to evaluate the association between ACE I/D gene polymorphism and MCNS susceptibility. Method. A predefined literature search and selection of eligible relevant studies were performed to collect the data from electronic databases. Results. Six articles were identified for the analysis of association between ACE I/D gene polymorphism and MCNS risk, including 4 for Asians, one in Caucasian population and one for Africans. There was a markedly positive association between D allele or DD genotype and MCNS susceptibility in Asians (D: P = .01, DD: P = .02), but not for Caucasians and Africans (Caucasians: D: P = .16, DD: P = .98; Africans: D: P = .81, DD: P = .49). Furthermore, the II genotype seemed not to play a protective role against MCNS risk for Asians, Caucasians and Africans (P = .12, P = .09, P = .76, resp.). Interestingly, there was also significant association between ACE I/D gene polymorphism and MCNS susceptibility in overall populations (D: P = .007, DD: P = .04, II: P = .03). Conclusion. D allele or DD genotype might be a significant genetic molecular marker for MCNS susceptibility in Asians and overall populations, but not for Caucasians and Africans. More larger and rigorous genetic epidemiological investigations are required to further explore this association.