Low Major Histocompatibility Complex Class II Variation in the Endangered Indo-Pacific Humpback Dolphin (Sousa chinensis): Inferences About the Role of Balancing Selection

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

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

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.

Comparing Class II MHC DRB3 Diversity in Colombian Simmental and Simbrah Cattle Across Worldwide Bovine Populations

The major histocompatibility complex (MHC) exerts great influence on responses to infectious diseases and vaccination due to its fundamental role in the adaptive immune system. Knowledge about MHC polymorphism distribution among breeds can provide insights into cattle evolution and diversification as well as population-based immune response variability, thus guiding further studies. Colombian Simmental and Simbrah cattle’s BoLA-DRB3 genetic diversity was compared to that of taurine and zebuine breeds worldwide to estimate functional diversity. High allele richness was observed for Simmental and Simbrah cattle; nevertheless, high homozygosity was associated with individual low sequence variability in both the β1 domain and the peptide binding region (PBR), thereby implying reduced MHC-presented peptide repertoire size. There were strong signals of positive selection acting on BoLA-DRB3 in all populations, some of which were poorly structured and displayed common alleles accounting for their high genetic similarity. PBR sequence correlation analysis suggested that, except for a few populations exhibiting some divergence at PBR, global diversity regarding potential MHC-presented peptide repertoire could be similar for the cattle populations analyzed here, which points to the retention of functional diversity in spite of the selective pressures imposed by breeding.

Population-level effects of polychlorinated biphenyl (PCB) exposure on highly vulnerable Indo-Pacific humpback dolphins from their largest habitat

While polychlorinated biphenyl (PCB)-related risks have been reported at the cellular, organ, and individual levels in some marine mammals, studies quantifying the PCB-associated population-level effects are limited. Here, we combined chemical analysis and individual-based model simulation to investigate the impact of PCBs on the Indo-Pacific humpback dolphin (sub)population from the Pearl River Estuary (PRE). An annual PCB accumulation rate of 0.29 ± 0.07 mg/kg lipid per year was estimated based on the measured age-specific male data as males continue to accumulate PCBs throughout their lifetime, without depurating contaminant loads. Using the Taiwan Strait dolphin population with low PCBs as a baseline, we compare our model simulations in PRE population to estimate relative population impacts of PCBs and other stressors. When using the current vital rates of the PRE dolphins which have been affected by PCBs and other stressors (e.g., underwater noise, prey limitation, etc.), our simulations revealed a substantial decline (8.1%) in the annual population growth rate (λ) of PRE metapopulation compared to baseline over the next 100 years. At the estimated PCB accumulation rate, the PCB-mediated effects on calf survival and immunity would cause a slight decline (0.9%) in λ relative to baseline. Our findings suggest a relatively limited impact of PCBs on the long-term survival of PRE dolphins among all stressors. However, it should be noted that even under model simulations where dietary PCBs were eliminated, humpback dolphins would still need a long time to reduce their PCB burdens to a relatively "safe" level through biological cycling. Considering that the baseline vital rates might also have been affected by PCBs and other stressors, our results are considered relative rather than absolute. This study provides a starting point for quantifying population-level consequences of contaminant exposure on humpback dolphins, although more efforts are needed to perfect this type of analysis.

Major Histocompatibility Complex Class II (DRB3) Genetic Diversity in Spanish Morucha and Colombian Normande Cattle Compared to Taurine and Zebu Populations

Bovine leukocyte antigens (BoLA) have been used as disease markers and immunological traits in cattle due to their primary role in pathogen recognition by the immune system. A higher MHC allele diversity in a population will allow presenting a broader peptide repertoire. However, loss of overall diversity due to domestication process can decrease a population's peptide repertoire. Within the context of zebu and taurine cattle populations, BoLA-DRB3 genetic diversity in Spanish Morucha and Colombian Normande cattle was analyzed and an approach to estimate functional diversity was performed. Sequence-based typing was used for identifying 29, 23, 27, and 28 alleles in Spanish Morucha, Nariño-, Boyacá-, and Cundinamarca-Normande cattle, respectively. These breeds had remarkably low heterozygosity levels and the Hardy-Weinberg principle revealed significant heterozygote deficiency. F_ST and D_A genetic distance showed that Colombian Normande populations had greater variability than other phenotypically homogeneous breeds, such as Holstein. It was also found that Spanish Morucha cattle were strongly differentiated from other cattle breeds. Spanish Morucha had greater divergence in the peptide-binding region regarding other cattle breeds. However, peptide-binding region covariation indicated that the potential peptide repertoire seemed equivalent among cattle breeds. Despite the genetic divergence observed, the extent of the potential peptide repertoire in the cattle populations studied appears to be similar and thus their pathogen recognition potential should be equivalent, suggesting that functional diversity might persist in the face of bottlenecks imposed by domestication and breeding.

Whole Genome Sequencing of Chinese White Dolphin (Sousa chinensis) for High-Throughput Screening of Antihypertensive Peptides

Chinese white dolphin (Sousa chinensis), also known as the Indo-Pacific humpback dolphin, has been classified as “Vulnerable” on the IUCN Red List of Threatened Species. It is a special cetacean species that lives in tropical and subtropical nearshore waters, with significant differences from other cetaceans. Here, we sequenced and assembled a draft genome of the Chinese white dolphin with a total length of 2.3 Gb and annotation of 18,387 protein-coding genes. Genes from certain expanded families are potentially involved in DNA replication and repairing, suggesting that they may be related to adaptation of this marine mammal to nearshore environments. We also discovered that its historical population had undergone a remarkable bottleneck incident before the Mindel glaciation. In addition, a comparative genomic survey on antihypertensive peptides (AHTPs) among five representative mammals with various residential habitats (such as remarkable differences in exogenous ion concentrations and sea depth) revealed that these small bioactive peptides were highly conserved among these examined mammals, and they had the most abundant hits in collagen subunit proteins, especially for two putative AHTP peptides Gly-Leu-Pro (GLP) and Leu-Gly-Pro (LGP). Our genome assembly will be a valuable resource for further genetic researches on adaptive ecology and conservation biology of cetaceans, and for in-depth investigations into bioactive peptides in aquatic and terrestrial mammals for development of peptide-based drugs to treat various human cardiovascular diseases.

The genome resources for conservation of Indo-Pacific humpback dolphin, Sousa chinensis

The Indo-Pacific humpback dolphin (Sousa chinensis), is a threatened marine mammal and belongs to the First Order of the National Key Protected Wild Aquatic Animals List in China. However, limited genomic information is available for studies of its population genetics and biological conservation. Here, we have assembled a genomic sequence of this species using a whole genome shotgun (WGS) sequencing strategy after a pilot low coverage genome survey. The total assembled genome size was 2.34 Gb: with a contig N50 of 67 kb and a scaffold N50 of 9 Mb (107.6-fold sequencing coverage). The S. chinensis genome contained 24,640 predicted protein-coding genes and had approximately 37% repeated sequences. The completeness of the genome assembly was evaluated by benchmarking universal single copy orthologous genes (BUSCOs): 94.3% of a total 4,104 expected mammalian genes were identified as complete, and 2.3% were identified as fragmented. This newly produced high-quality assembly and annotation of the genome will greatly promote the future studies of the genetic diversity, conservation and evolution.

Design Type(s)	sequence assembly objective • sequence annotation objective
Measurement Type(s)	whole genome sequencing assay
Technology Type(s)	DNA sequencing
Factor Type(s)
Sample Characteristic(s)	Sousa chinensis • skin of body

Machine-accessible metadata file describing the reported data (ISA-Tab format)

The Marine Mammal Class II Major Histocompatibility Complex Organization

Sirenians share with cetaceans and pinnipeds several convergent traits selected for the aquatic lifestyle. Living in water poses new challenges not only for locomotion and feeding but also for combating new pathogens, which may render the immune system one of the best tools aquatic mammals have for dealing with aquatic microbial threats. So far, only cetaceans have had their class II Major Histocompatibility Complex (MHC) organization characterized, despite the importance of MHC genes for adaptive immune responses. This study aims to characterize the organization of the marine mammal class II MHC using publicly available genomes. We located class II sequences in the genomes of one sirenian, four pinnipeds and eight cetaceans using NCBI-BLAST and reannotated the sequences using local BLAST search with exon and intron libraries. Scaffolds containing class II sequences were compared using dotplot analysis and introns were used for phylogenetic analysis. The manatee class II region shares overall synteny with other mammals, however most DR loci were translocated from the canonical location, past the extended class II region. Detailed analysis of the genomes of closely related taxa revealed that this presumed translocation is shared with all other living afrotherians. Other presumptive chromosome rearrangements in Afrotheria are the deletion of DQ loci in Afrosoricida and deletion of DP in E. telfairi. Pinnipeds share the main features of dog MHC: lack of a functional pair of DPA/DPB genes and inverted DRB locus between DQ and DO subregions. All cetaceans share the Cetartiodactyla inversion separating class II genes into two subregions: class IIa, with DR and DQ genes, and class IIb, with non-classic genes and a DRB pseudogene. These results point to three distinct and unheralded class II MHC structures in marine mammals: one canonical organization but lacking DP genes in pinnipeds; one bearing an inversion separating IIa and IIb subregions lacking DP genes found in cetaceans; and one with a translocation separating the most diverse class II gene from the MHC found in afrotherians and presumptive functional DR, DQ, and DP genes. Future functional research will reveal how these aquatic mammals cope with pathogen pressures with these divergent MHC organizations.

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.

Genomic organization and adaptive evolution of IGHC genes in marine mammals

•

The genomic organization of IGHC genes in marine mammal is similar to that of terrestrial relatives.

•

The number of the immunoglobulin heavy chain constant region genes vary among different mammals.

•

Different levels of selective pressures were detected between marine and terrestrial mammalian lineages.

Immunoglobulins are important elements of the adaptive immune system that bind to an immense variety of microbial antigens to neutralize infectivity and specify effector functions. In the present study, the immunoglobulin heavy chain constant region (IGHC) genes from marine mammals were identified and compared with those of their terrestrial relatives to explore their genomic organization and evolutionary characteristics. The genomic organization of marine mammal IGHC genes was shown to be conservative with other eutherian mammals. Stronger signals of positive selection on IGHC were revealed in terrestrial mammals than that in marine mammals with the branch-site model, displaying different selective pressure, which might suggest their divergent adaptations to contrasted environments.

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

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