A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies.

The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon 2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives.

Highly-contiguous bovine genomes underpin accurate functional analyses and updated nomenclature of MHC class I.

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle.

The Influence of Selection on MHC DQA and DQB Haplotypes in the Endemic New Zealand Hector's and Maui Dolphins.

Strong balancing selection on the major histocompatibility complex (MHC) can lead to different patterns in gene frequencies and neutral genomic variation within species. We investigated diversity and geographic structure of MHC genes DQA and DQB, as well as their inferred functional haplotypes, from 2 regional populations (East and West Coast) of the endangered Hector's dolphin (Cephalorhynchus hectori hectori) and the critically endangered Maui dolphin (Cephalorhynchus hectori maui) (West Coast, North Island), and contrasted these results with patterns from neutral microsatellites. The Maui had the lowest number of alleles for DQA (2) and DQB (3), consistent with strong genetic drift acting on this remnant population. However, the 2 retained DQA alleles are among the most divergent combinations of all 4 alleles found across the Hector's metapopulation, potentially reflecting the retention of divergent alleles due to balancing selection. The high frequency of the divergent DQB*04 allele also gave this population the highest nucleotide diversity for DQB. Strong differentiation was evident for DQA, DQB, and DQA-DQB haplotypes between the regional populations of Hector's dolphins (FST > 0.213) and both subspecies (FST > 0.311). Differentiation was generally greater than observed at neutral microsatellite loci, suggesting the influence of selection between geographically proximate East and West Coast populations. This might be the result of spatial differences in directional selection on those opposite coastlines. In addition, measures of the ratio of nonsynonymous to synonymous substitutions (dN/dS) were consistent with balancing selection over evolutionary time. Together, these results suggest a complex interplay of balancing selection, directional selection, local fidelity, and genetic drift.

The evolution of the natural killer complex; a comparison between mammals using new high-quality genome assemblies and targeted annotation.

Natural killer (NK) cells are a diverse population of lymphocytes with a range of biological roles including essential immune functions. NK cell diversity is in part created by the differential expression of cell surface receptors which modulate activation and function, including multiple subfamilies of C-type lectin receptors encoded within the NK complex (NKC). Little is known about the gene content of the NKC beyond rodent and primate lineages, other than it appears to be extremely variable between mammalian groups. We compared the NKC structure between mammalian species using new high-quality draft genome assemblies for cattle and goat; re-annotated sheep, pig, and horse genome assemblies; and the published human, rat, and mouse lemur NKC. The major NKC genes are largely in the equivalent positions in all eight species, with significant independent expansions and deletions between species, allowing us to propose a model for NKC evolution during mammalian radiation. The ruminant species, cattle and goats, have independently evolved a second KLRC locus flanked by KLRA and KLRJ, and a novel KLRH-like gene has acquired an activating tail. This novel gene has duplicated several times within cattle, while other activating receptor genes have been selectively disrupted. Targeted genome enrichment in cattle identified varying levels of allelic polymorphism between the NKC genes concentrated in the predicted extracellular ligand-binding domains. This novel recombination and allelic polymorphism is consistent with NKC evolution under balancing selection, suggesting that this diversity influences individual immune responses and may impact on differential outcomes of pathogen infection and vaccination.

Molecular species identification of Astrotoma agassizii from planktonic embryos: further evidence for a cryptic species complex.

The phrynophiurid brittle star Astrotoma agassizii is abundant in the cold temperate Magellanic region of South America and has a circumpolar Antarctic distribution. Three genetically distinct lineages were recently identified, with one in Antarctica geographically and genetically isolated from both South American lineages (Hunter R, Halanych KM. 2008. Evaluating connectivity in the brooding brittle star Astrotoma agassizii across the Drake Passage in the Southern Ocean. J Hered. 99:137-148.). Despite being an apparent brooding species, A. agassizii displayed a high genetic homogeneity at 2 mitochondrial markers (16s and COII) across a geographical range of more than 500 km along the Antarctic Peninsula. Here, using 16s ribosomal RNA sequences, we match a variety of early developmental stages (fertilized eggs, embryos; n = 12) collected from plankton samples in the Ross Sea to sequences of A. agassizii from the Antarctic Peninsula. The single 16s haplotype reported here is an identical match to one 16s haplotype found for A. agassizii from the Antarctic Peninsula, more than 5000 km away. Based on the regular occurrence of A. agassizii developmental stages in plankton samples, we propose that the Antarctic lineage of this species has a planktonic dispersive stage, with brooding restricted to the South American lineages. A different developmental mode would provide further evidence for cryptic speciation in this brittle star.

Using DNA barcoding and phylogenetics to identify Antarctic invertebrate larvae: Lessons from a large scale study.

Ecological studies of the diversity and distribution of marine planktonic larvae are increasingly depending on molecular methods for accurate taxonomic identification. The greater coverage of reference marine species on genetic databases such as GenBank and BoLD (Barcoding of Life Data Systems; www.boldystems.org); together with the decreasing costs for DNA sequencing have made large scale larval identification studies using molecular methods more feasible. Here, we present the development and implementation of a practical molecular approach to identify over 2000 individual marine invertebrate larvae that were collected in the Ross Sea, Antarctica, during the austral summer over five years (2002-2007) as part of the LGP (Latitudinal Gradient Project). Larvae for molecular ID were morphologically identified to belong to the Phyla Mollusca, Echinodermata, Nemertea and Annelida (Class Polychaeta), but also included unidentified early developmental stages which could not be assigned a specific taxon (e.g., eggs, blastulae). The use of a 100µm mesh plankton net makes this one of the first larval identification studies to simultaneously consider both embryos and larvae. Molecular identification methods included amplification of up to three molecular loci for each specimen, a pre-identification step using BLAST with GenBank, phylogenetic reconstructions and cross-validation of assigned Molecular Operational Taxonomic Units (MOTUs). This combined approach of morphological and molecular methods assigned about 700 individuals to 53 MOTUs, which were identified to the lowest possible taxonomic level. During the course of this long-term study we identified several procedural difficulties, including issues with the collection of larvae, locus amplification, contamination, assignment and validation of MOTUs. The practical guidelines that we describe here should greatly assist other researchers to conduct reliable molecular identification studies of larvae in the future.