Detection and characterization of small insertion and deletion genetic variants in modern layer chicken genomes
© Boschiero et al. 2015
Received: 17 December 2014
Accepted: 22 June 2015
Published: 31 July 2015
Small insertions and deletions (InDels) constitute the second most abundant class of genetic variants and have been found to be associated with many traits and diseases. The present study reports on the detection and characterisation of about 883 K high quality InDels from the whole-genome analysis of several modern layer chicken lines from diverse breeds.
To reduce the error rates seen in InDel detection, this study used the consensus set from two InDel-calling packages: SAMtools and Dindel, as well as stringent post-filtering criteria. By analysing sequence data from 163 chickens from 11 commercial and 5 experimental layer lines, this study detected about 883 K high quality consensus InDels with 93 % validation rate and an average density of 0.78 InDels/kb over the genome. Certain chromosomes, viz, GGAZ, 16, 22 and 25 showed very low densities of InDels whereas the highest rate was observed on GGA6. In spite of the higher recombination rates on microchromosomes, the InDel density on these chromosomes was generally lower relative to macrochromosomes possibly due to their higher gene density. About 43–87 % of the InDels were found to be fixed within each line. The majority of detected InDels (86 %) were 1–5 bases and about 63 % were non-repetitive in nature while the rest were tandem repeats of various motif types. Functional annotation identified 613 frameshift, 465 non-frameshift and 10 stop-gain/loss InDels. Apart from the frameshift and stopgain/loss InDels that are expected to affect the translation of protein sequences and their biological activity, 33 % of the non-frameshift were predicted as evolutionary intolerant with potential impact on protein functions. Moreover, about 2.5 % of the InDels coincided with the most-conserved elements previously mapped on the chicken genome and are likely to define functional elements. InDels potentially affecting protein function were found to be enriched for certain gene-classes e.g. those associated with cell proliferation, chromosome and Golgi organization, spermatogenesis, and muscle contraction.
The large catalogue of InDels presented in this study along with their associated information such as functional annotation, estimated allele frequency, etc. are expected to serve as a rich resource for application in future research and breeding in the chicken.
Small insertions and deletions (InDels) are the second most abundant kind of genetic variants in the genome after single nucleotide polymorphisms (SNPs) and are the most common type of structural variants (SV) [1, 2]. InDels have been implicated in many diseases and other traits. For instance, in human approximately one quarter of all known Mendelian diseases are associated with InDels . In other organisms as well, such as in chicken, InDels have been found to be associated with different phenotypes e.g. growth , plumage colour , egg production , performance , body weight  and retinal degeneration and embryonic mortality . Understanding InDels in greater detail is therefore important for profiling genetic variations within genomes , detecting causal mutations of genetic disorders [11, 12], studying the evolutionary relationship of species  and detecting footprints of selection [14, 15].
In spite of their importance and the massive advancements in high-throughput next-generation sequencing (NGS) technologies, the discovery of InDels has lagged behind that of SNPs for a number of reasons: (i) it is generally more difficult to map short reads by NGS covering InDel sites to the correct locations since this involves more complex gapped alignment or paired-end sequencing inference [16, 17]; (ii) since the read mappers align each fragment independently of other fragments, the InDels may appear as stretches of SNPs rather than as gaps ; (iii) the majority of InDels occur as short tandem repeats, which are difficult to map [19, 20]; (iv) higher coverage is required to detect InDels due to their relatively low frequency ; and (v) distinguishing true InDels from sequencing errors is difficult since there is no accurate sequencing error model for InDel . InDel calling, therefore, requires a suitable genome aligner that can perform gapped alignment and software such as BWA [22, 23] and NOVOalign [23, 24] have been used for this purpose. Moreover, to reduce the problems of misalignments, local realignment of reads with InDels have been suggested and to minimise the detection of false positives due to other factors, stringent criteria for post-alignment filtrations needs to be applied .
Unlike SNPs, which are generated through point mutations, SVs (including InDels) can arise through a number of mechanisms, such as replication slippage, recombination, unequal crossing over and tandem duplications caused by imperfect repair of double-strand breaks [26–29]. As a consequence, InDel sizes vary widely from as small as 1 base to over several kilobases (kb) [10, 12, 30, 31]. This variation in size means that robust detection of SVs requires the use of a range of detection methods. Small SVs, like InDels, are frequently detected by mapping the small NGS sequencing reads against the reference genome. However, the size of the reads dictates the maximum size of InDels that can be detected by this method. In the 1000 Genomes Project Consortium , for example, this approach was used to detect InDels with a size range of 1–50 nucleotides.
In the present study we aimed to detect small InDels (less than 100 nucleotides) in the chicken genome, which is an important farm animal and a key model organism for genomic and developmental biology studies. Although a number of studies [14, 33–35] have focused on detecting SNPs from the chicken resulting in the discovery of millions of these variants, only a few studies [36–38] have analysed InDels. As a result the number of InDels reported in the public databases is quite low; for example dbSNP (build 140) reports only about 439 K InDels from the chicken genome. Only in a recent study, Yan et al.  have reported about 1.3 M InDels by analysing 12 diverse chicken lines. While this number has been a great contribution to the InDel database, study by Yan et al. analysed only single bird from each population and as a result could not shed light on certain aspects such as whether these InDels were segregating or fixed within populations, their allele frequencies, etc.
We used NGS sequence data generated from 163 chickens from 11 commercial and 5 experimental layer lines (samples pooled within lines)  to define our InDel set. To reduce the error rates often seen in InDel detection, we took only the consensus set called by two software packages, namely SAMtools  and Dindel , followed by stringent filtration steps. The approach of taking consensus variants from multiple callers has been used in many recent studies on different species [38, 40–42] . Although several InDel callers are now available, we used these two for a number of reasons. First, a recent study found that Dindel and SAMtools mpileup have the highest sensitivity in InDel calling at low coverage (less than 30X) compared to the other callers such as VarScan and GATK . In our study the sequence coverage ranged between 7 and 17X for different populations and hence the choice of these packages appeared reasonable. Second, SAMtools is one of the most commonly used tools for detection of variants due to its simple workflow with many advanced features, such as its ability to perform local realignment . Nevertheless, while the package has been reported to have a low false discovery rate (FDR) for SNP calling [44, 45], for InDel calling it was reported to have a rather high rate of 4.8 % on real data . The Dindel package, on the other hand, has been modelled specifically for InDel calling and was found to have a much lower FDR (1.6 % on real data) compared to SAMtools [18, 16] but it has a complex workflow and long running time. Our study reports the discovery of about 883 K high quality consensus InDels using these two packages followed by stringent filtration criteria, and discusses the physical and functional characteristics of these InDels.
InDel calling using SAMtools and Dindel
Both the SAMtools and Dindel packages use Bayesian models for the detection of InDels. SAMtools takes into account mapping qualities, base qualities and error rates (from raw sequence quality scores) in approximating the posterior probabilities of consensus genotypes, which are then used for calling the variants . The SAMtools also performs local realignment around candidate InDel sites . The basic principle underlying the Dindel package is to realign all reads mapped to a region to a number of candidate haplotypes of at least 120-bp length, which represent alternatives to the reference genome sequence . The probabilistic realignment model of Dindel takes into account base qualities, map qualities of reads, insert size distribution and position-dependent rates of InDel error based on homopolymer run length. Upon realignment of the reads, InDels are called by comparing the posterior probabilities of pairs of haplotypes with and without InDels. Both tools provide Phred-based quality scores to the called variants that can be used for further filtration.
Even though the quality scores for individual variants were different from the two packages and the scores from SAMtools were 16–34 % higher compared to those from Dindel (the average quality scores were 108 and 83, respectively), there was a significant correlation between the scores generated by the two packages (r ≈ 0.36; P < 0.001).
Some InDels were detected exclusively by either SAMtools (n = 162,851) or Dindel (n = 245,005). One source of difference between these non-consensus sets was the size of the InDels. SAMtools was able to detect larger InDels and the maximum length detected was 81 bp, whereas the largest InDel detected by Dindel was 56 bp. Another difference was that many variants detected exclusively by Dindel had low quality scores, possibly originating from calling variants associated with poor base or map qualities for which no minimum thresholds were applied as explained above. Moreover, some InDels detected exclusively by SAMtools had low quality, and also low depth of coverage.
Stringent filtration of consensus InDels defines a high quality set
In order to improve the confidence in the detected InDels, only the consensus set from the two packages were chosen for further filtration and analyses. Among all the variants detected by the two packages, about 74 % (n = 1,110,149) were called by both methods and hence were considered consensus. About 11 % of the SAMtools detected variants and 16 % of the Dindel variants failed to be within the consensus set and were excluded from further analyses.
Even after obtaining the consensus variants, further filtrations were applied to extract a catalogue of high quality InDels. Major filtration criteria were: InDel quality ≥30, coverage ≥5 and ≤ mean coverage in a line + 3 SD, non-reference allele supported by both forward and reverse strands, and gap between consecutive InDels >1 base. This resulted in the retention of 883,411 variants (about 80 % of the consensus set) including 397,438 insertions, 476,793 deletions and 9,180 block substitutions (i.e. substitution of a stretch of consecutive nucleotides with a new block). The number of consensus filtered (CF) InDels varied widely across the chicken lines (Fig. 1; Additional file 1) with the minimum number detected from the inbred line, I5 (n = 171,680) and the maximum number from the white egg layer line, WEL6 (n = 391,796). The average number of InDels detected from brown egg layer (BEL) lines however, was higher (410 K ± 28 K), compared to those from WELs (395 K ±55 K) and inbred lines (357 K ±26 K). This InDel diversity was found to be highly correlated (r = 0.78, P < 0.001) with SNP diversity from these lines presented in a recent paper by Gheyas et al. .
A large proportion of our CF InDels were fixed within lines for the non-reference alleles (allele frequency ≥0.9) and the proportion of these InDel varied widely across the chicken lines ranging from 43 to 87 %. In general, the inbred lines showed the greatest proportion of fixed variants (on average about 82 %) followed by WEL (average 62 %) and BEL lines (average 47 %). The greater proportion of fixed variants in the inbred lines was expected as these lines have been developed by many generations of sib-mating for experimental purposes.
InDel validation and false discovery rates
False discovery rates (FDR) for the InDels identified by next generation sequencing in 24 validation regions
FP rate (1-specificity)
FN rate (1-sensitivity)
Physical properties of consensus filtered InDels
Apart from length, we also classified the insertions (n = 397,438) and deletions (n = 476,793) based on whether they consisted of non-repeat or repeat expansions. About 65 % of the insertions and 61 % of the deletions were non-repeat in nature consisting of single (47 % insertions and 38 % deletions), two (11 % insertions and 14 % deletions), three (3 % insertions and 5 % deletions) or four or more nucleotides (3 % of both the insertions and deletions). The rest of the InDels consisted of either perfect or imperfect repeat motifs. Those with perfect repeats included: monomeric (14 % for insertions and 12 % for deletions), dimeric (0.5 % for both insertions and deletions) and trimeric motifs (0.2 % for insertions and 0.1 % for deletions). About 21 % of the insertions and 27 % of the deletions consisted of imperfect repeats consisting of combinations of multiple motifs of different order. We did not remove the InDels within tandem repeat sequences as they constituted a substantial proportion.
Functional annotation of InDels and their predicted biological effects
Annotation of genetic variants against functional elements in the genome (e.g. genes, or non-coding functional elements) is a major step towards their characterisation. In the following sub-sections we describe the annotation of the CF InDels against known genes, predict the effects of coding variants, explore the distribution pattern of the coding InDels in relation of polypeptide chains and also annotate the InDels against most conserved elements (MCEs).
InDels within coding and non-coding regions
Summary of annotation of consensus filtered InDels based on Ensembl gene annotations (release 71)
Total number of InDels
1 kb downstream
1 kb upstream
Non-coding RNA (ncRNA)
Effects of non-frameshift InDels on protein function
Relative positions of coding InDels within polypeptide chains
InDels within conserved elements in the genome
Allele frequency pattern of different InDel categories
In the Additional file 3 we present a heatmap using the mean AAF of putative functional InDels (viz. the variants categorized as exonic, splicing, ncRNA and those coinciding with MCEs) that are fixed (AAF ≥ 0.9) in at least one of the chicken groups. Out of 14,033 InDels used for creating the heatmap, about 94 % overlapped with one or more QTLs specified in chicken QTLdb (http://www.animalgenome.org/cgi-bin/QTLdb/GG/index). For instance, about 80 % of these InDels coincided with egg-related traits (e.g. egg production, shell quality, age and weight at first egg, etc.), 56 % coincided with bone-related traits (e.g. mineral content, strength, weight, etc.), 36 % overlapped with antibody response to various diseases and 25 % overlapped with feed conversion ratio or feed efficiency traits. About 16 % of the InDels presented in the heatmap were fixed in all of the three groups and most probably represented old variants. On the other hand, there were variants that were fixed in only a specific group while segregating at much lower frequency (e.g. ≤0.5) or not detected at all in other groups; their proportion varied with 7.3 % in BEL, 10.3 % in WEL and 8.4 % in inbred lines. These group-specific variants may be associated with traits relevant to that particular chicken group.
Enrichment of functional classes of genes with potential loss-of-function mutations
In order to investigate if the potentially loss-of-function (LOF) mutations (viz. the FS, stop-gain/loss, and intolerant-NFS InDels) are enriched within certain groups of genes with functional similarity, we performed a clustered analysis of the genes using the DAVID Gene Functional Classification tool . Performing the analyses separately on the three chicken groups showed enrichment of different classes of genes. We selected only those classes with enrichment scores (ES) ≥1.3 and genes with kappa score ≥0.75 (as suggested by DAVID) that indicate strong agreement between genes. This resulted in the retention of one class from each group. The class from the WEL group was related to cell proliferation, chromosome and Golgi organization, and epithelial cell adhesion. The gene class from the BEL group was related to cell differentiation, spermatogenesis, muscle contraction and blood coagulation and the class from the inbred group was related to cell proliferation and respiration. These classes consisted of 12 genes in total (Additional file 4). Five of these genes were shared by at least two groups: the gene STK31 was common to all three groups; the genes - SMCHD1, TBRG4 and DDX31 were common between the inbred and the WEL groups; and the gene, DCLK3 was shared by the WEL and BEL groups. The rest of the genes were present exclusively in one particular group: the genes AQR and KIFC3 in the WEL group, and the genes MAK, PBK, SRMS, PRKCH and DES in the BEL group.
We observed that one FS InDel within significant gene classes was identified from multiple lines and was present at very high frequency. This variant (insertion of a single T at the position GGA2:46501278) was located within the first exon of the gene DCLK3 (doublecortin gene family) and was detected from all the eleven lines from the WEL and BEL groups and reached near fixation, indicating possible selective advantage for these chicken groups. On the other hand, three InDels from the significant gene classes were detected at very high frequency (AAF = 1) only within single lines with the possibility to have functional implications within those specific lines. These were: a FS insertion of a single C (at the position GGA11:483395) in the exon 9 of the KIFC3 gene from a single WEL line, an intolerant NFS deletion at the position GGA5: 31527889 from the AQR (or IBP160) gene from another WEL line, and a FS insertion of a single C (at the position GGA7:21670698) in the exon 7 of the DES (desmin) gene from a BEL line.
This study aimed to detect and characterise high quality InDels from the chicken genome by screening multiple commercial and experimental layer chicken lines. A major strength of our study was the use of multiple birds (10–15) from each of the populations analysed, which allowed better characterisation of the detected variants e.g. provided the opportunity to explore the allele-frequency patterns of InDels. Although in a recent study Yan et al.  reported about 1.3 M InDels, they only used single bird from each of the 12 breeds analysed and did not investigate to any depth the functional characteristics of the InDels. For instance, they made no attempt to predict the functional effects, if any, of non-frameshift InDels, or to estimate allele frequency patterns or to find InDels overlapping with conserved elements. Our study provides a more comprehensive functional characterization of the InDels using most of the currently available resources and genomic databases on chicken. This study is also the first detailed InDel characterization in any avian species, and hence can serve as an important resource for other birds.
Since InDel calling inherently suffers from multiple problems, we adopted a number of approaches to ensure a high degree of fidelity in our calls including taking only the consensus variants detected by two different callers, and applying a number of stringent filtration criteria. This approach was highly successful as we observed a very high rate of validation (93 %). This validation rate is better than that reported in the similar recent study (88 %) by Yan et al. . The number of InDels used for validation in our study, however, was small and as a result, the rates may vary on larger dataset.
Some of the shortfalls in our study originated from the use of a pooled-sequencing approach. Although the major strength of pooled sequencing is that it facilitates the screening of many individuals within limited budget and time, it has certain drawbacks. One major drawback is that this approach is prone to miss rare and low frequency variants when sequencing is performed at low coverage [52, 53]. This was observed in our study too by the large estimate of false negative rate (26 %) and by the absence of rare and low frequency variants in allele frequency graphs. The sequence coverage in our study was low- only 7-17X for pooled samples of 10–15 individuals per line. This coverage however, dropped further during downstream analyses when only good quality reads were recruited for variant calling. This factor was further compounded with the use of stringent filtration criteria, especially, the use of at least five reads covering an InDel site and the support of non-reference allele by at least one read from both the strands. Apart from these factors, it is generally more difficult to detect heterozygous InDels compared to homozygous ones when the sequencing coverage is low [54–56]. According to an estimate based on single sample sequencing, it would require coverage of at least 20X to detect 99 % of the heterozygous variants . These details explain why in our study we have seen mostly high frequency InDels in all groups of chickens and in all annotation categories. However, we concede that the allele frequency estimation from pooled sequencing may not be accurate as it is impossible to assess if all the samples had equal contribution in the reads covering InDel sites. Besides the frequency estimates are likely to be upwardly biased due to use of only 10–15 individuals per line for the estimations. In spite of all these issues, we argue that we were able to shed some light on the frequency spectrum of the InDels and the findings suggest that we have detected mostly common variants. We would also like to emphasize that it is the common variants that are most likely to have major implications for poultry breeding.
The main findings observed about our CF InDels were: (1) majority were small in size (1–5 nucleotides), (2) a substantial proportion was located within tandem repeats, (3) large proportions were fixed within lines, (4) the InDel densities varied widely across chromosomes, (5) higher number of frameshift than non-frameshift variants were detected but their ratio was much lower compared to that of non-triplet and triplet InDels in the non-coding regions and (6) FS mutations were located more frequently either near the N- or C-terminal part of the proteins compared with NFS mutations. Most of these observations are consistent with the findings of other studies. For instance, previous studies on chicken have also found majority of InDels to be within small size range of 1–10 nucleotides [33, 36, 38]. Similarly, several studies in humans and chickens have reported excess of InDels to be present within tandem duplicates [20, 36]. The proportion of fixed InDels (43 to 87 %) within lines is a function of level of inbreeding, although the large rate of fixation observed in our study is most probably inflated due to the use of limited number of individuals and also due to the difficulty in detecting heterozygous InDels as discussed above. During our InDel call we removed those InDels which were fixed in all or most of the lines (see Materials and Methods) to ensure that these are not errors in reference genome.
The present study revealed that micro-chromosomes had significantly lower density of InDels compared with macro- and intermediate sized chromosomes. This is contrary to the expectation based on the higher recombination rates in microchromosomes [57, 58] as studies on other species have revealed that recombination rate is positively correlated with polymorphism rate [59, 60]. When we compared the InDel densities with the SNP densities in chromosomes - using the SNPs described by Kranis et al. - we observed a different picture; the micro- and intermediate chromosomes had significantly higher (P < 0.05) densities of SNPs compared to that in macrochromosomes. This indicates that even though SNP density is affected by recombination rates of the chromosomes, for InDel density other factors, such as selection pressure probably play a more important role. It is likely that the high gene content in the microchromosomes promotes purifying selection against harmful mutations like InDels [57, 58]. Chromosome 16 is an example of rich harbour of many important genes or gene clusters viz. the major histocompatibility complex (MHC), nucleolus organiser region (NOR), olfactory receptor (OR), cysteine-rich domain scavenger receptor (SRCR) and putative immunoglobulin-like receptors . This might have prompted the purifying selection to keep the InDel rate in this chromosome as one of the lowest (0.52 per kb) among the autosomes. Besides, this chromosome is rich in duplicated regions, which possibly have affected the alignment of sequence reads leading to detection of very few InDels in the first place . Chromosome Z also showed a very low density of InDels. In contrast to the microchromosomes, the low density of InDels in GGAZ is a reflection of reduced overall genetic diversity observed in this chromosome due to a number of possible reasons such as selection on sex-linked characters and low male effective population size [62, 34]. Similar to our findings, several other studies have also reported lower densities of InDels in micro- and sex chromosomes [36, 38]. Besides, the sequenced birds consisted of both male and female samples and as a result the depth of coverage for sex chromosomes is expected to be less than that of autosomes, which may also had an impact on the number variants detected from these chromosomes.
The detection of relatively higher number of FS compared to NFS InDels is consistent with the results from previous studies on different species [38, 63]. Non-triplet InDels (such as FS) are at least twice more likely to be generated by chance alone compared to the triplet InDels (e.g. NFS). However, when present in coding-regions, non-triplet InDels can be highly disruptive as they cause complete change in the sequence of C-terminal part of the protein or may uncover a stop codon to truncate the translated product. The peptides produced by FS mutations have been found to be associated with many diseases [64–67]. As a consequence, a much lower proportion of non-triplet variant is expected in coding regions due to purifying selection compared to that in non-coding regions. Our finding conformed to this expectation as the ratio of non-triplet and triplet InDels was about 4 times less in coding regions compared with that in non-coding parts.
The observation that FS mutations appeared more frequently at the beginning or end of the proteins corroborates with the understanding that FS mutations in the middle of the proteins will have more harmful impacts than those at the end. However, one caveat of this finding is that the presence of many FS InDels at the extreme ends of proteins may reflect wrong annotation of genes and their corresponding translation start site in the current gene database. It is possible that these InDels are actually located outside the boundary of coding region. To assess this possibility, we investigated some randomly selected genes (ca. 10) with FS mutations that were located at the beginning of the protein (position ≤0.1 relative to polypeptide length). Blastp analyses of the proteins encoded by these genes showed that in most cases these are known proteins in other species. However, when we predicted the translation start sites of these genes with the NetStart 1.0 programme , in most cases (9 out of 10 genes checked) different initiation codons were predicted from the ones specified in Ensembl. In all these cases, the FS InDels were located before the NetStart predicted translation initiation site, indicating the possibility of inaccurate gene annotation. On the other hand, when we randomly checked 10 genes with NFS mutation as a control, in most cases the NetStart predicted initiation sites were same as those given by Ensembl. However, correcting the annotation of coding regions is beyond the scope of the present study and would require further work, such as a genome-wide proteomic investigation.
In the present study we identified a number of gene classes that were significantly enriched with potentially harmful mutations and some of the FS and intolerant-NFS InDels within these gene classes were present in high frequency in multiple or single lines. It is difficult to explain how genes associated with important physiological functions can harbour such harmful mutations. However, presence of paralogs for majority (67 %) of these genes can be a possible explanation as these can help to compensate for the loss-of-function caused by harmful mutations. We explored if the fixed FS and NFS-InDels overlapped with any chicken QTLs as this might help find association of these variants with important traits. We found that the FS insertion (GGA2:46501278) in the DCLK3 gene that was near-fixed in both WEL and BEL groups, coincided with QTLs associated with a wide range of traits such as tibia bone mineral content, egg production, blood quality, fatness, etc. as specified in Chicken QTLdb (http://www.animalgenome.org/cgi-bin/QTLdb/GG/index). The FS insertion (GGA11:483395) detected from the KIFC3 gene from an individual WEL line overlapped with QTLs associated with tibia area and stress, intestine length and carcass traits. The intolerant NFS mutation (GGA5: 31527889) detected in the AQR gene from a WEL line and the FS insertion (GGA7:21670698) from the DES gene overlapped with QTLs associated with fatness traits, egg production, carcass traits, antibody responses to diseases, etc. Overlap with QTLs associated with many traits suggests possible pleiotropic effects of these genes and their functional variants. An alternative possibility to selective advantage is that these InDels are moderately (or lowly) harmful under the environmental conditions where the birds are maintained but have attained fixation through hitchhiking by being linked to some beneficial mutations under strong positive selection .
The CF InDels detected in our study were compared with those available in the public domain (dbSNP build 140) and those recently described by Yan et al. . When the comparison was made based on coordinates alone, about 81 % of our CF variants were found to be shared with the already known InDels. However, when the allelic information was combined with the coordinates for comparison, the percentage of shared variants reduced to about 53 %, indicating that the non-reference allele at an InDel site may differ among populations. The major reason for the lower sharing of alleles between studies is that a substantial proportion of these InDels are tandem repeats that are prone to be highly variable. For instance, about 59 % of these non-shared variants due to allele were perfect tandem repeats of monomeric nature (i.e. repeats of either A, T, C or G), while many others were imperfect or complex repeats. Another possible reason can be the misalignment of the reads, particularly when majority of these inconsistent alleles consisted of tandem repeats. To investigate this possibility, we randomly checked about 50 InDels and most of them appeared true InDels with correct alignments and good coverage in our study. The observation of large proportion of common InDels across studies has served as a validation that most of these InDels are true variants.
This study provides a large catalogue of small insertion and deletion genetic variants and their detailed characterization by analysing many chickens from diverse commercial and experimental layer lines. Use of consensus InDels from two different bioinformatics packages followed by stringent filtration criteria provided confidence in the detected set and the FDR estimation suggested a high rate of validation. Moreover, overlap of a large proportion (~81 %) of the 883 K InDels detected in the present study with Yan et al.  is a further proof that majority of the InDels detected in our study and those present in the public domain are correct. This paper adds about 168 K novel InDels over what has already been detected from chicken genome by previous studies. Most importantly, this study provides an in-depth characterisation of the InDels by adopting different approaches and using available resources. To the best of our knowledge, this is the first study on chicken to shed any light on the frequency spectrum of detected InDels and it indicates that majority of the InDels detected so far from chickens by different studies are probably high frequency, common variants. The results of this study also suggest that sequencing at much higher coverage will be required to detect rare and low frequency InDels.
The resource created in this study is expected to have major implications in future studies not only in chicken but also in other avian species. For instance, the large catalogue of InDels along with their functional characteristics can help gain insights into the genetic variant profile of chicken genome and in the identification of causal variants underlying various diseases and other important traits. Particularly, potentially functional variants that have been found to be fixed or near fixed in different groups can be further explored to understand their possible effects on phenotypes.
Whole genome re-sequencing of layer chickens lines
The NGS sequence data used for InDel detection was generated under a previous study . In brief the data was generated by Illumina sequencing of 163 chickens originating from 11 commercial and 5 experimental layer lines (Additional file 1). The samples from commercial lines were supplied by Hy-line International and Lohmann in the Synbreed Consortium and consisted of 6 lines of white egg layers (WEL1-6) and 5 lines of brown egg layers (BEL1-5). Other samples originated from four inbred lines (Wellcome, N, 15 and 0 lines) from the Institute of Animal Health (now The National Avian Research Facility (NARF), Edinburgh) and an unselected brown leghorn line from the Roslin Institute (RI-J line). For each line, except WEL6, 10–15 individuals were sequenced as pooled DNA using a paired end protocol. For WEL6, however, three samples were sequenced individually. The length of paired reads varied from 76–101 nucleotides. Sequencing was performed with a mean observed coverage of 8-17X per line. Further details on library preparation, sequencing and alignment can be retrieved from Kranis et al. . We mapped detected InDels on the published chicken reference genome (Gallus_gallus_4.0).
InDel calling by SAMtools
The InDel calling with SAMtools (version 0.1.18)  was performed using its mpileup function and BCFtools. The mpileup function computes the likelihood of the data given each possible genotype and stores the likelihoods in BCF format (binary variant call). BCFtools then applies the prior and does the actual variant calling. We used the following commands for calling the variants using SAMtools.
samtools mpileup –q20 –Q20 -AB -ugf <referenceFile.fa> <bamFile.bam> | bcftools view -bvcg><var.raw.bcf>
bcftools view <var.raw.bcf> | vcfutils.pl varFilter -D99999> <var.flt.vcf>
In these commands, the –q and –Q options were given to specify the minimum thresholds for base and map qualities to be 20. The option –A was used to force the analysis of all the reads including anomalous read pairs to avoid issues in the analyses such as sudden stops. The option –B was used to disable the BAQ (Base Alignment Quality) calculation as it is computationally very demanding. Disabling BAQ may result in increased quality scores in some false positive SNP calls close to InDels, but it will not affect the alignment process in InDel calling. The maximum depth of coverage was limited to 99999 by the option –D.
InDel calling by Dindel
InDel calling using the Dindel package (version 1.01)  consisted of four stages:
Stage 1– extraction of all candidate InDels
dindel --analysis getCIGARInDels –bamFile < bamFile.bam > --outputFile < outFile-stage1 > --ref < referenceFile.fa>
Extra step in Stage 1– reduction of the number of candidate InDels:
python selectCandidates.py -i <outFile-stage1> −o <outFile-extrastage>
Stage 2– grouping the candidate InDels into windows of ~120 bp
python makeWindows.py --inputVarFile <outFile_extrastage> --windowFilePrefix <outFile-stage2> --numWindowsPerFile 20000
Stage 3– generation for every window of candidate haplotypes from candidate InDels
dindel --analysis indels --doPooled –bamFile <bamFile.bam> --ref <referenceFile.fa> --varFile <outFile-stage2> --libFile <librarieFile> --outputFile <outFile-stage3>
Stage 4 – production of a VCF4 file format with the InDels calls:
python mergeOutputPooled.py --inputFiles <outFile-stage3> --outputFile <outFile-stage4> --ref <referenceFile.fa> --numSamples 3–15 --numBamFiles 1
We used this procedure for calling InDels from a pool, and the default parameters which included the selection of candidate InDels seen at least twice (in the stage to reduce the number of candidate InDels).
Filtration of consensus InDels
We obtained the common or consensus set of InDels from the two InDel callers, SAMtools and Dindel, by comparing the co-ordinates of the InDels from these two packages. The following filtration criteria were then applied on the consensus set to reduce the number of false-positives: (i) InDel quality score ≥30 based on SAMtools provided score; (ii) coverage at the InDel position ≥5 and ≤ mean coverage in a lines + 3 SD (standard deviation); (iii) the non-reference allele supported by both forward and reverse strands; and (iv) gap between consecutive InDels >1 base. Finally, we removed those InDels which are fixed for the non-reference allele in all or most of the lines (14–16 lines) as these may actually represent possible sequencing errors in reference genome or alignment error. Most of the consensus filtered InDels (n = 865,597) have been submitted to dbSNP (NCBI) using the handle “DWBURT” with the submitter batch ID “Chicken_indel_dwburt” (http://www.ncbi.nlm.nih.gov/projects/SNP/snp_viewBatch.cgi?sbid=1062064).
The large InDels (>50 bp) and the majority of block substitution InDels have been submitted to EVA database (provisional accession number PRJEB9374; http://www.ebi.ac.uk/eva/).
Calculation of FDRs for InDels
We estimated the FDRs of InDel calling by resequencing several random regions with the Sanger method followed by comparing the detected InDels from the Sanger and NGS data from these regions. For Sanger sequencing, we initially selected 28 regions from different chromosomes (GGA1-6, 10–12, 15–17, 19, 21, 22, 24, 26, 28 and linkage group LGE22C19W28_E50C23). All the ten chickens from the RI-J line were individually sequenced using the Sanger method. The sequences were mapped to reference genome using the long read mapping algorithm in BWA and also separately with DNASTAR MegAlign™ (http://www.dnastar.com/t-megalign.aspx) programme with default parameters. The BWA alignment was used with SAMtools to obtain the InDel locations fast, whereas the trace files from MegAlignTM were used to manually check the correctness of alignments and for manual calling of the variants.
Low confidence Sanger InDels i.e. those with poor sequence quality, poor alignment or lacking support from both forward and reverse strands were excluded from the calculations. Three regions from GGA10, 15 and 17 had to be removed due to missing or poor quality sequence from one or more samples by Sanger method even after repeated attempts. Further, one region from GGA3 was removed due to alignment issue, leaving 24 regions for estimation of FDRs.
Based on the comparison of Sanger and NGS variants from the same region an InDel was called true positive (TP) when it was detected by both NGS and Sanger, a false positive (FP) when it was detected only by NGS, a false negative (FN) when it was detected only by Sanger. Any sequenced bases which were not called as InDels by either Sanger or NGS were considered to be true negatives (TN). Based on these, FDR calculations were defined as follows: Sensitivity (rate of true InDels correctly identified by NGS) = TP/(TP + FN); Specificity (rate of true non-InDels that were correctly recognized by NGS) = TN/(TN + FP); False positive rate (FPR) = (1-specificity); and False negative rate (FNR) = (1-sensitivity). The proportion of false positive InDels in the filtered NGS list was calculated as (total number of NGS InDels – total number of InDels detected by both NGS and Sanger)/total number of NGS InDels).
Functional annotation and effect prediction of InDels
The genomic positions of the InDels and their effect on protein coding regions were predicted by annotating them against the Ensembl gene annotation database (release 71) for chicken. The software ANNOVAR (version July 06, 2012)  was used for this purpose. In addition, we also annotated the InDels against 1,608 novel non-coding RNA transcripts (ncRNA), which have recently been characterized . To identify if any of the InDels coincided with these ncRNA transcripts we used the BEDtools (version 2.17.0, http://bedtools.readthedocs.org/en/latest/).
The InDels within the coding regions (non-frameshift and stop-gain/loss InDels) were further analysed using the Protein Variation Effect Analyser (PROVEAN, version 1.1) , which predicts whether an amino acid substitution, insertion or deletion are likely to affect protein function or not. PROVEAN was run with the default parameters on the NCBI non-redundant protein database. Delta score of −2.5 was used as the threshold below which any InDel was predicted “Deleterious” or “Intolerant” in evolutionary term.
The filtered list of InDels were also annotated using ANNOVAR against the PhastCons predicted “Most Conserved Elements” (MCE) for chicken. For this annotation, however, we first had to map the InDels against the previous version of chicken reference sequence (Gallus_gallus_2.1) as the MCE dataset contains co-ordinates in relation to this reference build. It was possible to map unambiguously only 585,154 (66 %) InDels on Gallus_gallus_2.1 and as a result we could only annotate these variants. The MCE data were downloaded from UCSC database: ftp://hgdownload.cse.ucsc.edu/goldenPath/galGal3/database/phastConsElements7way.txt.gz. This data contains PhastCons scores  for chicken (Gallus_gallus_2.1 or 3, May 2006).
The DAVID Gene Functional Classification tool  was used to find if any particular classes of genes were enriched for frameshift and intolerant non-frameshift mutations in the BEL, WEL and inbred chicken groups. For this analysis, we selected only those genes with kappa score ≥0.75 and clusters with enrichment scores ≥1.3.
Venn diagram, heatmap and allele frequency estimation
Venn diagrams of InDels shared between different lines were created by the BioVenn tool (available at http://www.cmbi.ru.nl/cdd/biovenn/). Estimation of allele frequency of InDels was based on the proportion of high quality reads supporting the non-reference alleles. Mean frequencies within groups (e.g. BEL, WEL and Inbred) were calculated based on the populations where the InDels were detected. Heatmap of fixed InDels (AAF ≥ 0.9) was generated by the programme, Genesis .
This work was funded by a DEFRA/BBSRC LINK grant in collaboration with The Roslin Institute (RI) of University of Edinburgh, Hy-line International, Affymetrix and the Synbreed Consortium. CB received fellowship from CNPq, National Council for Scientific and Technological Development – Brazil.
- Vali U, Brandstrom M, Johansson M, Ellegren H. Insertion-deletion polymorphisms (indels) as genetic markers in natural populations. BMC Genet. 2008;9:8.PubMedPubMed CentralView ArticleGoogle Scholar
- Mullaney JM, Mills RE, Pittard WS, Devine SE. Small insertions and deletions (INDELs) in human genomes. Hum Mol Genet. 2010;19:R131–R6.PubMedPubMed CentralView ArticleGoogle Scholar
- Stenson PD, Ball EV, Mort M, Phillips AD, Shiel JA, Thomas NST, et al. Human gene mutation database (HGMD (R)): 2003 update. Hum Mutat. 2003;21(6):577–81.PubMedView ArticleGoogle Scholar
- Fang MX, Nie QH, Luo CL, Zhang DX, Zhang XQ. An 8 bp indel in exon 1 of Ghrelin gene associated with chicken growth. Domest Anim Endocrin. 2007;32(3):216–25.View ArticleGoogle Scholar
- Kerje S, Sharma P, Gunnarsson U, Kim H, Bagchi S, Fredriksson R, et al. The Dominant white, Dun and Smoky color variants in chicken are associated with insertion/deletion polymorphisms in the PMEL17 gene. Genetics. 2004;168(3):1507–18.PubMedPubMed CentralView ArticleGoogle Scholar
- Cui JX, Du HL, Liang Y, Deng XM, Li N, Zhang XQ. Association of polymorphisms in the promoter region of chicken prolactin with egg production. Poult Sci. 2006;85(1):26–31.PubMedView ArticleGoogle Scholar
- Zhang S, Han RL, Gao ZY, Zhu SK, Tian YD, Sun GR, et al. A novel 31-bp indel in the paired box 7 (PAX7) gene is associated with chicken performance traits. Brit Poult Sci. 2014;55(1):31–6.View ArticleGoogle Scholar
- Tang S, Ou J, Sun D, Zhang Y, Xu G, Zhang Y. A novel 62-bp indel mutation in the promoter region of transforming growth factor-beta 2 (TGFB2) gene is associated with body weight in chickens. Anim Genet. 2011;42(1):108–12.PubMedView ArticleGoogle Scholar
- Tummala H, Ali M, Getty P, Hocking PM, Burt DW, Inglehearn CF, et al. Mutation in the guanine nucleotide-binding protein beta-3 causes retinal degeneration and embryonic mortality in chickens. Invest Ophth Vis Sci. 2006;47(11):4714–8.View ArticleGoogle Scholar
- Mills RE, Pittard WS, Mullaney JM, Farooq U, Creasy TH, Mahurkar AA, et al. Natural genetic variation caused by small insertions and deletions in the human genome. Genome Res. 2011;21(6):830–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Iengar P. An analysis of substitution, deletion and insertion mutations in cancer genes. Nucleic Acids Res. 2012;40(14):6401–13.PubMedPubMed CentralView ArticleGoogle Scholar
- Ball EV, Stenson PD, Abeysinghe SS, Krawczak M, Cooper DN, Chuzhanova NA. Microdeletions and microinsertions causing human genetic disease: Common mechanisms of mutagenesis and the role of local DNA sequence complexity. Hum Mutat. 2005;26(3):205–13.PubMedView ArticleGoogle Scholar
- Britten RJ. Divergence between samples of chimpanzee and human DNA sequences is 5 %, counting indels. Proc Natl Acad Sci U S A. 2002;99(21):13633–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Rubin CJ, Zody MC, Eriksson J, Meadows JRS, Sherwood E, Webster MT, et al. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010;464(7288):587–91.PubMedView ArticleGoogle Scholar
- Chen CH, Chuang TJ, Liao BY, Chen FC. Scanning for the signatures of positive selection for human-specific insertions and deletions. Genome Biol Evol. 2009;1:415–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Neuman JA, Isakov O, Shomron N. Analysis of insertion-deletion from deep-sequencing data: software evaluation for optimal detection. Brief Bioinform. 2013;14(1):46–55.PubMedView ArticleGoogle Scholar
- Li H, Ruan J, Durbin R. Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008;18(11):1851–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Albers CA, Lunter G, MacArthur DG, McVean G, Ouwehand WH, Durbin R. Dindel: Accurate indel calls from short-read data. Genome Res. 2011;21(6):961–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Messer PW, Arndt PF. The majority of recent short DNA insertions in the human genome are tandem duplications. Mol Biol Evol. 2007;24(5):1190–7.PubMedView ArticleGoogle Scholar
- Madsen BE, Villesen P, Wiuf C. Short tandem repeats in human exons: a target for disease mutations. BMC Genomics. 2008;9:410.PubMedPubMed CentralView ArticleGoogle Scholar
- Bansal V, Libiger O. A probabilistic method for the detection and genotyping of small indels from population-scale sequence data. Bioinformatics. 2011;27(15):2047–53.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Homer N. A survey of sequence alignment algorithms for next-generation sequencing. Brief Bioinform. 2010;11(5):473–83.PubMedPubMed CentralView ArticleGoogle Scholar
- Krawitz P, Rodelsperger C, Jager M, Jostins L, Bauer S, Robinson PN. Microindel detection in short-read sequence data. Bioinformatics. 2010;26(6):722–9.PubMedView ArticleGoogle Scholar
- Jia PL, Li F, Xia JF, Chen HQ, Ji HB, Pao W, et al. Consensus rules in variant detection from next-generation sequencing data. PLoS One. 2012;7(6):e38470.PubMedPubMed CentralView ArticleGoogle Scholar
- Gregory TR. Insertion-deletion biases and the evolution of genome size. Gene. 2004;324:15–34.PubMedView ArticleGoogle Scholar
- Lovett ST. Encoded errors: mutations and rearrangements mediated by misalignment at repetitive DNA sequences. Mol Microbiol. 2004;52(5):1243–53.PubMedView ArticleGoogle Scholar
- Tanay A, Siggia ED. Sequence context affects the rate of short insertions and deletions in flies and primates. Genome Biol. 2008;9(2):R37.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim R, Guo JT. Systematic analysis of short internal indels and their impact on protein folding. BMC Struct Biol. 2010;10:24.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen FC, Chen CJ, Li WH, Chuang TJ. Human-specific insertions and deletions inferred from mammalian genome sequences. Genome Res. 2007;17(1):16–22.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang ZDD, Du J, Lam H, Abyzov A, Urban AE, Snyder M, et al. Identification of genomic indels and structural variations using split reads. BMC Genomics. 2011;12:375.PubMedPubMed CentralView ArticleGoogle Scholar
- Consortium TGP. A map of human genome variation from population-scale sequencing. Nature. 2010;467(7319):1061–73.View ArticleGoogle Scholar
- Wong GKS, Liu B, Wang J, Zhang Y, Yang X, Zhang ZJ, et al. A genetic variation map for chicken with 2.8 million single-nucleotide polymorphisms. Nature. 2004;432(7018):717–22.PubMedView ArticleGoogle Scholar
- Kranis A, Gheyas AA, Boschiero C, Turner F, Yu L, Smith S, et al. Development of a high density 600 K SNP genotyping array for chicken. BMC Genomics. 2013;14:59.PubMedPubMed CentralView ArticleGoogle Scholar
- Gheyas AA, Boschiero C, Eory L, Ralph H, Kuo R, Woolliams JA, et al. Functional classification of 15 million SNPs detected from diverse chicken populations. DNA Res. 2015;22(3):205–17.PubMedPubMed CentralView ArticleGoogle Scholar
- Brandstrom M, Ellegren H. The genomic landscape of short insertion and deletion polymorphisms in the chicken (Gallus gallus) genome: A high frequency of deletions in tandem duplicates. Genetics. 2007;176(3):1691–701.PubMedPubMed CentralView ArticleGoogle Scholar
- Fan WL, Ng CS, Chen CF, Lu MYJ, Chen YH, Liu CJ, et al. Genome-wide patterns of genetic variation in two domestic chickens. Genome Biol Evol. 2013;5(7):1376–92.PubMedPubMed CentralView ArticleGoogle Scholar
- Yan Y, Yi G, Sun C, Qu L, Yang N. Genome-wide characterization of insertion and deletion variation in chicken using next generation sequencing. PLoS One. 2014;9(8):e104652.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhan BJ, Fadista J, Thomsen B, Hedegaard J, Panitz F, Bendixen C. Global assessment of genomic variation in cattle by genome resequencing and high-throughput genotyping. BMC Genomics. 2011;12:557.PubMedPubMed CentralView ArticleGoogle Scholar
- Sebastiani P, Riva A, Montano M, Pham P, Torkamani A, Scherba E, et al. Whole genome sequences of a male and female supercentenarian, ages greater than 114 years. Front Genet. 2011;2:90.PubMedGoogle Scholar
- Kettleborough RN, Busch-Nentwich EM, Harvey SA, Dooley CM, de Bruijn E, van Eeden F, et al. A systematic genome-wide analysis of zebrafish protein-coding gene function. Nature. 2013;496(7446):494–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Samtools. http://samtools.sourceforge.net/mpileup.shtml (site visited on 02/09/2014).
- Park MH, Rhee H, Park JH, Woo HM, Choi BO, Kim BY, et al. Comprehensive analysis to improve the validation rate for single nucleotide variants detected by next-generation sequencing. PLoS One. 2014;9(1):e86664.PubMedPubMed CentralView ArticleGoogle Scholar
- Raineri E, Ferretti L, Esteve-Codina A, Nevado B, Heath S, Perez-Enciso M. SNP calling by sequencing pooled samples. BMC Bioinformatics. 2012;13:239.PubMedPubMed CentralView ArticleGoogle Scholar
- Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PLoS One. 2012;7(10):e46688.PubMedPubMed CentralView ArticleGoogle Scholar
- Ng PC, Levy S, Huang J, Stockwell TB, Walenz BP, Li K, et al. Genetic variation in an individual human exome. PLoS Genet. 2008;4(8):e1000160.PubMedPubMed CentralView ArticleGoogle Scholar
- Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou MM, Rosenbloom K, et al. Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005;15(8):1034–50.PubMedPubMed CentralView ArticleGoogle Scholar
- UCSC. Multiple alignments of 6 vertebrate genomes with chicken. http://genome.ucsc.edu/goldenpath/help/ftp.html (site visited on 03/09/2014).
- Cooper DN. Functional intronic polymorphisms: Buried treasure awaiting discovery within our genes. Hum Genomics. 2010;4(5):284–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.View ArticleGoogle Scholar
- Cutler DJ, Jensen JD. To Pool, or Not to Pool? Genetics. 2010;186(1):41–3.PubMedPubMed CentralView ArticleGoogle Scholar
- Ingman M, Gyllensten U. SNP frequency estimation using massively parallel sequencing of pooled DNA. Eur J Hum Genet. 2009;17(3):383–6.PubMedView ArticleGoogle Scholar
- Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008;456(7218):53–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Eck SH, Benet-Pages A, Flisikowski K, Meitinger T, Fries R, Strom TM. Whole genome sequencing of a single Bos taurus animal for single nucleotide polymorphism discovery. Genome Biol. 2009;10(8):R82.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen K, McLellan MD, Ding L, Wendl MC, Kasai Y, Wilson RK, et al. PolyScan: An automatic indel and SNP detection approach to the analysis of human resequencing data. Genome Res. 2007;17(5):659–66.PubMedPubMed CentralView ArticleGoogle Scholar
- Burt DW. Chicken genome: current status and future opportunities. Genome Res. 2005;15(12):1692–8.PubMedView ArticleGoogle Scholar
- Hillier LW, Miller W, Birney E, Warren W, Hardison RC, Ponting CP, et al. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004;432(7018):695–716.View ArticleGoogle Scholar
- Begun DJ, Aquadro CF. Levels of naturally occurring DNA polymorphism correlate with recombination rates in D. melanogaster. Nature. 1992;356(6369):519–20.PubMedView ArticleGoogle Scholar
- Nachman MW. Single nucleotide polymorphisms and recombination rate in humans. Trends Genet. 2001;17(9):481–5.PubMedView ArticleGoogle Scholar
- Miller MM, Robinson CM, Abernathy J, Goto RM, Hamilton MK, Zhou H, et al. Mapping genes to chicken microchromosome 16 and discovery of olfactory and scavenger receptor genes near the major histocompatibility complex. J Hered. 2014;105(2):203–15.PubMedView ArticleGoogle Scholar
- Sundstrom H, Webster MT, Ellegren H. Reduced variation on the chicken Z chromosome. Genetics. 2004;167(1):377–85.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee KT, Chung WH, Lee SY, Choi JW, Kim J, Lim D, et al. Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity. BMC Genomics. 2013;14:519.PubMedPubMed CentralView ArticleGoogle Scholar
- Shiihara T, Watanabe M, Moriyama K, Uematsu M, Sameshima K. A novel PLP1 frameshift mutation causing a milder form of Pelizaeus-Merzbacher disease. Brain Dev. 2014;37(4):455–8.PubMedView ArticleGoogle Scholar
- Ogura Y, Bonen DK, Inohara N, Nicolae DL, Chen FF, Ramos R, et al. A frameshift mutation in NOD2 associated with susceptibility to Crohn’s disease. Nature. 2001;411(6837):603–6.PubMedView ArticleGoogle Scholar
- Myerowitz R. Tay-Sachs disease-causing mutations and neutral polymorphisms in the Hex A gene. Hum Mutat. 1997;9(3):195–208.PubMedView ArticleGoogle Scholar
- Iannuzzi MC, Stern RC, Collins FS, Hon CT, Hidaka N, Strong T, et al. Two frameshift mutations in the cystic fibrosis gene. Am J Hum Genet. 1991;48(2):227–31.PubMedPubMed CentralGoogle Scholar
- Pedersen AG, Nielsen H. Neural network prediction of translation initiation sites in eukaryotes: perspectives for EST and genome analysis. Proceedings/International Conference on Intelligent Systems for Molecular Biology; ISMB International Conference on Intelligent Systems for Molecular Biology. 1997;5:226–33.Google Scholar
- Chun S, Fay JC. Evidence for hitchhiking of deleterious mutations within the human genome. PLoS Genet. 2011;7(8):e1002240.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang K, Li M, Hakonarson H. ANNOVAR: Functional annotation of genetic variants from next-generation sequencing data. Nucleic Acids Res. 2010;38:e164.PubMedPubMed CentralView ArticleGoogle Scholar
- Gardner PP, Fasold M, Burge SW, Ninova M, Hertel J, Kehr S, et al. Conservation and losses of Non-Coding RNAs in Avian genomes. PLoS One. 2015;10(3):e0121797.PubMedPubMed CentralView ArticleGoogle Scholar
- Sturn A, Quackenbush J, Trajanoski Z. Genesis: cluster analysis of microarray data. Bioinformatics. 2002;18(1):207–8.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.