Pim van Hooft

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

BackgroundNext 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.ResultsMore 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.ConclusionWe 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.

Genetic monitoring detects an overlooked cryptic species and reveals the diversity and distribution of three invasive Rattus congeners in South Africa.

BackgroundSouth Africas long and extensive trade activity has ensured ample opportunities for exotic species introduction. Whereas the rich biodiversity of endemic southern African fauna has been the focus of many studies, invasive vertebrates are generally overlooked despite potential impacts on biodiversity, health and agriculture. Genetic monitoring of commensal rodents in South Africa which uncovered the presence of Rattus tanezumi, a South-East Asian endemic not previously known to occur in Africa, provided the impetus for expanded studies on all invasive Rattus species present.ResultsTo this end, intensified sampling at 28 South African localities and at one site in Swaziland, identified 149 Rattus specimens. Cytochrome b gene sequencing revealed the presence of two R. tanezumi, seven Rattus rattus and five Rattus norvegicus haplotypes in south Africa. Phylogenetic results were consistent with a single, recent R. tanezumi introduction and indicated that R. norvegicus and R. rattus probably became established following at least two and three independent introductions, respectively. Intra- and inter-specific diversity was highest in informal human settlements, with all three species occurring at a single metropolitan township site. Rattus norvegicus and R. rattus each occurred sympatrically with Rattus tanezumi at one and five sites, respectively. Karyotyping of selected R. rattus and R. tanezumi individuals identified diploid numbers consistent with those reported previously for these cryptic species. Ordination of bioclimatic variables and MaxEnt ecological niche modelling confirmed that the bioclimatic niche occupied by R. tanezumi in south Africa was distinct from that occupied in its naturalised range in south-east Asia suggesting that factors other than climate may influence the distribution of this species.ConclusionsThis study has highlighted the value of genetic typing for detecting cryptic invasive species, providing historical insights into introductions and for directing future sampling. The apparent ease with which a cryptic species can become established signals the need for broader implementation of genetic monitoring programmes. In addition to providing baseline data and potentially identifying high-risk introduction routes, the predictive power of ecological niche modelling is enhanced when species records are genetically verified.

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 worlds mallard population. These diverse approaches all support the conclusion that there is a lack of clear population structure, suggesting that the worlds mallards, perhaps with minor exceptions, form a single large, mainly interbreeding population.

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

BackgroundMain 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.ResultsPhylogenetic 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 (FST = 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).ConclusionsBased on this first complete assessment of the mallards 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.

Pan-African Genetic Structure in the African Buffalo (Syncerus caffer): Investigating Intraspecific Divergence

The African buffalo (Syncerus caffer) exhibits extreme morphological variability, which has led to controversies about the validity and taxonomic status of the various recognized subspecies. The present study aims to clarify these by inferring the pan-African spatial distribution of genetic diversity, using a comprehensive set of mitochondrial D-loop sequences from across the entire range of the species. All analyses converged on the existence of two distinct lineages, corresponding to a group encompassing West and Central African populations and a group encompassing East and Southern African populations. The former is currently assigned to two to three subspecies (S. c. nanus, S. c. brachyceros, S. c. aequinoctialis) and the latter to a separate subspecies (S. c. caffer). Forty-two per cent of the total amount of genetic diversity is explained by the between-lineage component, with one to seventeen female migrants per generation inferred as consistent with the isolation-with-migration model. The two lineages diverged between 145 000 to 449 000 years ago, with strong indications for a population expansion in both lineages, as revealed by coalescent-based analyses, summary statistics and a star-like topology of the haplotype network for the S. c. caffer lineage. A Bayesian analysis identified the most probable historical migration routes, with the Cape buffalo undertaking successive colonization events from Eastern toward Southern Africa. Furthermore, our analyses indicate that, in the West-Central African lineage, the forest ecophenotype may be a derived form of the savanna ecophenotype and not vice versa, as has previously been proposed. The African buffalo most likely expanded and diverged in the late to middle Pleistocene from an ancestral population located around the current-day Central African Republic, adapting morphologically to colonize new habitats, hence developing the variety of ecophenotypes observed today.

Reintroductions and genetic introgression from domestic pigs have shaped the genetic population structure of Northwest European wild boar

BackgroundPopulation genetic studies focus on natural dispersal and isolation by landscape barriers as the main drivers of genetic population structure. However, anthropogenic factors such as reintroductions, translocations and wild x domestic hybridization may also have strong effects on genetic population structure. In this study we genotyped 351 Single Nucleotide Polymorphism markers evenly spread across the genome in 645 wild boar (Sus scrofa) from Northwest Europe to evaluate determinants of genetic population structure.ResultsWe show that wild boar genetic population structure is influenced by historical reintroductions and by genetic introgression from domestic pigs. Six genetically distinct and geographically coherent wild boar clusters were identified in the Netherlands and Western Germany. The Dutch Veluwe cluster is known to be reintroduced, and three adjacent Dutch and German clusters are suspected to be a result of reintroduction, based on clustering results, low levels of heterozygosity and relatively high genetic distances to nearby populations. Recent wild x domestic hybrids were found geographically widespread across clusters and at low frequencies (average 3.9%). The relationship between pairwise kinship coefficients and geographic distance showed male-biased dispersal at the population genetic level.ConclusionsOur results demonstrate that wildlife and landscape management by humans are shaping the genetic diversity of an iconic wildlife species. Historical reintroductions, translocation and recent restocking activities with farmed wild boar have all influenced wild boar genetic population structure. The current trend of wild boar population growth and range expansion has recently led to a number of contact zones between clusters, and further admixture between the different wild boar clusters is to be expected.

The Development of a Genome Wide SNP Set for the Barnacle Goose Branta leucopsis

Migratory birds are of particular interest for population genetics because of the high connectivity between habitats and populations. A high degree of connectivity requires using many genetic markers to achieve the required statistical power, and a genome wide SNP set can fit this purpose. Here we present the development of a genome wide SNP set for the Barnacle Goose Branta leucopsis, a model species for the study of bird migration. We used the genome of a different waterfowl species, Mallard Anas platyrhynchos, as a reference to align Barnacle Goose second generation sequence reads from an RRL library and detected 2188 SNPs genome wide. Furthermore, we used chimeric flanking sequences, merged from both Mallard and Barnacle Goose DNA sequence information, to create primers for validation by genotyping. Validation with a 384 SNP genotyping set resulted in 374 (97%) successfully typed SNPs in the assay, of which 358 (96%) were polymorphic. Additionally, we validated our SNPs on relatively old (30 years) museum samples, which resulted in a success rate of at least 80%. This shows that museum samples could be used in standard SNP genotyping assays. Our study also shows that the genome of a related species can be used as reference to detect genome wide SNPs in birds, because genomes of birds are highly conserved. This is illustrated by the use of chimeric flanking sequences, which showed that the incorporation of flanking nucleotides from Mallard into Barnacle Goose sequences lead to equal genotyping performance when compared to flanking sequences solely composed of Barnacle Goose sequence.

Autosomal and mtDNA Markers Affirm the Distinctiveness of Lions in West and Central Africa

The evolutionary history of a species is key for understanding the taxonomy and for the design of effective management strategies for species conservation. The knowledge about the phylogenetic position of the lion (Panthera leo) in West/Central Africa is largely based on mitochondrial markers. Previous studies using mtDNA only have shown this region to hold a distinct evolutionary lineage. In addition, anthropogenic factors have led to a strong decline in West/Central African lion numbers, thus, the conservation value of these populations is particularly high. Here, we investigate whether autosomal markers are concordant with previously described phylogeographic patterns, and confirm the unique position of the West/Central African lion. Analysis of 20 microsatellites and 1,454 bp of the mitochondrial DNA in 16 lion populations representing the entire geographic range of the species found congruence in both types of markers, identifying four clusters: 1) West/Central Africa, 2) East Africa, 3) Southern Africa and 4) India. This is not in line with the current taxonomy, as defined by the IUCN, which only recognizes an African and an Asiatic subspecies. There are no indications that genetic diversity in West/Central Africa lions is lower than in either East or Southern Africa, however, given this genetic distinction and the recent declines of lion numbers in this region, we strongly recommend prioritization of conservation projects in West/Central Africa. As the current taxonomic nomenclature does not reflect the evolutionary history of the lion, we suggest that a taxonomic revision of the lion is warranted.

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.

Rainfall-driven sex-ratio genes in African buffalo suggested by correlations between Y-chromosomal haplotype frequencies and foetal sex ratio

BackgroundThe Y-chromosomal diversity in the African buffalo (Syncerus caffer) population of Kruger National Park (KNP) is characterized by rainfall-driven haplotype frequency shifts between year cohorts. Stable Y-chromosomal polymorphism is difficult to reconcile with haplotype frequency variations without assuming frequency-dependent selection or specific interactions in the population dynamics of X- and Y-chromosomal genes, since otherwise the fittest haplotype would inevitably sweep to fixation. Stable Y-chromosomal polymorphism due one of these factors only seems possible when there are Y-chromosomal distorters of an equal sex ratio, which act by negatively affecting X-gametes, or Y-chromosomal suppressors of a female-biased sex ratio. These sex-ratio (SR) genes modify (suppress) gamete transmission in their own favour at a fitness cost, allowing for stable polymorphism.ResultsHere we show temporal correlations between Y-chromosomal haplotype frequencies and foetal sex ratios in the KNP buffalo population, suggesting SR genes. Frequencies varied by a factor of five; too high to be alternatively explained by Y-chromosomal effects on pregnancy loss. Sex ratios were male-biased during wet and female-biased during dry periods (male proportion: 0.47-0.53), seasonally and annually. Both wet and dry periods were associated with a specific haplotype indicating a SR distorter and SR suppressor, respectively.ConclusionsThe distinctive properties suggested for explaining Y-chromosomal polymorphism in African buffalo may not be restricted to this species alone. SR genes may play a broader and largely overlooked role in mammalian sex-ratio variation.