Search

The animal gut microbiome can have a strong influence on the health, fitness, and behavior of its hosts. The composition of the gut microbial community can be influenced by factors such as diet, environment, and evolutionary history (phylosymbiosis). However, the relative influence of these factors is unknown in most bird species. Furthermore, phylosymbiosis studies have largely focused on clades that diverged tens of millions of years ago, and little is known about the degree of gut microbiome divergence in more recent species radiations. This study explores the drivers of microbiome variation across the unique and recent Hawaiian honeycreeper radiation (Fringillidae: Drepanidinae). Fecal samples were collected from 14 extant species spanning the main islands of the Hawaiian archipelago and were sequenced using three metabarcoding markers to characterize the gut microbiome, invertebrate diet, and plant diet of Hawaiian honeycreepers. We then used these metabarcoding data and the honeycreeper host phylogeny to evaluate their relative roles in shaping the gut microbiome. Microbiome variation across birds was highly individualized; however, source island had a small but significant effect on microbiome structure. The microbiomes did not recapitulate the host phylogenetic tree, indicating that evolutionary history does not strongly influence microbiome structure in the honeycreeper clade. These results expand our understanding of the roles of diet, geography, and phylogeny on avian microbiome structure, while also providing important ecological information about the diet and gut microbiota of wild Hawaiian honeycreepers., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.)

Island radiations, such as those of the Australo-Pacific, offer unique insight into diversification, extinction, and early speciation processes. Yet, their speciation and colonization histories are often obscured by conflicting genomic signals from incomplete lineage sorting or hybridization. Here, we integrated mitogenomes and genome-wide SNPs to unravel the evolutionary history of one of the world's most geographically widespread island radiations. The Australo-Pacific reed warblers (Acrocephalus luscinius complex) are a speciose lineage including five species that have become extinct since the 19th century and ten additional species of conservation concern. The radiation spans over 10,000 km across Australo-Papua, Micronesia and Polynesia, including the Mariana, Hawaii and Pitcairn Island archipelagos. Earlier mtDNA studies suggested a stepping-stone colonization process, resulting in archipelago-level secondary sympatry of divergent mtDNA lineages in the Mariana Islands and Marquesas. These studies hypothesised that morphologically similar species on neighbouring islands arose from ecological convergence. Using hDNA from historical museum specimens and modern genetic samples, we show that incomplete lineage sorting and/or gene flow have shaped the radiation of Australo-Pacific reed warblers rather than secondary sympatry. The nuclear genome reconstructs a simpler biogeographic history than mtDNA, showing close relationships between species in the Mariana Islands and Marquesas despite their paraphyletic mtDNA lineages. Gene flow likely involved early and late colonizing waves of the radiation before the loss of ancestral dispersive ability. Our results highlight how collection genomics can elucidate evolutionary history and inform conservation efforts for threatened species., (Published by Oxford University Press for The Society for the Study of Evolution (SSE) 2024. This work is written by (a) US Government employee(s) and is in the public domain in the US.)

Museum genomics provide an opportunity to investigate population demographics of extinct species, especially valuable when research prior to extinction was minimal. The Bachman's warbler (Vermivora bachmanii) is hypothesized to have gone extinct due to loss of its specialized habitat. However, little is known about other potential contributing factors such as natural rarity or changes to connectivity following habitat fragmentation. We examined mitochondrial DNA (mtDNA) and genome-wide SNPs using specimens collected from breeding and migration sites across the range of the Bachman's warbler. We found no signals of strong population structuring across the breeding range of Bachman's warblers in both mtDNA and genome-wide SNPs. Thus, long-term population isolation did not appear to be a significant contributor to the extinction of the Bachman's warbler. Instead, our findings support the theory that Bachman's warblers underwent a rapid decline likely driven by habitat destruction, which may have been exacerbated by the natural rarity, habitat specificity and low genetic diversity of the species., (© 2024. The Author(s).)

Genetic and genomic data are increasingly used to aid conservation management of endangered species by providing insights into evolutionary histories, factors associated with extinction risks, and potential for future adaptation. For the 'Alalā, or Hawaiian crow (Corvus hawaiiensis), genetic concerns include negative correlations between inbreeding and hatching success. However, it is unclear if low genetic diversity and inbreeding depression are consequences of a historical population bottleneck, or if 'Alalā had historically low genetic diversity that predated human influence, perhaps as a result of earlier declines or founding events. In this study, we applied a hybridization-based sequence capture to generate a genome-wide single nucleotide polymorphism (SNP) dataset for comparing historical specimens collected in the 1890s, when 'Alalā were more numerous, to samples taken between 1973 and 1998, when 'Alalā population densities were near the lowest documented levels in the wild, prior to all individuals being collected for captive rearing. We found low genome-wide diversity in both sample groups, however, the modern sample group (1973 to 1998 cohort) exhibited relatively fewer polymorphic alleles, a lower proportion of polymorphic loci, and lower observed heterozygosity, consistent with a population decline and potential bottleneck effects. These results combined with a current low population size highlight the importance of continued efforts by conservation managers to mitigate inbreeding and maintain founder representation to preserve what genetic diversity remains., (© The Author(s) 2023. Published by Oxford University Press on behalf of The American Genetic Association. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Plasmodium parasites infect thousands of species and provide an exceptional system for studying host-pathogen dynamics, especially for multi-host pathogens. However, understanding these interactions requires an accurate assay of infection. Assessing Plasmodium infections using microscopy on blood smears often misses infections with low parasitemias (the fractions of cells infected), and biases in malaria prevalence estimates will differ among hosts that differ in mean parasitemias. We examined Plasmodium relictum infection and parasitemia using both microscopy of blood smears and quantitative polymerase chain reaction (qPCR) on 299 samples from multiple bird species in Hawai'i and fit models to predict parasitemias from qPCR cycle threshold (Ct) values. We used these models to quantify the extent to which microscopy underestimated infection prevalence and to more accurately estimate infection patterns for each species for a large historical study done by microscopy. We found that most qPCR-positive wild-caught birds in Hawaii had low parasitemias (Ct scores ≥35), which were rarely detected by microscopy. The fraction of infections missed by microscopy differed substantially among eight species due to differences in species' parasitemia levels. Infection prevalence was likely 4-5-fold higher than previous microscopy estimates for three introduced species, including Zosterops japonicus, Hawaii's most abundant forest bird, which had low average parasitemias. In contrast, prevalence was likely only 1.5-2.3-fold higher than previous estimates for Himatione sanguinea and Chlorodrepanis virens, two native species with high average parasitemias. Our results indicate that relative patterns of infection among species differ substantially from those observed in previous microscopy studies, and that differences depend on variation in parasitemias among species. Although microscopy of blood smears is useful for estimating the frequency of different Plasmodium stages and host attributes, more sensitive quantitative methods, including qPCR, are needed to accurately estimate and compare infection prevalence among host species., (Copyright © 2023 Australian Society for Parasitology. All rights reserved.)

Of the estimated 55 Hawaiian honeycreepers (subfamily Carduelinae) only 17 species remain, nine of which the International Union for Conservation of Nature considers endangered. Among the most pressing threats to honeycreeper survival is avian malaria, caused by the introduced blood parasite Plasmodium relictum, which is increasing in distribution in Hawai'i as a result of climate change. Preventing further honeycreeper decline will require innovative conservation strategies that confront malaria from multiple angles. Research on mammals has revealed strong connections between gut microbiome composition and malaria susceptibility, illuminating a potential novel approach to malaria control through the manipulation of gut microbiota. One honeycreeper species, Hawai'i 'amakihi (Chlorodrepanis virens), persists in areas of high malaria prevalence, indicating they have acquired some level of immunity. To investigate if avian host-specific microbes may be associated with malaria survival, we characterized cloacal microbiomes and malaria infection for 174 'amakihi and 172 malaria-resistant warbling white-eyes (Zosterops japonicus) from Hawai'i Island using 16S rRNA gene metabarcoding and quantitative polymerase chain reaction. Neither microbial alpha nor beta diversity covaried with infection, but 149 microbes showed positive associations with malaria survivors. Among these were Escherichia and Lactobacillus spp., which appear to mitigate malaria severity in mammalian hosts, revealing promising candidates for future probiotic research for augmenting malaria immunity in sensitive endangered species., (© 2022 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

The unprecedented rise in the number of new and emerging infectious diseases in the last quarter century poses direct threats to human and wildlife health. The introduction to the Hawaiian archipelago of Plasmodium relictum and the mosquito vector that transmits the parasite has led to dramatic losses in endemic Hawaiian forest bird species. Understanding how mechanisms of disease immunity to avian malaria may evolve is critical as climate change facilitates increased disease transmission to high elevation habitats where malaria transmission has historically been low and the majority of the remaining extant Hawaiian forest bird species now reside. Here, we compare the transcriptomic profiles of highly susceptible Hawai'i 'amakihi (Chlorodrepanis virens) experimentally infected with P. relictum to those of uninfected control birds from a naïve high elevation population. We examined changes in gene expression profiles at different stages of infection to provide an in-depth characterization of the molecular pathways contributing to survival or mortality in these birds. We show that the timing and magnitude of the innate and adaptive immune response differed substantially between individuals that survived and those that succumbed to infection, and likely contributed to the observed variation in survival. These results lay the foundation for developing gene-based conservation strategies for Hawaiian honeycreepers by identifying candidate genes and cellular pathways involved in the pathogen response that correlate with a bird's ability to recover from malaria infection., (© The Author(s) 2023. Published by Oxford University Press on behalf of The American Genetic Association. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Dynamic interactions between host, pathogen and host-associated microbiome dictate infection outcomes. Pathogens including Batrachochytrium dendrobatidis (Bd) threaten global biodiversity, but conservation efforts are hindered by limited understanding of amphibian host, Bd and microbiome interactions. We conducted a vaccination and infection experiment using Eastern hellbender salamanders ( Cryptobranchus alleganiensis alleganiensis ) challenged with Bd to observe infection, skin microbial communities and gene expression of host skin, pathogen and microbiome throughout the experiment. Most animals survived high Bd loads regardless of their vaccination status and vaccination did not affect pathogen load, but host gene expression differed based on vaccination. Oral vaccination (exposure to killed Bd) stimulated immune gene upregulation while topically and sham-vaccinated animals did not significantly upregulate immune genes. In early infection, topically vaccinated animals upregulated immune genes but orally and sham-vaccinated animals downregulated immune genes. Bd increased pathogenicity-associated gene expression in late infection when Bd loads were highest. The microbiome was altered by Bd, but there was no correlation between anti-Bd microbe abundance or richness and pathogen burden. Our observations suggest that hellbenders initially generate a vigorous immune response to Bd, which is ineffective at controlling disease and is subsequently modulated. Interactions with antifungal skin microbiota did not influence disease progression., (© 2023 The Authors.)

Connectivity among wildlife populations facilitates exchange of genetic material between groups. Changes to historical connectivity patterns resulting from anthropogenic activities can therefore have negative consequences for genetic diversity, particularly for small or isolated populations. DNA obtained from museum specimens can enable direct comparison of temporal changes in connectivity among populations, which can aid in conservation planning and contribute to the understanding of population declines. However, museum DNA can be degraded and only available in low quantities, rendering it challenging for use in population genomic analyses. Applications of genomic methodologies such as targeted sequencing address this issue by enabling capture of shared variable sites, increasing quantity and quality of recovered genomic information. We used targeted sequencing of ultra-conserved Elements (UCEs) to evaluate potential changes in connectivity and genetic diversity of roseate terns (Sterna dougallii) with a breeding distribution in the northwestern Atlantic and the Caribbean. Both populations experienced range contractions and population declines due to anthropogenic activity in the 20th century, which has the potential to alter historical connectivity regimes. Instead, we found that the two populations were differentiated historically as well as contemporaneously, with little evidence of migration between them for either time period. We also found no evidence for temporal changes in genetic diversity, although these interpretations may have been limited due to sequencing artifacts caused by the degraded nature of the museum samples. Population structuring in migratory seabirds is typically reflective of low rates of divergence and high connectivity among geographically segregated subpopulations. Our contrasting results suggest the potential presence of ecological mechanisms driving population differentiation, and highlight the value of targeted sequencing on DNA derived from museum specimens to uncover long-term patterns of genetic differentiation in wildlife populations., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Society for Integrative and Comparative Biology.)

Robert Fleischer and colleagues introduce the unique songbird fauna of Hawaii and the threats it faces., (Published by Elsevier Inc.)

Background: Plasmodium parasites that cause bird malaria occur in all continents except Antarctica and are primarily transmitted by mosquitoes in the genus Culex. Culex quinquefasciatus, the mosquito vector of avian malaria in Hawai'i, became established in the islands in the 1820s. While the deadly effects of malaria on endemic bird species have been documented for many decades, vector-parasite interactions in avian malaria systems are relatively understudied., Methods: To evaluate the gene expression response of mosquitoes exposed to a Plasmodium infection intensity known to occur naturally in Hawai'i, offspring of wild-collected Hawaiian Cx. quinquefasciatus were fed on a domestic canary infected with a fresh isolate of Plasmodium relictum GRW4 from a wild-caught Hawaiian honeycreeper. Control mosquitoes were fed on an uninfected canary. Transcriptomes of five infected and three uninfected individual mosquitoes were sequenced at each of three stages of the parasite life cycle: 24 h post feeding (hpf) during ookinete invasion; 5 days post feeding (dpf) when oocysts are developing; 10 dpf when sporozoites are released and invade the salivary glands., Results: Differential gene expression analyses showed that during ookinete invasion (24 hpf), genes related to oxidoreductase activity and galactose catabolism had lower expression levels in infected mosquitoes compared to controls. Oocyst development (5 dpf) was associated with reduced expression of a gene with a predicted innate immune function. At 10 dpf, infected mosquitoes had reduced expression levels of a serine protease inhibitor, and further studies should assess its role as a Plasmodium agonist in C. quinquefasciatus. Overall, the differential gene expression response of Hawaiian Culex exposed to a Plasmodium infection intensity known to occur naturally in Hawai'i was low, but more pronounced during ookinete invasion., Conclusions: This is the first analysis of the transcriptional responses of vectors to malaria parasites in non-mammalian systems. Interestingly, few similarities were found between the response of Culex infected with a bird Plasmodium and those reported in Anopheles infected with human Plasmodium. The relatively small transcriptional changes observed in mosquito genes related to immune response and nutrient metabolism support conclusions of low fitness costs often documented in experimental challenges of Culex with avian Plasmodium., (© 2022. The Author(s).)

Understanding predator population dynamics is important for conservation management because of the critical roles predators play within ecosystems. Noninvasive genetic sampling methods are useful for the study of predators like canids that can be difficult to capture or directly observe. Here, we introduce the FAECES* method (Fast and Accurate Enrichment of Canid Excrement for Species* and other analyses) which expands the toolbox for canid researchers and conservationists by using in-solution hybridization sequence capture to produce single nucleotide polymorphism (SNP) genotypes for multiple canid species from scat-derived DNA using a single enrichment. We designed a set of hybridization probes to genotype both coyotes (Canis latrans) and kit foxes (Vulpes macrotis) at hundreds of polymorphic SNP loci and we tested the probes on both tissues and field-collected scat samples. We enriched and genotyped by sequencing 52 coyote and 70 kit fox scats collected in and around a conservation easement in the Nevada Mojave Desert. We demonstrate that the FAECES* method produces genotypes capable of differentiating coyotes and kit foxes, identifying individuals and their sex, and estimating genetic diversity and effective population sizes, even using highly degraded, low-quantity DNA extracted from scat. We found that the study area harbours a large and diverse population of kit foxes and a relatively smaller population of coyotes. By replicating our methods in the future, conservationists can assess the impacts of management decisions on canid populations. The method can also be adapted and applied more broadly to enrich and sequence multiple loci from any species of interest using scat or other noninvasive genetic samples., (© 2021 John Wiley & Sons Ltd.)

Conservation benefits from incorporating genomics to explore the impacts of population declines, inbreeding, loss of genetic variation and hybridization. Here we use the near-extinct Mariana Islands reedwarbler radiation to showcase how ancient DNA approaches can allow insights into the population dynamics of extinct species and threatened populations for which historical museum specimens or material with low DNA yield (e.g., scats, feathers) are the only sources for DNA. Despite their having paraphyletic mitochondrial DNA (mtDNA), nuclear single nucleotide polymorphisms (SNPs) support the distinctiveness of critically endangered Acrocephalus hiwae and the other three species in the radiation that went extinct between the 1960s and 1990s. Two extinct species, A. yamashinae and A. luscinius, were deeply divergent from each other and from a third less differentiated lineage containing A. hiwae and extinct A. nijoi. Both mtDNA and SNPs suggest that the two isolated populations of A. hiwae from Saipan and Alamagan Islands are sufficiently distinct to warrant subspecies recognition and separate conservation management. We detected no significant differences in genetic diversity or inbreeding between Saipan and Alamagan, nor strong signatures of geographical structuring within either island. However, the implications of possible signatures of inbreeding in both Saipan and Alamagan, and long-term population declines in A. hiwae that pre-date modern anthropogenic threats require further study with denser population sampling. Our study highlights the value that conservation genomics studies of island radiations have as windows onto the possible future for the world's biota as climate change and habitat destruction increasingly fragment their ranges and contribute to rapid declines in population abundances., (© 2022 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode., (Copyright © 2021 Elsevier Inc. All rights reserved.)

We designed two probiotic treatments to control chytridiomycosis caused by Batrachochytrium dendrobatidis (Bd) on infected Panamanian golden frogs (Atelopus zeteki), a species that is thought to be extinct in the wild due to Bd. The first approach disrupted the existing skin microbe community with antibiotics then exposed the frogs to a core golden frog skin microbe (Diaphorobacter sp.) that we genetically modified to produce high titers of violacein, a known antifungal compound. One day following probiotic treatment, the engineered Diaphorobacter and the violacein-producing pathway could be detected on the frogs but the treatment failed to improve frog survival when exposed to Bd. The second approach exposed frogs to the genetically modified bacterium mixed into a consortium with six other known anti-Bd bacteria isolated from captive A. zeteki, with no preliminary antibiotic treatment. The consortium treatment increased the frequency and abundance of three probiotic isolates (Janthinobacterium, Chryseobacterium, and Stenotrophomonas) and these persisted on the skin 4 weeks after probiotic treatment. There was a temporary increase in the frequency and abundance of three other probiotics isolates (Masillia, Serratia, and Pseudomonas) and the engineered Diaphorobacter isolate, but they subsequently disappeared from the skin. This treatment also failed to reduce frog mortality upon exposure., (© 2021. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

Introduced into Hawaii in the early 1900s, the Japanese white-eye or warbling white-eye (Zosterops japonicus) is now the most abundant land bird in the archipelago. Here, we present the first Z. japonicus genome, sequenced from an individual in its invasive range. This genome provides an important resource for future studies in invasion genomics. We annotated the genome using two workflows-standalone AUGUSTUS and BRAKER2. We found that AUGUSTUS was more conservative with gene predictions when compared with BRAKER2. The final number of annotated gene models was similar between the two workflows, but standalone AUGUSTUS had over 70% of gene predictions with Blast2GO annotations versus under 30% using BRAKER2. Additionally, we tested whether using RNA-seq data from 47 samples had a significant impact on annotation quality when compared with data from a single sample, as generating RNA-seq data for genome annotation can be expensive and requires well preserved tissue. We found that more data did not significantly change the number of annotated genes using AUGUSTUS but using BRAKER2 the number increased substantially. The results presented here will aid researchers in annotating draft genomes of nonmodel species as well as those studying invasion success., (Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution 2021. This work is written by US Government employees and is in the public domain in the US.)

The malaria parasite Plasmodium relictum (lineage GRW4) was introduced less than a century ago to the native avifauna of Hawai'i, where it has since caused major declines of endemic bird populations. One of the native bird species that is frequently infected with GRW4 is the Hawai'i 'amakihi ( Chlorodrepanis virens ). To achieve a better understanding of the transcriptional activities of this virulent parasite, we performed a controlled challenge experiment of 15 'amakihi that were infected with GRW4. Blood samples containing malaria parasites were collected at two time points (intermediate and peak infection stages) from host individuals that were either experimentally infected by mosquitoes or inoculated with infected blood. We then used RNA sequencing to assemble a high-quality blood transcriptome of P. relictum GRW4, allowing us to quantify parasite expression levels inside individual birds. We found few significant differences (one to two transcripts) in GRW4 expression levels between host infection stages and between inoculation methods. However, 36 transcripts showed differential expression levels among all host individuals, indicating a potential presence of host-specific gene regulation across hosts. To reduce the extinction risk of the remaining native bird species in Hawai'i, genetic resources of the local Plasmodium lineage are needed to enable further molecular characterization of this parasite. Our newly built Hawaiian GRW4 transcriptome assembly, together with analyses of the parasite's transcriptional activities inside the blood of Hawai'i 'amakihi, can provide us with important knowledge on how to combat this deadly avian disease in the future., Competing Interests: The authors declare that they have no conflicts of interest., (© 2021 Smithosnian Institution. Ecology and Evolution published by John Wiley & Sons Ltd. This article has been contributed to by US Government employees and their work is in the public domain in the USA.)

Very few animals habitually manufacture and use tools. It has been suggested that advanced tool behaviour co-evolves with a suite of behavioural, morphological and life history traits. In fact, there are indications for such an adaptive complex in tool-using crows (genus Corvus species). Here, we sequenced the genomes of two habitually tool-using and ten non-tool-using crow species to search for genomic signatures associated with a tool-using lifestyle. Using comparative genomic and population genetic approaches, we screened for signals of selection in protein-coding genes in the tool-using New Caledonian and Hawaiian crows. While we detected signals of recent selection in New Caledonian crows near genes associated with bill morphology, our data indicate that genetic changes in these two lineages are surprisingly subtle, with little evidence at present for convergence. We explore the biological explanations for these findings, such as the relative roles of gene regulation and protein-coding changes, as well as the possibility that statistical power to detect selection in recently diverged lineages may have been insufficient. Our study contributes to a growing body of literature aiming to decipher the genetic basis of recently evolved complex behaviour., (© 2020 John Wiley & Sons Ltd.)

Rapid advances in genomic tools for use in ecological contexts and non-model systems allow unprecedented insight into interactions that occur beyond direct observation. We developed an approach that couples microbial forensics with molecular dietary analysis to identify species interactions and scavenging by invasive rats on native and introduced birds in Hawaii. First, we characterized bacterial signatures of bird carcass decay by conducting 16S rRNA high-throughput sequencing on chicken ( Gallus gallus domesticus ) tissues collected over an 11-day decomposition study in natural Hawaiian habitats. Second, we determined if field-collected invasive black rats ( Rattus rattus ; n  = 51, stomach and fecal samples) had consumed birds using molecular diet analysis with two independent PCR assays (mitochondrial Cytochrome Oxidase I and Cytochrome b genes) and Sanger sequencing. Third, we characterized the gut microbiome of the same rats using 16S rRNA high-throughput sequencing and identified 15 bacterial taxa that were (a) detected only in rats that consumed birds ( n  = 20/51) and (b) were indicative of decaying tissue in the chicken decomposition experiment. We found that 18% of rats ( n  = 9/51) likely consumed birds as carrion by the presence of bacterial biomarkers of decayed tissue in their gut microbiome. One species of native bird ( Myadestes obscurus ) and three introduced bird species ( Lophura leucomelanos , Meleagris gallopavo , Zosterops japonicus ) were detected in the rats' diets, with individuals from these species (except L. nycthemera ) likely consumed through scavenging. Bacterial biomarkers of bird carcass decay can persist through rat digestion and may serve as biomarkers of scavenging. Our approach can be used to reveal trophic interactions that are challenging to measure through direct observation., Competing Interests: The authors have declared that no competing interests exist., (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity 1-4 . Sparse taxon sampling has previously been proposed to confound phylogenetic inference 5 , and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

The disease chytridiomycosis caused by the fungus Bd has devastated amphibian populations worldwide. Functional genomic contributions to host susceptibility remain enigmatic and vary between species and populations. We conducted experimental Bd infections in Rana yavapaiensis, a species with intraspecific variation in chytridiomycosis susceptibility, to assess the skin and spleen transcriptomic response to infection over time. We predicted that increased immune gene expression would be associated with a positive disease outcome, but we instead found that surviving frogs had significantly reduced immune gene expression compared to susceptible frogs and to uninfected controls. MHC class IIβ gene expression was also significantly higher in susceptible frogs compared to surviving frogs. Furthermore, susceptible frogs expressed a significantly larger number of distinct class IIβ alleles, demonstrating a negative correlation between class IIβ expression, functional diversity, and survival. Expression of the MHC class IIβ locus previously associated with Bd disease outcomes was a significant predictor of Bd infection intensity at early infection stages but not at late infection stages, suggesting initial MHC-linked immune processes are important for ultimate disease outcomes. We infer through disease association and phylogenetic analysis that certain MHC variants are linked to the immune expression that was negatively associated with survival, and we hypothesize that frogs that did not express these alleles could better survive infections. Our study finds that MHC expression at early and late infection stages predicts Bd infection intensity, and suggests that generating a sustained immune response against Bd may be counterproductive for surviving chytridiomycosis in this partially susceptible species., (© 2020 John Wiley & Sons Ltd.)

The introduction of nonnative species and reductions in native biodiversity have resulted in substantial changes in vector and host communities globally, but the consequences for pathogen transmission are poorly understood. In lowland Hawaii, bird communities are composed of primarily introduced species, with scattered populations of abundant native species. We examined the influence of avian host community composition, specifically the role of native and introduced species, as well as host diversity, on the prevalence of avian malaria (Plasmodium relictum) in the southern house mosquito (Culex quinquefasciatus). We also explored the reciprocal effect of malaria transmission on native host populations and demography. Avian malaria infection prevalence in mosquitoes increased with the density and relative abundance of native birds, as well as host community competence, but was uncorrelated with host diversity. Avian malaria transmission was estimated to reduce population growth rates of Hawai'i 'amakihi (Chlorodrepanis virens) by 7-14%, but mortality from malaria could not explain gaps in this species' distribution at our sites. Our results suggest that, in Hawaii, native host species increase pathogen transmission to mosquitoes, but introduced species can also support malaria transmission alone. The increase in pathogen transmission with native bird abundance leads to additional disease mortality in native birds, further increasing disease impacts in an ecological feedback cycle. In addition, vector abundance was higher at sites without native birds and this overwhelmed the effects of host community composition on transmission such that infected mosquito abundance was highest at sites without native birds. Higher disease risk at these sites due to higher vector abundance could inhibit recolonization and recovery of native species to these areas. More broadly, this work shows how differences in host competence for a pathogen among native and introduced taxa can influence transmission and highlights the need to examine this question in other systems to determine the generality of this result., (© 2020 by the Ecological Society of America.)

Coffea arabica is a highly traded commodity worldwide, and its plantations are habitat to a wide range of organisms. Coffee farmers are shifting away from traditional shade coffee farms in favor of sun-intensive, higher yield farms, which can impact local biodiversity. Using plant-associated microorganisms in biofertilizers, particularly fungi collected from local forests, to increase crop yields has gained traction among coffee producers. However, the taxonomic and spatial distribution of many fungi in coffee soil, nearby forests and biofertilizers is unknown. We collected soil samples from a sun coffee system, shade coffee system, and nearby forest from Izalco, Sonsonate, El Salvador. At each coffee system, we collected soil from the surface (upper) and 10 cm below the surface (lower), and from the coffee plant drip line (drip line) and the walkway between two plants (walkway). Forest soils were collected from the surface only. We used ITS metabarcoding to characterize fungal communities in soil and in the biofertilizer (applied in both coffee systems), and assigned fungal taxa to functional guilds using FUNGuild. In the sun and shade coffee systems, we found that drip line soil had higher richness in pathotrophs, symbiotrophs, and saprotrophs than walkway soil, suggesting that fungi select for microhabitats closer to coffee plants. Upper and lower soil depths did not differ in fungal richness or composition, which may reflect the shallow root system of Coffea arabica. Soil from shade, sun, and forest had similar numbers of fungal taxa, but differed dramatically in community composition, indicating that local habitat differences drive fungal species sorting among systems. Yet, some fungal taxa were shared among systems, including seven fungal taxa present in the biofertilizer. Understanding the distribution of coffee soil mycobiomes can be used to inform sustainable, ecologically friendly farming practices and identify candidate plant-growth promoting fungi for future studies., Competing Interests: The authors have declared that no competing interests exist.

Indochina and Sundaland are biologically diverse, interconnected regions of Southeast Asia with complex geographic histories. Few studies have examined phylogeography of bird species that span the two regions because of inadequate population sampling. To determine how geographic barriers/events and disparate dispersal potential have influenced the population structure, gene flow, and demographics of species that occupy the entire area, we studied five largely codistributed rainforest bird species: Arachnothera longirostra , Irena puella , Brachypodius atriceps , Niltava grandis , and Stachyris nigriceps . We accomplished relatively thorough sampling and data collection by sequencing ultraconserved elements (UCEs) using DNA extracted from modern and older (historical) specimens. We obtained a genome-wide set of 753-4,501 variable loci and 3,919-18,472 single nucleotide polymorphisms. The formation of major within-species lineages occurred within a similar span of time (0.5-1.5 mya). Major patterns in population genetic structure are largely consistent with the dispersal potential and habitat requirements of the study species. A population break across the Isthmus of Kra was shared only by the two hill/submontane insectivores ( N. grandis and S. nigriceps ). Across Sundaland, there is little structure in B. atriceps , which is a eurytopic and partially frugivorous species that often utilizes forest edges. Two other eurytopic species, A. longirostra and I. puella , possess highly divergent populations in peripheral Sunda Islands (Java and/or Palawan) and India. These species probably possess intermediate dispersal abilities that allowed them to colonize new areas, and then remained largely isolated subsequently. We also observed an east-west break in Indochina that was shared by B. atriceps and S. nigriceps , species with very different habitat requirements and dispersal potential. By analyzing high-throughput DNA data, our study provides an unprecedented comparative perspective on the process of avian population divergence across Southeast Asia, a process that is determined by geography, species characteristics, and the stochastic nature of dispersal and vicariance events., Competing Interests: None declared., (© 2020 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

The Hawai'ian honeycreepers (drepanids) are a classic example of adaptive radiation: they adapted to a variety of novel dietary niches, evolving a wide range of bill morphologies. Here we investigated genomic diversity, demographic history, and genes involved in bill morphology phenotypes in 2 honeycreepers: the 'akiapōlā'au (Hemignathus wilsoni) and the Hawai'i 'amakihi (Chlorodrepanis virens). The 'akiapōlā'au is an endangered island endemic, filling the "woodpecker" niche by using a unique bill morphology, while the Hawai'i 'amakihi is a dietary generalist common on the islands of Hawai'i and Maui. We de novo sequenced the 'akiapōlā'au genome and compared it to the previously sequenced 'amakihi genome. The 'akiapōlā'au is far less heterozygous and has a smaller effective population size than the 'amakihi, which matches expectations due to its smaller census population and restricted ecological niche. Our investigation revealed genomic islands of divergence, which may be involved in the honeycreeper radiation. Within these islands of divergence, we identified candidate genes (including DLK1, FOXB1, KIF6, MAML3, PHF20, RBP1, and TIMM17A) that may play a role in honeycreeper adaptations. The gene DLK1, previously shown to influence Darwin's finch bill size, may be related to honeycreeper bill morphology evolution, while the functions of the other candidates remain unknown., (Published by Oxford University Press on behalf of The American Genetic Association 2019.)

Adaptive radiation plays a fundamental role in our understanding of the evolutionary process. However, the concept has provoked strong and differing opinions concerning its definition and nature among researchers studying a wide diversity of systems. Here, we take a broad view of what constitutes an adaptive radiation, and seek to find commonalities among disparate examples, ranging from plants to invertebrate and vertebrate animals, and remote islands to lakes and continents, to better understand processes shared across adaptive radiations. We surveyed many groups to evaluate factors considered important in a large variety of species radiations. In each of these studies, ecological opportunity of some form is identified as a prerequisite for adaptive radiation. However, evolvability, which can be enhanced by hybridization between distantly related species, may play a role in seeding entire radiations. Within radiations, the processes that lead to speciation depend largely on (1) whether the primary drivers of ecological shifts are (a) external to the membership of the radiation itself (mostly divergent or disruptive ecological selection) or (b) due to competition within the radiation membership (interactions among members) subsequent to reproductive isolation in similar environments, and (2) the extent and timing of admixture. These differences translate into different patterns of species accumulation and subsequent patterns of diversity across an adaptive radiation. Adaptive radiations occur in an extraordinary diversity of different ways, and continue to provide rich data for a better understanding of the diversification of life., (© The American Genetic Association 2020. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens., (© 2019 John Wiley & Sons Ltd.)

Milk is inhabited by a community of bacteria and is one of the first postnatal sources of microbial exposure for mammalian young. Bacteria in breast milk may enhance immune development, improve intestinal health, and stimulate the gut-brain axis for infants. Variation in milk microbiome structure (e.g., operational taxonomic unit [OTU] diversity, community composition) may lead to different infant developmental outcomes. Milk microbiome structure may depend on evolutionary processes acting at the host species level and ecological processes occurring over lactation time, among others. We quantified milk microbiomes using 16S rRNA high-throughput sequencing for nine primate species and for six primate mothers sampled over lactation. Our data set included humans (Homo sapiens, Philippines and USA) and eight nonhuman primate species living in captivity (bonobo [Pan paniscus], chimpanzee [Pan troglodytes], western lowland gorilla [Gorilla gorilla gorilla], Bornean orangutan [Pongo pygmaeus], Sumatran orangutan [Pongo abelii], rhesus macaque [Macaca mulatta], owl monkey [Aotus nancymaae]) and in the wild (mantled howler monkey [Alouatta palliata]). For a subset of the data, we paired microbiome data with nutrient and hormone assay results to quantify the effect of milk chemistry on milk microbiomes. We detected a core primate milk microbiome of seven bacterial OTUs indicating a robust relationship between these bacteria and primate species. Milk microbiomes differed among primate species with rhesus macaques, humans and mantled howler monkeys having notably distinct milk microbiomes. Gross energy in milk from protein and fat explained some of the variations in microbiome composition among species. Microbiome composition changed in a predictable manner for three primate mothers over lactation time, suggesting that different bacterial communities may be selected for as the infant ages. Our results contribute to understanding ecological and evolutionary relationships between bacteria and primate hosts, which can have applied benefits for humans and endangered primates in our care., (© 2019 Wiley Periodicals, Inc.)

Reptiles show varying degrees of facultative parthenogenesis. Here we use genetic methods to determine that an isolated, captive female Asian water dragon produced at least nine offspring via parthenogenesis. We identified microsatellites for the species from shotgun genomic sequences, selected and optimized primer sets, and tested all of the offspring for a set of seven microsatellites that were heterozygous in the mother. We verified that the seven loci showed high levels of polymorphism in four wild Asian water dragons from Vietnam. In all cases, the offspring (unhatched, but developed eggs, or hatched young) had only a single allele at each locus, and contained only alleles present in the mother's genotype (i.e., were homozygous or hemizygous). The probability that our findings resulted from the female mating with one or more males is extremely small, indicating that the offspring were derived from a single female gamete (either alone or via duplication and/or fusion) and implicating parthenogenesis. This is the first documented case of parthenogenesis in the Squamate family Agamidae., Competing Interests: The authors have declared that no competing interests exist.

The first report of giant pandas (Ailuropoda melanoleuca) infected with a novel Hepatozoon species is presented. An intraleukocytic parasite was detected via routine blood smear from a zoo-housed giant panda at the National Zoological Park. Ribosomal DNA sequences indicated a previously undescribed Hepatozoon species. Phylogenetic and distance analyses of the sequences placed it within its own branch, clustered with Old World species with carnivore (primarily ursid and mustelid) hosts. Retrospective and opportunistic testing of other individuals produced additional positive detections (17/23, 73.9%), demonstrating 100% prevalence (14/14) across five institutions. All animals were asymptomatic at time of sampling, and health implications for giant pandas remain unknown.

Pathogens compete with host microbiomes for space and resources. Their shared environment impacts pathogen-microbiome-host interactions, which can lead to variation in disease outcome. The skin microbiome of red-backed salamanders (Plethodon cinereus) can reduce infection by the pathogen Batrachochytrium dendrobatidis (Bd) at moderate infection loads, with high species richness and high abundance of competitors as putative mechanisms. However, it is unclear if the skin microbiome can reduce epizootic Bd loads across temperatures. We conducted a laboratory experiment to quantify skin microbiome and host responses (P. cinereus: n = 87) to Bd at mimicked epizootic loads across temperatures (13, 17 and 21°C). We quantified skin microbiomes using 16S rRNA gene metabarcoding and identified operational taxonomic units (OTUs) taxonomically similar to culturable bacteria known to kill Bd (anti-Bd OTUs). Prior to pathogen exposure, temperature changed the microbiome (OTU richness decreased by 12% and the abundance of anti-Bd OTUs increased by 18% per degree increase in temperature), but these changes were not predictive of disease outcome. After exposure, Bd changed the microbiome (OTU richness decreased by 0.1% and the abundance of anti-Bd OTUs increased by 0.2% per 1% increase in Bd load) and caused high host mortality across temperatures (35/45: 78%). Temperature indirectly impacted microbiome change and mortality through its direct effect on pathogen load. We did not find support for the microbiome impacting Bd load or host survival. Our research reveals complex host, pathogen, microbiome and environmental interactions to demonstrate that during epizootic events the microbiome will be unlikely to reduce pathogen invasion, even for putatively Bd-resistant species., (© 2019 John Wiley & Sons Ltd.)

Species with narrow environmental tolerances are often distributed within fragmented patches of suitable habitat, and dispersal among these subpopulations can be difficult to directly observe. Genetic data can help quantify gene flow between localities, which is especially important for vulnerable species with a disjunct range. The Shenandoah salamander ( Plethodon shenandoah ) is a federally endangered species known only from three mountaintops in Virginia, USA. To reconstruct the evolutionary history and population connectivity of this species, we generated both mitochondrial and nuclear data using sequence capture from individuals collected across all three mountaintops. Applying population and landscape genetic methods, we found strong population structure that was independent of geographic distance. Both the nuclear markers and mitochondrial genomes indicated a deep split between the most southern population and the genetically similar central and northern populations. Although there was some mitochondrial haplotype-splitting between the central and northern populations, there was admixture in nuclear markers. This is indicative of either a recent split or current male-biased dispersal among mountain isolates. Models of landscape resistance found that dispersal across north-facing slopes at mid-elevation levels best explain the observed genetic structure among populations. These unexpected results highlight the importance of incorporating landscape features in understanding and predicting the movement and fragmentation of this range-restricted salamander species across space., Competing Interests: The authors have no conflicts of interest or competing interests to disclose.

Infectious diseases are causing catastrophic losses to global biodiversity. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian pathogen susceptibility. MHC genes encode two classes of polymorphic cell-surface molecules that can recognize and bind to diverse pathogen peptides. While MHC class I genes are the classic mediators of viral-acquired immunity, larval amphibians do not express them. Consequently, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Rana sylvatica), which are highly susceptible to Ranavirus, across 17 ponds and 2 years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIβ locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if supertypes or alleles influenced host responses to Ranavirus. Among 381 sampled individuals, 26% were infected with Ranavirus. We recovered 20 unique MHC class IIβ alleles that fell into two deeply diverged clades and seven supertypes. MHC genotypes were associated with Ranavirus infection intensity, but not prevalence. Specifically, MHC heterozygotes and supertype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and other supertypes. We conclude that MHC class IIβ functional genetic variation is an important component of Ranavirus susceptibility. Identifying immunogenetic signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

The Mariana Crow, or Åga ( Corvus kubaryi ), is a critically endangered species (IUCN -International Union for Conservation of Nature), endemic to the islands of Guam and Rota in the Mariana Archipelago. It is locally extinct on Guam, and numbers have declined dramatically on Rota to a historical low of less than 55 breeding pairs throughout the island in 2013. Because of its extirpation on Guam and population decline on Rota, it is of critical importance to assess the genetic variation among individuals to assist ongoing recovery efforts. We conducted a population genomics analysis comparing the Guam and Rota populations and studied the genetic structure of the Rota population. We used blood samples from five birds from Guam and 78 birds from Rota. We identified 145,552 candidate single nucleotide variants (SNVs) from a genome sequence of an individual from Rota and selected a subset of these to develop an oligonucleotide in-solution capture assay. The Guam and Rota populations were genetically differentiated from each other. Crow populations sampled broadly across their range on Rota showed significant genetic structuring ⁻ a surprising result given the small size of this island and the good flight capabilities of the species. Knowledge of its genetic structure will help improve management strategies to help with its recovery.

Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria, Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i 'amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations where P. relictum exists, and can sustain infection without major fitness consequences. High-elevation, unexposed populations of 'amakihi display little to no tolerance. To explore the genomic basis of adaptation to P. relictum in low-elevation 'amakihi, we genotyped 125 'amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containing SNPs and used the reference 'amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low- and high-elevation population pairs and identified loci with signatures of selection within low-elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta-defensin, glycoproteins and interleukin-related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations., (© 2018 John Wiley & Sons Ltd.)

Many species, including humans, have emerged via complex reticulate processes involving hybridisation. Under certain circumstances, hybridisation can cause distinct lineages to collapse into a single lineage with an admixed mosaic genome. Most known cases of such 'speciation reversal' or 'lineage fusion' involve recently diverged lineages and anthropogenic perturbation. Here, we show that in western North America, Common Ravens (Corvus corax) have admixed mosaic genomes formed by the fusion of non-sister lineages ('California' and 'Holarctic') that diverged ~1.5 million years ago. Phylogenomic analyses and concordant patterns of geographic structuring in mtDNA, genome-wide SNPs and nuclear introns demonstrate long-term admixture and random interbreeding between the non-sister lineages. In contrast, our genomic data support reproductive isolation between Common Ravens and Chihuahuan Ravens (C. cryptoleucus) despite extensive geographic overlap and a sister relationship between Chihuahuan Ravens and the California lineage. These data suggest that the Common Raven genome was formed by secondary lineage fusion and most likely represents a case of ancient speciation reversal that occurred without anthropogenic causes.

The amphibian skin microbiome is recognized for its role in defence against pathogens, including the deadly fungal pathogen Batrachochytrium dendrobatidis (Bd). Yet, we have little understanding of evolutionary and ecological processes that structure these communities, especially for salamanders and closely related species. We investigated patterns in the distribution of bacterial communities on Plethodon salamander skin across host species and environments. Quantifying salamander skin microbiome structure contributes to our understanding of how host-associated bacteria are distributed across the landscape, among host species, and their putative relationship with disease. We characterized skin microbiome structure (alpha-diversity, beta-diversity and bacterial operational taxonomic unit [OTU] abundances) using 16S rRNA gene sequencing for co-occurring Plethodon salamander species (35 Plethodon cinereus, 17 Plethodon glutinosus, 10 Plethodon cylindraceus) at three localities to differentiate the effects of host species from environmental factors on the microbiome. We sampled the microbiome of P. cinereus along an elevational gradient (n = 50, 700-1,000 m a.s.l.) at one locality to determine whether elevation predicts microbiome structure. Finally, we quantified prevalence and abundance of putatively anti-Bd bacteria to determine if Bd-inhibitory bacteria are dominant microbiome members. Co-occurring salamanders had similar microbiome structure, but among sites salamanders had dissimilar microbiome structure for beta-diversity and abundance of 28 bacterial OTUs. We found that alpha-diversity increased with elevation, beta-diversity and the abundance of 17 bacterial OTUs changed with elevation (16 OTUs decreasing, 1 OTU increasing). We detected 11 putatively anti-Bd bacterial OTUs that were present on 90% of salamanders and made up an average relative abundance of 83% (SD ± 8.5) per salamander. All salamanders tested negative for Bd. We conclude that environment is more influential in shaping skin microbiome structure than host differences in these congeneric species, and suggest that environmental characteristics that covary with elevation influence microbiome structure. High prevalence and abundance of anti-Bd bacteria may contribute to low Bd levels in these populations of Plethodon salamanders., (© 2017 The Authors. Journal of Animal Ecology © 2017 British Ecological Society.)

Panama is a major hub for commercial shipping between two oceans, making it an ideal location to examine parasite biogeography, potential invasions, and the spread of infectious agents. Our goals were to (i) characterise the diversity and genetic connectivity of Perkinsus spp. haplotypes across the Panamanian Isthmus and (ii) combine these data with sequences from around the world to evaluate the current phylogeography and genetic connectivity of these widespread molluscan parasites. We collected 752 bivalves from 12 locations along the coast of Panama including locations around the Bocas del Toro archipelago and the Caribbean and Pacific entrances to the Panama Canal, from December 2012 to February 2013. We used molecular genetic methods to screen for Perkinsus spp. and obtained internal transcribed spacer region (ITS) ribosomal DNA (rDNA) sequences for all positive samples. Our sequence data were used to evaluate regional haplotype diversity and distribution across both coasts of Panama, and were then combined with publicly available sequences to create global haplotype networks. We found 26 ITS haplotypes from four Perkinsus spp. (1-12 haplotypes per species) in Panama. Perkinsus beihaiensis haplotypes had the highest genetic diversity, were the most regionally widespread, and were associated with the greatest number of hosts. On a global scale, network analyses demonstrated that some haplotypes found in Panama were cosmopolitan (Perkinsus chesapeaki, Perkinsus marinus), while others were more geographically restricted (Perkinsus olseni, P. beihaiensis), indicating different levels of genetic connectivity and dispersal. We found some Perkinsus haplotypes were shared across the Isthmus of Panama and several regions around the world, including across ocean basins. We also found that haplotype diversity is currently underestimated and directly related to the number of sequences. Nevertheless, our results demonstrate long-range dispersal and global connectivity for many haplotypes, suggesting that dispersal through shipping probably contributes to these biogeographical patterns., (Published by Elsevier Ltd.)

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Aim: Lower species diversity, increased population densities and ecological niche enlargement are common characteristics of island faunas. However it remains to be determined if they extend to the parasite community. We tested if Haemosporidia parasite pressure varies between islands and the mainland with two different levels of analysis: i) at the host community level, and ii) with paired-species comparisons between islands and the mainland., Location: Gulf of Guinea, West Africa., Methods: We used molecular-based methods to identify avian Haemosporidian parasites ( Plasmodium , Haemoproteus and Leucocytozoon ) to describe their diversity, prevalence, host specificity and their phylogenetic relationships in five islands of the Gulf of Guinea and in nearby mainland areas., Results: We found reduced Haemosporidia diversity on islands for Haemoproteus and Leucocytozoon , but not for Plasmodium . In addition, lower parasite prevalence on islands was found using a paired-species approach. Although the mean host specificity of the parasite community on islands did not differ from the mainland, we found a very distinct parasite species assemblage on the islands, which was composed of both the most generalist and the most specialist lineages., Main Conclusions: This study supports the hypothesis that parasite pressure is reduced on islands. Colonization is made by generalists with high host switching capacities, with some subsequently evolving into highly specialised parasites. This suggests that 'taxon cycle' dynamics may explain the assemblage of insular parasite communities.

We engaged pet salamander owners in the United States to screen their animals for two amphibian chytrid fungal pathogens Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal). We provided pet owners with a sampling kit and instructional video to swab the skin of their animals. We received 639 salamander samples from 65 species by mail, and tested them for Bd and Bsal using qPCR. We detected Bd on 1.3% of salamanders (95% CI 0.0053-0.0267) and did not detect Bsal (95% CI 0.0000-0.0071). If Bsal is present in the U.S. population of pet salamanders, it occurs at a very low prevalence. The United States Fish and Wildlife Service listed 201 species of salamanders as "injurious wildlife" under the Lacey Act (18 U.S.C. § 42) on January 28, 2016, a precautionary action to prevent the introduction of Bsal to the U.S. through the importation of salamanders. This action reduced the number of salamanders imported to the U.S. from 2015 to 2016 by 98.4%. Our results indicate that continued precautions should be taken to prevent the introduction and establishment of Bsal in the U.S., which is a hotspot of salamander biodiversity.

Both large-wildlife loss and climatic changes can independently influence the prevalence and distribution of zoonotic disease. Given growing evidence that wildlife loss often has stronger community-level effects in low-productivity areas, we hypothesized that these perturbations would have interactive effects on disease risk. We experimentally tested this hypothesis by measuring tick abundance and the prevalence of tick-borne pathogens ( Coxiella burnetii and Rickettsia spp . ) within long-term, size-selective, large-herbivore exclosures replicated across a precipitation gradient in East Africa. Total wildlife exclusion increased total tick abundance by 130% (mesic sites) to 225% (dry, low-productivity sites), demonstrating a significant interaction of defaunation and aridity on tick abundance. When differing degrees of exclusion were tested for a subset of months, total tick abundance increased from 170% (only mega-herbivores excluded) to 360% (all large wildlife excluded). Wildlife exclusion differentially affected the abundance of the three dominant tick species, and this effect varied strongly over time, likely due to differences among species in their host associations, seasonality, and other ecological characteristics. Pathogen prevalence did not differ across wildlife exclusion treatments, rainfall levels, or tick species, suggesting that exposure risk will respond to defaunation and climate change in proportion to total tick abundance. These findings demonstrate interacting effects of defaunation and aridity that increase disease risk, and they highlight the need to incorporate ecological context when predicting effects of wildlife loss on zoonotic disease dynamics., (© 2017 The Author(s).)

White-nose syndrome, first diagnosed in North America in 2006, causes mass deaths among bats in North America. We found the causative fungus, Pseudogymnoascus destructans, in a 1918 sample collected in Europe, where bats have now adapted to the fungus. These results are consistent with a Eurasian origin of the pathogen.

Female preference for local cultural traits has been proposed as a barrier to breeding among animal populations. As such, several studies have found correlations between male bird song dialects and population genetics over relatively large distances. To investigate whether female choice for local dialects could act as a barrier to breeding between nearby and contiguous populations, we tested whether variation in male song dialects explains genetic structure among eight populations of rufous-collared sparrows (Zonotrichia capensis) in Ecuador. Our study sites lay along a transect, and adjacent study sites were separated by approximately 25 km, an order of magnitude less than previously examined for this and most other species. This transect crossed an Andean ridge and through the Quijos River Valley, both of which may be barriers to gene flow. Using a variance partitioning approach, we show that song dialect is important in explaining population genetics, independent of the geographic variables: distance, the river valley and the Andean Ridge. This result is consistent with the hypothesis that song acts as a barrier to breeding among populations in close proximity. In addition, songs of contiguous populations differed by the same degree or more than between two populations previously shown to exhibit female preference for local dialect, suggesting that birds from these populations would also breed preferentially with locals. As expected, all geographic variables (distance, the river valley and the Andean Ridge) also predicted population genetic structure. Our results have important implications for the understanding whether, and at what spatial scale, culture can affect population divergence., (© 2017 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2017 European Society For Evolutionary Biology.)

Diverse bacteria inhabit amphibian skin; some of those bacteria inhibit growth of the fungal pathogen Batrachochytrium dendrobatidis Yet there has been no systematic survey of anti- B. dendrobatidis bacteria across localities, species, and elevations. This is important given geographic and taxonomic variations in amphibian susceptibility to B. dendrobatidis Our collection sites were at locations within the Appalachian Mountains where previous sampling had indicated low B. dendrobatidis prevalence. We determined the numbers and identities of anti- B. dendrobatidis bacteria on 61 Plethodon salamanders (37 P. cinereus , 15 P. glutinosus , 9 P. cylindraceus ) via culturing methods and 16S rRNA gene sequencing. We sampled co-occurring species at three localities and sampled P. cinereus along an elevational gradient (700 to 1,000 meters above sea level [masl]) at one locality. We identified 50 anti- B. dendrobatidis bacterial operational taxonomic units (OTUs) and found that the degree of B. dendrobatidis inhibition was not correlated with relatedness. Five anti- B. dendrobatidis bacterial strains occurred on multiple amphibian species at multiple localities, but none were shared among all species and localities. The prevalence of anti- B. dendrobatidis bacteria was higher at Shenandoah National Park (NP), VA, with 96% (25/26) of salamanders hosting at least one anti- B. dendrobatidis bacterial species compared to 50% (7/14) at Catoctin Mountain Park (MP), MD, and 38% (8/21) at Mt. Rogers National Recreation Area (NRA), VA. At the individual level, salamanders at Shenandoah NP had more anti- B. dendrobatidis bacteria per individual (μ = 3.3) than those at Catoctin MP (μ = 0.8) and at Mt. Rogers NRA (μ = 0.4). All salamanders tested negative for B. dendrobatidis Anti- B. dendrobatidis bacterial species are diverse in central Appalachian Plethodon salamanders, and their distribution varied geographically. The antifungal bacterial species that we identified may play a protective role for these salamanders. IMPORTANCE Amphibians harbor skin bacteria that can kill an amphibian fungal pathogen, Batrachochytrium dendrobatidis Some amphibians die from B. dendrobatidis infection, whereas others do not. The bacteria that can kill B. dendrobatidis , called anti- B. dendrobatidis bacteria, are thought to influence the B. dendrobatidis infection outcome for the amphibian. Yet how anti- B. dendrobatidis bacterial species vary among amphibian species and populations is unknown. We determined the distribution of anti- B. dendrobatidis bacterial species among three salamander species ( n = 61) sampled at three localities. We identified 50 unique anti- B. dendrobatidis bacterial species and found that all of the tested salamanders were negative for B. dendrobatidis Five anti- B. dendrobatidis bacterial species were commonly detected, suggesting a stable, functional association with these salamanders. The number of anti- B. dendrobatidis bacteria per individual varied among localities but not among co-occurring salamander species, demonstrating that environment is more influential than host factors in structuring the anti- B. dendrobatidis bacterial community. These anti- B. dendrobatidis bacteria may serve a protective function for their salamander hosts., (Copyright © 2017 American Society for Microbiology.)