Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

Partial molecular characterization, expression pattern and polymorphism analysis of MHC I genes in Chinese domestic goose (Anser cygnoides)

Major histocompatibility complex (MHC) allelic polymorphism is critically important for mediating antigen presentation in vertebrates. Presently, there are insufficient studies of MHC genetic diversity in domestic Anseriform birds. In this study, we analyzed the expression profile of MHC I genes and screened for MHC I exon 2 polymorphism in one domestic goose population from China using Illumina MiSeq sequencing. The results showed that four MHC I alleles (Ancy-IE2*09/*11/*13/*21) in one goose were identified based on cDNA cloning and sequencing using four primer combinations, and the varying number of cDNA clones implied that these four classical sequences showed differential expression patterns. Through next-generation sequencing, 27 alleles were obtained from 68 geese with 3-10 putative alleles per individual, indicating at least the existence of 5 MHC I loci in the goose. The marked excess of the non-synonymous over the synonymous substitution in the peptide-binding region (PBR) along 27 alleles and five positively selected sites (PSSs) detected around the PBR indicated that balancing selection might be the major force in shaping high MHC variation in the goose. Additionally, IA alleles displaying lower polymorphism were subject to less positive selection pressure than non-IA alleles with a higher level of polymorphism.

Transcriptomic Insights and the Development of Microsatellite Markers to Assess Genetic Diversity in the Broodstock Management of Litopenaeus stylirostris

The Pacific blue shrimp (Litopenaeus stylirostris) is a premium product in the international seafood market. However, intensified farming has increased disease incidence and reduced genetic diversity. In this study, we developed a transcriptome database for L. stylirostris and mined microsatellite markers to analyze their genetic diversity. Using the Illumina HiSeq 4000 platform, we identified 53,263 unigenes from muscle, hepatopancreas, the intestine, and lymphoid tissues. Microsatellite analysis identified 36,415 markers from 18,657 unigenes, predominantly dinucleotide repeats. Functional annotation highlighted key disease resistance pathways and enriched categories. The screening and PCR testing of 42 transcriptome-based and 58 literature-based markers identified 40 with successful amplification. The genotyping of 200 broodstock samples revealed that Na, Ho, He, PIC, and FIS values were 3, 0.54 ± 0.05, 0.43 ± 0.09, 0.41 ± 0.22, and 0.17 ± 0.27, respectively, indicating moderate genetic variability and significant inbreeding. Four universal microsatellite markers (CL1472.Contig13, CL517.Contig2, Unigene5692, and Unigene7147) were identified for precise diversity analysis in Pacific blue, Pacific white (Litopenaeus vannamei), and black tiger shrimps (Penaeus monodon). The transcriptome database supports the development of markers and functional gene analysis for selective breeding programs. Our findings underscore the need for an appropriate genetic management system to mitigate inbreeding depression, reduce disease susceptibility, and preserve genetic diversity in farmed shrimp populations.

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

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

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

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

Major histocompatibility complex class II DR and DQ evolution and variation in wild capuchin monkey species (Cebinae)

The major histocompatibility complex (MHC) is an important gene complex contributing to adaptive immunity. Studies of platyrrhine MHC have focused on identifying experimental models of immune system function in the equivalent Human Leukocyte Antigen (HLA). These genes have thus been explored primarily in captive platyrrhine individuals from research colonies. However, investigations of standing MHC variation and evolution in wild populations are essential to understanding its role in immunity, sociality and ecology. Capuchins are a promising model group exhibiting the greatest habitat diversity, widest diet breadth and arguably the most social complexity among platyrrhines, together likely resulting in varied immunological challenges. We use high-throughput sequencing to characterize polymorphism in four Class II DR and DQ exons for the first time in seven capuchin species. We find evidence for at least three copies for DQ genes and at least five for DRB, with possible additional unrecovered diversity. Our data also reveal common genotypes that are inherited across our most widely sampled population, Cebus imitator in Sector Santa Rosa, Costa Rica. Notably, phylogenetic analyses reveal that platyrrhine DQA sequences form a monophyletic group to the exclusion of all Catarrhini sequences examined. This result is inconsistent with the trans-species hypothesis for MHC evolution across infraorders in Primates and provides further evidence for the independent origin of current MHC genetic diversity in Platyrrhini. Identical allele sharing across cebid species, and more rarely genera, however, does underscore the complexity of MHC gene evolution and the need for more comprehensive assessments of allelic diversity and genome structure.

Genotyping and De Novo Discovery of Allelic Variants at the Brassicaceae Self-Incompatibility Locus from Short-Read Sequencing Data

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinarily high levels of polymorphism, both in terms of the number of functionally distinct S-alleles that segregate in SI species and in terms of their nucleotide sequence divergence. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we develop a powerful methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain of SRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected from natural populations of Arabidopsis halleri. Furthermore, using a de novo assembly approach of the filtered short reads, we obtain full-length sequences of eSRK even when the initial sequence in the database was only partial, and we discover putative new SRK alleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

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.

Design, delivery and perception of condition-dependent chemical signals in strepsirrhine primates: implications for human olfactory communication

The study of human chemical communication benefits from comparative perspectives that relate humans, conceptually and empirically, to other primates. All major primate groups rely on intraspecific chemosignals, but strepsirrhines present the greatest diversity and specialization, providing a rich framework for examining design, delivery and perception. Strepsirrhines actively scent mark, possess a functional vomeronasal organ, investigate scents via olfactory and gustatory means, and are exquisitely sensitive to chemically encoded messages. Variation in delivery, scent mixing and multimodality alters signal detection, longevity and intended audience. Based on an integrative, 19-species review, the main scent source used (excretory versus glandular) differentiates nocturnal from diurnal or cathemeral species, reflecting differing socioecological demands and evolutionary trajectories. Condition-dependent signals reflect immutable (species, sex, identity, genetic diversity, immunity and kinship) and transient (health, social status, reproductive state and breeding history) traits, consistent with socio-reproductive functions. Sex reversals in glandular elaboration, marking rates or chemical richness in female-dominant species implicate sexual selection of olfactory ornaments in both sexes. Whereas some compounds may be endogenously produced and modified (e.g. via hormones), microbial analyses of different odorants support the fermentation hypothesis of bacterial contribution. The intimate contexts of information transfer and varied functions provide important parallels applicable to olfactory communication in humans. This article is part of the Theo Murphy meeting issue 'Olfactory communication in humans'.

Lack of statistical power as a major limitation in understanding MHC-mediated immunocompetence in wild vertebrate populations

Disentangling the sources of variation in developing an effective immune response against pathogens is of major interest to immunoecology and evolutionary biology. To date, the link between immunocompetence and genetic variation at the major histocompatibility complex (MHC) has received little attention in wild animals, despite the key role of MHC genes in activating the adaptive immune system. Although several studies point to a link between MHC and immunocompetence, negative findings have also been reported. Such disparate findings suggest that limited statistical power might be affecting studies on this topic, owing to insufficient sample sizes and/or a generally small effect of MHC on the immunocompetence of wild vertebrates. To clarify this issue, we investigated the link between MHC variation and seven immunocompetence proxies in a large sample of barn owls and estimated the effect sizes and statistical power of this and published studies on this topic. We found that MHC poorly explained variation in immunocompetence of barn owls, with small-to-moderate associations between MHC and immunocompetence in owls (effect size: .1 ≥ r ≤ .3) similar to other vertebrates studied to date. Such small-to-moderate effects were largely associated with insufficient power, which was only sufficient (>0.8) to detect moderate-to-large effect sizes (r ≥ .3). Thus, studies linking MHC variation with immunocompetence in wild populations are underpowered to detect MHC effects, which are likely to be of generally small magnitude. Larger sample sizes (>200) will be required to achieve sufficient power in future studies aiming to robustly test for a link between MHC variation and immunocompetence.

Genetic variation at MHC class II loci influences both olfactory signals and scent discrimination in ring-tailed lemurs

BACKGROUND

Diversity at the Major Histocompatibility Complex (MHC) is critical to health and fitness, such that MHC genotype may predict an individual's quality or compatibility as a competitor, ally, or mate. Moreover, because MHC products can influence the components of bodily secretions, an individual's body odors may signal its MHC composition and influence partner identification or mate choice. Here, we investigated MHC-based signaling and recipient sensitivity by testing for odor-gene covariance and behavioral discrimination of MHC diversity and pairwise dissimilarity in a strepsirrhine primate, the ring-tailed lemur (Lemur catta).

METHODS

First, we coupled genotyping of the MHC class II gene, DRB, with gas chromatography-mass spectrometry of genital gland secretions to investigate if functional genetic diversity is signaled by the chemical diversity of lemur scent secretions. We also assessed if the chemical similarity between individuals correlated with their MHC-DRB similarity. Next, we assessed if lemurs discriminated this chemically encoded, genetic information in opposite-sex conspecifics.

RESULTS

We found that both sexes signaled overall MHC-DRB diversity and pairwise MHC-DRB similarity via genital secretions, but in a sex- and season-dependent manner. Additionally, the sexes discriminated absolute and relative MHC-DRB diversity in the genital odors of opposite-sex conspecifics, suggesting that lemur genital odors function to advertise genetic quality.

CONCLUSIONS

In summary, genital odors of ring-tailed lemurs provide honest information about an individual's absolute and relative MHC quality. Complementing evidence in humans and Old World monkeys, we suggest that reliance on scent signals to communicate MHC quality may be important across the primate lineage.

Challenges of next-generation sequencing in conservation management: Insights from long-term monitoring of corridor effects on the genetic diversity of mouse lemurs in a fragmented landscape

Long-term genetic monitoring of populations is essential for efforts aimed at preserving genetic diversity of endangered species. Here, we employ a framework of long-term genetic monitoring to evaluate the effects of fragmentation and the effectiveness of the establishment of corridors in restoring population connectivity and genetic diversity of mouse lemurs Microcebus ganzhorni. To this end, we supplement estimates of neutral genetic diversity with the assessment of adaptive genetic variability of the major histocompatibility complex (MHC). In addition, we address the challenges of long-term genetic monitoring of functional diversity by comparing the genotyping performance and estimates of MHC variability generated by single-stranded conformation polymorphism (SSCP)/Sanger sequencing with those obtained by high-throughput sequencing (next-generation sequencing [NGS], Illumina), an issue that is particularly relevant when previous work serves as a baseline for planning management strategies that aim to ensure the viability of a population. We report that SSCP greatly underestimates individual diversity and that discrepancies in estimates of MHC diversity attributable to the comparisons of traditional and NGS genotyping techniques can influence the conclusions drawn from conservation management scenarios. Evidence of migration among fragments in Mandena suggests that mouse lemurs are robust to the process of fragmentation and that the effect of corridors is masked by ongoing gene flow. Nonetheless, results based on a larger number of shared private alleles at neutral loci between fragment pairs found after the establishment of corridors in Mandena suggest that gene flow is augmented as a result of enhanced connectivity. Our data point out that despite low effective population size, M. ganzhorni maintains high individual heterozygosity at neutral loci and at MHC II DRB gene and that selection plays a predominant role in maintaining MHC diversity. These findings highlight the importance of long-term genetic monitoring in order to disentangle between the processes of drift and selection maintaining adaptive genetic diversity in small populations.

Determining Mhc-DRB profiles in wild populations of three congeneric true lemur species by noninvasive methods

The major histocompatibility complex (MHC) is a highly polymorphic and polygenic genomic region that plays a crucial role in immune-related diseases. Given the need for comparative studies on the variability of immunologically important genes among wild populations and species, we investigated the allelic variation of MHC class II DRB among three congeneric true lemur species: the red-fronted lemur (Eulemur rufifrons), red-bellied lemur (Eulemur rubriventer), and black lemur (Eulemur macaco). We noninvasively collected hair and faecal samples from these species across different regions in Madagascar. We assessed DRB exon 2 polymorphism with a newly developed primer set, amplifying nearly all non-synonymous codons of the antigen-binding sites. We defined 26 DRB alleles from 45 individuals (17 alleles from E. rufifrons (N = 18); 5 from E. rubriventer (N = 7); and 4 from E. macaco (N = 20). All detected alleles are novel and show high levels of nucleotide (26.8%) and non-synonymous codon polymorphism (39.4%). In these lemur species, we found neither evidence of a duplication of DRB genes nor a sharing of alleles among sympatric groups or allopatric populations of the same species. The non-sharing of alleles may be the result of a geographical separation over a long time span and/or different pathogen selection pressures. We found dN/dS rates > 1 in the functionally important antigen recognition sites, providing evidence for balancing selection. Especially for small and isolated populations, quantifying and monitoring DRB variation are recommended to establish successful conservation plans that mitigate the possible loss of immunogenetic diversity in lemurs.

Genotyping strategy matters when analyzing hypervariable major histocompatibility complex-Experience from a passerine bird

Genotyping of classical major histocompatibility complex (MHC) genes is challenging when they are hypervariable and occur in multiple copies. In this study, we used several different approaches to genotype the moderately variable MHC class I exon 3 (MHCIe3) and the highly polymorphic MHC class II exon 2 (MHCIIβe2) in the bluethroat (Luscinia svecica). Two family groups (eight individuals) were sequenced in replicates at both markers using Ion Torrent technology with both a single- and a dual-indexed primer structure. Additionally, MHCIIβe2 was sequenced on Illumina MiSeq. Allele calling was conducted by modifications of the pipeline developed by Sommer et al. (BMC Genomics, 14, 2013, 542) and the software AmpliSAS. While the different genotyping strategies gave largely consistent results for MHCIe3, with a maximum of eight alleles per individual, MHCIIβe2 was remarkably complex with a maximum of 56 MHCIIβe2 alleles called for one individual. Each genotyping strategy detected on average 50%-82% of all MHCIIβe2 alleles per individual, but dropouts were largely allele-specific and consistent within families for each strategy. The discrepancies among approaches indicate PCR biases caused by the platform-specific primer tails. Further, AmpliSAS called fewer alleles than the modified Sommer pipeline. Our results demonstrate that allelic dropout is a significant problem when genotyping the hypervariable MHCIIβe2. As these genotyping errors are largely nonrandom and method-specific, we caution against comparing genotypes across different genotyping strategies. Nevertheless, we conclude that high-throughput approaches provide a major advance in the challenging task of genotyping hypervariable MHC loci, even though they may not reveal the complete allelic repertoire.

Genetic wealth, population health: Major histocompatibility complex variation in captive and wild ring-tailed lemurs (Lemur catta)

Across species, diversity at the major histocompatibility complex (MHC) is critical to individual disease resistance and, hence, to population health; however, MHC diversity can be reduced in small, fragmented, or isolated populations. Given the need for comparative studies of functional genetic diversity, we investigated whether MHC diversity differs between populations which are open, that is experiencing gene flow, versus populations which are closed, that is isolated from other populations. Using the endangered ring-tailed lemur (Lemur catta) as a model, we compared two populations under long-term study: a relatively "open," wild population (n = 180) derived from Bezà Mahafaly Special Reserve, Madagascar (2003-2013) and a "closed," captive population (n = 121) derived from the Duke Lemur Center (DLC, 1980-2013) and from the Indianapolis and Cincinnati Zoos (2012). For all animals, we assessed MHC-DRB diversity and, across populations, we compared the number of unique MHC-DRB alleles and their distributions. Wild individuals possessed more MHC-DRB alleles than did captive individuals, and overall, the wild population had more unique MHC-DRB alleles that were more evenly distributed than did the captive population. Despite management efforts to maintain or increase genetic diversity in the DLC population, MHC diversity remained static from 1980 to 2010. Since 2010, however, captive-breeding efforts resulted in the MHC diversity of offspring increasing to a level commensurate with that found in wild individuals. Therefore, loss of genetic diversity in lemurs, owing to small founder populations or reduced gene flow, can be mitigated by managed breeding efforts. Quantifying MHC diversity within individuals and between populations is the necessary first step to identifying potential improvements to captive management and conservation plans.

Shared evolutionary origin of major histocompatibility complex polymorphism in sympatric lemurs

Genes of the major histocompatibility complex (MHC) play a central role in adaptive immune responses of vertebrates. They exhibit remarkable polymorphism, often crossing species boundaries with similar alleles or allelic motifs shared across species. This pattern may reflect parallel parasite-mediated selective pressures, either favouring the long maintenance of ancestral MHC allelic lineages across successive speciation events by balancing selection ("trans-species polymorphism"), or alternatively favouring the independent emergence of functionally similar alleles post-speciation via convergent evolution. Here, we investigate the origins of MHC similarity across several species of dwarf and mouse lemurs (Cheirogaleidae). We examined MHC class II variation in two highly polymorphic loci (DRB, DQB) and evaluated the overlap of gut-parasite communities in four sympatric lemurs. We tested for parasite-MHC associations across species to determine whether similar parasite pressures may select for similar MHC alleles in different species. Next, we integrated our MHC data with those previously obtained from other Cheirogaleidae to investigate the relative contribution of convergent evolution and co-ancestry to shared MHC polymorphism by contrasting patterns of codon usage at functional vs. neutral sites. Our results indicate that parasites shared across species may select for functionally similar MHC alleles, implying that the dynamics of MHC-parasite co-evolution should be envisaged at the community level. We further show that balancing selection maintaining trans-species polymorphism, rather than convergent evolution, is the primary mechanism explaining shared MHC sequence motifs between species that diverged up to 30 million years ago.

Integrating Genomic Data Sets for Knowledge Discovery: An Informed Approach to Management of Captive Endangered Species

Many endangered captive populations exhibit reduced genetic diversity resulting in health issues that impact reproductive fitness and quality of life. Numerous cost effective genomic sequencing and genotyping technologies provide unparalleled opportunity for incorporating genomics knowledge in management of endangered species. Genomic data, such as sequence data, transcriptome data, and genotyping data, provide critical information about a captive population that, when leveraged correctly, can be utilized to maximize population genetic variation while simultaneously reducing unintended introduction or propagation of undesirable phenotypes. Current approaches aimed at managing endangered captive populations utilize species survival plans (SSPs) that rely upon mean kinship estimates to maximize genetic diversity while simultaneously avoiding artificial selection in the breeding program. However, as genomic resources increase for each endangered species, the potential knowledge available for management also increases. Unlike model organisms in which considerable scientific resources are used to experimentally validate genotype-phenotype relationships, endangered species typically lack the necessary sample sizes and economic resources required for such studies. Even so, in the absence of experimentally verified genetic discoveries, genomics data still provides value. In fact, bioinformatics and comparative genomics approaches offer mechanisms for translating these raw genomics data sets into integrated knowledge that enable an informed approach to endangered species management.