Comparative Genomics of the Waterfowl Innate Immune System.

Animal species differ considerably in their ability to fight off infections. Finding the genetic basis of these differences is not easy, as the immune response is comprised of a complex network of proteins that interact with one another to defend the body against infection. Here, we used population- and comparative genomics to study the evolutionary forces acting on the innate immune system in natural hosts of the avian influenza virus (AIV). For this purpose, we used a combination of hybrid capture, next- generation sequencing and published genomes to examine genetic diversity, divergence, and signatures of selection in 127 innate immune genes at a micro- and macroevolutionary time scale in 26 species of waterfowl. We show across multiple immune pathways (AIV-, toll-like-, and RIG-I -like receptors signalling pathways) that genes involved genes in pathogen detection (i.e., toll-like receptors) and direct pathogen inhibition (i.e., antimicrobial peptides and interferon-stimulated genes), as well as host proteins targeted by viral antagonist proteins (i.e., mitochondrial antiviral-signaling protein, [MAVS]) are more likely to be polymorphic, genetically divergent, and under positive selection than other innate immune genes. Our results demonstrate that selective forces vary across innate immune signaling signalling pathways in waterfowl, and we present candidate genes that may contribute to differences in susceptibility and resistance to infectious diseases in wild birds, and that may be manipulated by viruses. Our findings improve our understanding of the interplay between host genetics and pathogens, and offer the opportunity for new insights into pathogenesis and potential drug targets.

Avian Immunome DB: an example of a user-friendly interface for extracting genetic information.

BACKGROUND: Genomic and genetic studies often require a target list of genes before conducting any hypothesis testing or experimental verification. With the ever-growing number of sequenced genomes and a variety of different annotation strategies, comes the potential for ambiguous gene symbols, making it cumbersome to capture the "correct" set of genes. In this article, we present and describe the Avian Immunome DB (AVIMM) for easy gene property extraction as exemplified by avian immune genes. The avian immune system is characterised by a cascade of complex biological processes underlaid by more than 1000 different genes. It is a vital trait to study particularly in birds considering that they are a significant driver in spreading zoonotic diseases. With the completion of phase II of the B10K ("Bird 10,000 Genomes") consortium's whole-genome sequencing effort, we have included 363 annotated bird genomes in addition to other publicly available bird genome data which serve as a valuable foundation for AVIMM. CONSTRUCTION AND CONTENT: A relational database with avian immune gene evidence from Gene Ontology, Ensembl, UniProt and the B10K consortium has been designed and set up. The foundation stone or the "seed" for the initial set of avian immune genes is based on the well-studied model organism chicken (Gallus gallus). Gene annotations, different transcript isoforms, nucleotide sequences and protein information, including amino acid sequences, are included. Ambiguous gene names (symbols) are resolved within the database and linked to their canonical gene symbol. AVIMM is supplemented by a command-line interface and a web front-end to query the database. UTILITY AND DISCUSSION: The internal mapping of unique gene symbol identifiers to canonical gene symbols allows for an ambiguous gene property search. The database is organised within core and feature tables, which makes it straightforward to extend for future purposes. The database design is ready to be applied to other taxa or biological processes. Currently, the database contains 1170 distinct avian immune genes with canonical gene symbols and 612 synonyms across 363 bird species. While the command-line interface readily integrates into bioinformatics pipelines, the intuitive web front-end with download functionality offers sophisticated search functionalities and tracks the origin for each record. AVIMM is publicly accessible at https://avimm.ab.mpg.de .

Contrasting effects of host species and phylogenetic diversity on the occurrence of HPAI H5N1 in European wild birds.

Studies on the highly pathogenic avian influenza (HPAI) H5N1 suggest that wild bird migration may facilitate its long-distance spread, yet the role of wild bird community composition in its transmission risk remains poorly understood. Furthermore, most studies on the diversity-disease relationship focused on host species diversity without considering hosts' phylogenetic relationships, which may lead to rejecting a species diversity effect when the community has host species that are only distantly related. Here, we explored the influence of waterbird community composition for determining HPAI H5N1 occurrence in wild birds in a continental-scale study across Europe. In particular, we tested the diversity-disease relationship using both host species diversity and host phylogenetic diversity. Our results provide the first demonstration that host community composition-compared with previously identified environmental risk factors-can also effectively explain the spatial pattern of H5N1 occurrence in wild birds. We further show that communities with more higher risk host species and more closely related species have a higher risk of H5N1 outbreaks. Thus, both host species diversity and community phylogenetic structure, in addition to environmental factors, jointly influence H5N1 occurrence. Our work not only extends the current theory on the diversity-disease relationship, but also has important implications for future monitoring of H5N1 and other HPAI subtypes.

A history of hybrids? Genomic patterns of introgression in the True Geese.

BACKGROUND: The impacts of hybridization on the process of speciation are manifold, leading to distinct patterns across the genome. Genetic differentiation accumulates in certain genomic regions, while divergence is hampered in other regions by homogenizing gene flow, resulting in a heterogeneous genomic landscape. A consequence of this heterogeneity is that genomes are mosaics of different gene histories that can be compared to unravel complex speciation and hybridization events. However, incomplete lineage sorting (often the outcome of rapid speciation) can result in similar patterns. New statistical techniques, such as the D-statistic and hybridization networks, can be applied to disentangle the contributions of hybridization and incomplete lineage sorting. We unravel patterns of hybridization and incomplete lineage sorting during and after the diversification of the True Geese (family Anatidae, tribe Anserini, genera Anser and Branta) using an exon-based hybridization network approach and taking advantage of discordant gene tree histories by re-sequencing all taxa of this clade. In addition, we determine the timing of introgression and reconstruct historical effective population sizes for all goose species to infer which demographic or biogeographic factors might explain the observed patterns of introgression. RESULTS: We find indications for ancient interspecific gene flow during the diversification of the True Geese and were able to pinpoint several putative hybridization events. Specifically, in the genus Branta, both the ancestor of the White-cheeked Geese (Hawaiian Goose, Canada Goose, Cackling Goose and Barnacle Goose) and the ancestor of the Brent Goose hybridized with Red-breasted Goose. One hybridization network suggests a hybrid origin for the Red-breasted Goose, but this scenario seems unlikely and it not supported by the D-statistic analysis. The complex, highly reticulated evolutionary history of the genus Anser hampered the estimation of ancient hybridization events by means of hybridization networks. The reconstruction of historical effective population sizes shows that most species showed a steady increase during the Pliocene and Pleistocene. These large effective population sizes might have facilitated contact between diverging goose species, resulting in the establishment of hybrid zones and consequent gene flow. CONCLUSIONS: Our analyses suggest that the evolutionary history of the True Geese is influenced by introgressive hybridization. The approach that we have used, based on genome-wide phylogenetic incongruence and network analyses, will be a useful procedure to reconstruct the complex evolutionary histories of many naturally hybridizing species groups.

The Evolution of Innate Immune Genes: Purifying and Balancing Selection on ß-Defensins in Waterfowl.

In disease dynamics, high immune gene diversity can confer a selective advantage to hosts in the face of a rapidly evolving and diverse pathogen fauna. This is supported empirically for genes involved in pathogen recognition and signalling. In contrast, effector genes involved in pathogen clearance may be more constrained. ß-Defensins are innate immune effector genes; their main mode of action is via disruption of microbial membranes. Here, five ß-defensin genes were characterized in mallards (Anas platyrhynchos) and other waterfowl; key reservoir species for many zoonotic diseases. All five genes showed remarkably low diversity at the individual-, population-, and species-level. Furthermore, there was widespread sharing of identical alleles across species divides. Thus, specific ß-defensin alleles were maintained not only spatially but also over long temporal scales, with many amino acid residues being fixed across all species investigated. Purifying selection to maintain individual, highly efficacious alleles was the primary evolutionary driver of these genes in waterfowl. However, we also found evidence for balancing selection acting on the most recently duplicated ß-defensin gene (AvBD3b). For this gene, we found that amino acid replacements were more likely to be radical changes, suggesting that duplication of ß-defensin genes allows exploration of wider functional space. Structural conservation to maintain function appears to be crucial for avian ß-defensin effector molecules, resulting in low tolerance for new allelic variants. This contrasts with other types of innate immune genes, such as receptor and signalling molecules, where balancing selection to maintain allelic diversity has been shown to be a strong evolutionary force.

Genetic isolation between coastal and fishery-impacted, offshore bottlenose dolphin (Tursiops spp.) populations.

The identification of species and population boundaries is important in both evolutionary and conservation biology. In recent years, new population genetic and computational methods for estimating population parameters and testing hypotheses in a quantitative manner have emerged. Using a Bayesian framework and a quantitative model-testing approach, we evaluated the species status and genetic connectedness of bottlenose dolphin (Tursiops spp.) populations off remote northwestern Australia, with a focus on pelagic 'offshore' dolphins subject to incidental capture in a trawl fishery. We analysed 71 dolphin samples from three sites beyond the 50 m depth contour (the inshore boundary of the fishery) and up to 170 km offshore, including incidentally caught and free-ranging individuals associating with trawl vessels, and 273 dolphins sampled at 12 coastal sites inshore of the 50 m depth contour and within 10 km of the coast. Results from 19 nuclear microsatellite markers showed significant population structure between dolphins from within the fishery and coastal sites, but also among dolphins from coastal sites, identifying three coastal populations. Moreover, we found no current or historic gene flow into the offshore population in the region of the fishery, indicating a complete lack of recruitment from coastal sites. Mitochondrial DNA corroborated our findings of genetic isolation between dolphins from the offshore population and coastal sites. Most offshore individuals formed a monophyletic clade with common bottlenose dolphins (T. truncatus), while all 273 individuals sampled coastally formed a well-supported clade of Indo-Pacific bottlenose dolphins (T. aduncus). By including a quantitative modelling approach, our study explicitly took evolutionary processes into account for informing the conservation and management of protected species. As such, it may serve as a template for other, similarly inaccessible study populations.

A tree of geese: A phylogenomic perspective on the evolutionary history of True Geese.

Phylogenetic incongruence can be caused by analytical shortcomings or can be the result of biological processes, such as hybridization, incomplete lineage sorting and gene duplication. Differentiation between these causes of incongruence is essential to unravel complex speciation and diversification events. The phylogeny of the True Geese (tribe Anserini, Anatidae, Anseriformes) was, until now, contentious, i.e., the phylogenetic relationships and the timing of divergence between the different goose species could not be fully resolved. We sequenced nineteen goose genomes (representing seventeen species of which three subspecies of the Brent Goose, Branta bernicla) and used an exon-based phylogenomic approach (41,736 exons, representing 5887 genes) to unravel the evolutionary history of this bird group. We thereby provide general guidance on the combination of whole genome evolutionary analyses and analytical tools for such cases where previous attempts to resolve the phylogenetic history of several taxa could not be unravelled. Identical topologies were obtained using either a concatenation (based upon an alignment of 6,630,626 base pairs) or a coalescent-based consensus method. Two major lineages, corresponding to the genera Anser and Branta, were strongly supported. Within the Branta lineage, the White-cheeked Geese form a well-supported sub-lineage that is sister to the Red-breasted Goose (Branta ruficollis). In addition, two main clades of Anser species could be identified, the White Geese and the Grey Geese. The results from the consensus method suggest that the diversification of the genus Anser is heavily influenced by rapid speciation and by hybridization, which may explain the failure of previous studies to resolve the phylogenetic relationships within this genus. The majority of speciation events took place in the late Pliocene and early Pleistocene (between 4 and 2millionyears ago), conceivably driven by a global cooling trend that led to the establishment of a circumpolar tundra belt and the emergence of temperate grasslands. Our approach will be a fruitful strategy for resolving many other complex evolutionary histories at the level of genera, species, and subspecies.

Population-level consequences of complementary sex determination in a solitary parasitoid.

BACKGROUND: Sex determination mechanisms are known to be evolutionarily labile but the factors driving transitions in sex determination mechanisms are poorly understood. All insects of the Hymenoptera are haplodiploid, with males normally developing from unfertilized haploid eggs. Under complementary sex determination (CSD), diploid males can be produced from fertilized eggs that are homozygous at the sex locus. Diploid males have near-zero fitness and thus represent a genetic load, which is especially severe under inbreeding. Here, we study mating structure and sex determination in the parasitoid Cotesia vestalis to investigate what may have driven the evolution of two complementary sex determination loci in this species. RESULTS: We genotyped Cotesia vestalis females collected from eight fields in four townships in Western Taiwan. 98 SNP markers were developed by aligning Illumina sequence reads of pooled DNA of eight different females against a de novo assembled genome of C. vestalis. This proved to be an efficient method for this non-model species and provides a resource for future use in related species. We found significant genetic differentiation within the sampled population but variation could not be attributed to sampling locations by AMOVA. Non-random mating was detected, with 8.1% of matings between siblings. Diploid males, detected by flow cytometry, were produced at a rate of 1.4% among diploids. CONCLUSIONS: We think that the low rate of diploid male production is best explained by a CSD system with two independent sex loci, supporting laboratory findings on the same species. Fitness costs of diploid males in C. vestalis are high because diploid males can mate with females and produce infertile triploid offspring. This severe fitness cost of diploid males combined with non-random mating may have resulted in evolution from single locus CSD to CSD with two independent loci.

Generation of SNP datasets for orangutan population genomics using improved reduced-representation sequencing and direct comparisons of SNP calling algorithms.

BACKGROUND: High-throughput sequencing has opened up exciting possibilities in population and conservation genetics by enabling the assessment of genetic variation at genome-wide scales. One approach to reduce genome complexity, i.e. investigating only parts of the genome, is reduced-representation library (RRL) sequencing. Like similar approaches, RRL sequencing reduces ascertainment bias due to simultaneous discovery and genotyping of single-nucleotide polymorphisms (SNPs) and does not require reference genomes. Yet, generating such datasets remains challenging due to laboratory and bioinformatical issues. In the laboratory, current protocols require improvements with regards to sequencing homologous fragments to reduce the number of missing genotypes. From the bioinformatical perspective, the reliance of most studies on a single SNP caller disregards the possibility that different algorithms may produce disparate SNP datasets. RESULTS: We present an improved RRL (iRRL) protocol that maximizes the generation of homologous DNA sequences, thus achieving improved genotyping-by-sequencing efficiency. Our modifications facilitate generation of single-sample libraries, enabling individual genotype assignments instead of pooled-sample analysis. We sequenced ~1% of the orangutan genome with 41-fold median coverage in 31 wild-born individuals from two populations. SNPs and genotypes were called using three different algorithms. We obtained substantially different SNP datasets depending on the SNP caller. Genotype validations revealed that the Unified Genotyper of the Genome Analysis Toolkit and SAMtools performed significantly better than a caller from CLC Genomics Workbench (CLC). Of all conflicting genotype calls, CLC was only correct in 17% of the cases. Furthermore, conflicting genotypes between two algorithms showed a systematic bias in that one caller almost exclusively assigned heterozygotes, while the other one almost exclusively assigned homozygotes. CONCLUSIONS: Our enhanced iRRL approach greatly facilitates genotyping-by-sequencing and thus direct estimates of allele frequencies. Our direct comparison of three commonly used SNP callers emphasizes the need to question the accuracy of SNP and genotype calling, as we obtained considerably different SNP datasets depending on caller algorithms, sequencing depths and filtering criteria. These differences affected scans for signatures of natural selection, but will also exert undue influences on demographic inferences. This study presents the first effort to generate a population genomic dataset for wild-born orangutans with known population provenance.

Global lack of flyway structure in a cosmopolitan bird revealed by a genome wide survey of single nucleotide polymorphisms.

Knowledge about population structure and connectivity of waterfowl species, especially mallards (Anas platyrhynchos), is a priority because of recent outbreaks of avian influenza. Ringing studies that trace large-scale movement patterns have to date been unable to detect clearly delineated mallard populations. We employed 363 single nucleotide polymorphism markers in combination with population genetics and phylogeographical approaches to conduct a population genomic test of panmixia in 801 mallards from 45 locations worldwide. Basic population genetic and phylogenetic methods suggest no or very little population structure on continental scales. Nor could individual-based structuring algorithms discern geographical structuring. Model-based coalescent analyses for testing models of population structure pointed to strong genetic connectivity among the world's mallard population. These diverse approaches all support the conclusion that there is a lack of clear population structure, suggesting that the world's mallards, perhaps with minor exceptions, form a single large, mainly interbreeding population.

The meaning of wild: Genetic and adaptive consequences from large-scale releases of domestic mallards.

The translocation of individuals around the world is leading to rising incidences of anthropogenic hybridization, particularly between domestic and wild congeners. We apply a landscape genomics approach for thousands of mallard (Anas platyrhynchos) samples across continental and island populations to determine the result of over a century of supplementation practices. We establish that a single domestic game-farm mallard breed is the source for contemporary release programs in Eurasia and North America, as well as for established feral populations in New Zealand and Hawaii. In particular, we identify central Europe and eastern North America as epicenters of ongoing anthropogenic hybridization, and conclude that the release of game-farm mallards continues to affect the genetic integrity of wild mallards. Conversely, self-sustaining feral populations in New Zealand and Hawaii not only show strong differentiation from their original stock, but also signatures of local adaptation occurring in less than a half-century since game-farm mallard releases have ceased. We conclude that 'wild' is not singular, and that even feral populations are capable of responding to natural processes. Although considered paradoxical to biological conservation, understanding the capacity for wildness among feral and feral admixed populations in human landscapes is critical as such interactions increase in the Anthropocene.

Evaluating Effects of AIV Infection Status on Ducks Using a Flow Cytometry-Based Differential Blood Count.

Ducks have recently received a lot of attention from the research community due to their importance as natural reservoirs of avian influenza virus (AIV). Still, there is a lack of tools to efficiently determine the immune status of ducks. The purpose of this work was to develop an automated differential blood count for the mallard duck (Anas platyrhynchos), to assess reference values of white blood cell (WBC) counts in this species, and to apply the protocol in an AIV field study. We established a flow cytometry-based duck WBC differential based on a no-lyse no-wash single-step one-tube technique, applying a combination of newly generated monoclonal antibodies with available duck-specific as well as cross-reacting chicken markers. The blood cell count enables quantification of mallard thrombocytes, granulocytes, monocytes, B cells, CD4+ T cells (T helper) and CD8+ cytotoxic T cells. The technique is reproducible, accurate, and much faster than traditional evaluations of blood smears. Stabilization of blood samples enables analysis up to 1 week after sampling, thus allowing for evaluation of blood samples collected in the field. We used the new technique to investigate a possible influence of sex, age, and AIV infection status on WBC counts in wild mallards. We show that age has an effect on the WBC counts in mallards, as does sex in juvenile mallards. Interestingly, males naturally infected with low pathogenic AIV showed a reduction of lymphocytes (lymphocytopenia) and thrombocytes (thrombocytopenia), which are both common in influenza A infection in humans. IMPORTANCE Outbreaks of avian influenza in poultry and humans are a global public health concern. Aquatic birds are the primary natural reservoir of avian influenza viruses (AIVs), and strikingly, AIVs mainly cause asymptomatic or mild infection in these species. Hence, immunological studies in aquatic birds are important for investigating variation in disease outcome of different hosts to AIV and may aid in early recognition and a better understanding of zoonotic events. Unfortunately, immunological studies in these species were so far hampered by the lack of diagnostic tools. Here, we present a technique that enables high-throughput white blood cell (WBC) analysis in the mallard and report changes in WBC counts in wild mallards naturally infected with AIV. Our protocol permits large-scale immune status monitoring in a widespread wild and domesticated duck species and provides a tool to further investigate the immune response in an important reservoir host of zoonotic viruses.

A Comparison of Host Responses to Infection with Wild-Type Avian Influenza Viruses in Chickens and Tufted Ducks.

Cross-species transmission of influenza A virus (IAV) from wild waterfowl to poultry is the first step in a chain of events that can ultimately lead to exposure and infection of humans. Herein, we study the outcome of infection with eight different mallard-origin IAV subtypes in two different avian hosts: tufted ducks and chickens. We found that infection and shedding patterns as well as innate immune responses were highly dependent on viral subtypes, host species, and inoculation routes. For example, intraoesophageal inoculation, commonly used in mallard infection experiments, resulted in no infections in contrast to oculonasal inoculation, suggesting a difference in transmission routes. Despite H9N2 being endemic in chickens, inoculation of mallard-origin H9N2 failed to cause viable infection beyond 1 day postinfection in our study design. The innate immune responses were markedly different in chickens and tufted ducks, and despite the presence of retinoic acid-inducible gene-I (RIG-I) in tufted duck transcriptomes, it was neither up nor downregulated in response to infection. Overall, we have revealed the heterogeneity of infection patterns and responses in two markedly different avian hosts following a challenge with mallard-origin IAV. These virus-host interactions provide new insights into important aspects of interspecies transmission of IAV. IMPORTANCE Our current findings highlight important aspects of IAV infection in birds that have implications for our understanding of its zoonotic ecology. In contrast to mallards where the intestinal tract is the main site of IAV replication, chickens and tufted ducks show limited or no signs of intestinal infection suggesting that the fecal-oral transmission route might not apply to all bird IAV host species. Our results indicate that mallard-origin IAVs undergo genetic changes upon introduction into new hosts, suggesting rapid adaptation to a new environment. However, similar to the mallard, chickens and tufted ducks show a limited immune response to infection with low pathogenic avian influenza viruses. These findings and future studies in different IAV hosts are important for our understanding of barriers to IAV transmission between species and ultimately from the wild reservoir to humans.

Widespread horizontal genomic exchange does not erode species barriers among sympatric ducks.

BACKGROUND: The study of speciation and maintenance of species barriers is at the core of evolutionary biology. During speciation the genome of one population becomes separated from other populations of the same species, which may lead to genomic incompatibility with time. This separation is complete when no fertile offspring is produced from inter-population matings, which is the basis of the biological species concept. Birds, in particular ducks, are recognised as a challenging and illustrative group of higher vertebrates for speciation studies. There are many sympatric and ecologically similar duck species, among which fertile hybrids occur relatively frequently in nature, yet these species remain distinct. RESULTS: We show that the degree of shared single nucleotide polymorphisms (SNPs) between five species of dabbling ducks (genus Anas) is an order of magnitude higher than that previously reported between any pair of eukaryotic species with comparable evolutionary distances. We demonstrate that hybridisation has led to sustained exchange of genetic material between duck species on an evolutionary time scale without disintegrating species boundaries. Even though behavioural, genetic and ecological factors uphold species boundaries in ducks, we detect opposing forces allowing for viable interspecific hybrids, with long-term evolutionary implications. Based on the superspecies concept we here introduce the novel term "supra-population" to explain the persistence of SNPs identical by descent within the studied ducks despite their history as distinct species dating back millions of years. CONCLUSIONS: By reviewing evidence from speciation theory, palaeogeography and palaeontology we propose a fundamentally new model of speciation to accommodate our genetic findings in dabbling ducks. This model, we argue, may also shed light on longstanding unresolved general speciation and hybridisation patterns in higher organisms, e.g. in other bird groups with unusually high hybridisation rates. Observed parallels to horizontal gene transfer in bacteria facilitate the understanding of why ducks have been such an evolutionarily successful group of animals. There is large evolutionary potential in the ability to exchange genes among species and the resulting dramatic increase of effective population size to counter selective constraints.

Reference genomes for conservation.

High-quality reference genomes for non-model species can benefit conservation.

Genome wide SNP discovery, analysis and evaluation in mallard (Anas platyrhynchos).

BACKGROUND: Next generation sequencing technologies allow to obtain at low cost the genomic sequence information that currently lacks for most economically and ecologically important organisms. For the mallard duck genomic data is limited. The mallard is, besides a species of large agricultural and societal importance, also the focal species when it comes to long distance dispersal of Avian Influenza. For large scale identification of SNPs we performed Illumina sequencing of wild mallard DNA and compared our data with ongoing genome and EST sequencing of domesticated conspecifics. This is the first study of its kind for waterfowl. RESULTS: More than one billion base pairs of sequence information were generated resulting in a 16× coverage of a reduced representation library of the mallard genome. Sequence reads were aligned to a draft domesticated duck reference genome and allowed for the detection of over 122,000 SNPs within our mallard sequence dataset. In addition, almost 62,000 nucleotide positions on the domesticated duck reference showed a different nucleotide compared to wild mallard. Approximately 20,000 SNPs identified within our data were shared with SNPs identified in the sequenced domestic duck or in EST sequencing projects. The shared SNPs were considered to be highly reliable and were used to benchmark non-shared SNPs for quality. Genotyping of a representative sample of 364 SNPs resulted in a SNP conversion rate of 99.7%. The correlation of the minor allele count and observed minor allele frequency in the SNP discovery pool was 0.72. CONCLUSION: We identified almost 150,000 SNPs in wild mallards that will likely yield good results in genotyping. Of these, ~101,000 SNPs were detected within our wild mallard sequences and ~49,000 were detected between wild and domesticated duck data. In the ~101,000 SNPs we found a subset of ~20,000 SNPs shared between wild mallards and the sequenced domesticated duck suggesting a low genetic divergence. Comparison of quality metrics between the total SNP set (122,000 + 62,000 = 184,000 SNPs) and the validated subset shows similar characteristics for both sets. This indicates that we have detected a large amount (~150,000) of accurately inferred mallard SNPs, which will benefit bird evolutionary studies, ecological studies (e.g. disentangling migratory connectivity) and industrial breeding programs.

Evolution and connectivity in the world-wide migration system of the mallard: inferences from mitochondrial DNA.

BACKGROUND: Main waterfowl migration systems are well understood through ringing activities. However, in mallards (Anas platyrhynchos) ringing studies suggest deviations from general migratory trends and traditions in waterfowl. Furthermore, surprisingly little is known about the population genetic structure of mallards, and studying it may yield insight into the spread of diseases such as Avian Influenza, and in management and conservation of wetlands. The study of evolution of genetic diversity and subsequent partitioning thereof during the last glaciation adds to ongoing discussions on the general evolution of waterfowl populations and flyway evolution. Hypothesised mallard flyways are tested explicitly by analysing mitochondrial mallard DNA from the whole northern hemisphere. RESULTS: Phylogenetic analyses confirm two mitochondrial mallard clades. Genetic differentiation within Eurasia and North-America is low, on a continental scale, but large differences occur between these two land masses (F(ST) = 0.51). Half the genetic variance lies within sampling locations, and a negligible portion between currently recognised waterfowl flyways, within Eurasia and North-America. Analysis of molecular variance (AMOVA) at continent scale, incorporating sampling localities as smallest units, also shows the absence of population structure on the flyway level. Finally, demographic modelling by coalescence simulation proposes a split between Eurasia and North-America 43,000 to 74,000 years ago and strong population growth (~100fold) since then and little migration (not statistically different from zero). CONCLUSIONS: Based on this first complete assessment of the mallard's world-wide population genetic structure we confirm that no more than two mtDNA clades exist. Clade A is characteristic for Eurasia, and clade B for North-America although some representatives of clade A are also found in North-America. We explain this pattern by evaluating competing hypotheses and conclude that a complex mix of historical, recent and anthropogenic factors shaped the current mallard populations. We refute population classification based on flyways proposed by ornithologists and managers, because they seem to have little biological meaning. Our results have implications for wetland management and conservation, with special regard to the release of farmed mallards for hunting, as well as for the possible transmission of Avian Influenza by mallards due to migration.

Health monitoring in birds using bio-loggers and whole blood transcriptomics.

Monitoring and early detection of emerging infectious diseases in wild animals is of crucial global importance, yet reliable ways to measure immune status and responses are lacking for animals in the wild. Here we assess the usefulness of bio-loggers for detecting disease outbreaks in free-living birds and confirm detailed responses using leukocyte composition and large-scale transcriptomics. We simulated natural infections by viral and bacterial pathogens in captive mallards (Anas platyrhynchos), an important natural vector for avian influenza virus. We show that body temperature, heart rate and leukocyte composition change reliably during an acute phase immune response. Using genome-wide gene expression profiling of whole blood across time points we confirm that immunostimulants activate pathogen-specific gene regulatory networks. By reporting immune response related changes in physiological and behavioural traits that can be studied in free-ranging populations, we provide baseline information with importance to the global monitoring of zoonotic diseases.