Dataset Information

SNP detection and prediction of variability between chicken lines using genome resequencing of DNA pools.

ABSTRACT:

Background

Next-generation sequencing technologies are widely used for detection of millions of Single Nucleotide Polymorphisms (SNPs) and also provide a means of assessing their variation. This information is useful for composing subsets of highly informative SNPs for region-specific or genome-wide analysis and to identify mutations regulating phenotypic differences within or between populations. In this study, we investigated the sensitivity of SNP detection and introduced the flanking SNPs value (FSV) as a novel measure for predicting SNP-variability using ~5X genome resequencing with ABI SOLID and DNA pools from two chicken lines divergently selected for juvenile bodyweight.

Results

Genotyping with a 60 K SNP chip revealed polymorphisms within or between two divergently selected chicken lines for 31 363 SNPs, 48% of which were also detected using resequencing of DNA pools. SNP detection using resequencing was more powerful for positions with larger differences in allele frequency between the lines. About 50% of the SNPs with non-reference allele frequencies in the range 0.5-0.6 and 67% of those with frequencies > 0.9 could be detected. On average, ~3.7 SNPs/kb were detected by resequencing, with about 5% lower density on microchromosomes than on macrochromosomes. There was a positive correlation between the observed between-line SNP variation from the 60 K chip analysis and our proposed FSV score computed from the genome resequencing data. The strongest correlations on macrochromosomes and microchromosomes were observed when the FSV was calculated with total flanking regions of 62 kb (correlation 0.55) and 38 kb (correlation 0.45), respectively.

Conclusions

Genome resequencing with limited coverage (~5X) using pooled DNA samples and three non-reference reads as a threshold for SNP detection, identified 50 - 67% of the 60 K SNPs with a non-reference allele frequency larger than 0.5. The SNP density was around 5% lower on the microchromosomes, most likely because of their higher gene content. Our proposed method to estimate the SNP variation (FSV) uses additional sequence information to better predict SNP informativity. The FSV scores showed higher correlations for SNPs with a larger difference in allele frequency between the populations. The correlation was strongest on macrochromosomes, probably due to a lower recombination rate.

SUBMITTER: Marklund S

PROVIDER: S-EPMC3091782 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundThe Smyth line (SL) chicken is the only animal model for autoimmune vitiligo that spontaneously displays all clinical and biological manifestations of the human disorder. To understand the genetic components underlying the susceptibility to develop SL vitiligo (SLV), whole genome resequencing analysis was performed in SLV chickens compared with non-vitiliginous parental Brown line (BL) chickens, which maintain a very low incidence rate of vitiligo.ResultsIllumina sequencing technology and reference based assembly on Red Jungle Fowl genome sequences were used. Results of genome resequencing of pooled DNA of each 10 BL and SL chickens reached 5.1x and 7.0x coverage, respectively. The total number of SNPs was 4.8 and 5.5 million in BL and SL genome, respectively. Through a series of filtering processes, a total of ~1 million unique SNPs were found in the SL alone. Eventually of the 156 reliable marker SNPs, which can induce non-synonymous-, frameshift-, nonsense-, and no-start mutations in amino acid sequences in proteins, 139 genes were chosen for further analysis. Of these, 14 randomly chosen SNPs were examined for SNP verification by PCR and Sanger sequencing to detect SNP positions in 20 BL and 70 SL chickens. The results of the analysis of the 14 SNPs clearly showed differential frequencies of nucleotide bases in the SNP positions between BL and SL chickens. Bioinformatic analysis showed that the 156 most reliable marker SNPs included genes involved in dermatological diseases/conditions such as ADAMTS13, ASPM, ATP6V0A2, BRCA2, COL12A1, GRM5, LRP2, OBSCN, PLAU, RNF168, STAB2, and XIRP1. Intermolecular gene network analysis revealed that candidate genes identified in SLV play a role in networks centered on protein kinases (MAPK, ERK1/2, PKC, PRKDC), phosphatase (PPP1CA), ubiquitinylation (UBC) and amyloid production (APP).ConclusionsVarious potential genetic markers showing amino acid changes and potential roles in vitiligo development were identified in the SLV chicken through genome resequencing. The genetic markers and bioinformatic interpretations of amino acid mutations found in SLV chickens may provide insight into the genetic component responsible for the onset and the progression of autoimmune vitiligo and serve as valuable markers to develop diagnostic tools to detect vitiligo susceptibility.

Project description:DNA cytosine methylation is an important epigenetic modification that has significant effects on a variety of biological processes in animals. Avian species hold a crucial position in evolutionary history. In this study, we used whole-genome bisulfite sequencing (MethylC-seq) to generate single base methylation profiles of lungs in two genetically distinct and highly inbred chicken lines (Fayoumi and Leghorn) that differ in genetic resistance to multiple pathogens, and we explored the potential regulatory role of DNA methylation associated with immune response differences between the two chicken lines.The MethylC-seq was used to generate single base DNA methylation profiles of Fayoumi and Leghorn birds. In addition, transcriptome profiling using RNA-seq from the same chickens and tissues were obtained to interrogate how DNA methylation regulates gene transcription on a genome-wide scale.The general DNA methylation pattern across different regions of genes was conserved compared to other species except for hyper-methylation of repeat elements, which was not observed in chicken. The methylation level of miRNA and pseudogene promoters was high, which indicates that silencing of these genes may be partially due to promoter hyper-methylation. Interestingly, the promoter regions of more recently evolved genes tended to be more highly methylated, whereas the gene body regions of evolutionarily conserved genes were more highly methylated than those of more recently evolved genes. Immune-related GO (Gene Ontology) terms were significantly enriched from genes within the differentially methylated regions (DMR) between Fayoumi and Leghorn, which implicates DNA methylation as one of the regulatory mechanisms modulating immune response differences between these lines.This study establishes a single-base resolution DNA methylation profile of chicken lung and suggests a regulatory role of DNA methylation in controlling gene expression and maintaining genome transcription stability. Furthermore, profiling the DNA methylomes of two genetic lines that differ in disease resistance provides a unique opportunity to investigate the potential role of DNA methylation in host disease resistance. Our study provides a foundation for future studies on epigenetic modulation of host immune response to pathogens in chickens.

Project description:Genotyping-by-sequencing (GBS) represents a highly cost-effective high-throughput genotyping approach. By nature, however, GBS is subject to generating sizeable amounts of missing data and these will need to be imputed for many downstream analyses. The extent to which such missing data can be tolerated in calling SNPs has not been explored widely. In this work, we first explore the use of imputation to fill in missing genotypes in GBS datasets. Importantly, we use whole genome resequencing data to assess the accuracy of the imputed data. Using a panel of 301 soybean accessions, we show that over 62,000 SNPs could be called when tolerating up to 80% missing data, a five-fold increase over the number called when tolerating up to 20% missing data. At all levels of missing data examined (between 20% and 80%), the resulting SNP datasets were of uniformly high accuracy (96-98%). We then used imputation to combine complementary SNP datasets derived from GBS and a SNP array (SoySNP50K). We thus produced an enhanced dataset of >100,000 SNPs and the genotypes at the previously untyped loci were again imputed with a high level of accuracy (95%). Of the >4,000,000 SNPs identified through resequencing 23 accessions (among the 301 used in the GBS analysis), 1.4 million tag SNPs were used as a reference to impute this large set of SNPs on the entire panel of 301 accessions. These previously untyped loci could be imputed with around 90% accuracy. Finally, we used the 100K SNP dataset (GBS + SoySNP50K) to perform a GWAS on seed oil content within this collection of soybean accessions. Both the number of significant marker-trait associations and the peak significance levels were improved considerably using this enhanced catalog of SNPs relative to a smaller catalog resulting from GBS alone at ?20% missing data. Our results demonstrate that imputation can be used to fill in both missing genotypes and untyped loci with very high accuracy and that this leads to more powerful genetic analyses.

Project description:BACKGROUND: The chicken (Gallus gallus) is an important model organism that bridges the evolutionary gap between mammals and other vertebrates. Copy number variations (CNVs) are a form of genomic structural variation widely distributed in the genome. CNV analysis has recently gained greater attention and momentum, as the identification of CNVs can contribute to a better understanding of traits important to both humans and other animals. To detect chicken CNVs, we genotyped 475 animals derived from two broiler chicken lines divergently selected for abdominal fat content using chicken 60 K SNP array, which is a high-throughput method widely used in chicken genomics studies. RESULTS: Using PennCNV algorithm, we detected 438 and 291 CNVs in the lean and fat lines, respectively, corresponding to 271 and 188 CNV regions (CNVRs), which were obtained by merging overlapping CNVs. Out of these CNVRs, 99% were confirmed also by the CNVPartition program. These CNVRs covered 40.26 and 30.60 Mb of the chicken genome in the lean and fat lines, respectively. Moreover, CNVRs included 176 loss, 68 gain and 27 both (i.e. loss and gain within the same region) events in the lean line, and 143 loss, 25 gain and 20 both events in the fat line. Ten CNVRs were chosen for the validation experiment using qPCR method, and all of them were confirmed in at least one qPCR assay. We found a total of 886 genes located within these CNVRs, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed they could play various roles in a number of biological processes. Integrating the results of CNVRs, known quantitative trait loci (QTL) and selective sweeps for abdominal fat content suggested that some genes (including SLC9A3, GNAL, SPOCK3, ANXA10, HELIOS, MYLK, CCDC14, SPAG9, SOX5, VSNL1, SMC6, GEN1, MSGN1 and ZPAX) may be important for abdominal fat deposition in the chicken. CONCLUSIONS: Our study provided a genome-wide CNVR map of the chicken genome, thereby contributing to our understanding of genomic structural variations and their potential roles in abdominal fat content in the chicken.