|Fine mapping of complex traits in non-model species: using next generation sequencing and advanced intercross lines in Japanese quail.|
|Jump to Full Text|
|PMID: 23066875 Owner: NLM Status: MEDLINE|
|BACKGROUND: As for other non-model species, genetic analyses in quail will benefit greatly from a higher marker density, now attainable thanks to the evolution of sequencing and genotyping technologies. Our objective was to obtain the first genome wide panel of Japanese quail SNP (Single Nucleotide Polymorphism) and to use it for the fine mapping of a QTL for a fear-related behaviour, namely tonic immobility, previously localized on Coturnix japonica chromosome 1. To this aim, two reduced representations of the genome were analysed through high-throughput 454 sequencing: AFLP (Amplified Fragment Length Polymorphism) fragments as representatives of genomic DNA, and EST (Expressed Sequence Tag) as representatives of the transcriptome.
RESULTS: The sequencing runs produced 399,189 and 1,106,762 sequence reads from cDNA and genomic fragments, respectively. They covered over 434 Mb of sequence in total and allowed us to detect 17,433 putative SNP. Among them, 384 were used to genotype two Advanced Intercross Lines (AIL) obtained from three quail lines differing for duration of tonic immobility. Despite the absence of genotyping for founder individuals in the analysis, the previously identified candidate region on chromosome 1 was refined and led to the identification of a candidate gene.
CONCLUSIONS: These data confirm the efficiency of transcript and AFLP-sequencing for SNP discovery in a non-model species, and its application to the fine mapping of a complex trait. Our results reveal a significant association of duration of tonic immobility with a genomic region comprising the DMD (dystrophin) gene. Further characterization of this candidate gene is needed to decipher its putative role in tonic immobility in Coturnix.
|Laure Frésard; Sophie Leroux; Patrice Dehais; Bertrand Servin; Hélène Gilbert; Olivier Bouchez; Christophe Klopp; Cédric Cabau; Florence Vignoles; Katia Feve; Amélie Ricros; David Gourichon; Christian Diot; Sabine Richard; Christine Leterrier; Catherine Beaumont; Alain Vignal; Francis Minvielle; Frédérique Pitel|
Related Documents :
|16005165 - Asymmetry of intronic pre-mirna structures in functional risc assembly.
17267435 - De novo svm classification of precursor micrornas from genomic pseudo hairpins using gl...
24202755 - Comparative studies of the structure of chloroplast dna from four species of oryza: clo...
12509205 - Influence of stimulus distance in implicit learning of spatial and nonspatial event seq...
9573215 - Identification and enzymatic characterization of the maltose-inducible alpha-glucosidas...
12906855 - Characterization of the transmembrane channel-like (tmc) gene family: functional clues ...
|Type: Journal Article; Research Support, Non-U.S. Gov't Date: 2012-10-15|
|Title: BMC genomics Volume: 13 ISSN: 1471-2164 ISO Abbreviation: BMC Genomics Publication Date: 2012|
|Created Date: 2013-01-03 Completed Date: 2013-05-14 Revised Date: 2013-07-11|
Medline Journal Info:
|Nlm Unique ID: 100965258 Medline TA: BMC Genomics Country: England|
|Languages: eng Pagination: 551 Citation Subset: IM|
|INRA, UMR444 Laboratoire de Génétique Cellulaire, Castanet-Tolosan, F-31326, France.|
|APA/MLA Format Download EndNote Download BibTex|
Amplified Fragment Length Polymorphism Analysis
Avian Proteins / genetics*
Chickens / genetics
Coturnix / genetics*
Dystrophin / genetics*
Expressed Sequence Tags
Genetic Association Studies*
High-Throughput Nucleotide Sequencing
Immobility Response, Tonic*
Polymorphism, Single Nucleotide
Quantitative Trait Loci
|0/Avian Proteins; 0/Dystrophin|
Journal ID (nlm-ta): BMC Genomics
Journal ID (iso-abbrev): BMC Genomics
Publisher: BioMed Central
Copyright ©2012 Frésard et al.; licensee BioMed Central Ltd.
Received Day: 24 Month: 5 Year: 2012
Accepted Day: 8 Month: 10 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 15 Month: 10 Year: 2012
Volume: 13First Page: 551 Last Page: 551
PubMed Id: 23066875
Publisher Id: 1471-2164-13-551
|Fine mapping of complex traits in non-model species: using next generation sequencing and advanced intercross lines in Japanese quail|
|Laure Frésard12||Email: email@example.com|
|Sophie Leroux12||Email: firstname.lastname@example.org|
|Patrice Dehais12||Email: email@example.com|
|Bertrand Servin12||Email: firstname.lastname@example.org|
|Hélène Gilbert12||Email: email@example.com|
|Olivier Bouchez123||Email: firstname.lastname@example.org|
|Christophe Klopp4||Email: email@example.com|
|Cédric Cabau5||Email: firstname.lastname@example.org|
|Florence Vignoles12||Email: email@example.com|
|Katia Feve12||Email: firstname.lastname@example.org|
|Amélie Ricros12||Email: email@example.com|
|David Gourichon6||Email: firstname.lastname@example.org|
|Christian Diot78||Email: email@example.com|
|Sabine Richard91011||Email: firstname.lastname@example.org|
|Christine Leterrier91011||Email: email@example.com|
|Catherine Beaumont12||Email: firstname.lastname@example.org|
|Alain Vignal12||Email: email@example.com|
|Francis Minvielle1314||Email: firstname.lastname@example.org|
|Frédérique Pitel12||Email: email@example.com|
1INRA, UMR444 Laboratoire de Génétique Cellulaire, Castanet-Tolosan, F-31326, France
2ENVT, UMR444 Laboratoire de Génétique Cellulaire, Toulouse, F-31076, France
3INRA, GeT-PlaGe Genotoul, Castanet-Tolosan, F-31326, France
4INRA, Sigenae UR875 Biométrie et Intelligence Artificielle, Castanet-Tolosan, F-31326, France
5INRA, Sigenae, Nouzilly, F-37380, France
6INRA, PEAT Pôle d'Expérimentation Avicole de Tours, Nouzilly, F- 37380, France
7INRA, UMR1348 PEGASE, Saint-Gilles, F-35590, France
8Agrocampus Ouest, UMR1348 PEGASE, Rennes, F-35000, France
9INRA, UMR85 Physiologie de la Reproduction et des Comportements, Nouzilly, F-37380, France
10CNRS, UMR6175 Physiologie de la Reproduction et des Comportements, Nouzilly, F-37380, France
11Université François Rabelais Tours, UMR Physiologie de la Reproduction et des Comportements, Nouzilly, F-37380, France
12INRA, UR83 Recherche Avicoles, Nouzilly, F- 37380, France
13INRA, UMR1313 Génétique animale et biologie intégrative, Jouy en Josas, F-78350, France
14AgroParisTech, UMR Génétique animale et biologie intégrative, Jouy en Josas, F-78350, France
The two avian Coturnix species - Japanese quail (Coturnix japonica) and Common quail (Coturnix coturnix) - belong to the order Galliformes and to the family Phasianidae, which also includes the chicken (Gallus gallus domesticus). The distributions of the two quail species overlap in the lake Baikal region, to and from where wild quail migrate every year, toward Asia for C. japonica and toward Europe and Africa for C. coturnix. However, hybridization is possible between the two species , and hybrid males have a distinct crow .
In Asia, the Japanese quail was domesticated for its singing abilities in the 15th century, and selected lines were developed in Japan for egg production at the beginning of the 20th century . Japanese quail lines were then introduced successively in North America, Europe and the Near and Middle East in the first half of the 20th century . Its small size, short generation interval, good laying abilities and ease of reproduction interested scientists looking for a bird easy to maintain and to experiment with in the laboratory, or willing to develop alternative ways of producing eggs and meat for human consumption . However, for unknown reasons, the Common quail was not domesticated. Research on Japanese quail has led to the development and use of commercial quail strains for meat and egg production in Asia, Europe and other parts of the world  and to a better knowledge of its genetics .
The first modern genetic tools developed specifically for quail analyses were panels of markers [8,9] and genetic maps using AFLP (Amplified Fragment Length Polymorphism)  and microsatellites . Despite their limited extent, they have already been used extensively. At the population level, genome scans were carried out to detect QTL (Quantitive Trait Loci) for traits associated to fearfulness , growth, egg production and duration of tonic immobility  and egg laying curve , and markers were used to discriminate the two species and their hybrids [15,16]. At the genome level, genes carrying causal mutations were found for plumage colours, like TYRP1 for the roux plumage  and ASIP for the yellow  and recessive black  plumages. More recently, the existing quail maps were integrated and aligned to the assembled chicken sequence , thereby confirming the high degree of synteny conservation expected because of the high level of chromosome homology observed between both species , and the relevance of quail models for trait analyses in poultry.
Further developments in quail genomics are held back by the small number of available markers which delay the search for other causal mutations. It is especially critical for traits with no known homologous phenotypes in chicken, like the rusty plumage colour and the curly-shaped feathers , and for the fine mapping of QTL. Our objective was to develop a dedicated quail SNP (Single Nucleotide Polymorphism) panel, in order to fine map a QTL for duration of tonic immobility. Tonic immobility (TI) is an innate fear-related behaviour which occurs throughout the Animal Kingdom [23,24] and has been studied extensively in poultry (see ). Its duration (dTI) is considered to be a measure of fearfulness, representing an antipredatory behaviour. In birds, it is characterized by a strong motor inhibition, elicited by a brief period of physical restraint. The animal first struggles and attempts to escape, then gets immobile even after termination of restraint, and this freezing-like posture may last from a few seconds to several hours . Duration of tonic immobility is a quantitative polygenic trait with low to moderate heritability in poultry . It has been successfully selected upwards and downwards in experimental quail lines over a large number of generations . Despite its adaptive nature, fear may have deleterious effects on domestic animal welfare and health , and the relationship between duration of tonic immobility and productivity has been shown to vary from a positive to a negative one depending on the avian population and on the performance traits [28-30].
Deciphering the genetic architecture of tonic immobility will be challenging, as for other complex traits . We first aim to fine map a QTL associated with dTI previously mapped to quail chromosome 1 by using two F2 crosses between a Japanese quail line selected for its long tonic immobility (LTI) and two lines with shorter tonic immobility, STI (Short Tonic Immobility, selected for lower tonic immobility [12,13]) and DD (selected for early egg production [12,13]). To this end, two Advanced Intercross Lines (AIL) were built. Provided the density of information is sufficient, using advanced generations is expected to restrict the genome region involved in the variability of the trait , with the ultimate objective to identify a causative gene. We chose to obtain a high number of SNP directly from our quail populations by high-throughput sequencing (Titanium 454 GS-FLX, Roche) of two main types of reduced representations of the genome: restriction digested fractions of genomic DNA, and EST (Expressed Sequence Tag), extending a first experiment  by increasing the number of sequences obtained and detecting reliable SNP. A subset of these markers was selected to analyse the AIL populations and to fine map the QTL.
The sequences obtained in this study are available at NCBI (http://www.ncbi.nlm.nih.gov/, accession number SRP002189).
Despite the rapid reduction of sequencing costs, identifying SNP by whole-genome sequencing of many individuals in a non-model species remains expensive, notably because the reliability of the identified SNP depends on the sequence depth . Therefore, another approach was chosen  and reduced representation libraries (RRL) were built: a genomic library with tagged samples (AFLP fragments) and a transcriptomic library (cDNA fragments).
AFLP fragments were sequenced in one and a half run which produced 1,106,762 reads, totalizing about 280 Mb of sequence. Reads were mapped to the chicken genome (Figure 1) to determine if the chicken genome coverage was compatible with our aim to select quail markers evenly distributed across the genome, through the use of the chicken sequence assembly. The average number of reads per chromosome is consistent with the chicken chromosome lengths, with the exception of the Z chromosome, where more reads are observed than expected. This is probably due to the repetitive nature of a large part of the sequence of the Z chromosome . In total, 72% of the sequences mapped to the chicken genome. This result is comparable with that obtained in turkey, using a different alignment algorithm , confirming the possibility to obtain satisfactory results from alignment of heterologous sequences on a reference genome assembly. This illustrates the fact that both species have a similar level of divergence (around 35 MYA) with the chicken [20,38]. Except for some regions with high depth due to repetitive elements, the AFLP coverage was uniform enough (Figure 1b) to develop genome-wide SNP assays with confidence, given the proximity of the chicken and quail genomes .
From these data, 847,452 reads (76.6%) with identity tag and length higher than 120 bases were assembled into 42,860 contigs (N50 = 331 bp, consistent with the fraction of the AFLP fragments sequenced: between 200 and 400 bp), leaving 95,812 singletons. In the absence of a reference quail genome sequence, the SNP detection was performed from the obtained contigs. This led to the identification of 10,545 SNP.
Regarding the SNP type, a transition/transversion rate of 2.2/1 was observed. This value is similar to those reported for chicken (2.2/1), duck (2.3/1), and slightly higher than that for the great tit (1.7/1) [39,40]. Overall, our results confirm that reduced representation libraries from AFLP-derived fragments are successful tools for SNP detection. They also have been effectively used for phylogeographic studies [41-46].
In addition to the AFLP library, a second reduced representation of the genome was obtained from brain and embryo transcripts. The brain transcriptome may reflect a higher part of the genome than other tissues , and brain is probably the location of differential expression of behavioural genes between our quail lines. The embryo transcriptome was added to find SNP in genomic regions not covered by the brain transcriptome.
As previously described , a total of 399,189 cDNA reads have been obtained, of which 374,939 filtered sequences were assembled into 31,010 EST contigs, with an average length of 733 bases (from 100 to 6393, N50 = 786) and an average depth of 11 (from 2 to 3049), leaving 36,572 singletons. Among these 31,010 contigs, 27,919 (90%) mapped to the chicken genome.
Out of the 67,582 contigs or singletons built from cDNA data, 54% had significant matches (evalue < 10-5) on Swissprot, RefSeq_RNA or RefSeq_Protein databases. Sequences with more than 50% coverage with an identity level higher than 70% (n=18,416) corresponded to 8,550 unique referenced genes.
As already established for several species including Brassica rapa, Drosophila melanogaster or Homo sapiens[48-50], the transition/transversion ratio is higher in coding sequences than in untranscribed regions, and the 3.1:1 ratio in our EST data versus 2.2:1 from AFLP data confirmed this assessment. SNP were mainly discovered in non-coding region (Additional file 1: Figure S1), and 30% of the observed SNP in the coding region were non-synonymous SNP, like in chicken .
A total of 17,433 SNP were detected: 10,545 came from the genomic fragments while 6,888 were identified from the cDNA sequences, 29 of which being common between the two analyses based on their mapping on the chicken genome. SNP with at least 25 bp of flanking sequence each side have been submitted to dbSNP. This panel constitutes the first genome-wide collection of SNP markers for genetic studies in Coturnix.
The distribution of dTI was similar for the two AIL populations (Additional file 2: Figure S2). Yet, the distribution for DD-LTI was shifted slightly toward higher values of dTI, and was therefore more affected by the censoring of dTI at 300 sec. This difference probably reflects the distinct origin and history of selection of the two lines with lower dTI, DD (not selected on dTI) and STI (selected for short dTI) which were used to start the two AIL populations in crosses with Line LTI (selected for long dTI).
As previous experiments had shown the presence of a suggestive QTL on quail chromosome 1 [12,13], SNP were selected from sequences mapping to the homologous chicken chromosome 1. A total of 1,098 SNP selected from AFLP and 1,567 observed in cDNA were obtained, and 384 SNP putatively covering the entire quail chromosome 1 were chosen for genotyping. Twenty-five of them, detected from quail cDNA putatively homologous to several chicken genes, appeared to map on different chromosomes when considering their best localization (Additional file 3: Table S1). Out of the 384 markers, 352 (92%) gave reliable genotypes in our populations (AIL_DD: 326 SNP, AIL_STI: 336 SNP), and a total of 320 (91%) were polymorphic. After testing for non deviation from Hardy-Weinberg equilibrium, respectively 305 (AIL_DD) and 332 (AIL_STI) markers were used for further analyses.
A framework map of quail chromosome 1 was built, comprising 60 markers with an average number of informative meioses of 505 (from 79 to 835), for a genetic length of 287.4 cM and with an average between-marker distance of 4.9 ± 3.9 cM (Figure 2). The length of our genetic map is comparable to that (274.8 cM) obtained previously for the same chromosome by Kayang and coworkers , but is shorter than the one (402.9 cM) obtained in a backcross study . Several factors could explain this difference in length, like the structure of the populations and the number of individuals genotyped in each study, the type of markers used or the method employed to calculate the distance.
Assuming quail and chicken chromosomes 1 have similar physical length, the relative coverage for our quail genetic map was 1.5 cM/Mb, which is much lower than the 2.4 cM/Mb ratio observed in chicken . This result, probably partially related to the variability of recombination rates between populations, is consistent with the smaller number of chiasmata reported in quail lampbrush chromosomes .
AIL populations are expected to offer a very powerful tool to fine map QTL, especially when inbred lines are used as founders , and the final AIL generation is analyzed using classical QTL methods, taking advantage of the recombination events that have occurred over the course of multiple successive generations. We used this strategy with our lines to fine map a QTL previously detected in two distinct F2 crosses [12,13]. Two AIL populations were used: the AIL_DD families based upon F2 crosses of DD and LTI quail, and the AIL_STI population resulting from F2 crosses between STI and LTI birds.
First, QTL detection was applied within AIL sire families by interval mapping with QTLMap . Unexpectedly, no QTL for dTI could be detected in either of the final AIL generations, even after performing a bulk analysis using the two AIL populations. We suspected this result was the consequence of a lower detection power because the origin of alleles was not accounted for. Indeed, recent AIL studies included genotyping and phenotyping of all individuals in the AIL line , or were set up from at least one inbred line , which made it possible to trace the origin of alleles back to the F0 lines. To overcome the problem, an association analysis was performed in our AIL, with the assumption that enough power would be gained by adding information from dams within families of sires homozygous for the QTL, that were uninformative in the QTL mapping analysis based on familial linkage. This strategy was successful in the AIL_DD population (Figure 3): 4 markers, in a genomic region orthologous to position 121.54-121.78 Mb in chicken, showed a significant association with dTI (Bonferroni corrected for multiple testing p < 0.01): FSIEL0M02C8AXD_286, FSIEL0M02D00IT_258, FSIEL0M02C7RBX_180, and FSIEL0M02DMFKD_281 (Additional file 3: Table S1). Non-significant markers encompassing the region were respectively at 119.07 (F3SDQUK02B73C1_105) and 123.28 Mb (FSIEL0M01A1W48_163). As the AIL approach allows us to reasonably consider that linkage disequilibrium is limited, our investigation could be first restricted around this 4.21 Mb area. This result confirms the first mapping of a suggestive QTL in the original F2 population, as the region detected here belongs to the confidence interval estimated by the two-LOD drop-off method  from the original LTI-DD QTL position [12,13,58]. The same region also contained a significant QTL in the original LTI-STI F2 population , but had not been previously shown to be associated with dTI in chicken .
The association study showed also that the most significant marker (FSIEL0M02C8AXD_286, orthologous position 121.54 Mb on the chicken genome) had a substitution effect of 0.45 phenotypic standard deviations and that only 3 F6 sires were heterozygous at the QTL. By using these estimates, the power of detection of the QTL in the AIL_DD familial analysis was estimated to reach only 40%, which accounts for the fact that the QTL was not detected using this approach.
The same analysis in the AIL_STI population did not reveal any significant association between dTI and a marker. This may be due to the lack of linkage disequilibrium between the markers and the causative mutation in the STI population: whereas the most significant markers from the AIL_DD were fixed in the DD line (all F0 parents homozygous for the same allele), the STI population showed allele sharing, indicating that a complete association between a putative fixed mutation and a SNP marker could not be observed in the AIL_STI. The facts that the AIL were not phenotyped under the same setting (housing, technician) and that the genetic origins of the F0 were a same base line population for LTI and STI, and a different unrelated one for DD  may have also contributed to the diverging results of the association analyses.
Most of the genes recorded in the investigated region have no known association with tonic immobility (Additional file 4: Table S2), with the notable exception of the Dystrophin gene (Q9PTS7_CHICK or DMD). In the candidate region, one gene (IL1RAPL1, IL-1-receptor accessory protein-like 1, 119.54-120.16 Mb) is involved in mental retardation or autism [61,62]; structural variations including another gene (PTCHD1, patched domain containing 1, 121.89-121.94 Mb) have been shown to be involved in mental disorders like autism , but no clear connection could be made with tonic immobility. The 4 significant markers are included in 4 different genes (Additional file 4: Table S2): Q5ZLT2, or PDK3 (pyruvate dehydrogenase kinase, isoenzyme 3), a kinase involved in the regulation of glucose metabolism, has been showed to be associated with drug resistance and colon cancer ; EIF2S3 (eukaryotic translation initiation factor 2, subunit 3 gamma) is a gene coding for a subunit of eukaryotic translation initiation factor 2, involved in protein synthesis by playing a role in translation start site selection (see ); KLHL15 (kelch-like protein 15) encodes for a protein belonging to the Kelch-like proteins family, implicated in embryogenesis and carcinogenesis mainly through cytoskeleton organization ; all these three genes do not seem to have a function directly linked with the duration of tonic immobility. The fourth one, SAT1 (spermidine/spermine N1-acetyltransferase 1), catalyzing the acetylation of polyamines, has been associated with antineoplastic activity of several antitumor  and is essential for the polyamine homeostasis (see ). SAT1 has been shown in transgenic mice to be involved in several behavioural traits as activity and spatial learning in females, but not in defensive reactivity . It may constitute a correct candidate gene. But the DMD gene, localized between 118.07 and 119.07 Mb, seems to be a very interesting one even if only the 3’end of the gene is included in the candidate region. Sekiguchi et al.  showed that DMD is involved in defensive behaviors in mouse, including a strong freezing response, reminiscent of the strong motor inhibition exhibited by quail in the tonic immobility response. This gene, well known for its implication in Duchenne muscular dystrophy in Human through mutations leading to loss-of-function of the gene products in striated muscles , is also expressed in the brain, where it might account for the mental deficiency reported in Duchenne muscular dystrophy patients (see  for a review). However, further work is needed to investigate the precise involvement of brain dystrophin in the control of emotions and cognition.
Using Japanese quail as an example, this work confirms the efficiency of AFLP fragments as a reduced representation of the genome to develop SNP markers for species for which no whole-genome draft sequence is available. Transcript sequencing brings a high number of quail EST that can be used to help annotate a future quail genome assembly. This study provides also a high number of Japanese quail SNP, for genetic analyses in Coturnix. It underlines the importance of modern technologies as Next Generation Sequencing in the renewed ability to perform association analyses, even in non-model species.
The association study carried out with a SNP subpanel for chromosome 1 led to a candidate gene, DMD, as being partly responsible for the variability of tonic immobility in Japanese quail. It should be the focus of further investigations to help unravel this complex trait.
Tonic immobility was the selection criterion for the divergent selection of two quail lines: the LTI and STI lines . These two lines and a third one, line DD unselected on tonic immobility and of a different origin , were used for the study. Two AIL populations were produced: the "AIL_DD" population was obtained by crossing 6 DD male to 6 LTI females and 6 LTI males to 6 DD females to produce the F1 generation. The offspring were then intermated within each generation until the F7 generation (603 offspring, Table 1). The "AIL_STI" population was founded by crossing LTI and STI quails and carried out until the F6 generation (679 offspring, Table 1). Between 9 and 11 days of age, duration of tonic immobility was measured on each quail of the final AIL generation. TI duration is the time during which a bird remains immobile after the fear reaction has been induced by keeping it on its back in a U-shaped cradle for 10 seconds (Additional file 5: Figure S3). dTI was set to 0 second if tonic immobility was not induced after five attempts. If the quail was still in a tonic immobility state after 5 min, the observation was censored and dTI set to 300 seconds. Statistical analyses were run on the transformed variable LOG(dTI+1).
Transcripts and AFLP fragments libraries' preparation, and 454 sequencing, have already been described . Three quail lines were used: LTI, STI, DD [27,73], and total RNA was extracted from about 500 mg of adult quail brains (4 samples from each line) and total embryos (Embryonic day 8, 2 samples from each line) for library preparation. Blood samples from two individuals from each line were used for the AFLP library preparation.
Trimming of low-quality sequence and 454 adaptors were performed using the Newbler software (Roche/454).
Identification of tags within each sequence was performed by using crossmatch (from the PhredPhrap suite, http://www.phrap.org/phredphrapconsed.html, Phil Green) with tagged AFLP primers as query sequences. The crossmatch results were used to split the initial sequences when primers were observed within the fragment sequence (and not only at its ends). We then selected the sequences long enough for the clustering step (>120 bp). Sequences with 1 or 2 tags were used to generate sequence clusters, based on more than 95% identity (over a length of 80% of the shorter sequence) between cluster fragments, as processed by Megablast . Next, contigs were assembled from the clusters with the CAP3 sequence assembly program .
SNP detection was performed with an in-house script: we selected SNP with a minimum depth of 5 sequences, a minus allele count of 2 and having complete similarity between the sequences in the ten base pairs located on the left and right sides of the SNP.
Alignment of the reads to the chicken genomic sequence (Chicken May 2006 (WUGSC 2.1/galGal3) assembly) was performed with the Burrows-Wheeler Alignment algorithm bwa, version 0.6.1-r104, command bwasw, at default parameters (mismatch penalty of 3, maximum one gap opened, gap open penalty of 11, gap extension penalty of 4…). This algorithm, more permissive about alignment gaps, has been developed for long queries but was used here to deal with heterologous alignment . Genome coverage was calculated from the SAMtools v0.1.18 (mpileup tool) utilities  (samtools mpileup -A -B) and visualised through a custom script.
Artefactual sequence duplicates due to the 454 technology (98% identity, same start, same strand, same 454 run) were removed .
Because of the large number of sequences, we used a two-step process to assemble these sequences into contigs. The first step built clusters of sequences sharing at least 75 bp with an identity rate of 96% using MegaBlast . The second step constructed coherent contigs from the previous clusters using CAP3 , at the recommended stringency of 40 bp overlap with 90% sequence identity.
As no annotated quail reference genome is available, the contigs or singletons were annotated by BLAST  on the chicken genome using several databases (SwissProt, Refseq_RNA, RefSeq_Prot). Sequences with more than 50% coverage and with an identity level higher than 70% were considered as homologous to a reference gene, and their GO annotation were retrieved from the databases. All data were loaded in a locally-adapted Ensembl database.
The SNP detection was performed with the same script and parameters as for the AFLP fragments. Functional annotation of the SNP was performed with custom Perl scripts, using the Ensembl APIs (Application Programme Interface) .
In order to select SNP from contigs uniquely mapped to GGA01, a Megablast analysis  was performed to map the AFLP contigs on the chicken reference genome (Chicken May 2006 (WUGSC 2.1/galGal3) assembly). Similarly, cDNA SNP were selected from genes localized on chicken chromosome 1. Then, by using their orthologous positions on the chicken genome, 384 SNP were chosen as putatively covering the entire quail chromosome 1 (Additional file 3: Table S1). They were analysed on 1399 animals (AIL F6 or F7, their parents, and 26 F0 individuals, corresponding to the 24 AIL_DD founders and 2 LTI individuals) using the Illumina GoldenGate assay. Genotype calling was performed with GenomeStudio V2010.1 (Illumina). Genotypes were filtered for call rate (> 0.85), call frequency (> 0.50) and cluster dispersion. Markers with an unexpected deviation from Hardy-Weinberg equilibrium (p < 10-6) were removed with PLINK (option -- hardy, http://pngu.mgh.harvard.edu/purcell/plink/ ) .
A genetic map of quail chromosome 1 was built by using the CRI-MAP program  version 2.503 (revised version of Phil Green's CRI-MAP v2.4, modified by Jill Maddox and Ian Evans) with the AIL_STI population, which had a higher number of reliable markers. The "twopoint" option allowed to group markers in a unique linkage group (twopoint LOD score >= 3). The “build” option was used to order markers within the linkage group, while the “flips” option enabled to confirm the order of the markers. We build a framework map with a LOD threshold of 3, using a stepwise locus-adding strategy, starting from the markers with the higher number of informative meioses.
QTL interval mapping was performed with the QTLMap software using likelihood computation at every centimorgan, ignoring the maternal chromosomes . The linearized likelihood under the alternative hypothesis was computed as:
with i the sire family in the design; ik the progeny within sire family i; fyik/hsi, Mi the likelihood of individual ik phenotype computed conditionnally to hsi the sire i phase and Mi the marker information in the sire family i, considered as normally distributed with mean comprising some nuisance fixed effects αp, the sire family mean of the phenotypes µi, and the QTL substitution effect in sire family i aix, and standard deviation of the sire family i σi (heteroscedastic model). Pikx is the transmission probability of one of the sire chromosome to its progeny ik at position x. The technician having recorded the phenotypes, the sex and the hatch were set as nuisance fixed effects in the model. Chromosome-wide significance thresholds at 5% were calculated through 1,000 phenotype permutations. Power of detection was tested on AIL_DD modulating its two principal components, QTL frequency (ie proportion of heterozygote sires in the parental F5 generation) and QTL effect, by simulating QTL using QTLMap.
Marker-trait association was tested by using a linear mixed model:
β : fixed effects (intercept, hatch, sex, technician, genotype at marker)
a : random sire effect on trait mean a~N(0,Iσ2a)
b : random dam effect within sire on trait mean b~N(0,Iσ2b)
e : residual e~N(0,Iσ2e)
where X,Z1 and Z2 are design matrices and I denotes the identity matrix.
The model was implemented using the R « nlme » package (Linear and Nonlinear Mixed Effects Models) . The Bonferroni correction was applied to correct for multiple testing.
A search for candidate genes was performed through AnnotQTL  (http://annotqtl.genouest.org). This software gathers the functional annotation of genes mapped to a specific chromosomal region, from several public databases (Gene Ontology, Mammalian Phenotype, HGNC and Pubmed), in order to help find the best candidate genes.
The authors declare that they have no competing interests.
CB, AV, FM and FP conceived the study. FM, CD and FP participated in its design and coordination. CB, FM, SR and DG provided the AIL pedigrees. CL, FM, SR and CB supervised phenotyping. SL, KF, FV, OB and LF carried out the molecular studies and analysis. PD, CK, CC and AR performed the bioinformatic studies. SL and LF conducted SNP analyses. LF, BS and HG performed QTL and association studies. LF, FM and FP drafted the manuscript. All authors read and approved the final manuscript.
Figure S1. Distribution of the functional annotation of the detected SNP (log scale). Annotations are from the Ensembl APIs (Application Programme Interface): the functional consequence of each SNP in each transcript has been predicted using the Variant Effect Predictor (VEP). All the non intergenic consequences are represented.
Additional file 2
Figure S2. Distribution of the trait in both AIL populations. TI : Tonic Immobility.
Additional file 3
Table S1. List of the 384 SNP used in this study.
Click here for additional data file (1471-2164-13-551-S3.xlsx)
Additional file 4
Table S2. List of the genes localized in the candidate interval. Source : Biomart, Ensembl Genes 61 (Sanger UK), Gallus Gallus Genes (WASHUCS52).
Click here for additional data file (1471-2164-13-551-S4.xlsx)
Additional file 5
Figure S3. Quail in a tonic immobility state. Immobilization is induced by keeping the animal on its back in a U-shaped cradle for 10 seconds.
The authors thank Jean-Marie Hervouet and Michel Couty for their technical assistance in animal rearing, and Jill Maddox for providing a revised version of CRI-MAP. Genotyping and sequencing were performed at the GeT-PlaGe Genotoul platform, Toulouse Midi-Pyrenees. The work was supported by grants from INRA AIP Bioressources and Animal Genetics Division.
|Guyomarc'h J-C,Combreau O,Puigcerver M,Fontoura P,Aebischer N,Coturnix coturnix quailBWP UpdateYear: 199822746|
|Deregnaucourt S,Guyomarc'h J-C,Aebischer N,Hybridization between European Quail Coturnix coturnix and Japanese Quail Coturnix japonicaArdeaYear: 2002901521|
|Deregnaucourt S,Guyomarc'h J,Richard V,Classification of hybrid crows in quail using artificial neural networksBehav ProcessesYear: 20015610311210.1016/S0376-6357(01)00188-711672936|
|Wakasugi N,Japanese quailEvolution of Domesticated AnimalsYear: 1984Longman, London319321|
|Woodard AE,Abplanalp H,Wilson WO,Vohra P,Japanese quail husbandry in the laboratory (Coturnix coturnix japonica)Year: 1973Department of Avian Sciences, University of California, Davis, Ca|
|Minvielle F,Genetics and breeding of Japanese quail for production around the worldProceedings of the 6th Asian Pacific Poultry CongressYear: 1998Nagoya, Japan122127|
|Tsudzuki M,Mutations of Japanese Quail (Coturnix japonica) and Recent Advances of Molecular Genetics for This SpeciesJ Poult SciYear: 20084515917910.2141/jpsa.45.159|
|Mannen H,Murata K,Kikuchi S,Fujita D,Sasazaki S,Fujiwara A,Tsuji S,Development and Mapping of Microsatellite Markers Derived from cDNA in Japanese Quail (Coturnix japonica)J Poult SciYear: 20054226327110.2141/jpsa.42.263|
|Kayang BB,Inoue-Murayama M,Hoshi T,Matsuo K,Takahashi H,Minezawa M,Mizutani M,Ito S,Microsatellite loci in Japanese quail and cross-species amplification in chicken and guinea fowlGenet Sel EvolYear: 20023423325310.1186/1297-9686-34-2-23312081810|
|Roussot O,Feve K,Plisson-Petit F,Pitel F,Faure JM,Beaumont C,Vignal A,AFLP linkage map of the Japanese quail Coturnix japonicaGenet Sel EvolYear: 20033555957210.1186/1297-9686-35-6-55912939205|
|Kayang BB,Vignal A,Inoue-Murayama M,Miwa M,Monvoisin JL,Ito S,Minvielle F,A first-generation microsatellite linkage map of the Japanese quailAnim GenetYear: 20043519520010.1111/j.1365-2052.2004.01135.x15147390|
|Beaumont C,Roussot O,Feve K,Vignoles F,Leroux S,Pitel F,Faure JM,Mills AD,Guemene D,Sellier N,Mignon-Grasteau S,Le Roy P,Vignal A,A genome scan with AFLP markers to detect fearfulness-related QTLs in Japanese quailAnim GenetYear: 20053640140710.1111/j.1365-2052.2005.01336.x16167983|
|Minvielle F,Kayang BB,Inoue-Murayama M,Miwa M,Vignal A,Gourichon D,Neau A,Monvoisin JL,Ito S,Microsatellite mapping of QTL affecting growth, feed consumption, egg production, tonic immobility and body temperature of Japanese quailBMC GenomicsYear: 200568710.1186/1471-2164-6-8715941487|
|Minvielle F,Kayang BB,Inoue-Murayama M,Miwa M,Vignal A,Gourichon D,Neau A,Monvoisin JL,Ito S,Search for QTL affecting the shape of the egg laying curve of the Japanese quailBMC GenetYear: 200672616677378|
|Chazara O,Minvielle F,Roux D,Bed'hom B,Feve K,Coville J-L,Kayang BB,Lumineau S,Vignal A,Boutin J-M,Rognon X,Evidence for introgressive hybridization of wild common quail (Coturnix coturnix) by domesticated Japanese quail (Coturnix japonica) in FranceConserv GenetYear: 2010111051106210.1007/s10592-009-9951-8|
|Barilani M,Deregnaucourt S,Gallego S,Galli L,Mucci N,Piombo R,Puigcerver M,Rimondi S,Rodríguez-Teijeiro JD,Spanò S,Randi E,Detecting hybridization in wild (Coturnix c. coturnix) and domesticated (Coturnix c. japonica) quail populationsBiol ConservYear: 200512644545510.1016/j.biocon.2005.06.027|
|Nadeau NJ,Mundy NI,Gourichon D,Minvielle F,Association of a single-nucleotide substitution in TYRP1 with roux in Japanese quail (Coturnix japonica)Anim GenetYear: 20073860961310.1111/j.1365-2052.2007.01667.x18028514|
|Nadeau NJ,Minvielle F,Ito S,Inoue-Murayama M,Gourichon D,Follett SA,Burke T,Mundy NI,Characterization of Japanese quail yellow as a genomic deletion upstream of the avian homolog of the mammalian ASIP (agouti) geneGeneticsYear: 200817877778610.1534/genetics.107.07707318287407|
|Hiragaki T,Inoue-Murayama M,Miwa M,Fujiwara A,Mizutani M,Minvielle F,Ito S,Recessive black is allelic to the yellow plumage locus in Japanese quail and associated with a frameshift deletion in the ASIP geneGeneticsYear: 200817877177510.1534/genetics.107.07704018287406|
|Kayang BB,Fillon V,Inoue-Murayama M,Miwa M,Leroux S,Feve K,Monvoisin JL,Pitel F,Vignoles M,Mouilhayrat C,Beaumont C,Ito S,Minvielle F,Vignal A,Integrated maps in quail (Coturnix japonica) confirm the high degree of synteny conservation with chicken (Gallus gallus) despite 35 million years of divergenceBMC GenomicsYear: 2006710110.1186/1471-2164-7-10116669996|
|Shibusawa M,Minai S,Nishida-Umehara C,Suzuki T,Mano T,Yamada K,Namikawa T,Matsuda Y,A comparative cytogenetic study of chromosome homology between chicken and Japanese quailCytogenet Cell GenetYear: 20019510310910.1159/00005702611978979|
|Minvielle F,Gourichon D,Moussu C,Two new plumage mutations in the Japanese quail: "curly" feather and "rusty" plumageBMC GenetYear: 200561415760477|
|Quadros AF,Bugs PS,Araujo PB,Tonic immobility in terrestrial isopods: intraspecific and interspecific variabilityZookeysYear: 201217615517010.3897/zookeys.176.235522536106|
|Gallup GGJ,Animal hypnosis: factual status of a fictional conceptPsychol BullYear: 1974818368534612575|
|Forkman B,Boissy A,Meunier-Salaun MC,Canali E,Jones RB,A critical review of fear tests used on cattle, pigs, sheep, poultry and horsesPhysiol BehavYear: 20079234037410.1016/j.physbeh.2007.03.01618046784|
|Altan O,Settar P,Unver Y,Cabuk M,Heritabilities of Tonic Immobility and Leucocytic Response in Sire and Dam Layer LinesTurk J Vet Anim SciYear: 20052938|
|Mills AD,Faure JM,Divergent selection for duration of tonic immobility and social reinstatement behavior in Japanese quail (Coturnix coturnix japonica) chicksJ Comp PsycholYear: 199110525382032452|
|Minvielle F,Coville J-L,Krupa A,Monvoisin J-L,Maeda Y,Okamoto S,Genetic similarity and relationships of DNA fingerprints with performance and with heterosis in Japanese quail lines from two origins and under reciprocal recurrent or within-line selection for early egg productionGenet Sel EvolYear: 20003228930210.1186/1297-9686-32-3-28914736393|
|Mignon-Grasteau S,Roussot O,Delaby C,Faure JM,Mills A,Leterrier C,Guemene D,Constantin P,Mills M,Lepape G,Beaumont C,Factorial correspondence analysis of fear-related behaviour traits in Japanese quailBehav ProcessesYear: 200361697510.1016/S0376-6357(02)00162-612543484|
|Mignon-Grasteau S,Minvielle F,Relation between tonic immobility and production estimated by factorial correspondence analysis in Japanese quailPoult SciYear: 2003821839184414717540|
|Frazer KA,Murray SS,Schork NJ,Topol EJ,Human genetic variation and its contribution to complex traitsNat Rev GenetYear: 20091024125119293820|
|Darvasi A,Soller M,Advanced intercross lines, an experimental population for fine genetic mappingGeneticsYear: 1995141119912078582624|
|Leroux S,Feve K,Vignoles F,Bouchez O,Klopp C,Noirot C,Gourichon D,Richard S,Leterrier C,Beaumont C,Minvielle F,Vignal A,Pitel F,Non PCR-amplified Transcripts and AFLP fragments as reduced representations of the quail genome for 454 Titanium sequencingBMC Res NotesYear: 2010321410.1186/1756-0500-3-21420667075|
|Zhan B,Fadista J,Thomsen B,Hedegaard J,Panitz F,Bendixen C,Global assessment of genomic variation in cattle by genome resequencing and high-throughput genotypingBMC GenomicsYear: 20111255710.1186/1471-2164-12-55722082336|
|Van Tassell CP,Smith TP,Matukumalli LK,Taylor JF,Schnabel RD,Lawley CT,Haudenschild CD,Moore SS,Warren WC,Sonstegard TS,SNP discovery and allele frequency estimation by deep sequencing of reduced representation librariesNat MethodsYear: 2008524725210.1038/nmeth.118518297082|
|Bellott DW,Skaletsky H,Pyntikova T,Mardis ER,Graves T,Kremitzki C,Brown LG,Rozen S,Warren WC,Wilson RK,Page DC,Convergent evolution of chicken Z and human X chromosomes by expansion and gene acquisitionNatureYear: 201046661261610.1038/nature0917220622855|
|Kerstens HH,Crooijmans RP,Veenendaal A,Dibbits BW,Chin AWTF,den Dunnen JT,Groenen MA,Large scale single nucleotide polymorphism discovery in unsequenced genomes using second generation high throughput sequencing technology: applied to turkeyBMC GenomicsYear: 20091047910.1186/1471-2164-10-47919835600|
|Dalloul RA,Long JA,Zimin AV,Aslam L,Beal K,Blomberg Le A,Bouffard P,Burt DW,Crasta O,Crooijmans RP,Cooper K,Coulombe RA,De S,Delany ME,Dodgson JB,Dong JJ,Evans C,Frederickson KM,Flicek P,Florea L,Folkerts O,Groenen MA,Harkins TT,Herrero J,Hoffmann S,Megens HJ,Jiang A,de Jong P,Kaiser P,Kim H,Kim KW,Kim S,Langenberger D,Lee MK,Lee T,Mane S,Marcais G,Marz M,McElroy AP,Modise T,Nefedov M,Notredame C,Paton IR,Payne WS,Pertea G,Prickett D,Puiu D,Qioa D,Raineri E,Ruffier M,Salzberg SL,Schatz MC,Scheuring C,Schmidt CJ,Schroeder S,Searle SM,Smith EJ,Smith J,Sonstegard TS,Stadler PF,Tafer H,Tu ZJ,Van Tassell CP,Vilella AJ,Williams KP,Yorke JA,Zhang L,Zhang HB,Zhang X,Zhang Y,Reed KM,Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysisPLoS BiolYear: 20108|
|van Bers NE,van Oers K,Kerstens HH,Dibbits BW,Crooijmans RP,Visser ME,Groenen MA,Genome-wide SNP detection in the great tit Parus major using high throughput sequencingMol EcolYear: 201019Suppl 1899920331773|
|Kraus RH,Kerstens HH,Van Hooft P,Crooijmans RP,Van Der Poel JJ,Elmberg J,Vignal A,Huang Y,Li N,Prins HH,Groenen MA,Genome wide SNP discovery, analysis and evaluation in mallard (Anas platyrhynchos)BMC GenomicsYear: 20111215010.1186/1471-2164-12-15021410945|
|van Orsouw NJ,Hogers RC,Janssen A,Yalcin F,Snoeijers S,Verstege E,Schneiders H,van der Poel H,van Oeveren J,Verstegen H,van Eijk MJ,Complexity reduction of polymorphic sequences (CRoPS): a novel approach for large-scale polymorphism discovery in complex genomesPLoS OneYear: 20072e117210.1371/journal.pone.000117218000544|
|Gompert Z,Forister ML,Fordyce JA,Nice CC,Williamson RJ,Buerkle CA,Bayesian analysis of molecular variance in pyrosequences quantifies population genetic structure across the genome of Lycaeides butterfliesMol EcolYear: 2010192455247320497324|
|Williams LM,Ma X,Boyko AR,Bustamante CD,Oleksiak MF,SNP identification, verification, and utility for population genetics in a non-model genusBMC GenetYear: 2010113220433726|
|Baird NA,Etter PD,Atwood TS,Currey MC,Shiver AL,Lewis ZA,Selker EU,Cresko WA,Johnson EA,Rapid SNP discovery and genetic mapping using sequenced RAD markersPLoS OneYear: 20083e337610.1371/journal.pone.000337618852878|
|Mammadov JA,Chen W,Ren R,Pai R,Marchione W,Yalcin F,Witsenboer H,Greene TW,Thompson SA,Kumpatla SP,Development of highly polymorphic SNP markers from the complexity reduced portion of maize [Zea mays L.] genome for use in marker-assisted breedingTheor Appl GenetYear: 201012157758810.1007/s00122-010-1331-820401646|
|Trebbi D,Maccaferri M,de Heer P,Sorensen A,Giuliani S,Salvi S,Sanguineti MC,Massi A,van der Vossen EA,Tuberosa R,High-throughput SNP discovery and genotyping in durum wheat (Triticum durum Desf.)Theor Appl GenetYear: 201112355556910.1007/s00122-011-1607-721611761|
|Ramskold D,Wang ET,Burge CB,Sandberg R,An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence dataPLoS Comput BiolYear: 20095e100059810.1371/journal.pcbi.100059820011106|
|Moriyama EN,Powell JR,Intraspecific nuclear DNA variation in DrosophilaMol Biol EvolYear: 19961326127710.1093/oxfordjournals.molbev.a0255638583899|
|Park L,Relative mutation rates of each nucleotide for another estimated from allele frequency spectra at human gene lociGenet Res (Camb)Year: 20099129330310.1017/S001667230999016419640324|
|Park S,Yu HJ,Mun JH,Lee SC,Genome-wide discovery of DNA polymorphism in Brassica rapaMol Genet GenomicsYear: 201028313514510.1007/s00438-009-0504-020024583|
|Burt DW,White SJ,Avian genomics in the 21st centuryCytogenet Genome ResYear: 200711761310.1159/00010315917675839|
|Sasazaki S,Hinenoya T,Lin B,Fujiwara A,Mannen H,A comparative map of macrochromosomes between chicken and Japanese quail based on orthologous genesAnim GenetYear: 20063731632010.1111/j.1365-2052.2006.01454.x16879339|
|Groenen MA,Wahlberg P,Foglio M,Cheng HH,Megens HJ,Crooijmans RP,Besnier F,Lathrop M,Muir WM,Wong GK,Gut I,Andersson L,A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rateGenome ResYear: 20091951051919088305|
|Rodionov AV,Chechik MS,[Lampbrush chromosomes in the japanese quail Coturnix coturnix japonica: cytological maps of macro chromosomes and meiotic crossover frequency in females]GenetikaYear: 2002381246125112391886|
|Elsen JM,Mangin B,Goffinet B,Boichard D,Le Roy P,Alternative models for QTL detection in livestock. I. General introductionGenet Sel EvolYear: 19993121322410.1186/1297-9686-31-3-213|
|Besnier F,Wahlberg P,Ronnegard L,Ek W,Andersson L,Siegel PB,Carlborg O,Fine mapping and replication of QTL in outbred chicken advanced intercross linesGenet Sel EvolYear: 201143310.1186/1297-9686-43-321241486|
|Redmond SB,Chuammitri P,Andreasen CB,Palic D,Lamont SJ,Genetic control of chicken heterophil function in advanced intercross lines: associations with novel and with known Salmonella resistance loci and a likely mechanism for cell death in extracellular trap productionImmunogeneticsYear: 20116344945810.1007/s00251-011-0523-y21455609|
|Mangin B,Goffinet B,Comparison of several confidence intervals for QTL locationHeredityYear: 19977834535310.1038/hdy.1997.57|
|Schutz KE,Kerje S,Jacobsson L,Forkman B,Carlborg O,Andersson L,Jensen P,Major growth QTLs in fowl are related to fearful behavior: possible genetic links between fear responses and production traits in a red junglefowl x white leghorn intercrossBehav GenetYear: 20043412113014739702|
|Minvielle F,Mills AD,Faure JM,Monvoisin JL,Gourichon D,Fearfulness and performance related traits in selected lines of Japanese quail (Coturnix japonica)Poult SciYear: 20028132132611902406|
|Carrie A,Jun L,Bienvenu T,Vinet MC,McDonell N,Couvert P,Zemni R,Cardona A,Van Buggenhout G,Frints S,Hamel B,Moraine C,Ropers HH,Strom T,Howell GR,Whittaker A,Ross MT,Kahn A,Fryns JP,Beldjord C,Marynen P,Chelly J,A new member of the IL-1 receptor family highly expressed in hippocampus and involved in X-linked mental retardationNat GenetYear: 199923253110471494|
|Piton A,Michaud JL,Peng H,Aradhya S,Gauthier J,Mottron L,Champagne N,Lafreniere RG,Hamdan FF,Joober R,Fombonne E,Marineau C,Cossette P,Dube MP,Haghighi P,Drapeau P,Barker PA,Carbonetto S,Rouleau GA,Mutations in the calcium-related gene IL1RAPL1 are associated with autismHum Mol GenetYear: 2008173965397410.1093/hmg/ddn30018801879|
|Noor A,Whibley A,Marshall CR,Gianakopoulos PJ,Piton A,Carson AR,Orlic-Milacic M,Lionel AC,Sato D,Pinto D,Drmic I,Noakes C,Senman L,Zhang X,Mo R,Gauthier J,Crosbie J,Pagnamenta AT,Munson J,Estes AM,Fiebig A,Franke A,Schreiber S,Stewart AF,Roberts R,McPherson R,Guter SJ,Cook EH Jr,Dawson G,Schellenberg GD,Battaglia A,Maestrini E,Jeng L,Hutchison T,Rajcan-Separovic E,Chudley AE,Lewis SM,Liu X,Holden JJ,Fernandez B,Zwaigenbaum L,Bryson SE,Roberts W,Szatmari P,Gallagher L,Stratton MR,Gecz J,Brady AF,Schwartz CE,Schachar RJ,Monaco AP,Rouleau GA,Hui CC,Lucy Raymond F,Scherer SW,Vincent JB,Disruption at the PTCHD1 Locus on Xp22.11 in Autism spectrum disorder and intellectual disabilitySci Transl MedYear: 2010249ra6810.1126/scitranslmed.3001267|
|Lu CW,Lin SC,Chien CW,Lee CT,Lin BW,Lee JC,Tsai SJ,Overexpression of pyruvate dehydrogenase kinase 3 increases drug resistance and early recurrence in colon cancerAm J PatholYear: 20111791405141410.1016/j.ajpath.2011.05.05021763680|
|Alone PV,Cao C,Dever TE,Translation initiation factor 2gamma mutant alters start codon selection independent of Met-tRNA bindingMol Cell BiolYear: 2008286877688810.1128/MCB.01147-0818794367|
|Yoshida K,Identification and characterization of a novel kelch-like gene KLHL15 in silicoOncol RepYear: 2005131133113715870933|
|Wang Y,Devereux W,Stewart TM,Casero RA Jr,Cloning and characterization of human polyamine-modulated factor-1, a transcriptional cofactor that regulates the transcription of the spermidine/spermine N(1)-acetyltransferase geneJ Biol ChemYear: 1999274220952210110.1074/jbc.274.31.2209510419538|
|Pegg AE,Spermidine/spermine-N(1)-acetyltransferase: a key metabolic regulatorAm J Physiol Endocrinol MetabYear: 2008294E995E101010.1152/ajpendo.90217.200818349109|
|Kaasinen SK,Oksman M,Alhonen L,Tanila H,Janne J,Spermidine/spermine N1-acetyltransferase overexpression in mice induces hypoactivity and spatial learning impairmentPharmacol Biochem BehavYear: 200478354510.1016/j.pbb.2004.02.00115159132|
|Sekiguchi M,Zushida K,Yoshida M,Maekawa M,Kamichi S,Yoshida M,Sahara Y,Yuasa S,Takeda S,Wada K,A deficit of brain dystrophin impairs specific amygdala GABAergic transmission and enhances defensive behaviour in miceBrainYear: 200913212413510.1093/brain/awp04018927146|
|Hoffman EP,Brown RH Jr,Kunkel LM,Dystrophin: the protein product of the Duchenne muscular dystrophy locusCellYear: 19875191992810.1016/0092-8674(87)90579-43319190|
|Perronnet C,Vaillend C,Dystrophins, utrophins, and associated scaffolding complexes: role in mammalian brain and implications for therapeutic strategiesJ Biomed BiotechnolYear: 2010201084942620625423|
|Minvielle F,Monvoisin JL,Costa J,Frenot A,Maeda Y,Changes in heterosis under within-line selection or reciprocal recurrent selection: an experiment on early egg production in Japanese quailJ Anim Breed GenetYear: 199911636337710.1046/j.1439-0388.1999.00218.x|
|Zhang Z,Schwartz S,Wagner L,Miller W,A greedy algorithm for aligning DNA sequencesJ Comput BiolYear: 2000720321410.1089/1066527005008147810890397|
|Huang X,Madan A,CAP3: A DNA sequence assembly programGenome ResYear: 1999986887710.1101/gr.9.9.86810508846|
|Li H,Durbin R,Fast and accurate long-read alignment with Burrows-Wheeler transformBioinformaticsYear: 20102658959510.1093/bioinformatics/btp69820080505|
|Li H,Handsaker B,Wysoker A,Fennell T,Ruan J,Homer N,Marth G,Abecasis G,Durbin R,The Sequence Alignment/Map format and SAMtoolsBioinformaticsYear: 2009252078207910.1093/bioinformatics/btp35219505943|
|Mariette J,Noirot C,Klopp C,Assessment of replicate bias in 454 pyrosequencing and a multi-purpose read-filtering toolBMC Res NotesYear: 2011414910.1186/1756-0500-4-14921615897|
|Altschul SF,Gish W,Miller W,Myers EW,Lipman DJ,Basic local alignment search toolJ Mol BiolYear: 19902154034102231712|
|McLaren W,Pritchard B,Rios D,Chen Y,Flicek P,Cunningham F,Deriving the consequences of genomic variants with the Ensembl API and SNP Effect PredictorBioinformaticsYear: 2010262069207010.1093/bioinformatics/btq33020562413|
|Purcell S,Neale B,Todd-Brown K,Thomas L,Ferreira MA,Bender D,Maller J,Sklar P,de Bakker PI,Daly MJ,Sham PC,PLINK: a tool set for whole-genome association and population-based linkage analysesAm J Hum GenetYear: 20078155957510.1086/51979517701901|
|Green P,Falls K,Crooks S,Documentation for CRI-MAP, version 2.4 (3/26/90)Year: 1990Washington University School of Medicine, St Louis, MO|
|Pinheiro J,Bates D,DebRoy S,Sarkar D,Team RDC,nlme: Linear and Nonlinear Mixed Effects ModelsYear: 2011|
|Lecerf F,Bretaudeau A,Sallou O,Desert C,Blum Y,Lagarrigue S,Demeure O,AnnotQTL: a new tool to gather functional and comparative information on a genomic regionNucleic Acids ResYear: 201139W328W33310.1093/nar/gkr36121596783|
Number of animals in each generation of the AIL pedigrees
Keywords: Quail, Tonic immobility, Sequencing, AFLP, Transcripts, SNP, AIL.
Previous Document: Internet access and online cancer information seeking among Latino immigrants from safety net clinic...
Next Document: Interaction of Intraluminal Tissue and Coronary Sinus Lead Stabilized with Stent Placement.