Spatial transcriptomics reveals molecular dysfunction associated with cortical Lewy pathology.

A key hallmark of Parkinson's disease (PD) is Lewy pathology. Composed of α-synuclein, Lewy pathology is found both in dopaminergic neurons that modulate motor function, and cortical regions that control cognitive function. Recent work has established the molecular identity of dopaminergic neurons susceptible to death, but little is known about cortical neurons susceptible to Lewy pathology or molecular changes induced by aggregates. In the current study, we use spatial transcriptomics to capture whole transcriptome signatures from cortical neurons with α-synuclein pathology compared to neurons without pathology. We find, both in PD and related PD dementia, dementia with Lewy bodies and in the pre-formed fibril α-synucleinopathy mouse model, that specific classes of excitatory neurons are vulnerable to developing Lewy pathology. Further, we identify conserved gene expression changes in aggregate-bearing neurons that we designate the Lewy-associated molecular dysfunction from aggregates (LAMDA) signature. Neurons with aggregates downregulate synaptic, mitochondrial, ubiquitin-proteasome, endo-lysosomal, and cytoskeletal genes and upregulate DNA repair and complement/cytokine genes. Our results identify neurons vulnerable to Lewy pathology in the PD cortex and describe a conserved signature of molecular dysfunction in both mice and humans.

Disruption of myelin structure and oligodendrocyte maturation in a pigtail macaque model of congenital Zika infection.

Zika virus (ZikV) infection during pregnancy can cause congenital Zika syndrome (CZS) and neurodevelopmental delay in non-microcephalic infants, of which the pathogenesis remains poorly understood. We utilized an established pigtail macaque maternal-to-fetal ZikV infection/exposure model to study fetal brain pathophysiology of CZS manifesting from ZikV exposure in utero. We found prenatal ZikV exposure led to profound disruption of fetal myelin, with extensive downregulation in gene expression for key components of oligodendrocyte maturation and myelin production. Immunohistochemical analyses revealed marked decreases in myelin basic protein intensity and myelinated fiber density in ZikV-exposed animals. At the ultrastructural level, the myelin sheath in ZikV-exposed animals showed multi-focal decompaction consistent with perturbation or remodeling of previously formed myelin, occurring concomitant with dysregulation of oligodendrocyte gene expression and maturation. These findings define fetal neuropathological profiles of ZikV-linked brain injury underlying CZS resulting from ZikV exposure in utero. Because myelin is critical for cortical development, ZikV-related perturbations in oligodendrocyte function may have long-term consequences on childhood neurodevelopment, even in the absence of overt microcephaly.

GSK-3 Inhibitor Elraglusib Enhances Tumor-Infiltrating Immune Cell Activation in Tumor Biopsies and Synergizes with Anti-PD-L1 in a Murine Model of Colorectal Cancer.

Glycogen synthase kinase-3 (GSK-3) is a serine/threonine kinase that has been implicated in numerous oncogenic processes. GSK-3 inhibitor elraglusib (9-ING-41) has shown promising preclinical and clinical antitumor activity across multiple tumor types. Despite promising early-phase clinical trial results, there have been limited efforts to characterize the potential immunomodulatory properties of elraglusib. We report that elraglusib promotes immune cell-mediated tumor cell killing of microsatellite stable colorectal cancer (CRC) cells. Mechanistically, elraglusib sensitized CRC cells to immune-mediated cytotoxicity and enhanced immune cell effector function. Using western blots, we found that elraglusib decreased CRC cell expression of NF-κB p65 and several survival proteins. Using microarrays, we discovered that elraglusib upregulated the expression of proapoptotic and antiproliferative genes and downregulated the expression of cell proliferation, cell cycle progression, metastasis, TGFß signaling, and anti-apoptotic genes in CRC cells. Elraglusib reduced CRC cell production of immunosuppressive molecules such as VEGF, GDF-15, and sPD-L1. Elraglusib increased immune cell IFN-γ secretion, which upregulated CRC cell gasdermin B expression to potentially enhance pyroptosis. Elraglusib enhanced immune effector function resulting in augmented granzyme B, IFN-γ, TNF-α, and TRAIL production. Using a syngeneic, immunocompetent murine model of microsatellite stable CRC, we evaluated elraglusib as a single agent or combined with immune checkpoint blockade (anti-PD-1/L1) and observed improved survival in the elraglusib and anti-PD-L1 group. Murine responders had increased tumor-infiltrating T cells, augmented granzyme B expression, and fewer regulatory T cells. Murine responders had reduced immunosuppressive (VEGF, VEGFR2) and elevated immunostimulatory (GM-CSF, IL-12p70) cytokine plasma concentrations. To determine the clinical significance, we then utilized elraglusib-treated patient plasma samples and found that reduced VEGF and BAFF and elevated IL-1 beta, CCL22, and CCL4 concentrations correlated with improved survival. Using paired tumor biopsies, we found that tumor-infiltrating immune cells had a reduced expression of inhibitory immune checkpoints (VISTA, PD-1, PD-L2) and an elevated expression of T-cell activation markers (CTLA-4, OX40L) after elraglusib treatment. These results address a significant gap in knowledge concerning the immunomodulatory mechanisms of GSK-3 inhibitor elraglusib, provide a rationale for the clinical evaluation of elraglusib in combination with immune checkpoint blockade, and are expected to have an impact on additional tumor types, besides CRC.

Ultra High-plex Spatial Proteogenomic Investigation of Giant Cell Glioblastoma Multiforme Immune Infiltrates Reveals Distinct Protein and RNA Expression Profiles.

A deeper understanding of complex biological processes, including tumor development and immune response, requires ultra high-plex, spatial interrogation of multiple "omes". Here we present the development and implementation of a novel spatial proteogenomic (SPG) assay on the GeoMx Digital Spatial Profiler platform with next-generation sequencing readout that enables ultra high-plex digital quantitation of proteins (>100-plex) and RNA (whole transcriptome, >18,000-plex) from a single formalin-fixed paraffin-embedded (FFPE) sample. This study highlighted the high concordance, R > 0.85 and <15% change in sensitivity between the SPG assay and the single-analyte assays on various cell lines and tissues from human and mouse. Furthermore, we demonstrate that the SPG assay was reproducible across multiple users. When used in conjunction with advanced cellular neighborhood segmentation, distinct immune or tumor RNA and protein targets were spatially resolved within individual cell subpopulations in human colorectal cancer and non-small cell lung cancer. We used the SPG assay to interrogate 23 different glioblastoma multiforme (GBM) samples across four pathologies. The study revealed distinct clustering of both RNA and protein based on pathology and anatomic location. The in-depth investigation of giant cell glioblastoma multiforme (gcGBM) revealed distinct protein and RNA expression profiles compared with that of the more common GBM. More importantly, the use of spatial proteogenomics allowed simultaneous interrogation of critical protein posttranslational modifications alongside whole transcriptomic profiles within the same distinct cellular neighborhoods. Significance: We describe ultra high-plex spatial proteogenomics; profiling whole transcriptome and high-plex proteomics on a single FFPE tissue section with spatial resolution. Investigation of gcGBM versus GBM revealed distinct protein and RNA expression profiles.

Spatial transcriptomics reveals molecular dysfunction associated with Lewy pathology.

Lewy pathology composed of α-synuclein is the key pathological hallmark of Parkinson's disease (PD), found both in dopaminergic neurons that control motor function, and throughout cortical regions that control cognitive function. Recent work has investigated which dopaminergic neurons are most susceptible to death, but little is known about which neurons are vulnerable to developing Lewy pathology and what molecular changes an aggregate induces. In the current study, we use spatial transcriptomics to selectively capture whole transcriptome signatures from cortical neurons with Lewy pathology compared to those without pathology in the same brains. We find, both in PD and in a mouse model of PD, that there are specific classes of excitatory neurons that are vulnerable to developing Lewy pathology in the cortex. Further, we identify conserved gene expression changes in aggregate-bearing neurons that we designate the Lewy-associated molecular dysfunction from aggregates (LAMDA) signature. This gene signature indicates that neurons with aggregates downregulate synaptic, mitochondrial, ubiquitin-proteasome, endo-lysosomal, and cytoskeletal genes and upregulate DNA repair and complement/cytokine genes. However, beyond DNA repair gene upregulation, we find that neurons also activate apoptotic pathways, suggesting that if DNA repair fails, neurons undergo programmed cell death. Our results identify neurons vulnerable to Lewy pathology in the PD cortex and identify a conserved signature of molecular dysfunction in both mice and humans.

Viral Protein Accumulation of Zika Virus Variants Links with Regulation of Innate Immunity for Differential Control of Viral Replication, Spread, and Response to Interferon.

Asian lineage Zika virus (ZIKV) strains emerged globally, causing outbreaks linked with critical clinical disease outcomes unless the virus is effectively restricted by host immunity. We have previously shown that retinoic acid-inducible gene-I (RIG-I) senses ZIKV to trigger innate immunity to direct interferon (IFN) production and antiviral responses that can control ZIKV infection. However, ZIKV proteins have been demonstrated to antagonize IFN. Here, we conducted in vitro analyses to assess how divergent prototypic ZIKV variants differ in virologic properties, innate immune regulation, and infection outcome. We comparatively assessed African lineage ZIKV/Dakar/1984/ArD41519 (ZIKV/Dakar) and Asian lineage ZIKV/Malaysia/1966/P6740 (ZIKV/Malaysia) in a human epithelial cell infection model. De novo viral sequence determination identified amino acid changes within the ZIKV/Dakar genome compared to ZIKV/Malaysia. Viral growth analyses revealed that ZIKV/Malaysia accumulated viral proteins and genome copies earlier and to higher levels than ZIKV/Dakar. Both ZIKV strains activated RIG-I/IFN regulatory factor (IRF3) and NF-κB pathways to induce inflammatory cytokine expression and types I and III IFNs. However, ZIKV/Malaysia, but not ZIKV/Dakar, potently blocked downstream IFN signaling. Remarkably, ZIKV/Dakar protein accumulation and genome replication were rescued in RIG-I knockout (KO) cells late in acute infection, resulting in ZIKV/Dakar-mediated blockade of IFN signaling. We found that RIG-I signaling specifically restricts viral protein accumulation late in acute infection where early accumulation of viral proteins in infected cells confers enhanced ability to limit IFN signaling, promoting viral replication and spread. Our results demonstrate that RIG-I-mediated innate immune signaling imparts restriction of ZIKV protein accumulation, which permits IFN signaling and antiviral actions controlling ZIKV infection. IMPORTANCE ZIKV isolates are classified under African or Asian lineages. Infection with emerging Asian lineage-derived ZIKV strains is associated with increased incidence of neurological symptoms that were not previously reported during infection with African or preemergent Asian lineage viruses. In this study, we utilized in vitro models to compare the virologic properties of and innate immune responses to two prototypic ZIKV strains from distinct lineages: African lineage ZIKV/Dakar and Asian lineage ZIKV/Malaysia. Compared to ZIKV/Dakar, ZIKV/Malaysia accumulates viral proteins earlier, replicates to higher levels, and robustly blocks IFN signaling during acute infection. Early accumulation of ZIKV/Malaysia NS5 protein confers enhanced ability to antagonize IFN signaling, dampening innate immune responses to promote viral spread. Our data identify the kinetics of viral protein accumulation as a major regulator of host innate immunity, influencing host-mediated control of ZIKV replication and spread. Importantly, these findings provide a novel framework for evaluating the virulence of emerging variants.

GSK-3 inhibitor elraglusib enhances tumor-infiltrating immune cell activation in tumor biopsies and synergizes with anti-PD-L1 in a murine model of colorectal cancer.

Inhibition of GSK-3 using small-molecule elraglusib has shown promising preclinical antitumor activity. Using in vitro systems, we found that elraglusib promotes immune cell-mediated tumor cell killing, enhances tumor cell pyroptosis, decreases tumor cell NF-κB-regulated survival protein expression, and increases immune cell effector molecule secretion. Using in vivo systems, we observed synergy between elraglusib and anti-PD-L1 in an immunocompetent murine model of colorectal cancer. Murine responders had more tumor-infiltrating T-cells, fewer tumor-infiltrating Tregs, lower tumorigenic circulating cytokine concentrations, and higher immunostimulatory circulating cytokine concentrations. To determine the clinical significance, we utilized human plasma samples from patients treated with elraglusib and correlated cytokine profiles with survival. Using paired tumor biopsies, we found that CD45+ tumor-infiltrating immune cells had lower expression of inhibitory immune checkpoints and higher expression of T-cell activation markers in post-elraglusib patient biopsies. These results introduce several immunomodulatory mechanisms of GSK-3 inhibition using elraglusib, providing a rationale for the clinical evaluation of elraglusib in combination with immunotherapy. Statement of significance: Pharmacologic inhibition of GSK-3 using elraglusib sensitizes tumor cells, activates immune cells for increased anti-tumor immunity, and synergizes with anti-PD-L1 immune checkpoint blockade. These results introduce novel biomarkers for correlations with response to therapy which could provide significant clinical utility and suggest that elraglusib, and other GSK-3 inhibitors, should be evaluated in combination with immune checkpoint blockade.

Cytomegalovirus-vaccine-induced unconventional T cell priming and control of SIV replication is conserved between primate species.

Strain 68-1 rhesus cytomegalovirus expressing simian immunodeficiency virus (SIV) antigens (RhCMV/SIV) primes MHC-E-restricted CD8+ T cells that control SIV replication in 50%-60% of the vaccinated rhesus macaques. Whether this unconventional SIV-specific immunity and protection is unique to rhesus macaques or RhCMV or is intrinsic to CMV remains unknown. Here, using cynomolgus CMV vectors expressing SIV antigens (CyCMV/SIV) and Mauritian cynomolgus macaques, we demonstrate that the induction of MHC-E-restricted CD8+ T cells requires matching CMV to its host species. RhCMV does not elicit MHC-E-restricted CD8+ T cells in cynomolgus macaques. However, cynomolgus macaques vaccinated with species-matched 68-1-like CyCMV/SIV mounted MHC-E-restricted CD8+ T cells, and half of the vaccinees stringently controlled SIV post-challenge. Protected animals manifested a vaccine-induced IL-15 transcriptomic signature that is associated with efficacy in rhesus macaques. These findings demonstrate that the ability of species-matched CMV vectors to elicit MHC-E-restricted CD8+ T cells that are required for anti-SIV efficacy is conserved in nonhuman primates, and these data support the development of HCMV/HIV for a prophylactic HIV vaccine.

ZIKV can infect human term placentas in the absence of maternal factors.

Zika virus infection can result in devastating pregnancy outcomes when it crosses the placental barrier. For human pregnancies, the mechanisms of vertical transmission remain enigmatic. Utilizing a human placenta-cotyledon perfusion model, we examined Zika virus exposure in the absence of maternal factors. To distinguish responses related to viral infection vs. recognition, we evaluated cotyledons perfused with either active or inactivated Zika virus. Active Zika virus exposure resulted in infection, cell death and syncytium injury. Pathology corresponded with transcriptional changes related to inflammation and innate immunity. Inactive Zika virus exposure also led to syncytium injury and related changes in gene expression but not cell death. Our observations reveal pathologies and innate immune responses that are dependent on infection or virus placenta interactions independent of productive infection. Importantly, our findings indicate that Zika virus can infect and compromise placentas in the absence of maternal humoral factors that may be protective.

Interleukin-15 response signature predicts RhCMV/SIV vaccine efficacy.

Simian immunodeficiency virus (SIV) challenge of rhesus macaques (RMs) vaccinated with strain 68-1 Rhesus Cytomegalovirus (RhCMV) vectors expressing SIV proteins (RhCMV/SIV) results in a binary outcome: stringent control and subsequent clearance of highly pathogenic SIV in ~55% of vaccinated RMs with no protection in the remaining 45%. Although previous work indicates that unconventionally restricted, SIV-specific, effector-memory (EM)-biased CD8+ T cell responses are necessary for efficacy, the magnitude of these responses does not predict efficacy, and the basis of protection vs. non-protection in 68-1 RhCMV/SIV vector-vaccinated RMs has not been elucidated. Here, we report that 68-1 RhCMV/SIV vector administration strikingly alters the whole blood transcriptome of vaccinated RMs, with the sustained induction of specific immune-related pathways, including immune cell, toll-like receptor (TLR), inflammasome/cell death, and interleukin-15 (IL-15) signaling, significantly correlating with subsequent vaccine efficacy. Treatment of a separate RM cohort with IL-15 confirmed the central involvement of this cytokine in the protection signature, linking the major innate and adaptive immune gene expression networks that correlate with RhCMV/SIV vaccine efficacy. This change-from-baseline IL-15 response signature was also demonstrated to significantly correlate with vaccine efficacy in an independent validation cohort of vaccinated and challenged RMs. The differential IL-15 gene set response to vaccination strongly correlated with the pre-vaccination activity of this pathway, with reduced baseline expression of IL-15 response genes significantly correlating with higher vaccine-induced induction of IL-15 signaling and subsequent vaccine protection, suggesting that a robust de novo vaccine-induced IL-15 signaling response is needed to program vaccine efficacy. Thus, the RhCMV/SIV vaccine imparts a coordinated and persistent induction of innate and adaptive immune pathways featuring IL-15, a known regulator of CD8+ T cell function, that support the ability of vaccine-elicited unconventionally restricted CD8+ T cells to mediate protection against SIV challenge.

Developing a transcriptomic framework for testing testosterone-mediated handicap hypotheses.

Sexually selected traits are hypothesized to be honest signals of individual quality due to the costs associated with their maintenance, development, and/or production. Testosterone, a sex steroid associated with the development and/or production of sexually selected traits, has been proposed to enforce the honesty of sexually selected traits via its immunosuppressive effects (i.e., the Immunocompetence Handicap Hypothesis) and/or by influencing an individual's exposure/susceptibility to oxidative stress (i.e., the Oxidation Handicap Hypothesis). Previous work testing these hypotheses has primarily focused on physiological measurements of immunity or oxidative stress, but little is known about the molecular pathways by which testosterone could influence immunity and/or oxidative stress pathways. To further understand the transcriptomic consequences of experimentally elevated testosterone in the context of handicap hypotheses, we used previously published RNA-seq data from studies that measured the transcriptome of individuals treated with either a testosterone-filled or an empty (i.e., control) implant. Two studies encompassing three species of bird and three tissue types fit our selection criteria and we reanalyzed the data using weighted gene co-expression network analysis. Testosterone-treated individuals exhibited signatures of immunosuppression and our results describe the molecular pathways underlying this effect. We also provide some evidence to suggest that the transcriptomic signature of immunosuppression is evolutionarily conserved between the three species. While our results provide no evidence to suggest testosterone mediates handicaps via pathways associated with oxidative stress, they do support the hypothesis that testosterone enforces the honesty of sexually-selected traits by influencing an individual's immunocompetence. Overall, this study develops a framework for testing testosterone-mediated handicap hypotheses and provides guidelines for future integrative and comparative studies focused on the proximate mechanisms mediating sexually selected traits.

Parent and offspring genotypes influence gene expression in early life.

Parents can have profound effects on offspring fitness. Little, however, is known about the mechanisms through which parental genetic variation influences offspring physiology in natural systems. White-throated sparrows (Zonotrichia albicollis, WTSP) exist in two genetic morphs, tan and white, controlled by a large polymorphic supergene. Morphs mate disassortatively, resulting in two pair types: tan male × white female (T × W) pairs, which provide biparental care and white male × tan female (W × T) pairs, which provide female-biased care. To investigate how parental composition impacts offspring, we performed RNA-seq on whole blood of WTSP nestlings sampled from nests of both pair types. Parental pair type had a large effect on nestling gene expression, with 881 genes differentially expressed (DE) and seven correlated gene coexpression modules. The DE genes and modules expressed at higher levels in W × T nests with female-biased parental care function in metabolism and stress-related pathways resulting from the overrepresentation of proteolysis and stress-response genes (e.g., SOD2, NR3C1). These results show that parental genotypes and/or associated behaviours influence nestling physiology, and highlight avenues of further research investigating the ultimate implications for the maintenance of this polymorphism. Nestlings also exhibited morph-specific gene expression, with 92 differentially expressed genes, comprising immunity genes and genes encompassed by the supergene. Remarkably, we identified the same regulatory hub genes in these blood-derived expression networks as were previously identified in adult WTSP brains (EPM2A, BPNT1, TAF5L). These hub genes were located within the supergene, highlighting the importance of this gene complex in structuring regulatory networks across diverse tissues.

Light pollution increases West Nile virus competence of a ubiquitous passerine reservoir species.

Among the many anthropogenic changes that impact humans and wildlife, one of the most pervasive but least understood is light pollution. Although detrimental physiological and behavioural effects resulting from exposure to light at night are widely appreciated, the impacts of light pollution on infectious disease risk have not been studied. Here, we demonstrate that artificial light at night (ALAN) extends the infectious-to-vector period of the house sparrow (Passer domesticus), an urban-dwelling avian reservoir host of West Nile virus (WNV). Sparrows exposed to ALAN maintained transmissible viral titres for 2 days longer than controls but did not experience greater WNV-induced mortality during this window. Transcriptionally, ALAN altered the expression of gene regulatory networks including key hubs (OASL, PLBD1 and TRAP1) and effector genes known to affect WNV dissemination (SOCS). Despite mounting anti-viral immune responses earlier, transcriptomic signatures indicated that ALAN-exposed individuals probably experienced pathogen-induced damage and immunopathology, potentially due to evasion of immune effectors. A simple mathematical modelling exercise indicated that ALAN-induced increases of host infectious-to-vector period could increase WNV outbreak potential by approximately 41%. ALAN probably affects other host and vector traits relevant to transmission, and additional research is needed to advise the management of zoonotic diseases in light-polluted areas.

Transcriptional response to West Nile virus infection in the zebra finch (Taeniopygia guttata).

West Nile virus (WNV) is a widespread arbovirus that imposes a significant cost to both human and wildlife health. WNV exists in a bird-mosquito transmission cycle in which passerine birds act as the primary reservoir host. As a public health concern, the mammalian immune response to WNV has been studied in detail. Little, however, is known about the avian immune response to WNV. Avian taxa show variable susceptibility to WNV and what drives this variation is unknown. Thus, to study the immune response to WNV in birds, we experimentally infected captive zebra finches (Taeniopygia guttata). Zebra finches provide a useful model, as like many natural avian hosts they are moderately susceptible to WNV and thus provide sufficient viremia to infect mosquitoes. We performed RNAseq in spleen tissue during peak viremia to provide an overview of the transcriptional response. In general, we find strong parallels with the mammalian immune response to WNV, including upregulation of five genes in the Rig-I-like receptor signalling pathway, and offer insights into avian-specific responses. Together with complementary immunological assays, we provide a model of the avian immune response to WNV and set the stage for future comparative studies among variably susceptible populations and species.

Divergence and Functional Degradation of a Sex Chromosome-like Supergene.

A major challenge in biology is to understand the genetic basis of adaptation. One compelling idea is that groups of tightly linked genes (i.e., "supergenes" [1, 2]) facilitate adaptation in suites of traits that determine fitness. Despite their likely importance, little is known about how alternate supergene alleles arise and become differentiated, nor their ultimate fate within species. Herein we address these questions by investigating the evolutionary history of a supergene in white-throated sparrows, Zonotrichia albicollis. This species comprises two morphs, tan and white, that differ in pigmentation and components of social behavior [3-5]. Morph is determined by alternative alleles at a balanced >100-Mb inversion-based supergene, providing a unique system for studying gene-behavior relationships. Using over two decades of field data, we document near-perfect disassortative mating among morphs, as well as the fitness consequences of rare assortative mating. We use de novo whole-genome sequencing coupled with population- and phylogenomic data to show that alternate supergene alleles are highly divergent at over 1,000 genes and that these alleles originated prior to the split of Z. albicollis from its sister species and may be polymorphic in Z. albicollis due to a past hybridization event. We provide evidence that the "white" allele may be degrading, similar to neo-Y/W sex chromosomes. We further show that the "tan" allele has surprisingly low levels of genetic diversity yet does not show several canonical signatures of recurrent positive selection. We discuss these results in the context of the origin, molecular evolution, and possible fate of this remarkable polymorphism.

High major histocompatibility complex class I polymorphism despite bottlenecks in wild and domesticated populations of the zebra finch (Taeniopygia guttata).

BACKGROUND: Two subspecies of zebra finch, Taeniopygia guttata castanotis and T. g. guttata are native to Australia and the Lesser Sunda Islands, respectively. The Australian subspecies has been domesticated and is now an important model system for research. Both the Lesser Sundan subspecies and domesticated Australian zebra finches have undergone population bottlenecks in their history, and previous analyses using neutral markers have reported reduced neutral genetic diversity in these populations. Here we characterize patterns of variation in the third exon of the highly variable major histocompatibility complex (MHC) class I α chain. As a benchmark for neutral divergence, we also report the first mitochondrial NADH dehydrogenase 2 (ND2) sequences in this important model system. RESULTS: Despite natural and human-mediated population bottlenecks, we find that high MHC class I polymorphism persists across all populations. As expected, we find higher levels of nucleotide diversity in the MHC locus relative to neutral loci, and strong evidence of positive selection acting on important residues forming the peptide-binding region (PBR). Clear population differentiation of MHC allele frequencies is also evident, and this may be due to adaptation to new habitats and associated pathogens and/or genetic drift. Whereas the MHC Class I locus shows broad haplotype sharing across populations, ND2 is the first locus surveyed to date to show reciprocal monophyly of the two subspecies. CONCLUSIONS: Despite genetic bottlenecks and genetic drift, all surveyed zebra finch populations have maintained high MHC Class I diversity. The diversity at the MHC Class I locus in the Lesser Sundan subspecies contrasts sharply with the lack of diversity in previously examined neutral loci, and may thus be a result of selection acting to maintain polymorphism. Given uncertainty in historical population demography, however, it is difficult to rule out neutral processes in maintaining the observed diversity. The surveyed populations also differ in MHC Class I allele frequencies, and future studies are needed to assess whether these changes result in functional immune differences.