Reduction of genetic diversity in 'Alalā (Hawaiian crow; Corvus hawaiiensis) between the late 1800s and the late 1900s

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.

Linking avian malaria parasitemia estimates from quantitative PCR and microscopy reveals new infection patterns in Hawai'i

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.

Sequence capture identifies fastidious chytrid fungi directly from host tissue

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.

Microbiomes associated with avian malaria survival differ between susceptible Hawaiian honeycreepers and sympatric malaria-resistant introduced birds

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.

Gene expression reveals immune response strategies of naïve Hawaiian honeycreepers experimentally infected with introduced avian malaria

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.

Immune priming prior to pathogen exposure sheds light on the relationship between host, microbiome and pathogen in disease

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.

Museum Genomics Provide Evidence for Persistent Genetic Differentiation in a Threatened Seabird Species in the Western Atlantic

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.

Transcriptional response of individual Hawaiian Culex quinquefasciatus mosquitoes to the avian malaria parasite Plasmodium relictum

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04271-x.

A comparative study of RNA yields from museum specimens, including an optimized protocol for extracting RNA from formalin-fixed specimens

Animal specimens in natural history collections are invaluable resources in examining the historical context of pathogen dynamics in wildlife and spillovers to humans. For example, natural history specimens may reveal new associations between bat species and coronaviruses. However, RNA viruses are difficult to study in historical specimens because protocols for extracting RNA from these specimens have not been optimized. Advances have been made in our ability to recover nucleic acids from formalin-fixed paraffin-embedded samples (FFPE) commonly used in human clinical studies, yet other types of formalin preserved samples have received less attention. Here, we optimize the recovery of RNA from formalin-fixed ethanol-preserved museum specimens in order to improve the usability of these specimens in surveys for zoonotic diseases. We provide RNA quality and quantity measures for replicate tissues subsamples of 22 bat specimens from five bat genera (Rhinolophus, Hipposideros, Megareops, Cynopterus, and Nyctalus) collected in China and Myanmar from 1886 to 2003. As tissues from a single bat specimen were preserved in a variety of ways, including formalin-fixed (8 bats), ethanol-preserved and frozen (13 bats), and flash frozen (2 bats), we were able to compare RNA quality and yield across different preservation methods. RNA extracted from historical museum specimens is highly fragmented, but usable for short-read sequencing and targeted amplification. Incubation of formalin-fixed samples with Proteinase-K following thorough homogenization improves RNA yield. This optimized protocol extends the types of data that can be derived from existing museum specimens and facilitates future examinations of host and pathogen RNA from specimens.

Conservation genomics and systematics of a near-extinct island radiation

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 mtDNA, nuclear 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 geographic 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 predate modern anthropogenic threats require further study with denser population sampling. Our study highlights the value 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 fragments their ranges and contributes to rapid declines in population abundances.

An efficient method for simultaneous species, individual, and sex identification via in-solution single nucleotide polymorphism capture from low-quality scat samples

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.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra

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.

Genetically modifying skin microbe to produce violacein and augmenting microbiome did not defend Panamanian golden frogs from disease

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.

Adaptive Radiation Genomics of Two Ecologically Divergent Hawai'ian Honeycreepers: The 'akiapōlā'au and the Hawai'i 'amakihi

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.

Transcriptome assembly and differential gene expression of the invasive avian malaria parasite Plasmodium relictum in Hawai'i

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.

Comparative Analysis of Annotation Pipelines Using the First Japanese White-Eye (Zosterops japonicus) Genome

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.

Identification of novel bacterial biomarkers to detect bird scavenging by invasive rats

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.

A genome-wide investigation of adaptive signatures in protein-coding genes related to tool behaviour in New Caledonian and Hawaiian crows

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.

Sustained immune activation is associated with susceptibility to the amphibian chytrid fungus

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.

Soil fungal communities differ between shaded and sun-intensive coffee plantations in El Salvador

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.

The role of native and introduced birds in transmission of avian malaria in Hawaii

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.

Opening the door to greater phylogeographic inference in Southeast Asia: Comparative genomic study of five codistributed rainforest bird species using target capture and historical DNA

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.

Comparing Adaptive Radiations Across Space, Time, and Taxa

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.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks

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.

Diversity and temporal dynamics of primate milk microbiomes

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.

Parthenogenesis in a captive Asian water dragon (Physignathus cocincinus) identified with novel microsatellites

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.

Fungal disease and temperature alter skin microbiome structure in an experimental salamander system

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.

First Report of a Novel Hepatozoon sp. in Giant Pandas (Ailuropoda melanoleuca)

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.

North-facing slopes and elevation shape asymmetric genetic structure in the range-restricted salamander Plethodon shenandoah

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.

Parallel evolution of gene classes, but not genes: Evidence from Hawai'ian honeycreeper populations exposed to avian malaria

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.

White-Nose Syndrome Fungus in a 1918 Bat Specimen from France

White-nose syndrome, first diagnosed in North America in 2006, causes mass deaths among bats in North America. We found the causative fungus, Pseudogymnoascusdestructans, 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.

Genomic evidence of speciation reversal in ravens

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.

Interacting effects of wildlife loss and climate on ticks and tick-borne disease

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.

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

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.

Insularity effects on the assemblage of the blood parasite community of the birds from the Gulf of Guinea

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.

Batrachochytrium salamandrivorans not detected in U.S. survey of pet salamanders

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.

KNOWN SECONDARY CONTACT AND RAPID GENE FLOW AMONG SUBSPECIES AND DIALECTS IN THE BROWN-HEADED COWBIRD

We analyze geographic variation in morphology for Sierra Nevadan populations of Brown-headed Cowbirds (Molothrus ater) in relation to two levels of geographic structure of populations: 1) subspecies and 2) vocal dialects. Two morphologically distinct subspecies occupy opposite slopes of the Sierra Nevada: M. a artemisiae on the east slope is larger than M. a. obscurus on the west slope, and its juveniles have predominantly white versus yellow rictal (bill) flanges. Populations of obscurus moved into California from the lower Colorado River around 1900 and invaded the western Sierra during the 1930's. A relatively constant type of flight whistle occurs in obscurus populations up to 700 km apart, but east-slope artemisiae are divided into many distinct dialects. The means of seven morphological measurements and of principal component 1 (i.e., body size) for 2,287 individuals show similar clinal patterns for adult males, yearling males, and females over a 60-km north-south transect along the eastern Sierra Nevada: large sizes occur at the southernmost study site at Bishop, decrease clinally to the north to Mammoth Lakes, and then increase clinally to the north from Mammoth Lakes to the northernmost site at Lee Vining. This reversal in cline slope is reiterated for variation in the frequency of white-flanged juveniles. Our data do not support adaptive explanations for the clines. Rather, the most parsimonious explanation is that there has been flow of obscurus genes into artemisiae from the west slope across the crest of the Sierra Nevada into the Mammoth Lakes area. Concordantly, the Sierran crest at Mammoth Lakes is considerably lower than any other point within about 80 km to the north or 180 km to the south. Two other findings strongly support this interpretation. First, wing lengths of our longest-winged samples (at Bishop and Lee Vining) are significantly shorter than those of cowbirds collected in the same region between 1912 and 1922. Second, the wing lengths of Mammoth Lakes adult males and females during 1978-1981 are significantly larger than those of the 1982-1985 period. This putative hybridization must have been rapid, as there have been fewer than 50 years since secondary contact could have occurred. Calculation of the gene-flow parameter Ι of Endler (1977) indicates that the gene-flow rate reported here is higher than for nearly all of the species he surveyed. This gene flow occurs between two subspecies with different flight whistles and across several cowbird flight-whistle dialects within artemisiae. Therefore, vocal differences among these dialects do not appear to be a strong deterrent to gene flow.

Genetic population structure in an equatorial sparrow: roles for culture and geography

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.