High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single-island endemic bird, but no evidence of MHC-based mate choice

Avian Island Radiations Shed Light on the Dynamics of Adaptive and Nonadaptive Radiation

Understanding the mechanisms underlying species formation and differentiation is a central goal of evolutionary biology and a formidable challenge. This understanding can provide valuable insights into the origins of the astonishing diversity of organisms living on our planet. Avian evolutionary radiations on islands have long fascinated biologists as they provide the ideal variation to study the ecological and evolutionary forces operating on the continuum between incipient lineages to complete speciation. In this review, we summarize the key insights gained from decades of research on adaptive and nonadaptive radiations of both extant and extinct insular bird species. We present a new comprehensive global list of potential avian radiations on oceanic islands, based on published island species checklists, taxonomic studies, and phylogenetic analyses. We demonstrate that our understanding of evolutionary processes is being greatly enhanced through the use of genomic tools. However, to advance the field, it is critical to complement this information with a solid understanding of the ecological and behavioral traits of both extinct and extant avian island species.

Using bioinformatics to investigate functional diversity: a case study of MHC diversity in koalas

Conservation genomics can greatly improve conservation outcomes of threatened populations, including those impacted by disease. Understanding diversity within immune gene families, including the major histocompatibility complex (MHC) and toll-like receptors (TLR), is important due to the role they play in disease resilience and susceptibility. With recent advancements in sequencing technologies and bioinformatic tools, the cost of generating high-quality sequence data has significantly decreased and made it possible to investigate diversity across entire gene families in large numbers of individuals compared to investigating only a few genes or a few populations previously. Here, we use the koala as a case study for investigating functional diversity across populations. We utilised previous target enrichment data and 438 whole genomes to firstly, determine the level of sequencing depth required to investigate MHC diversity and, secondly, determine the current level of diversity in MHC genes in koala populations. We determined for low complexity, conserved genes such as TLR genes 10 × sequencing depth is sufficient to reliably genotype more than 90% of variants, whereas for complex genes such as the MHC greater than 20 × and preferably 30 × sequencing depth is required. We used whole genome data to identify 270 biallelic SNPs across 24 MHC genes as well as copy number variation (CNV) within class I and class II genes and conduct supertype analysis. Overall, we have provided a bioinformatic workflow for investigating variation in a complex immune gene family from whole genome sequencing data and determined current levels of diversity within koala MHC genes.

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.

MHCtools 1.5: Analysis of MHC Sequencing Data in R

The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.

Structural variations in crocodile lizard populations provide novel insights into genomic variations in endangered animals

Population-scale genome resequencing of endangered animals may contribute to gaining an understanding of how genomes vary as population sizes become smaller, as well as the functional implications of such variation. In this study, we analysed structural variations and gene presence and absence variations in the genomes of population of the endangered crocodile lizards. We found that the frequencies of some genes showed significant differences between crocodile lizards in different regions, indicating the influence of environmental selection, as well as potential contributions from demography and isolation, in shaping gene presence and absence variations. The haplotype diversity of major histocompatibility complex (MHC) genes was also found to differ among crocodile lizards inhabiting different regions. These findings indicate that well-designed interbreeding of crocodile lizards from different regions may facilitate the exchange of genes between different lizard populations and increase the haplotype diversity of MHC genes, which may be beneficial for the survival of these lizards. Our findings in this study, based on differences in gene structural variation, provide new insights into genomic variation and may contribute to the conservation of endangered animals.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Genetic Diversity and Differentiation of MHC Class I Genes in Red-Crowned Crane Populations

The red-crowned crane (Grus japonensis) has been demoted to “vulnerable species” because its populations have apparently stabilized in Japan and Korea. Low variation and genetic drift may cause damage to the nascent recovery of the G. japonensis population. The major histocompatibility complex (MHC) is one of the most polymorphic gene families in the vertebrate genome and can reflect information on the adaptive evolution of endangered species. In this study, variations in MHC I exon 3 of captive G. japonensis in China were assessed and compared with those in cranes from Japan. Forty MHC alleles of 274 base pairs were isolated from 32 individuals from two captive populations in China. There was high variability in the nucleotide and amino acid composition, showing the proportion of polymorphic sites of 18.98 and 32.97%, respectively. Comparative analyses of the Chinese and Japanese populations based on 222 base pair sequences revealed more alleles and higher variation in the Chinese population. The lack of significant geographical differentiation of G. japonensis was supported by the genetic differentiation coefficient (0.04506) between the Chinese and Japanese populations. Positive selection of antigen-binding sites was observed, which contributed to maintaining the diversity of MHC class I genes. Phylogenetic analysis suggested the persistence of trans-species polymorphisms among MHC class I genes in Gruidae species. Our results may contribute to optimizing the management of G. japonensis populations and population recovery of this threatened species.

MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms

The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the R package MHCtools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed MHCtools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belonged to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of MHCtools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility

Whole‐genome sequencing of non‐model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data‐handling. Three such regions are the major histocompatibility complex (MHC), sex‐determining regions (SDRs) and the plant self‐incompatibility locus (S‐locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short‐read sequences against multiple alternative references improves genotyping comprehensiveness at the S‐locus thereby contributing to more accurate assessments of allelic frequencies. Long‐read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S‐locus haplotypes. These advances are adding value to short‐read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

Female preferences for male golden snub-nosed monkeys vary with male age and social context

Age is a key factor affecting sexual selection, as many physical and social traits are age-related. Although studies of primate mate choice often consider particular age-related traits, few consider the collective effects of male age. We tested the hypothesis that female golden snub-nosed monkeys Rhinopithecus roxellana prefer prime aged males (10–15 years) over younger and older males. We examined a habituated, provisioned troop during a 3-year study in the Qinling Mountains, China. Prime age males were more likely to be resident males of 1-male units (OMUs) than males of other ages. Since females are free to transfer between OMUs, the number of females per OMU can be indicative of female preferences. We examined the number of females per OMU, and found that it increased with resident male age up to 7–8 years, and declined after 12 years, such that prime age resident males had more females than other resident males. Females also initiated extra-unit copulations with high-ranking prime age males at significantly higher rates than with other males. Nevertheless, females tended to transfer from OMUs with high-ranking, older resident males to those with low-ranking, younger resident males. Thus, females appear to use different strategies when choosing social mates and extra-unit mates (i.e., different social contexts). We speculate that females may perceive early signs of aging in males and trade off the benefits and costs of high rank versus male senescence. This study lays the groundwork for future studies that examine possible direct and indirect benefits of such strategies.

Many loci make light work: High individual diversity despite low population diversity and random mating at class I MHC in a Critically Endangered island songbird

Multi-cellular organisms are under constant attack from parasites, making immune defence a critical aspect of fitness. In vertebrate animals, genes of the major histocompatibility complex (MHC) determine the breadth of pathogens to which individuals can respond. Having many MHC alleles can confer better protection against infectious disease, and balancing selection at MHC is widespread. Indeed, MHC loci are famously variable, with some populations harbouring thousands of alleles (Biedrzycka et al., 2018; Robinson, Soormally, Hayhurst, & Marsh, 2016). MHC has also long fascinated behavioural ecologists because mate choice-for example, preferring MHC-dissimilar partners-may amplify the effects of natural selection (Penn & Potts, 1999). But despite keen interest in the evolutionary ecology of MHC, extensive duplication (Minias, Pikus, Whittingham, & Dunn, 2019) has made these genes challenging to study. In a From the Cover article in this issue of Molecular Ecology, Stervander, Dierickx, Thorley, Brooke, and Westerdahl (2020) characterizes class I MHC in a Critically Endangered songbird, relating genotype to mate choice and survivorship. By inferring copy number and patterns of allelic co-segregation, the authors pave the way to elucidating the genomic architecture of MHC in this bottlenecked population. These insights help reconcile apparently counterintuitive findings: no effect of MHC genotype on mate choice or survival, and high MHC diversity within individuals despite low diversity at the population level. The latter finding is cause for optimism regarding conservation prospects. Moreover, these results suggest that ancient duplication events can have longstanding effects on the adaptive landscapes of natural and sexual selection.