Allele-Specific HLA Loss and Immune Escape in Lung Cancer Evolution.

Immune evasion is a hallmark of cancer. Losing the ability to present neoantigens through human leukocyte antigen (HLA) loss may facilitate immune evasion. However, the polymorphic nature of the locus has precluded accurate HLA copy-number analysis. Here, we present loss of heterozygosity in human leukocyte antigen (LOHHLA), a computational tool to determine HLA allele-specific copy number from sequencing data. Using LOHHLA, we find that HLA LOH occurs in 40% of non-small-cell lung cancers (NSCLCs) and is associated with a high subclonal neoantigen burden, APOBEC-mediated mutagenesis, upregulation of cytolytic activity, and PD-L1 positivity. The focal nature of HLA LOH alterations, their subclonal frequencies, enrichment in metastatic sites, and occurrence as parallel events suggests that HLA LOH is an immune escape mechanism that is subject to strong microenvironmental selection pressures later in tumor evolution. Characterizing HLA LOH with LOHHLA refines neoantigen prediction and may have implications for our understanding of resistance mechanisms and immunotherapeutic approaches targeting neoantigens. VIDEO ABSTRACT.

Germline HLA landscape does not predict efficacy of pembrolizumab monotherapy across solid tumor types.

We evaluated the impact of class I and class II human leukocyte antigen (HLA) genotypes, heterozygosity, and diversity on the efficacy of pembrolizumab. Seventeen pembrolizumab clinical trials across eight tumor types and one basket trial in patients with advanced solid tumors were included (n > 3,500 analyzed). Germline DNA was genotyped using a custom genotyping array. HLA diversity (measured by heterozygosity and evolutionary divergence) across class I loci was not associated with improved response to pembrolizumab, either within each tumor type evaluated or across all patients. Similarly, HLA heterozygosity at each class I and class II gene was not associated with response to pembrolizumab after accounting for the number of tests conducted. No conclusive association between HLA genotype and response to pembrolizumab was identified in this dataset. Germline HLA genotype or diversity alone is not an important independent determinant of response to pembrolizumab and should not be used for clinical decision-making in patients treated with pembrolizumab.

Demography and environment modulate the effects of genetic diversity on extinction risk in a butterfly metapopulation.

Linking genetic diversity to extinction is a common goal in genomic studies. Recently, a debate has arisen regarding the importance of genetic variation in conservation as some studies have failed to find associations between genome-wide genetic diversity and extinction risk. However, only rarely are genetic diversity and fitness measured together in the wild, and typically demographic history and environment are ignored. It is therefore difficult to infer whether a lack of an association is real or obscured by confounding factors. To address these shortcomings, we analyzed genetic data from 7,501 individuals with extinction data from 279 meadows and mortality of 1,742 larval nests in a butterfly metapopulation. We found a strong negative association between genetic diversity and extinction when considering only heterozygosity in models. However, this association disappeared when accounting for ecological covariates, suggesting a confounding between demography and genetics and a more complex role for heterozygosity in extinction risk. Modeling interactions between heterozygosity and demographic variables revealed that associations between extinction and heterozygosity were context-dependent. For example, extinction declined with increasing heterozygosity in large, but not currently small populations, although negative associations between heterozygosity, extinction, and mortality were detected in small populations with a recent history of decline. We conclude that low genetic diversity is an important predictor of extinction, predicting >25% increase in extinction beyond ecological factors in certain contexts. These results highlight that inferences about the importance of genetic diversity for population viability should not rely on genomic data alone but require investments in obtaining demographic and environmental data from natural populations.

Variants in ATRIP are associated with breast cancer susceptibility in the Polish population and UK Biobank.

Several breast cancer susceptibility genes have been discovered, but more are likely to exist. To identify additional breast cancer susceptibility genes, we used the founder population of Poland and performed whole-exome sequencing on 510 women with familial breast cancer and 308 control subjects. We identified a rare mutation in ATRIP (GenBank: NM_130384.3: c.1152_1155del [p.Gly385Ter]) in two women with breast cancer. At the validation phase, we found this variant in 42/16,085 unselected Polish breast cancer-affected individuals and in 11/9,285 control subjects (OR = 2.14, 95% CI = 1.13-4.28, p = 0.02). By analyzing the sequence data of the UK Biobank study participants (450,000 individuals), we identified ATRIP loss-of-function variants among 13/15,643 breast cancer-affected individuals versus 40/157,943 control subjects (OR = 3.28, 95% CI = 1.76-6.14, p < 0.001). Immunohistochemistry and functional studies showed the ATRIP c.1152_1155del variant allele is weakly expressed compared to the wild-type allele, and truncated ATRIP fails to perform its normal function to prevent replicative stress. We showed that tumors of women with breast cancer who have a germline ATRIP mutation have loss of heterozygosity at the site of ATRIP mutation and genomic homologous recombination deficiency. ATRIP is a critical partner of ATR that binds to RPA coating single-stranded DNA at sites of stalled DNA replication forks. Proper activation of ATR-ATRIP elicits a DNA damage checkpoint crucial in regulating cellular responses to DNA replication stress. Based on our observations, we conclude ATRIP is a breast cancer susceptibility gene candidate linking DNA replication stress to breast cancer.

A practical assembly guideline for genomes with various levels of heterozygosity.

Although current long-read sequencing technologies have a long-read length that facilitates assembly for genome reconstruction, they have high sequence errors. While various assemblers with different perspectives have been developed, no systematic evaluation of assemblers with long reads for diploid genomes with varying heterozygosity has been performed. Here, we evaluated a series of processes, including the estimation of genome characteristics such as genome size and heterozygosity, de novo assembly, polishing, and removal of allelic contigs, using six genomes with various heterozygosity levels. We evaluated five long-read-only assemblers (Canu, Flye, miniasm, NextDenovo and Redbean) and five hybrid assemblers that combine short and long reads (HASLR, MaSuRCA, Platanus-allee, SPAdes and WENGAN) and proposed a concrete guideline for the construction of haplotype representation according to the degree of heterozygosity, followed by polishing and purging haplotigs, using stable and high-performance assemblers: Redbean, Flye and MaSuRCA.

Global genomic instability caused by reduced expression of DNA polymerase ε in yeast.

SignificanceAlthough most studies of the genetic regulation of genome stability involve an analysis of mutations within the coding sequences of genes required for DNA replication or DNA repair, recent studies in yeast show that reduced levels of wild-type enzymes can also produce a mutator phenotype. By whole-genome sequencing and other methods, we find that reduced levels of the wild-type DNA polymerase ε in yeast greatly increase the rates of mitotic recombination, aneuploidy, and single-base mutations. The observed pattern of genome instability is different from those observed in yeast strains with reduced levels of the other replicative DNA polymerases, Pol α and Pol δ. These observations are relevant to our understanding of cancer and other diseases associated with genetic instability.

Abnormal biomarkers predict complex FAS or FADD defects missed by exome sequencing.

BACKGROUND: Elevated TCRαß+CD4-CD8- double-negative T cells (DNT) and serum biomarkers help identify FAS mutant patients with autoimmune lymphoproliferative syndrome (ALPS). However, in some patients with clinical features and biomarkers consistent with ALPS, germline or somatic FAS mutations cannot be identified on standard exon sequencing (ALPS-undetermined: ALPS-U). OBJECTIVE: We sought to explore whether complex genetic alterations in the FAS gene escaping standard sequencing or mutations in other FAS pathway-related genes could explain these cases. METHODS: Genetic analysis included whole FAS gene sequencing, copy number variation analysis, and sequencing of FAS cDNA and other FAS pathway-related genes. It was guided by FAS expression analysis on CD57+DNT, which can predict somatic loss of heterozygosity (sLOH). RESULTS: Nine of 16 patients with ALPS-U lacked FAS expression on CD57+DNT predicting heterozygous "loss-of-expression" FAS mutations plus acquired somatic second hits in the FAS gene, enriched in DNT. Indeed, 7 of 9 analyzed patients carried deep intronic mutations or large deletions in the FAS gene combined with sLOH detectable in DNT; 1 patient showed a FAS exon duplication. Three patients had reduced FAS expression, and 2 of them harbored mutations in the FAS promoter, which reduced FAS expression in reporter assays. Three of the 4 ALPS-U patients with normal FAS expression carried heterozygous FADD mutations with sLOH. CONCLUSION: A combination of serum biomarkers and DNT phenotyping is an accurate means to identify patients with ALPS who are missed by routine exome sequencing.

CRISPR/Cas gene editing in the human germline.

The ease and efficacy of CRISPR/Cas9 germline gene editing in animal models paved the way to human germline gene editing (HGGE), by which permanent changes can be introduced into the embryo. Distinct genes can be knocked out to examine their function during embryonic development. Alternatively, specific sequences can be introduced which can be applied to correct disease-causing mutations. To date, it has been shown that the success of HGGE is dependent on various experimental parameters and that various hurdles (i.e. loss-of-heterozygosity and mosaicism) need to be overcome before clinical applications should be considered. Due to the shortage of human germline material and the ethical constraints concerning HGGE, alternative models such as stem cells have been evaluated as well, in terms of their predictive value on the genetic outcome for HGGE approaches. This review will give an overview of the state of the art of HGGE in oocytes and embryos, and its accompanying challenges.

Phase plane dynamics of ERK phosphorylation.

The extracellular signal-regulated kinase (ERK) controls multiple critical processes in the cell and is deregulated in human cancers, congenital abnormalities, immune diseases, and neurodevelopmental syndromes. Catalytic activity of ERK requires dual phosphorylation by an upstream kinase, in a mechanism that can be described by two sequential Michaelis-Menten steps. The estimation of individual reaction rate constants from kinetic data in the full mechanism has proved challenging. Here, we present an analytically tractable approach to parameter estimation that is based on the phase plane representation of ERK activation and yields two combinations of six reaction rate constants in the detailed mechanism. These combinations correspond to the ratio of the specificities of two consecutive phosphorylations and the probability that monophosphorylated substrate does not dissociate from the enzyme before the second phosphorylation. The presented approach offers a language for comparing the effects of mutations that disrupt ERK activation and function in vivo. As an illustration, we use phase plane representation to analyze dual phosphorylation under heterozygous conditions, when two enzyme variants compete for the same substrate.

Genomic heterozygosity is associated with a lower risk of osteoarthritis.

BACKGROUND: Genomic heterozygosity has been shown to confer a health advantage in humans and play a protective role in complex diseases. Given osteoarthritis (OA) is a highly polygenic disease, we set out to determine if an association exists between OA and genomic heterozygosity. RESULTS: End-stage knee and hip OA patients and healthy controls were recruited from the Newfoundland and Labrador (NL) population. The Arthritis Research UK Osteoarthritis Genetics (arcOGEN) consortium database was utilized as a replication cohort. DNA was extracted from blood samples and genotyped. Individual rates of observed heterozygosity (HetRate) and heterozygosity excess (HetExcess) relative to the expected were mathematically derived, and standardized to a z-score. Logistic regression modeling was used to examine the association between OA and HetRate or HetExcess. A total of 559 knee and hip OA patients (mean age 66.5 years, body mass index (BMI) 33.7 kg/m2, and 55% females) and 118 healthy controls (mean age 56.4 years, BMI 29.5 kg/m2, and 59% female) were included in the NL cohort analysis. We found that OA had an inverse relationship with HetRate and HetExcess with odds ratios of 0.64 (95% CI: 0.45-0.91) and 0.65 (95% CI: 0.45-0.93) per standard deviation (SD), respectively. The arcOGEN data included 2,019 end-stage knee and hip OA patients and 2,029 healthy controls, validating our findings with HetRate and HetExcess odds ratios of 0.60 (95% CI: 0.56-0.64) and 0.44 (95% CI: 0.40-0.47) per SD, respectively. CONCLUSIONS: Our results are the first to clearly show evidence, from two separate cohorts, that reduced genomic heterozygosity confers a risk for the future development of OA.

A high heterozygosity genome assembly of Aedes albopictus enables the discovery of the association of PGANT3 with blood-feeding behavior.

The Asian tiger mosquito, Aedes albopictus, is a global invasive species, notorious for its role in transmitting dangerous human arboviruses such as dengue and Chikungunya. Although hematophagous behavior is repulsive, it is an effective strategy for mosquitoes like Aedes albopictus to transmit viruses, posing a significant risk to human health. However, the fragmented nature of the Ae. albopictus genome assembly has been a significant challenge, hindering in-depth biological and genetic studies of this mosquito. In this research, we have harnessed a variety of technologies and implemented a novel strategy to create a significantly improved genome assembly for Ae. albopictus, designated as AealbF3. This assembly boasts a completeness rate of up to 98.1%, and the duplication rate has been minimized to 1.2%. Furthermore, the fragmented contigs or scaffolds of AealbF3 have been organized into three distinct chromosomes, an arrangement corroborated through syntenic plot analysis, which compared the genetic structure of Ae. albopictus with that of Ae. aegypti. Additionally, the study has revealed a phylogenetic relationship suggesting that the PGANT3 gene is implicated in the hematophagous behavior of Ae. albopictus. This involvement was preliminarily substantiated through RNA interference (RNAi) techniques and behavioral experiment. In summary, the AealbF3 genome assembly will facilitate new biological insights and intervention strategies for combating this formidable vector of disease. The innovative assembly process employed in this study could also serve as a valuable template for the assembly of genomes in other insects characterized by high levels of heterozygosity.

Genetic diversity and signatures of selection in Icelandic horses and Exmoor ponies.

BACKGROUND: The Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds. RESULTS: HO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20). CONCLUSIONS: Evaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

Unraveling the Unusual Subgenomic Organization in the Neopolyploid Free-Living Flatworm Macrostomum lignano.

Whole genome duplication (WGD) is an evolutionary event resulting in a redundancy of genetic material. Different mechanisms of WGD, allo- or autopolyploidization, lead to distinct evolutionary trajectories of newly formed polyploids. Genome studies on such species are important for understanding the early stages of genome evolution. However, assembling neopolyploid is a challenging task due to the presence of 2 homologous (or homeologous) chromosome sets and therefore the existence of the extended paralogous regions in its genome. Post-WGD evolution of polyploids includes cytogenetic diploidization leading to the formation of species, whose polyploid origin might be hidden by disomic inheritance. Earlier we uncovered the hidden polyploid origin of the free-living flatworms of the genus Macrostomum (Macrostomum lignano, M. janickei, and M. mirumnovem). Cytogenetic diploidization in these species is accompanied by intensive chromosomal rearrangements including chromosomes fusions. In this study, we unravel the M. lignano genome organization through generation and sequencing of 2 sublines of the commonly used inbred line of M. lignano (called DV1) differing only in a copy number of the largest chromosome (MLI1). Using nontrivial assembly free comparative analysis of their genomes, we deciphered DNA sequences belonging to MLI1 and validated them by sequencing the pool of microdissected MLI1. Here we presented the uncommon mechanism of genome rediplodization of M. lignano, which consists of (i) presence of 3 subgenomes, which emerged via formation of large fused chromosomes and its variants, and (ii) sustaining their heterozygosity through inter- and intrachromosomal rearrangements.

Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus.

Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.

Telomere length and clonal chromosomal alterations in peripheral blood of patients with severe aplastic anaemia.

Severe aplastic anaemia (SAA) is a rare and life-threatening bone marrow failure disorder. We used data from the transplant outcomes in aplastic anaemia study to characterize mosaic chromosomal alterations (mCAs) in the peripheral blood of 738 patients with acquired SAA and evaluate their associations with telomere length (TL) and survival post-haematopoietic cell transplant (HCT). The median age at HCT was 20.4 years (range = 0.2-77.4). Patients with SAA had shorter TL than expected for their age (median TL percentile for age: 35.7th; range <1-99.99). mCAs were detected in 211 patients (28.6%), with chr6p copy-neutral loss of heterozygosity (6p-CNLOH) in 15.9% and chr7 loss in 3.0% of the patients; chrX loss was detected in 4.1% of female patients. Negative correlations between mCA cell fraction and measured TL (r = -0.14, p = 0.0002), and possibly genetically predicted TL (r = -0.07, p = 0.06) were noted. The post-HCT 3-year survival probability was low in patients with chr7 loss (39% vs. 72% in patients with chr6-CNLOH, 60% in patients with other mCAs and 70% in patients with no mCAs; p-log rank = 0.001). In multivariable analysis, short TL (p = 0.01), but not chr7 loss (p = 0.29), was associated with worse post-HCT survival. TL may guide clinical decisions in patients with SAA.

A case of hepatitis-associated aplastic anaemia following living-donor liver transplantation for fulminant hepatitis showing loss of heterozygosity in the 6p chromosome in the affected liver.

The mechanisms underlying hepatitis-associated aplastic anaemia (HAAA) that occurs several weeks after the development of acute hepatitis are unknown. A 20-year-old male developed HAAA following living-donor liver transplantation for fulminant hepatitis. The patient's leucocytes lacked HLA-class I due to loss of heterozygosity in the short arm of chromosome 6p (6pLOH). Interestingly, the patient's liver cells resected during the transplantation also exhibited 6pLOH that affected the same HLA haplotype as the leucocytes, suggesting that CD8+ T cells recognizing antigens presented by specific HLA molecules on liver cells may have attacked the haematopoietic stem cells of the patient, leading to the HAAA development.

Genetic and Phenotypic Consequences of Local Transitions between Sexual and Parthenogenetic Reproduction in the Wild.

AbstractTransitions from sexual to asexual reproduction have occurred in numerous lineages, but it remains unclear why asexual populations rarely persist. In facultatively parthenogenetic animals, all-female populations can arise when males are absent or become extinct, and such populations could help to understand the genetic and phenotypic changes that occur in the initial stages of transitions to asexuality. We investigated a naturally occurring spatial mosaic of mixed-sex and all-female populations of the facultatively parthenogenetic Australian phasmid Megacrania batesii. Analysis of single-nucleotide polymorphisms indicated multiple independent transitions between reproductive modes. All-female populations had much lower heterozygosity and allelic diversity than mixed-sex populations, but we found few consistent differences in fitness-related traits between population types. All-female populations exhibited more frequent and severe deformities in their (flight-incapable) wings but did not show higher rates of appendage loss. All-female populations also harbored more ectoparasites in swamp (but not beach) habitats. Reproductive mode explained little variation in female body size, fecundity, or egg hatch rate. Our results suggest that transitions to parthenogenetic reproduction can lead to dramatic genetic changes with little immediate effect on performance. All-female M. batesii populations appear to consist of high-fitness genotypes that might be able to thrive for many generations in relatively constant and benign environments but could be vulnerable to environmental challenges, such as increased parasite abundance.

Characteristics of Chinese breast cancer patients with double heterozygosity for BRCA1 and BRCA2 germline pathogenic variants.

PURPOSE: Despite of very rare, breast cancer patients with double heterozygosity (DH) variants in BRCA1 and BRCA2 genes have been identified in other ethnic groups and seem to be associated with distinctive phenotypes. However, little is known about the frequency and clinical characteristics of Chinese breast cancer patients with BRCA1/2 DH variants. METHODS: Four hundred and eleven unrelated patients with BRCA1 or BRCA2 pathogenic variants (PVs) were identified in a large series of unselected breast cancer patients. Another two siblings with metachronous bilateral breast cancer were referred for genetic counseling, after which BRCA1/2 DH variants were detected. RESULTS: Four unrelated breast cancer patients with BRCA1/2 DH were identified in the cohort of 411 patients with BRCA1 or BRCA2 PVs, the frequency of BRCA1/2 DH was 0.97%. In total, six BRCA1/2 DH patients from five families were found in this study. In two families, the hereditary pattern of DH was speculated to have originated from both sides of the family. BRCA1/2 DH patients were more likely to have a family history of breast cancer than patients with a BRCA1 (100% vs. 29.2%, P = 0.004) or BRCA2 (100% vs. 29.6%, P = 0.004) single PV. BRCA1/2 DH patients were more likely to be triple-negative breast tumors than patients with single BRCA2 PVs (66.7% vs. 14.1%, P = 0.020), which was comparable to the findings in patients with single BRCA1 PVs (66.7% vs. 56.9%, P = 1.00). CONCLUSION: Chinese patients with BRCA1/2 DH exhibit a high percentage of family history of breast cancer. The tumor pathological features of BRCA1/2 DH carriers are similar to those of BRCA1 PV carriers.

Multi-locus sequence typing and phylogenetics of Cryptococcus neoformans AD hybrids.

Hybrid AD strains of the human pathogenic Cryptococcus neoformans species complex have been reported from many parts of the world. However, their origin, diversity, and evolution are incompletely understood. In this study, we analyzed 102 AD hybrid strains representing 21 countries on five continents. For each strain, we obtained its mating type and its allelic sequences at each of the seven loci that have been used for genotyping haploid serotypes A and D strains of the species complex by the Cryptococcus research community. Our results showed that most AD hybrids exhibited loss of heterozygosity at one or more of the seven analyzed loci. Phylogenetic and population genetic analyses of the allelic sequences revealed multiple origins of the hybrids within each continent, dating back to one million years ago in Africa and up to the present in other continents. We found evidence for clonal reproduction and long-distance dispersal of these hybrids in nature. Comparisons with the global haploid serotypes A and D strains identified new alleles and new haploid multi-locus genotypes in AD hybrids, consistent with the presence of yet-to-be discovered genetic diversity in haploid populations of this species complex in nature. Together, our results indicate that AD hybrids can be effectively genotyped using the same multi-locus sequencing type approach as that established for serotypes A and D strains. Our comparisons of the AD hybrids among each other as well as with the global haploid serotypes A and D strains revealed novel genetic diversity as well as evidence for multiple origins and dynamic evolution of these hybrids in nature.

No evidence for a role of MHC class II genotype in the chemical encoding of heterozygosity and relatedness in Antarctic fur seals.

Despite decades of research, surprisingly little is known about the mechanism(s) by which an individual's genotype is encoded in odour. Many studies have focused on the role of the major histocompatibility complex (MHC) owing to its importance for survival and mate choice. However, the salience of MHC-mediated odours compared to chemicals influenced by the rest of the genome remains unclear, especially in wild populations where it is challenging to quantify and control for the effects of the genomic background. We addressed this issue in Antarctic fur seals by analysing skin swabs together with full-length MHC DQB II exon 2 sequences and data from 41 genome-wide distributed microsatellites. We did not find any effects of MHC relatedness on chemical similarity and there was also no relationship between MHC heterozygosity and chemical diversity. However, multilocus heterozygosity showed a significant positive association with chemical diversity, even after controlling for MHC heterozygosity. Our results appear to rule out a dominant role of the MHC in the chemical encoding of genetic information in a wild vertebrate population and highlight the need for genome-wide approaches to elucidate the mechanism(s) and specific genes underlying genotype-odour associations.