Low-coverage sequencing in a deep intercross of the Virginia body weight lines provides insight to the polygenic genetic architecture of growth: novel loci revealed by increased power and improved genome-coverage.

Genetic dissection of highly polygenic traits is a challenge, in part due to the power necessary to confidently identify loci with minor effects. Experimental crosses are valuable resources for mapping such traits. Traditionally, genome-wide analyses of experimental crosses have targeted major loci using data from a single generation (often the F2) with individuals from later generations being generated for replication and fine-mapping. Here, we aim to confidently identify minor-effect loci contributing to the highly polygenic basis of the long-term, bi-directional selection responses for 56-d body weight in the Virginia body weight chicken lines. To achieve this, a strategy was developed to make use of data from all generations (F2-F18) of the advanced intercross line, developed by crossing the low and high selected lines after 40 generations of selection. A cost-efficient low-coverage sequencing based approach was used to obtain high-confidence genotypes in 1Mb bins across 99.3% of the chicken genome for >3,300 intercross individuals. In total, 12 genome-wide significant, and 30 additional suggestive QTL reaching a 10% FDR threshold, were mapped for 56-d body weight. Only 2 of these QTL reached genome-wide significance in earlier analyses of the F2 generation. The minor-effect QTL mapped here were generally due to an overall increase in power by integrating data across generations, with contributions from increased genome-coverage and improved marker information content. The 12 significant QTL explain >37% of the difference between the parental lines, three times more than 2 previously reported significant QTL. The 42 significant and suggestive QTL together explain >80%. Making integrated use of all available samples from multiple generations in experimental crosses are economically feasible using the low-cost, sequencing-based genotyping strategies outlined here. Our empirical results illustrate the value of this strategy for mapping novel minor-effect loci contributing to complex traits to provide a more confident, comprehensive view of the individual loci that form the genetic basis of the highly polygenic, long-term selection responses for 56-d body weight in the Virginia body weight chicken lines.

A genomic inference of the White Plymouth Rock genealogy.

Crossing of populations has been, and still is, a central component in domestication and breed and variety formation. It is a way for breeders to utilize heterosis and to introduce new genetic variation into existing plant and livestock populations. During the mid-19th century, several chicken breeds that had been introduced to America from Europe and Asia became the founders for those formed in the USA. Historical records about the genealogy of these populations are often unclear and inconsistent. Here, we used genomics in an attempt to describe the ancestry of the White Plymouth Rock (WPR) chicken. In total, 150 chickens from the WPR and 8 other stocks that historical records suggested contributed to its formation were whole-genome re-sequenced. The admixture analyses of the autosomal and sex chromosomes showed that the WPR was likely founded as a cross between a paternal lineage that was primarily Dominique, and a maternal lineage where Black Java and Cochin contributed in essentially equal proportions. These results were consistent and provided quantification with the historical records that they were the main contributors to the WPR. The genomic analyses also revealed genome-wide contributions (<10% each) by Brahma, Langshan, and Black Minorca. When viewed on an individual chromosomal basis, contributions varied considerably among stocks.

Genomic signatures of 60 years of bidirectional selection for 8-week body weight in chickens.

Sixty years, constituting 60 generations, have passed since the founding of the Virginia body weight lines, an experimental population of White Plymouth Rock chickens. Using a stringent breeding scheme for divergent 8-week body weight, the lines, which originated from a common founder population, have responded to bidirectional selection with an approximate 15-fold difference in the selected trait. They provide a model system to study the genetics of complex traits in general and the influences of artificial selection on quantitative genetic architectures in particular. As we reflect on the 60th anniversary of the initiation of the Virginia body weight lines, there is opportunity to discuss the findings obtained using different analytical and experimental genetic and genomic strategies and integrate them with a recent pooled genome resequencing dataset. Hundreds of regions across the genome show differentiation between the 2 lines, reinforcing previous findings that response to selection relied on standing variation across many genes and giving insights into the haplotype complexity underlying regions associated with body weight.

Targeted resequencing of a genomic region influencing tameness and aggression reveals multiple signals of positive selection.

The identification of the causative genetic variants in quantitative trait loci (QTL) influencing phenotypic traits is challenging, especially in crosses between outbred strains. We have previously identified several QTL influencing tameness and aggression in a cross between two lines of wild-derived, outbred rats (Rattus norvegicus) selected for their behavior towards humans. Here, we use targeted sequence capture and massively parallel sequencing of all genes in the strongest QTL in the founder animals of the cross. We identify many novel sequence variants, several of which are potentially functionally relevant. The QTL contains several regions where either the tame or the aggressive founders contain no sequence variation, and two regions where alternative haplotypes are fixed between the founders. A re-analysis of the QTL signal showed that the causative site is likely to be fixed among the tame founder animals, but that several causative alleles may segregate among the aggressive founder animals. Using a formal test for the detection of positive selection, we find 10 putative positively selected regions, some of which are close to genes known to influence behavior. Together, these results show that the QTL is probably not caused by a single selected site, but may instead represent the joint effects of several sites that were targets of polygenic selection.

Methodological aspects of the genetic dissection of gene expression.

MOTIVATION: Dissection of the genetics underlying gene expression utilizes techniques from microarray analyses as well as quantitative trait loci (QTL) mapping. Available QLT mapping methods are not tailored for the highly automated analyses required to deal with the thousand of gene transcripts encountered in the mapping of QTL affecting gene expression (sometimes referred to as eQTL). This report focuses on the adaptation of QTL mapping methodology to perform automated mapping of QTL affecting gene expression. RESULTS: The analyses of expression data on > 12,000 gene transcripts in BXD recombinant inbred mice found, on average, 629 QTL exceeding the genome-wide 5% threshold. Using additional information on trait repeatabilities and QTL location, 168 of these were classified as 'high confidence' QTL. Current sample sizes of genetical genomics studies make it possible to detect a reasonable number of QTL using simple genetic models, but considerably larger studies are needed to evaluate more complex genetic models. After extensive analyses of real data and additional simulated data (altogether > 300,000 genome scans) we make the following recommendations for detection of QTL for gene expression: (1) For populations with an unbalanced number of replicates on each genotype, weighted least squares should be preferred above ordinary least squares. Weights can be based on repeatability of the trait and the number of replicates. (2) A genome scan based on multiple marker information but analysing only at marker locations is a good approximation to a full interval mapping procedure. (3) Significance testing should be based on empirical genome-wide significance thresholds that are derived for each trait separately. (4) The significant QTL can be separated into high and low confidence QTL using a false discovery rate that incorporates prior information such as transcript repeatabilities and co-localization of gene-transcripts and QTL. (5) Including observations on the founder lines in the QTL analysis should be avoided as it inflates the test statistic and increases the Type I error. (6) To increase the computational efficiency of the study, use of parallel computing is advised. These recommendations are summarized in a possible strategy for mapping of QTL in a least squares framework. AVAILABILITY: The software used for this study is available on request from the authors.

Simultaneous search for multiple QTL using the global optimization algorithm DIRECT.

MOTIVATION: A simultaneous search is necessary for maximizing the power to detect epistatic quantitative trait loci (QTL). The computational complexity demands that the traditional exhaustive search be replaced by a more efficient global optimization algorithm. RESULTS: We have the previously known algorithm adapted DIRECT, to the problem of simultaneous mapping of multiple QTL. We have compared DIRECT with standard exhaustive search and a genetic algorithm previously used for QTL mapping in two dimensions. In all two- and three-QTL test cases, DIRECT accurately finds the global optimum two to four orders of magnitude faster than when using an exhaustive search, and one order of magnitude faster than when using the genetic algorithm. Thus, randomization testing for determining empirical significance thresholds for at least three QTL is made feasible by the use of DIRECT. AVAILABILITY: The code of the prototype implementation is available at http://user.it.uu.se/~kl/qtl_software.html

The twofold difference in adult size between the red junglefowl and White Leghorn chickens is largely explained by a limited number of QTLs.

A large intercross between the domestic White Leghorn chicken and the wild ancestor, the red junglefowl, has been used in a Quantitative Trait Loci (QTL) study of growth and egg production. The linkage map based on 105 marker loci was in good agreement with the chicken consensus map. The growth of the 851 F2 individuals was lower than both parental lines prior to 46 days of age and intermediate to the two parental lines thereafter. The QTL analysis of growth traits revealed 13 loci that showed genome-wide significance. The four major growth QTLs explained 50 and 80% of the difference in adult body weight between the founder populations for females and males, respectively. A major QTL for growth, located on chromosome 1 appears to have pleiotropic effects on feed consumption, egg production and behaviour. There was a strong positive correlation between adult body weight and average egg weight. However, three QTLs affecting average egg weight but not body weight were identified. An interesting observation was that the estimated effects for the four major growth QTLs all indicated a codominant inheritance.

Melanocortin-4 receptor (MC4R) genotypes have no major effect on fatness in a Large White x Wild Boar intercross.

The melanocortin-4 receptor (MC4R), a G-protein coupled receptor, is implicated in mediating the effect of leptin on food intake and energy balance. A previous candidate gene study reported an association between an MC4R missense mutation (Asp298Asn) and fatness, growth and feed intake in pigs. To assess this association further, we analysed the segregation of this missense mutation in relation to variation in fatness traits using a Wild Boar x Large White intercross. The Wild Boar and Large White founders were homozygous for different MC4R alleles. The MC4R was assigned to the expected region on pig chromosome 1. The statistical evaluation did not reveal any indication of a significant effect on fatness related traits in this pedigree.

Parallel computing in interval mapping of quantitative trait loci.

Linear regression analysis is considered the least computationally demanding method for mapping quantitative trait loci (QTL). However, simultaneous search for multiple QTL, the use of permutations to obtain empirical significance thresholds, and larger experimental studies significantly increase the computational demand. This report describes an easily implemented parallel algorithm, which significantly reduces the computing time in both QTL mapping and permutation testing. In the example provided, the analysis time was decreased to less than 15% of a single processor system by the use of 18 processors. We indicate how the efficiency of the analysis could be improved by distributing the computations more evenly to the processors and how other ways of distributing the data facilitate the use of more processors. The use of parallel computing in QTL mapping makes it possible to routinely use permutations to obtain empirical significance thresholds for multiple traits and multiple QTL models. It could also be of use to improve the computational efficiency of the more computationally demanding QTL analysis methods.

The use of a genetic algorithm for simultaneous mapping of multiple interacting quantitative trait loci.

Here we describe a general method for improving computational efficiency in simultaneous mapping of multiple interacting quantitative trait loci (QTL). The method uses a genetic algorithm to search for QTL in the genome instead of an exhaustive enumerative ("step-by-step") search. It can be used together with any method of QTL mapping based on a genomic search, since it only provides a more efficient way to search the genome for QTL. The computational demand decreases by a factor of approximately 130 when using genetic algorithm-based mapping instead of an exhaustive enumerative search for two QTL in a genome size of 2000 cM using a resolution of 1 cM. The advantage of using a genetic algorithm increases further for larger genomes, higher resolutions, and searches for more QTL. We show that a genetic algorithm-based search has efficiency higher than or equal to a search method conditioned on previously identified QTL for all epistatic models tested and that this efficiency is comparable to that of an exhaustive search for multiple QTL. The genetic algorithm is thus a powerful and computationally tractable alternative to the exhaustive enumerative search for simultaneous mapping of multiple interacting QTL. The use of genetic algorithms for simultaneous mapping of more than two QTL and for determining empirical significance thresholds using permutation tests is also discussed.

The origin of the domestic pig: independent domestication and subsequent introgression.

The domestic pig originates from the Eurasian wild boar (Sus scrofa). We have sequenced mitochondrial DNA and nuclear genes from wild and domestic pigs from Asia and Europe. Clear evidence was obtained for domestication to have occurred independently from wild boar subspecies in Europe and Asia. The time since divergence of the ancestral forms was estimated at approximately 500,000 years, well before domestication approximately 9,000 years ago. Historical records indicate that Asian pigs were introduced into Europe during the 18th and early 19th centuries. We found molecular evidence for this introgression and the data indicated a hybrid origin of some major "European" pig breeds. The study is an advance in pig genetics and has important implications for the maintenance and utilization of genetic diversity in this livestock species.

High in situ rat intestinal permeability of artemisinin unaffected by multiple dosing and with no evidence of P-glycoprotein involvement.

The objective of this study was to investigate whether the decrease in artemisinin bioavailability after repeated oral dosing in humans can be a result of increased efflux of artemisinin by P-glycoprotein or decreased membrane transport at the intestinal barrier. The effective jejunal permeability (Peff) of artemisinin was investigated using an in situ rat perfusion model. Fifty-four rats were randomized to one of three treatment arms: no pretreatment, pretreatment with artemisinin emulsion for 5 days (60 mg/kg/day, p.o. ), or pretreatment with emulsion vehicle for 5 days. The rats within each treatment arm were randomized further to be jejunally perfused with either low (500 ng/ml) or high (5000 ng/ml) artemisinin concentration or low artemisinin concentration plus the P-glycoprotein inhibitor R,S-verapamil (400 microg/ml). Perfusate samples were assayed for content of artemisinin, R,S-verapamil, and perfusion viability markers. Artemisinin Peff was 1.44 +/- 0.38, 1. 17 +/- 0.32, and 1.71 +/- 0.29 (.10(-4), cm/s) in rats receiving no pretreatment and perfused with low, high, or low artemisinin concentration plus verapamil, respectively. Multiple oral dosing of artemisinin did not affect the jejunal permeability of artemisinin. R,S-verapamil Peff was similar in artemisinin-pretreated rats (1.09 +/- 0.54. 10(-4), cm/s) and rats pretreated with only vehicle (1.07 +/- 0.37. 10(-4), cm/s). The decrease in artemisinin bioavailability after multiple oral dosing in human is probably not a result of changes in P-glycoprotein expression or general intestinal transport. It seems more likely attributed to increased hepatocellular activity. Furthermore, artemisinin exhibits high jejunal permeability and is neither a substrate nor inducer of P-glycoprotein.