Frances Turner

Development of a high density 600K SNP genotyping array for chicken

BackgroundHigh density (HD) SNP genotyping arrays are an important tool for genetic analyses of animals and plants. Although the chicken is one of the most important farm animals, no HD array is yet available for high resolution genetic analysis of this species.ResultsWe report here the development of a 600 K Affymetrix® Axiom® HD genotyping array designed using SNPs segregating in a wide variety of chicken populations. In order to generate a large catalogue of segregating SNPs, we re-sequenced 243 chickens from 24 chicken lines derived from diverse sources (experimental, commercial broiler and layer lines) by pooling 10–15 samples within each line. About 139 million (M) putative SNPs were detected by mapping sequence reads to the new reference genome (Gallus_gallus_4.0) of which ~78 M appeared to be segregating in different lines. Using criteria such as high SNP-quality score, acceptable design scores predicting high conversion performance in the final array and uniformity of distribution across the genome, we selected ~1.8 M SNPs for validation through genotyping on an independent set of samples (n = 282). About 64% of the SNPs were polymorphic with high call rates (>98%), good cluster separation and stable Mendelian inheritance. Polymorphic SNPs were further analysed for their population characteristics and genomic effects. SNPs with extreme breach of Hardy-Weinberg equilibrium (P < 0.00001) were excluded from the panel. The final array, designed on the basis of these analyses, consists of 580,954 SNPs and includes 21,534 coding variants. SNPs were selected to achieve an essentially uniform distribution based on genetic map distance for both broiler and layer lines. Due to a lower extent of LD in broilers compared to layers, as reported in previous studies, the ratio of broiler and layer SNPs in the array was kept as 3:2. The final panel was shown to genotype a wide range of samples including broilers and layers with over 100 K to 450 K informative SNPs per line. A principal component analysis was used to demonstrate the ability of the array to detect the expected population structure which is an important pre-investigation step for many genome-wide analyses.ConclusionsThis Affymetrix® Axiom® array is the first SNP genotyping array for chicken that has been made commercially available to the public as a product. This array is expected to find widespread usage both in research and commercial application such as in genomic selection, genome-wide association studies, selection signature analyses, fine mapping of QTLs and detection of copy number variants.

A High Resolution Genome-Wide Scan for Significant Selective Sweeps: An Application to Pooled Sequence Data in Laying Chickens

In most studies aimed at localizing footprints of past selection, outliers at tails of the empirical distribution of a given test statistic are assumed to reflect locus-specific selective forces. Significance cutoffs are subjectively determined, rather than being related to a clear set of hypotheses. Here, we define an empirical p-value for the summary statistic by means of a permutation method that uses the observed SNP structure in the real data. To illustrate the methodology, we applied our approach to a panel of 2.9 million autosomal SNPs identified from re-sequencing a pool of 15 individuals from a brown egg layer line. We scanned the genome for local reductions in heterozygosity, suggestive of selective sweeps. We also employed a modified sliding window approach that accounts for gaps in the sequence and increases scanning resolution by moving the overlapping windows by steps of one SNP only, and suggest to call this a “creeping window” strategy. The approach confirmed selective sweeps in the region of previously described candidate genes, i.e. TSHR, PRL, PRLHR, INSR, LEPR, IGF1, and NRAMP1 when used as positive controls. The genome scan revealed 82 distinct regions with strong evidence of selection (genome-wide p-value<0.001), including genes known to be associated with eggshell structure and immune system such as CALB1 and GAL cluster, respectively. A substantial proportion of signals was found in poor gene content regions including the most extreme signal on chromosome 1. The observation of multiple signals in a highly selected layer line of chicken is consistent with the hypothesis that egg production is a complex trait controlled by many genes.

Design and development of exome capture sequencing for the domestic pig (Sus scrofa)

BackgroundThe domestic pig (Sus scrofa) is both an important livestock species and a model for biomedical research. Exome sequencing has accelerated identification of protein-coding variants underlying phenotypic traits in human and mouse. We aimed to develop and validate a similar resource for the pig.ResultsWe developed probe sets to capture pig exonic sequences based upon the current Ensembl pig gene annotation supplemented with mapped expressed sequence tags (ESTs) and demonstrated proof-of-principle capture and sequencing of the pig exome in 96 pigs, encompassing 24 capture experiments. For most of the samples at least 10x sequence coverage was achieved for more than 90% of the target bases. Bioinformatic analysis of the data revealed over 236,000 high confidence predicted SNPs and over 28,000 predicted indels.ConclusionsWe have achieved coverage statistics similar to those seen with commercially available human and mouse exome kits. Exome capture in pigs provides a tool to identify coding region variation associated with production traits, including loss of function mutations which may explain embryonic and neonatal losses, and to improve genomic assemblies in the vicinity of protein coding genes in the pig.

Genomic variation in macrophage-cultured European porcine reproductive and respiratory syndrome virus Olot/91 revealed using ultra-deep next generation sequencing

BackgroundPorcine Reproductive and Respiratory Syndrome (PRRS) is a disease of major economic impact worldwide. The etiologic agent of this disease is the PRRS virus (PRRSV). Increasing evidence suggest that microevolution within a coexisting quasispecies population can give rise to high sequence heterogeneity in PRRSV.FindingsWe developed a pipeline based on the ultra-deep next generation sequencing approach to first construct the complete genome of a European PRRSV, strain Olot/9, cultured on macrophages and then capture the rare variants representative of the mixed quasispecies population. Olot/91 differs from the reference Lelystad strain by about 5% and a total of 88 variants, with frequencies as low as 1%, were detected in the mixed population. These variants included 16 non-synonymous variants concentrated in the genes encoding structural and nonstructural proteins; including Glycoprotein 2a and 5.ConclusionUsing an ultra-deep sequencing methodology, the complete genome of Olot/91 was constructed without any prior knowledge of the sequence. Rare variants that constitute minor fractions of the heterogeneous PRRSV population could successfully be detected to allow further exploration of microevolutionary events.

Development of a Medium Density Combined-Species SNP Array for Pacific and European Oysters (Crassostrea gigas and Ostrea edulis)

SNP arrays are enabling tools for high-resolution studies of the genetic basis of complex traits in farmed and wild animals. Oysters are of critical importance in many regions from both an ecological and economic perspective, and oyster aquaculture forms a key component of global food security. The aim of our study was to design a combined-species, medium density SNP array for Pacific oyster (Crassostrea gigas) and European flat oyster (Ostrea edulis), and to test the performance of this array on farmed and wild populations from multiple locations, with a focus on European populations. SNP discovery was carried out by whole-genome sequencing (WGS) of pooled genomic DNA samples from eight C. gigas populations, and restriction site-associated DNA sequencing (RAD-Seq) of 11 geographically diverse O. edulis populations. Nearly 12 million candidate SNPs were discovered and filtered based on several criteria, including preference for SNPs segregating in multiple populations and SNPs with monomorphic flanking regions. An Affymetrix Axiom Custom Array was created and tested on a diverse set of samples (n = 219) showing ∼27 K high quality SNPs for C. gigas and ∼11 K high quality SNPs for O. edulis segregating in these populations. A high proportion of SNPs were segregating in each of the populations, and the array was used to detect population structure and levels of linkage disequilibrium (LD). Further testing of the array on three C. gigas nuclear families (n = 165) revealed that the array can be used to clearly distinguish between both families based on identity-by-state (IBS) clustering parental assignment software. This medium density, combined-species array will be publicly available through Affymetrix, and will be applied for genome-wide association and evolutionary genetic studies, and for genomic selection in oyster breeding programs.

Quantifying the Survival of Multiple Salmonella enterica Serovars In Vivo via Massively Parallel Whole-Genome Sequencing To Predict Zoonotic Risk

ABSTRACT Salmonella enterica is an animal and zoonotic pathogen of worldwide importance. Salmonella serovars that differ in their host and tissue tropisms exist. Cattle are an important reservoir of human nontyphoidal salmonellosis, and contaminated bovine peripheral lymph nodes enter the food chain via ground beef. The relative abilities of different serovars to survive within the bovine lymphatic system are poorly understood and constrain the development of control strategies. This problem was addressed by developing a massively parallel whole-genome sequencing method to study mixed-serovar infections in vivo. Salmonella serovars differ genetically by naturally occurring single nucleotide polymorphisms (SNPs) in certain genes. It was hypothesized that these SNPs could be used as markers to simultaneously identify serovars in mixed populations and quantify the abundance of each member in a population. The performance of the method was validated in vitro using simulated pools containing up to 11 serovars in various proportions. It was then applied to study serovar survival in vivo in cattle challenged orally with the same 11 serovars. All the serovars successfully colonized the bovine lymphatic system, including the peripheral lymph nodes, and thus pose similar risks of zoonosis. This method enables the fates of multiple genetically unmodified strains to be evaluated simultaneously in a single animal. It could be useful in reducing the number of animals required to study mixed-strain infections and in testing the cross-protective efficacy of vaccines and treatments. It also has the potential to be applied to diverse bacterial species which possess shared but polymorphic alleles. IMPORTANCE While some Salmonella serovars are more frequently isolated from lymph nodes rather than the feces and environment of cattle, the relative abilities of serovars to survive within the lymphatic system of cattle remain ill defined. A sequencing-based method which used available information from sequenced Salmonella genomes to study the dynamics of mixed-serovar infections in vivo was developed. The main advantages of the method include the simultaneous identification and quantification of multiple strains without any genetic modification and minimal animal use. This approach could be used in vaccination trials or in epidemiological surveys where an understanding of the dynamics of closely related strains of a pathogen in mixed populations could inform the prediction of zoonotic risk and the development of intervention strategies.

Genome-wide disruption of DNA methylation by 5-aza-2'-deoxycytidine in the parasitioid wasp Nasonia vitripennis

DNA methylation of cytosine residues across the genome influences how many genes and phenotypes are regulated. As such, understanding the role of DNA methylation and other epigenetic mechanisms has become very much a part of mapping genotype to phenotype, a major question in evolutionary biology. Ideally, we would like to manipulate DNA methylation patterns on a genome-wide scale, to elucidate the role of epigenetic modifications in phenotypic expression. Recently, the demethylating agent 5-aza-2’-deoxycytidine (5-aza-dC; commonly used in the epigenetic treatment of certain cancers), has been deployed to explore the epigenetic regulation of a number of traits of interest to evolutionary ecologists. Recently, we showed that treatment with 5-aza-dC shifted patterns of sex allocation as predicted by genomic conflict theory in the parasitoid wasp Nasonia vitripennis. This was the first (albeit indirect) experimental evidence for genomic conflict over sex allocation facilitated by DNA methylation. However, this work lacked confirmation of the effects of 5-aza-dC on DNA methylation, drawing commentary on the efficacy of 5-aza-dC in a novel system. Here, using whole-genome bisulphite sequencing, we demonstrate unequivocally that 5-aza-dC disrupts methylation across the Nasonia vitripennis genome. We show that disruption leads to both hypo- and hyper-methylation, may vary across tissues and time of sampling, and that the effects of 5-aza-dC are context- and sequence specific. We conclude that 5-aza-dC has the potential to be repurposed as a tool in evolutionary ecology for studying the role of DNA methylation.

Maintaining their genetic distance; limited gene flow between widely hybridising species of Geum with contrasting mating systems

Mating system transition from outcrossing to selfing frequently gives rise to sister lineages with contrasting outcrossing rates. The evolutionary fate of such lineages depends on the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing G. urbanum, two sister species derived by mating system transition, which frequently hybridise. A draft genome was generated for G. urbanum and used to develop dd-RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the two species diverged 2-3 Mya, and that long term gene flow between them has been very low (M=0.04). G. rivale showed greater genetic diversity in sympatry than allopatry, but genetic divergence between species was no lower in sympatry than allopatry, providing little evidence for recent introgression. Clustering of genotypes revealed that, apart from four early generation hybrids, individuals in sympatric populations fell into two genetically distinct groups with <1% admixture that corresponded exactly to their morphological species classification. Although our data suggest limited gene flow, we observed joint segregation of two putatively introgressed SNPs in G. urbanum populations that was associated with significant morphological variation; this provides tentative evidence for rare introduction of novel genetic diversity by interspecific gene flow. Our results indicate that despite frequent hybridisation, genetic exchange between G. rivale and G. urbanum has been very limited throughout their evolutionary history.

Maintaining their genetic distance: little evidence for introgression between widely hybridising species of Geum with contrasting mating systems

Abstract Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.