Polyglucosan bodies are accumulations of insoluble glucose polymers and proteins that form intracytoplasmic inclusions in the brain, large numbers of which can be indicative of neurodegenerative diseases such as Lafora disease. Montserrat orioles ( Icterus oberi) are an icterid passerine endemic to Montserrat with conservation populations maintained in captivity abroad. We demonstrate that polyglucosan bodies are unusually abundant in the cerebellar molecular and Purkinje cell layers and cerebellar peduncles of captive-bred and wild-caught Montserrat orioles. The bodies are periodic acid-Schiff positive and diastase resistant and label with concanavalin A and for ubiquitin, consistent with those seen in humans. We found no association of the polyglucosan bodies with concurrent neurological lesions or clinical signs, nor with EPM2A and EPM2B gene mutations associated with Lafora disease. We conclude that an abundance of cerebellar polyglucosan bodies may be a normal finding in aged Montserrat orioles and not a threat to the captive breeding population.

Avian influenza viruses (AIVs) are a major economic burden to the poultry industry and pose serious zoonotic risks, with human infections being reported every year. To date, the vaccination of birds remains the most important method for the prevention and control of AIV outbreaks. Most national vaccination strategies against AIV infection use whole virus-inactivated vaccines, which predominantly trigger a systemic antibody-mediated immune response. There are currently no studies that have examined the antibody repertoire of birds that were infected with and/or vaccinated against AIV. To this end, we evaluate the changes in the H9N2-specific IgM and IgY repertoires in chickens subjected to vaccination(s) and/or infectious challenge. We show that a large proportion of the IgM and IgY clones were shared across multiple individuals, and these public clonal responses are dependent on both the immunisation status of the birds and the specific tissue that was examined. Furthermore, the analysis revealed specific clonal expansions that are restricted to particular H9N2 immunisation regimes. These results indicate that both the nature and number of immunisations are important drivers of the antibody responses and repertoire profiles in chickens following H9N2 antigenic stimulation. We discuss how the repertoire biology of avian B-cell responses may affect the success of AIV vaccination in chickens, in particular the implications of public versus private clonal selection.

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

The pronounced growth in livestock populations since the 1950s has altered the epidemiological and evolutionary trajectory of their associated pathogens. For example, Marek's disease virus (MDV), which causes lymphoid tumors in chickens, has experienced a marked increase in virulence over the past century. Today, MDV infections kill >90% of unvaccinated birds, and controlling it costs more than US$1 billion annually. By sequencing MDV genomes derived from archeological chickens, we demonstrate that it has been circulating for at least 1000 years. We functionally tested the Meq oncogene, one of 49 viral genes positively selected in modern strains, demonstrating that ancient MDV was likely incapable of driving tumor formation. Our results demonstrate the power of ancient DNA approaches to trace the molecular basis of virulence in economically relevant pathogens.

BACKGROUND:The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome.METHODS:We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines.RESULTS:We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference.CONCLUSIONS:We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

Microbial colonisation is paramount to the normal development of the immune system, particularly at mucosal sites. However, the relationships between the microbiome and the adaptive immune repertoire have mostly been explored in rodents and humans. Here, we report a high-throughput sequencing analysis of the chicken TCRβ repertoire and the influences of microbial colonisation on tissue-resident TCRβ+ cells. The results reveal that the microbiome is an important driver of TCRβ diversity in both intestinal tissues and the bursa of Fabricius, but not in the spleen. Of note, public TCRβ sequences (shared across individuals) make a substantial contribution to the repertoire. Additionally, different tissues exhibit biases in terms of their V family and J gene usage, and these effects were influenced by the gut-associated microbiome. TCRβ clonal expansions were identified in both colonised and germ-free birds, but differences between the groups were indicative of an influence of the microbiota. Together, these findings provide an insight into the avian adaptive immune system and the influence of the microbiota on the TCRβ repertoire.

Penguins lost the ability to fly more than 60 million years ago, subsequently evolving a hyper-specialized marine body plan. Within the framework of a genome-scale, fossil-inclusive phylogeny, we identify key geological events that shaped penguin diversification and genomic signatures consistent with widespread refugia/recolonization during major climate oscillations. We further identify a suite of genes potentially underpinning adaptations related to thermoregulation, oxygenation, diving, vision, diet, immunity and body size, which might have facilitated their remarkable secondary transition to an aquatic ecology. Our analyses indicate that penguins and their sister group (Procellariiformes) have the lowest evolutionary rates yet detected in birds. Together, these findings help improve our understanding of how penguins have transitioned to the marine environment, successfully colonizing some of the most extreme environments on Earth.

The ability of zinc finger antiviral protein (ZAP) to recognize and respond to RNA virus sequences with elevated frequencies of CpG dinucleotides has been proposed as a functional part of the vertebrate innate immune antiviral response. It has been further proposed that ZAP activity shapes compositions of cytoplasmic mRNA sequences to avoid self-recognition, particularly mRNAs for interferons (IFNs) and IFN-stimulated genes (ISGs) expressed during the antiviral state. We investigated whether restriction of the replication of mutants of influenza A virus (IAV) and the echovirus 7 (E7) replicon with high CpG and UpA frequencies varied in different species of mammals and birds. Cell lines from different bird orders showed substantial variability in restriction of CpG-high mutants of IAV and E7 replicons, whereas none restricted UpA-high mutants, in marked contrast to universal restriction of both mutants in mammalian cells. Dinucleotide representation in ISGs and IFN genes was compared with those of cellular transcriptomes to determine whether potential differences in inferred ZAP activity between species shaped dinucleotide compositions of highly expressed genes during the antiviral state. While mammalian type 1 IFN genes typically showed often profound suppression of CpG and UpA frequencies, there was no oversuppression of either in ISGs in any species, irrespective of their ability to restrict CpG- or UpA-high mutants. Similarly, genome sequences of mammalian and avian RNA viruses were compositionally equivalent, as were IAV strains recovered from ducks, chickens and humans. Overall, we found no evidence for host variability in inferred ZAP function shaping host or viral transcriptome compositions.

Despite increasing interest in γδ T cells and their non-classical behaviour, most studies focus on animals with low numbers of circulating γδ T cells, such as mice and humans. Arguably, γδ T cell functions might be more prominent in chickens where these cells form a higher proportion of the circulatory T cell compartment. The TCR repertoire defines different subsets of γδ T cells, and such analysis is facilitated by well-annotated TCR loci. γδ T cells are considered at the cusp of innate and adaptive immunity but most functions have been identified in γδ low species. A deeper understanding of TCR repertoire biology in γδ high and γδ low animals is critical for defining the evolution of the function of γδ T cells. Repertoire dynamics will reveal populations that can be classified as innate-like or adaptive-like as well as those that straddle this definition. Here, a recent discrepancy in the structure of the chicken TCR gamma locus is resolved, demonstrating that tandem duplication events have shaped the evolution of this locus. Importantly, repertoire sequencing revealed large differences in the usage of individual TRGV genes, a pattern conserved across multiple tissues, including thymus, spleen and the gut. A single TRGV gene, TRGV3.3, with a highly diverse private CDR3 repertoire dominated every tissue in all birds. TRGV usage patterns were partly explained by the TRGV-associated recombination signal sequences. Public CDR3 clonotypes represented varying proportions of the repertoire of TCRs utilising different TRGVs, with one TRGV dominated by super-public clones present in all birds. The application of repertoire analysis enabled functional annotation of the TCRG locus in a species with a high circulating γδ phenotype. This revealed variable usage of TCRGV genes across multiple tissues, a pattern quite different to that found in γδ low species (human and mouse). Defining the repertoire biology of avian γδ T cells will be key to understanding the evolution and functional diversity of these enigmatic lymphocytes in an animal that is numerically more reliant on them. Practically, this will reveal novel ways in which these cells can be exploited to improve health in medical and veterinary contexts.

Penguins (Sphenisciformes) are an iconic order of flightless, diving seabirds distributed across a large latitudinal range in the Southern Hemisphere. The extensive area over which penguins are endemic is likely to have fostered variation in pathogen pressure, which in turn will have imposed differential selective pressures on the penguin immune system. At the front line of pathogen detection and response, the Toll-like receptors (TLRs) provide insight into host evolution in the face of microbial challenge. TLRs respond to conserved pathogen-associated molecular patterns and are frequently found to be under positive selection, despite retaining specificity for defined agonist classes. We undertook a comparative immunogenetics analysis of TLRs for all penguin species and found evidence of adaptive evolution that was largely restricted to the cell surface-expressed TLRs, with evidence of positive selection at, or near, key agonist-binding sites in TLR1B, TLR4, and TLR5. Intriguingly, TLR15, which is activated by fungal products, appeared to have been pseudogenized multiple times in the Eudyptes spp., but a full-length form was present as a rare haplotype at the population level. However, in vitro analysis revealed that even the full-length form of Eudyptes TLR15 was nonfunctional, indicating an ancestral cryptic pseudogenization prior to its eventual disruption multiple times in the Eudyptes lineage. This unusual pseudogenization event could provide an insight into immune adaptation to fungal pathogens such as Aspergillus, which is responsible for significant mortality in wild and captive bird populations.

ABSTRACT The ability of zinc finger antiviral protein (ZAP) to recognise and respond to RNA virus sequences with elevated frequencies of CpG dinucleotides has been proposed as a functional part of the vertebrate innate immune antiviral response. It has been further proposed that ZAP activity shapes compositions of cytoplasmic mRNA sequences to avoid self-recognition, particularly mRNAs for interferons (IFNs) and IFN-stimulated genes highly expressed when ZAP is upregulated during the antiviral state. We investigated the ZAP functional activity in different species of mammals and birds, and potential downstream effects of differences in CpG and UpA dinucleotide representations in host transcriptomes and in RNA viruses that infect them. Cell lines from different bird orders showed variability in restriction of influenza A virus and echovirus 7 replicons with elevated CpG frequencies and none restricted UpA-high mutants, in marked contrast to mammalian cell lines. Given this variability, we compared CpG and UpA representation in coding regions of ISGs and IFNs with the total cellular transcriptome to determine whether differences in ZAP activity shaped dinucleotide compositions of highly expressed genes during the antiviral state. While type 1 IFN genes typically showed often profound suppression of CpG and UpA frequencies, there was no over-suppression of CpGs or UpAs in ISGs in any species, irrespective of underlying ZAP activity. Similarly, mammalian and avian RNA virus genome sequences were compositionally equivalent as were IAV serotypes recovered from ducks, chickens and humans. Overall, we found no evidence for host variability in ZAP function impacting compositions of antiviral genes.

BACKGROUND:Penguins (Sphenisciformes) are a remarkable order of flightless wing-propelled diving seabirds distributed widely across the southern hemisphere. They share a volant common ancestor with Procellariiformes close to the Cretaceous-Paleogene boundary (66 million years ago) and subsequently lost the ability to fly but enhanced their diving capabilities. With ∼20 species among 6 genera, penguins range from the tropical Galápagos Islands to the oceanic temperate forests of New Zealand, the rocky coastlines of the sub-Antarctic islands, and the sea ice around Antarctica. To inhabit such diverse and extreme environments, penguins evolved many physiological and morphological adaptations. However, they are also highly sensitive to climate change. Therefore, penguins provide an exciting target system for understanding the evolutionary processes of speciation, adaptation, and demography. Genomic data are an emerging resource for addressing questions about such processes. RESULTS:Here we present a novel dataset of 19 high-coverage genomes that, together with 2 previously published genomes, encompass all extant penguin species. We also present a well-supported phylogeny to clarify the relationships among penguins. In contrast to recent studies, our results demonstrate that the genus Aptenodytes is basal and sister to all other extant penguin genera, providing intriguing new insights into the adaptation of penguins to Antarctica. As such, our dataset provides a novel resource for understanding the evolutionary history of penguins as a clade, as well as the fine-scale relationships of individual penguin lineages. Against this background, we introduce a major consortium of international scientists dedicated to studying these genomes. Moreover, we highlight emerging issues regarding ensuring legal and respectful indigenous consultation, particularly for genomic data originating from New Zealand Taonga species. CONCLUSIONS:We believe that our dataset and project will be important for understanding evolution, increasing cultural heritage and guiding the conservation of this iconic southern hemisphere species assemblage.

Over evolutionary time, pathogen challenge shapes the immune phenotype of the host to better respond to an incipient threat. The extent and direction of this selection pressure depend on the local pathogen composition, which is in turn determined by biotic and abiotic features of the environment. However, little is known about adaptation to local pathogen threats in wild animals. The Gentoo penguin (Pygoscelis papua) is a species complex that lends itself to the study of immune adaptation because of its circumpolar distribution over a large latitudinal range, with little or no admixture between different clades. In this study, we examine the diversity in a key family of innate immune genes—the Toll-like receptors (TLRs)—across the range of the Gentoo penguin. The three TLRs that we investigated present varying levels of diversity, with TLR4 and TLR5 greatly exceeding the diversity of TLR7. We present evidence of positive selection in TLR4 and TLR5, which points to pathogen-driven adaptation to the local pathogen milieu. Finally, we demonstrate that two positively selected cosegregating sites in TLR5 are sufficient to alter the responsiveness of the receptor to its bacterial ligand, flagellin. Taken together, these results suggest that Gentoo penguins have experienced distinct pathogen-driven selection pressures in different environments, which may be important given the role of the Gentoo penguin as a sentinel species in some of the world’s most rapidly changing environments.

Circoviruses infect a variety of animal species and have small (~1.8–2.2 kb) circular single-stranded DNA genomes. Recently a penguin circovirus (PenCV) was identified associated with an Adélie Penguin (Pygoscelis adeliae) with feather disorder and in the cloacal swabs of three asymptomatic Adélie Penguins at Cape Crozier, Antarctica. A total of 75 cloacal swab samples obtained from adults and chicks of three species of penguin (genus: Pygoscelis) from seven Antarctic breeding colonies (South Shetland Islands and Western Antarctic Peninsula) in the 2015−2016 breeding season were screened for PenCV. We identified new variants of PenCV in one Adélie Penguin and one Chinstrap Penguin (Pygoscelis antarcticus) from Port Charcot, Booth Island, Western Antarctic Peninsula, a site home to all three species of Pygoscelid penguins. These two PenCV genomes (length of 1986 nucleotides) share > 99% genome-wide nucleotide identity with each other and share ~87% genome-wide nucleotide identity with the PenCV sequences described from Adélie Penguins at Cape Crozier ~4400 km away in East Antarctica. We did not find any evidence of recombination among PenCV sequences. This is the first report of PenCV in Chinstrap Penguins and the first detection outside of Ross Island, East Antarctica. Given the limited knowledge on Antarctic animal viral diversity, future samples from Antarctic wildlife should be screened for these and other viruses to determine the prevalence and potential impact of viral infections.

The emergence of islands has been linked to spectacular radiations of diverse organisms. Although penguins spend much of their lives at sea, they rely on land for nesting, and a high proportion of extant species are endemic to geologically young islands. Islands may thus have been crucial to the evolutionary diversification of penguins. We test this hypothesis using a fossil-calibrated phylogeny of mitochondrial genomes (mitogenomes) from all extant and recently extinct penguin taxa. Our temporal analysis demonstrates that numerous recent island-endemic penguin taxa diverged following the formation of their islands during the Plio-Pleistocene, including the Galápagos (Galápagos Islands), northern rockhopper (Gough Island), erect-crested (Antipodes Islands), Snares crested (Snares) and royal (Macquarie Island) penguins. Our analysis also reveals two new recently extinct island-endemic penguin taxa from New Zealand's Chatham Islands: Eudyptes warhami sp. nov. and a dwarf subspecies of the yellow-eyed penguin, Megadyptes antipodes richdalei ssp. nov. Eudyptes warhami diverged from the Antipodes Islands erect-crested penguin between 1.1 and 2.5 Ma, shortly after the emergence of the Chatham Islands (∼3 Ma). This new finding of recently evolved taxa on this young archipelago provides further evidence that the radiation of penguins over the last 5 Ma has been linked to island emergence. Mitogenomic analyses of all penguin species, and the discovery of two new extinct penguin taxa, highlight the importance of island formation in the diversification of penguins, as well as the extent to which anthropogenic extinctions have affected island-endemic taxa across the Southern Hemisphere's isolated archipelagos.

Climate shifts are key drivers of ecosystem change. Despite the critical importance of Antarctica and the Southern Ocean for global climate, the extent of climate-driven ecological change in this region remains controversial. In particular, the biological effects of changing sea ice conditions are poorly understood. We hypothesize that rapid postglacial reductions in sea ice drove biological shifts across multiple widespread Southern Ocean species. We test for demographic shifts driven by climate events over recent millennia by analyzing population genomic datasets spanning 3 penguin genera (Eudyptes, Pygoscelis, and Aptenodytes). Demographic analyses for multiple species (macaroni/royal, eastern rockhopper, Adélie, gentoo, king, and emperor) currently inhabiting southern coastlines affected by heavy sea ice conditions during the Last Glacial Maximum (LGM) yielded genetic signatures of near-simultaneous population expansions associated with postglacial warming. Populations of the ice-adapted emperor penguin are inferred to have expanded slightly earlier than those of species requiring ice-free terrain. These concerted high-latitude expansion events contrast with relatively stable or declining demographic histories inferred for 4 penguin species (northern rockhopper, western rockhopper, Fiordland crested, and Snares crested) that apparently persisted throughout the LGM in ice-free habitats. Limited genetic structure detected in all ice-affected species across the vast Southern Ocean may reflect both rapid postglacial colonization of subantarctic and Antarctic shores, in addition to recent genetic exchange among populations. Together, these analyses highlight dramatic, ecosystem-wide responses to past Southern Ocean climate change and suggest potential for further shifts as warming continues.

Understanding the epidemiological dynamics of highly pathogenic avian influenza virus (HPAIV) in wild birds is crucial for guiding effective surveillance and control measures. The spread of H5 HPAIV has been well characterized over large geographical and temporal scales. However, information about the detailed dynamics and demographics of individual outbreaks in wild birds is rare and important epidemiological parameters remain unknown. We present data from a wild population of long-lived birds (mute swans; Cygnus olor) that has experienced three outbreaks of related H5 HPAIVs in the past decade, specifically, H5N1 (2007), H5N8 (2016) and H5N6 (2017). Detailed demographic data were available and intense sampling was conducted before and after the outbreaks; hence the population is unusually suitable for exploring the natural epidemiology, evolution and ecology of HPAIV in wild birds. We show that key epidemiological features remain remarkably consistent across multiple outbreaks, including the timing of virus incursion and outbreak duration, and the presence of a strong age-structure in morbidity that likely arises from an equivalent age-structure in immunological responses. The predictability of these features across a series of outbreaks in a complex natural population is striking and contributes to our understanding of HPAIV in wild birds. This article is part of the theme issue ‘Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes’. This issue is linked with the subsequent theme issue ‘Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control’.

a Department of Zoology, University of Otago, PO Box 56, Dunedin 9054, New Zealand. 8 b Manaaki Whenua Landcare Research, PO Box 69040, Lincoln, Canterbury 7640, New 9 Zealand. 10 c Bruce Museum, Greenwich, Connecticut, 06830, USA. 11 d Australian Centre for Ancient DNA, School of Biological Sciences, University of Adelaide, 12 Adelaide, SA, 5005, Australia. 13 e Museum of New Zealand Te Papa Tongarewa, PO Box 467, Wellington 6011, New 14 Zealand. 15 f Institute of Natural and Mathematical Sciences, Massey University, Auckland, New 16 Zealand. 17 g State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, 18 Chinese Academy of Sciences, 650223, Kunming, China. 19 h China National Genebank, BGI-Shenzhen, 518083, Shenzhen, Guangdong, China. 20 i Centre for Social Evolution, Department of Biology, Universitetsparken 15, University of 21 Copenhagen, DK-2100, Copenhagen, Denmark. 22 j ARAID Foundation. IUCA-Grupo Aragosaurus. Universidad de Zaragoza. C/ Pedro 23 Cerbuna, 12. 50009 Zaragoza, Spain. 24 k Bowling Green State University, Department of Biological Sciences, Bowling Green, Ohio, 25 43402, USA. 26 l Department of Zoology, University of Oxford, South Parks Road, Oxford, OX13PS, UK. 27 m Division of Pathology and Laboratory Medicine, University of Western Australia, Crawley, 28 Western Australia, 6009, Australia. 29 n Institute for Marine and Antarctic Studies, University of Tasmania, Hobart, Tasmania, 7001, 30 Australia. 31 FitzPatrick Institute of African Ornithology, DST-NRF Centre of Excellence, University of 32 Cape Town, Rondebosch 7701, South Africa. 33 p School of Environment, University of Auckland, Auckland, New Zealand. 34 35