A quick method for computing approximate thresholds for quantitative trait loci detection

Mapping QTL controlling count traits with excess zeros and ones using a zero-and-one-inflated generalized Poisson regression model

The research on the quantitative trait locus (QTL) mapping of count data has aroused the wide attention of researchers. There are frequent problems in applied research that limit the application of the conventional Poisson model in the analysis of count phenotypes, which include the overdispersion and excess zeros and ones. In this article, a novel model, that is, the zero-and-one-inflated generalized Poisson (ZOIGP) model, is proposed to deal with these problems. Based on the proposed model, a score test is performed for the inflation parameter, in which the ZOIGP model with a constant proportion of excess zeros and ones is compared with a standard generalized Poisson model. To illustrate the practicability of the ZOIGP model, we extend it to the QTL interval mapping application that underpins count phenotype with excess zeros and excess ones. The genetic effects are estimated utilizing the expectation-maximization algorithm embedded with the Newton-Raphson algorithm, and the genome-wide scan and likelihood ratio test is performed to map and test the potential QTLs. The statistical properties exhibited by the proposed method are investigated through simulation. Finally, a real data analysis example is used to illustrate the utility of the proposed method for QTL mapping.

Interpretation of Manhattan Plots and Other Outputs of Genome-Wide Association Studies

With increasing marker density, estimation of recombination rate between a marker and a causal mutation using linkage analysis becomes less important. Instead, linkage disequilibrium (LD) becomes the major indicator for gene mapping through genome-wide association studies (GWAS). In addition to the linkage between the marker and the causal mutation, many other factors may contribute to the LD, including population structure and cryptic relationships among individuals. As statistical methods and software evolve to improve statistical power and computing speed in GWAS, the corresponding outputs must also evolve to facilitate the interpretation of input data, the analytical process, and final association results. In this chapter, our descriptions focus on (1) considerations in creating a Manhattan plot displaying the strength of LD and locations of markers across a genome; (2) criteria for genome-wide significance threshold and the different appearance of Manhattan plots in single-locus and multiple-locus models; (3) exploration of population structure and kinship among individuals; (4) quantile-quantile (QQ) plot; (5) LD decay across the genome and LD between the associated markers and their neighbors; (6) exploration of individual and marker information on Manhattan and QQ plots via interactive visualization using HTML. The ultimate objective of this chapter is to help users to connect input data to GWAS outputs to balance power and false positives, and connect GWAS outputs to the selection of candidate genes using LD extent.

Overt and concealed genetic loads revealed by QTL mapping of genotype-dependent viability in the Pacific oyster Crassostrea gigas

Understanding the genetic bases of inbreeding depression, heterosis, and genetic load is integral to understanding how genetic diversity is maintained in natural populations. The Pacific oyster Crassostrea gigas, like many long-lived plants, has high fecundity and high early mortality (type-III survivorship), manifesting a large, overt, genetic load; the oyster harbors an even greater concealed genetic load revealed by inbreeding. Here, we map viability QTL (vQTL) in six interrelated F2 oyster families, using high-density linkage maps of single nucleotide polymorphisms generated by genotyping-by-sequencing (GBS) methods. Altogether, we detect 70 vQTL and provisionally infer 89 causal mutations, 11 to 20 per family. Genetic mortality caused by independent (unlinked) vQTL ranges from 94.2% to 97.8% across families, consistent with previous reports. High-density maps provide better resolution of genetic mechanisms, however. Models of one causal mutation present in both identical-by-descent (IBD) homozygotes and heterozygotes fit genotype frequencies at 37 vQTL; consistent with the mutation-selection balance theory of genetic load, 20 are highly deleterious, completely recessive mutations and 17 are less deleterious, partially dominant mutations. Another 22 vQTL require pairs of recessive or partially dominant causal mutations, half showing selection against recessive mutations linked in repulsion, producing pseudo-overdominance. Only eight vQTL appear to support the overdominance theory of genetic load, with deficiencies of both IBD homozygotes, but at least four of these are likely caused by pseudo-overdominance. Evidence for epistasis is absent. A high mutation rate, random genetic drift, and pseudo-overdominance may explain both the oyster's extremely high genetic diversity and a high genetic load maintained primarily by mutation-selection balance.

Genome survey and high-resolution genetic map provide valuable genetic resources for Fenneropenaeus chinensis

Fenneropenaeus chinensis is one of the most important aquaculture species in China. Research on its genomic and genetic structure not only helps us comprehend the genetic basis of complex economic traits, but also offers theoretical guidance in selective breeding. In the present study, a genome survey sequencing was performed to generate a rough reference genome utilized for groping preliminary genome characteristics and facilitate linkage and quantitative trait locus (QTL) mapping. Linkage mapping was conducted using a reduced-representation sequencing method 2b-RAD. In total, 36,762 SNPs were genotyped from 273 progenies in a mapping family, and a high-resolution linkage map was constructed. The consensus map contained 12,884 markers and spanned 5257.81 cM with an average marker interval of 0.41 cM, which was the first high-resolution genetic map in F. chinensis to our knowledge. QTL mapping and association analysis were carried out in 29 characters including body size, sex and disease resistance. 87 significant QTLs were detected in several traits and they were also evaluated by association analysis. Results of this study provide us valuable suggestions in genetic improvement and breeding of new varieties and also lay a basic foundation for further application of cloning of economic genes in selective breeding program and marker-assisted selection.

Genetic Determinism Exists for the Global DNA Methylation Rate in Sheep

Recent studies showed that epigenetic marks, including DNA methylation, influence production and adaptive traits in plants and animals. So far, most studies dealing with genetics and epigenetics considered DNA methylation sites independently. However, the genetic basis of the global DNA methylation rate (GDMR) remains unknown. The main objective of the present study was to investigate genetic determinism of GDMR in sheep. The experiment was conducted on 1,047 Romane sheep allocated into 10 half-sib families. After weaning, all the lambs were phenotyped for global GDMR in blood as well as for production and adaptive traits. GDMR was measured by LUminometric Methylation Analysis (LUMA) using a pyrosequencing approach. Association analyses were conducted on some of the lambs (n = 775) genotyped by using the Illumina OvineSNP50 BeadChip. Blood GDMR varied among the animals (average 70.7 ± 6.0%). Female lambs had significantly higher GDMR than male lambs. Inter-individual variability of blood GDMR had an additive genetic component and heritability was moderate (h² = 0.20 ± 0.05). No significant genetic correlation was found between GDMR and growth or carcass traits, birthcoat, or social behaviors. Association analyses revealed 28 QTLs associated with blood GDMR. Seven genomic regions on chromosomes 1, 5, 11, 17, 24, and 26 were of most interest due to either high significant associations with GDMR or to the relevance of genes located close to the QTLs. QTL effects were moderate. Genomic regions associated with GDMR harbored several genes not yet described as being involved in DNA methylation, but some are already known to play an active role in gene expression. In addition, some candidate genes, CHD1, NCO3A, KDM8, KAT7, and KAT6A have previously been described to be involved in epigenetic modifications. In conclusion, the results of the present study indicate that blood GDMR in domestic sheep is under polygenic influence and provide new insights into DNA methylation genetic determinism.

Quick approximation of threshold values for genome-wide association studies

Standard normal statistics, chi-squared statistics, Student's t statistics and F statistics are used to map quantitative trait nucleotides for both small and large sample sizes. In genome-wide association studies (GWASs) of single-nucleotide polymorphisms (SNPs), the statistical distributions depend on both genetic effects and SNPs but are independent of SNPs under the null hypothesis of no genetic effects. Therefore, hypothesis testing when a nuisance parameter is present only under the alternative was introduced to quickly approximate the critical thresholds of these test statistics for GWASs. When only the statistical probabilities are available for high-throughput SNPs, the approximate critical thresholds can be estimated with chi-squared statistics, formulated by statistical probabilities with a degree of freedom of two. High similarities in the critical thresholds between the accurate and approximate estimations were demonstrated by extensive simulations and real data analysis.

Mapping Stem Rust (Puccinia graminis f. sp. secalis) Resistance in Self-Fertile Winter Rye Populations

Rye stem rust caused by Puccinia graminis f. sp. secalis can be found in all European rye growing regions. When the summers are warm and dry, the disease can cause severe yield losses over large areas. To date only little research was done in Europe to trigger resistance breeding. To our knowledge, all varieties currently registered in Germany are susceptible. In this study, three biparental populations of inbred lines and one testcross population developed for mapping resistance were investigated. Over 2 years, 68-70 genotypes per population were tested, each in three locations. Combining the phenotypic data with genotyping results of a custom 10k Infinium iSelect single nucleotide polymorphism (SNP) array, we identified both quantitatively inherited adult plant resistance and monogenic all-stage resistance. A single resistance gene, tentatively named Pgs1, located at the distal end of chromosome 7R, could be identified in two independently developed populations. With high probability, it is closely linked to a nucleotide-binding leucine-rich repeat (NB-LRR) resistance gene homolog. A marker for a competitive allele-specific polymerase chain reaction (KASP) genotyping assay was designed that could explain 73 and 97% of the genetic variance in each of both populations, respectively. Additional investigation of naturally occurring rye leaf rust (caused by Puccinia recondita ROEBERGE) revealed a gene complex on chromosome 7R. The gene Pgs1 and further identified quantitative trait loci (QTL) have high potential to be used for breeding stem rust resistant rye.

Fine Mapping of the High-pH Tolerance and Growth Trait-Related Quantitative Trait Loci (QTLs) and Identification of the Candidate Genes in Pacific White Shrimp (Litopenaeus vannamei)

High-pH tolerance and growth are important traits for the shrimp culture industry in areas with saline-alkali water. In the present study, an F1 full-sib family of Pacific white shrimp (Litopenaeus vannamei) was generated with a new "semidirectional cross" method, and double-digest restriction site-associated DNA sequencing (ddRAD-Seq) technology was applied to genotype the 2 parents and 148 progenies. A total of 3567 high-quality markers were constructed for the genetic linkage map, and the total map length was 4161.555 centimorgans (cM), showing 48 linkage groups (LGs) with an average interlocus length of 1.167 cM. With a constrained logarithm of odds (LOD) score ≥ 2.50, 12 high-pH tolerance and 2 growth (body weight) QTLs were located. L. vannamei genomic scaffolds were used to assist with the detection of 21 stress- and 5 growth-related scaffold genes. According to the high-pH transcriptome data of our previous study, 6 candidate high-pH response genes were discovered, and 5 of these 6 genes were consistently expressed with the high-pH transcriptome data, validating the locations of the high-pH tolerance trait-related QTLs in this study. This paper is the first report of fine-mapping high-pH tolerance and growth (body weight) trait QTLs in one L. vannamei genetic map. Our results will further benefit marker-assisted selection work and might be useful for promoting genomic research on the shrimp L. vannamei.

New locus reveals the genetic architecture of sex reversal in the Chinese tongue sole (Cynoglossus semilaevis)

Sex reversal in insects, amphibians, reptiles, and fishes is a complicated and interesting biological phenomenon. Sex reversal changes the sex ratio of populations and may complicate breeding schemes. In the Chinese tongue sole (Cynoglossus semilaevis), genetic females may change into pseudomales, thereby increasing aquaculture costs because of the lower growth rate of the males than that of the females. Here we identify a new locus associated with sex reversal; this single nucleotide polymorphism (SNP) is located in the third intron of the doublesex and mab-3 related transcription factor 1 (Dmrt1) gene on the Z chromosome (named Cyn_Z_8564889) and has two alleles, A and G. Cyn_Z_8564889 regulates sex reversal interactively with our previously detected SNP (Cyn_Z_6676874), with the genetic females simultaneously carrying the T allele of Cyn_Z_6676874 and the A allele of Cyn_Z_8564889 changing into pseudomales. Other Dmrt1 polymorphisms were detected, which formed two haplotypes. Two SNPs in the second exon of Dmrt1 result in amino acid changes, suggesting that Dmrt1 is essential in sex reversal. We also verified that pseudomales produce no or little W sperm. The interaction and linkage between Cyn_Z_6676874 and Cyn_Z_8564889 and the absence of W sperm from pseudomales unravel the genetic architecture of sex reversal in C. semilaevis.

Identification of quantitative trait loci for the fatty acid composition in Korean native chicken

OBJECTIVE

Fatty acid composition is one of the most important meat quality traits because it can contribute to functional, sensorial, and nutritional factors. In this study, quantitative trait locus (QTL) analyses for fatty acid composition traits were investigated in thigh and breast meat of Korean native chicken (KNC).

METHODS

In total, 18 fatty acid composition traits were investigated from each meat sample using 83 parents, and 595 F1 chicks of 20 week old. Genotype assessment was performed using 171 informative DNA markers on 26 autosomes. The KNC linkage map was constructed by CRI-MAP software, which calculated genetic distances, with map orders between markers. The half-sib and full-sib QTL analyses were performed using GridQTL and SOLAR programs, respectively.

RESULTS

In total, 30 QTLs (12 in the thigh and 18 in the breast meat) were detected by the half-sib analysis and 7 QTLs (3 in the thigh and 4 in the breast meat) were identified by the full-sib analysis.

CONCLUSION

With further verification of the QTL regions using additional markers and positional candidate gene studies, these results can provide valuable information for determining causative mutations affecting the fatty acid composition of KNC meat. Moreover, these findings may aid in the selection of birds with favorable fatty acid composition traits.

Single Locus Maintains Large Variation of Sex Reversal in Half-Smooth Tongue Sole (Cynoglossus semilaevis)

Sex determination is a fundamental biological process for individual sex development and population sex ratios. However, for some species, the primary sex might be altered during development, and individuals can develop into the opposite sex. Sex reversal may happen in insects, reptiles, amphibians, and fishes. In half-smooth tongue sole (Cynoglossus semilaevis), some genetically female fish irreversibly reverse to pseudomales, resulting in higher costs in aquaculture owing to a lower growth rate of male fish during a 2-yr growth period. Here, we identified a locus with large controlling effect on sex reversal in the half-smooth tongue sole through genome-wide association study with high-density single nucleotide polymorphisms (SNPs). This SNP is located at the third intron of the F-box and leucine rich repeat protein 17 (FBXL17) gene on the Z chromosome, and it has two alleles, A and T. Genetic females with Z^AW genotypes will never reverse into phenotypic males, but those with Z^TW genotypes can sometimes undergo sex reversal. This SNP explains 82.7% of the genetic variation, or 58.4% of the phenotypic variation. Based on our results, a reproductive management program could be developed to improve the phenotypic female ratio in aquaculture, and elucidate the mechanism of sex reversal in half-smooth tongue sole. We expect that these findings will have a substantial impact on the population management in many harvested species where sex reversal occurs.

A Major Locus for Quantitatively Measured Shank Skin Color Traits in Korean Native Chicken

Shank skin color of Korean native chicken (KNC) shows large color variations. It varies from white, yellow, green, bluish or grey to black, whilst in the majority of European breeds the shanks are typically yellow-colored. Three shank skin color-related traits (i.e., lightness [L*], redness [a*], and yellowness [b*]) were measured by a spectrophotometer in 585 progeny from 68 nuclear families in the KNC resource population. We performed genome scan linkage analysis to identify loci that affect quantitatively measured shank skin color traits in KNC. All these birds were genotyped with 167 DNA markers located throughout the 26 autosomes. The SOLAR program was used to conduct multipoint variance-component quantitative trait locus (QTL) analyses. We detected a major QTL that affects b* value (logarithm of odds [LOD] = 47.5, p = 1.60×10⁻⁴⁹) on GGA24 (GGA for Gallus gallus). At the same location, we also detected a QTL that influences a* value (LOD = 14.2, p = 6.14×10⁻¹⁶). Additionally, beta-carotene dioxygenase 2 (BCDO2), the obvious positional candidate gene under the linkage peaks on GGA24, was investigated by the two association tests: i.e., measured genotype association (MGA) and quantitative transmission disequilibrium test (QTDT). Significant associations were detected between BCDO2 g.9367 A>C and a* (P_MGA = 1.69×10⁻²⁸; P_QTDT = 2.40×10⁻²⁵). The strongest associations were between BCDO2 g.9367 A>C and b* (P_MGA = 3.56×10⁻⁶⁶; P_QTDT = 1.68×10⁻⁶⁵). However, linkage analyses conditional on the single nucleotide polymorphism indicated that other functional variants should exist. Taken together, we demonstrate for the first time the linkage and association between the BCDO2 locus on GGA24 and quantitatively measured shank skin color traits in KNC.

High-resolution genetic linkage mapping, high-temperature tolerance and growth-related quantitative trait locus (QTL) identification in Marsupenaeus japonicus

The Kuruma prawn, Marsupenaeus japonicus, is one of the most promising marine invertebrates in the industry in Asia, Europe and Australia. However, the increasing global temperatures result in considerable economic losses in M. japonicus farming. In the present study, to select genetically improved animals for the sustainable development of the Kuruma prawn industry, a high-resolution genetic linkage map and quantitative trait locus (QTL) identification were performed using the RAD technology. The maternal map contained 5849 SNP markers and spanned 3127.23 cM, with an average marker interval of 0.535 cM. Instead, the paternal map contained 3927 SNP markers and spanned 3326.19 cM, with an average marker interval of 0.847 cM. The consensus map contained 9289 SNP markers and spanned 3610.90 cM, with an average marker interval of 0.388 cM and coverage of 99.06 % of the genome. The markers were grouped into 41 linkage groups in the maps. Significantly, negative correlation was detected between high-temperature tolerance (UTT) and body weight (BW). The QTL mapping revealed 129 significant QTL loci for UTT and four significant QTL loci for BW at the genome-wide significance threshold. Among these QTLs, 129 overlapped with linked SNPs, and the remaining four were located in regions between contiguous SNPs. They explained the total phenotypic variance ranging from 8.9 to 12.4 %. Because of a significantly negative correlation between growth and high-temperature tolerance, we demonstrate that this high-resolution linkage map and QTLs would be useful for further marker-assisted selection in the genetic improvement of M. japonicus.

Genetic inviability is a major driver of type III survivorship in experimental families of a highly fecund marine bivalve

The offspring of most highly fecund marine fish and shellfish suffer substantial mortality early in the life cycle, complicating prediction of recruitment and fisheries management. Early mortality has long been attributed to environmental factors and almost never to genetic sources. Previous work on a variety of marine bivalve species uncovered substantial genetic inviability among the offspring of inbred crosses, suggesting a large load of early-acting deleterious recessive mutations. However, genetic inviability of randomly bred offspring has not been addressed. Here, genome-wide surveys reveal widespread, genotype-dependent mortality in randomly bred, full-sib progenies of wild-caught Pacific oysters (Crassostrea gigas). Using gene-mapping methods, we infer that 11-19 detrimental alleles per family render 97.9-99.8% of progeny inviable. The variable genomic positions of viability loci among families imply a surprisingly large load of partially dominant or additive detrimental mutations in wild adult oysters. Although caution is required in interpreting the relevance of experimental results for natural field environments, we argue that the observed genetic inviability corresponds with type III survivorship, which is characteristic of both hatchery and field environments and that our results, therefore, suggest the need for additional experiments under the near-natural conditions of mesocosms. We explore the population genetic implications of our results, calculating a detrimental mutation rate that is comparable to that estimated for conifers and other highly fecund perennial plants. Genetic inviability ought to be considered as a potential major source of low and variable recruitment in highly fecund marine animals.

Variance Component Quantitative Trait Locus Analysis for Body Weight Traits in Purebred Korean Native Chicken

Quantitative trait locus (QTL) is a particular region of the genome containing one or more genes associated with economically important quantitative traits. This study was conducted to identify QTL regions for body weight and growth traits in purebred Korean native chicken (KNC). F₁ samples (n = 595) were genotyped using 127 microsatellite markers and 8 single nucleotide polymorphisms that covered 2,616.1 centi Morgan (cM) of map length for 26 autosomal linkage groups. Body weight traits were measured every 2 weeks from hatch to 20 weeks of age. Weight of half carcass was also collected together with growth rate. A multipoint variance component linkage approach was used to identify QTLs for the body weight traits. Two significant QTLs for growth were identified on chicken chromosome 3 (GGA3) for growth 16 to18 weeks (logarithm of the odds [LOD] = 3.24, Nominal p value = 0.0001) and GGA4 for growth 6 to 8 weeks (LOD = 2.88, Nominal p value = 0.0003). Additionally, one significant QTL and three suggestive QTLs were detected for body weight traits in KNC; significant QTL for body weight at 4 weeks (LOD = 2.52, nominal p value = 0.0007) and suggestive QTL for 8 weeks (LOD = 1.96, Nominal p value = 0.0027) were detected on GGA4; QTLs were also detected for two different body weight traits: body weight at 16 weeks on GGA3 and body weight at 18 weeks on GGA19. Additionally, two suggestive QTLs for carcass weight were detected at 0 and 70 cM on GGA19. In conclusion, the current study identified several significant and suggestive QTLs that affect growth related traits in a unique resource pedigree in purebred KNC. This information will contribute to improving the body weight traits in native chicken breeds, especially for the Asian native chicken breeds.

Inclusive Composite Interval Mapping of QTL by Environment Interactions in Biparental Populations

Identification of environment-specific QTL and stable QTL having consistent genetic effects across a wide range of environments is of great importance in plant breeding. Inclusive Composite Interval Mapping (ICIM) has been proposed for additive, dominant and epistatic QTL mapping in biparental populations for single environment. In this study, ICIM was extended to QTL by environment interaction (QEI) mapping for multi-environmental trials, where the QTL average effect and QEI effects could be properly estimated. Stepwise regression was firstly applied in each environment to identify the most significant marker variables which were then used to adjust the phenotypic values. One-dimensional scanning was then conducted on the adjusted phenotypic values across the environments in order to detect QTL with either average effect or QEI effects, or both average effect and QEI effects. In this way, the genetic background could be well controlled while the conventional interval mapping was applied. An empirical method to determine the threshold of logarithm of odds was developed, and the efficiency of the ICIM QEI mapping was demonstrated in simulated populations under different genetic models. One actual recombinant inbred line population was used to compare mapping results between QEI mapping and single-environment analysis.

High-density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis

The sex determination system in crabs is believed to be XY-XX from karyotypy, but centromeres could not be identified in some chromosomes and their morphology is not completely clear. Using quantitative trait locus mapping of the gender phenotype, we revealed a ZW-ZZ sex determination system in Eriocheir sinensis and presented a high-density linkage map covering ~98.5% of the genome, with 73 linkage groups corresponding to the haploid chromosome number. All sex-linked markers in the family we used were located on a single linkage group, LG60, and sex linkage was confirmed by genome-wide association studies (GWAS). Forty-six markers detected by GWAS were heterozygous and segregated only in the female parent. The female LG60 was thus the putative W chromosome, with the homologous male LG60 as the Z chromosome. The putative Z and W sex chromosomes were identical in size and carried many homologous loci. Sex ratio (5:1) skewing towards females in induced triploids using unrelated animals also supported a ZW-ZZ system. Transcriptome data were used to search for candidate sex-determining loci, but only one LG60 gene was identified as an ankyrin-2 gene. Double sex- and mab3-related transcription factor 1 (Dmrt1), a Z-linked gene in birds, was located on a putative autosome. With complete genome sequencing and transcriptomic data, more genes on putative sex chromosomes will be characterised, thus leading towards a comprehensive understanding of the sex determination and differentiation mechanisms of E. sinensis, and decapod crustaceans in general.

Identification of QTLs for behavioral reactivity to social separation and humans in sheep using the OvineSNP50 BeadChip

Background

Current trends in sheep farming practices rely on animals with a greater level of behavioral autonomy than before, a phenotype that actively contributes to the sustainability of animal production. Social reactivity and reactivity to humans are relevant behavioral traits in sheep, known for their strong gregariousness and weak tolerance to handling, which have previously been reported with moderate to high heritabilities. To identify loci underlying such behaviors, we performed a genome study in Romane lambs.

Results

The experiment was carried out on 934 male and female lambs allocated into 9 half-sib families (average of 103 lambs per family) and reared outside. After weaning, all the lambs were individually exposed to 4 standardized behavioral tests combining social isolation, exposure to humans or handling, confinement and novelty (i.e. arena test, corridor test, isolation box test, shearing test). A broad range of behaviors including vocalizations, locomotion, vigilance and flight distance, as well as the cortisol response to handling, were collected. All lambs were genotyped using the Illumina OvineSNP50 BeadChip. QTL detection was performed by linkage, association and joint linkage and association analyses using the QTLmap software. Five main QTL regions were identified on sheep chromosomes (Ovis Aries Region, OAR) 12, 16, 19, 21 and 23 among many other QTLs with small to moderate effects. The QTLs on OAR12, 16 and 21 showed significant associations with social reactivity. The QTLs on OAR19 and 23 were found to be associated with reactivity to humans. No overlapping QTLs were identified for the different traits measured in the behavioral tests, supporting the hypothesis that different genetic factors influence social reactivity and tolerance to humans.

Conclusion

The results of this study using ovine SNP data suggest that in domestic sheep the behavioral responses to social separation and exposure to humans are under polygenic influence. The most relevant QTLs reported in the present study contain interesting candidate genes previously described to be associated with various emotional and social behaviors in mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-778) contains supplementary material, which is available to authorized users.

RWGAIM: an efficient high-dimensional random whole genome average (QTL) interval mapping approach

Mapping of quantitative trait loci (QTLs) underlying variation in quantitative traits continues to be a powerful tool in genetic study of plants and other organisms. Whole genome average interval mapping (WGAIM), a mixed model QTL mapping approach using all intervals or markers simultaneously, has been demonstrated to outperform composite interval mapping, a common approach for QTL analysis. However, the advent of high-throughput high-dimensional marker platforms provides a challenge. To overcome this, a dimension reduction technique is proposed for WGAIM for efficient analysis of a large number of markers. This approach results in reduced computing time as it is dependent on the number of genetic lines (or individuals) rather than the number of intervals (or markers). The approach allows for the full set of potential QTL effects to be recovered. A proposed random effects version of WGAIM aims to reduce bias in the estimated size of QTL effects. Lastly, the two-stage outlier procedure used in WGAIM is replaced by a single stage approach to reduce possible bias in the selection of putative QTL in both WGAIM and the random effects version. Simulation is used to demonstrate the efficiency of the dimension reduction approach as well as demonstrate that while the approaches are very similar, the random WGAIM performs better than the original and modified fixed WGAIM by reducing bias and in terms of mean square error of prediction of estimated QTL effects. Finally, an analysis of a doubled haploid population is used to illustrate the three approaches.

High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc

Genetic linkage maps are indispensable tools in genetic and genomic studies. Recent development of genotyping-by-sequencing (GBS) methods holds great promise for constructing high-resolution linkage maps in organisms lacking extensive genomic resources. In the present study, linkage mapping was conducted for a bivalve mollusc (Chlamys farreri) using a newly developed GBS method-2b-restriction site-associated DNA (2b-RAD). Genome survey sequencing was performed to generate a preliminary reference genome that was utilized to facilitate linkage and quantitative trait locus (QTL) mapping in C. farreri. A high-resolution linkage map was constructed with a marker density (3806) that has, to our knowledge, never been achieved in any other molluscs. The linkage map covered nearly the whole genome (99.5%) with a resolution of 0.41 cM. QTL mapping and association analysis congruously revealed two growth-related QTLs and one potential sex-determination region. An important candidate QTL gene named PROP1, which functions in the regulation of growth hormone production in vertebrates, was identified from the growth-related QTL region detected on the linkage group LG3. We demonstrate that this linkage map can serve as an important platform for improving genome assembly and unifying multiple genomic resources. Our study, therefore, exemplifies how to build up an integrative genomic framework in a non-model organism.

Selecting informative traits for multivariate quantitative trait locus mapping helps to gain optimal power

A major consideration in multitrait analysis is which traits should be jointly analyzed. As a common strategy, multitrait analysis is performed either on pairs of traits or on all of traits. To fully exploit the power of multitrait analysis, we propose variable selection to choose a subset of informative traits for multitrait quantitative trait locus (QTL) mapping. The proposed method is very useful for achieving optimal statistical power for QTL identification and for disclosing the most relevant traits. It is also a practical strategy to effectively take advantage of multitrait analysis when the number of traits under consideration is too large, making the usual multivariate analysis of all traits challenging. We study the impact of selection bias and the usage of permutation tests in the context of variable selection and develop a powerful implementation procedure of variable selection for genome scanning. We demonstrate the proposed method and selection procedure in a backcross population, using both simulated and real data. The extension to other experimental mapping populations is straightforward.

Postzygotic isolation involves strong mitochondrial and sex-specific effects in Tigriopus californicus, a species lacking heteromorphic sex chromosomes

Detailed studies of the genetics of speciation have focused on a few model systems, particularly Drosophila. The copepod Tigriopus californicus offers an alternative that differs from standard animal models in that it lacks heteromorphic chromosomes (instead, sex determination is polygenic) and has reduced opportunities for sexual conflict, because females mate only once. Quantitative trait loci (QTL) mapping was conducted on reciprocal F2 hybrids between two strongly differentiated populations, using a saturated linkage map spanning all 12 autosomes and the mitochondrion. By comparing sexes, a possible sex ratio distorter was found but no sex chromosomes. Although studies of standard models often find an excess of hybrid male sterility factors, we found no QTL for sterility and multiple QTL for hybrid viability (indicated by non-Mendelian adult ratios) and other characters. Viability problems were found to be stronger in males, but the usual explanations for weaker hybrid males (sex chromosomes, sensitivity of spermatogenesis, sexual selection) cannot fully account for these male viability problems. Instead, higher metabolic rates may amplify deleterious effects in males. Although many studies of standard speciation models find the strongest genetic incompatibilities to be nuclear-nuclear (specifically X chromosome-autosome), we found the strongest deleterious interaction in this system was mito-nuclear. Consistent with the snowball theory of incompatibility accumulation, we found that trigenic interactions in this highly divergent cross were substantially more frequent (>6×) than digenic interactions. This alternative system thus allows important comparisons to studies of the genetics of reproductive isolation in more standard model systems.

A simulation study of permutation, bootstrap, and gene dropping for assessing statistical significance in the case of unequal relatedness

We used simulations to evaluate methods for assessing statistical significance in association studies. When the statistical model appropriately accounted for relatedness among individuals, unrestricted permutation tests and a few other simulation-based methods effectively controlled type I error rates; otherwise, only gene dropping controlled type I error but at the expense of statistical power.

Verification of QTL for grain starch content and its genetic correlation with oil content using two connected RIL populations in high-oil maize

Grain oil content is negatively correlated with starch content in maize in general. In this study, 282 and 263 recombinant inbred lines (RIL) developed from two crosses between one high-oil maize inbred and two normal dent maize inbreds were evaluated for grain starch content and its correlation with oil content under four environments. Single-trait QTL for starch content in single-population and joint-population analysis, and multiple-trait QTL for both starch and oil content were detected, and compared with the result obtained in the two related F(2∶3) populations. Totally, 20 single-population QTL for grain starch content were detected. No QTL was simultaneously detected across all ten cases. QTL at bins 5.03 and 9.03 were all detected in both populations and in 4 and 5 cases, respectively. Only 2 of the 16 joint-population QTL had significant effects in both populations. Three single-population QTL and 8 joint-population QTL at bins 1.03, 1.04-1.05, 3.05, 8.04-8.05, 9.03, and 9.05 could be considered as fine-mapped. Common QTL across F(2∶3) and RIL generations were observed at bins 5.04, 8.04 and 8.05 in population 1 (Pop.1), and at bin 5.03 in population 2 (Pop.2). QTL at bins 3.02-3.03, 3.05, 8.04-8.05 and 9.03 should be focused in high-starch maize breeding. In multiple-trait QTL analysis, 17 starch-oil QTL were detected, 10 in Pop.1 and 7 in Pop.2. And 22 single-trait QTL failed to show significance in multiple-trait analysis, 13 QTL for starch content and 9 QTL for oil content. However, QTL at bins 1.03, 6.03-6.04 and 8.03-8.04 might increase grain starch content and/or grain oil content without reduction in another trait. Further research should be conducted to validate the effect of these QTL in the simultaneous improvement of grain starch and oil content in maize.

Verification of QTL for grain starch content and its genetic correlation with oil content using two connected RIL populations in high-oil maize

Grain oil content is negatively correlated with starch content in maize in general. In this study, 282 and 263 recombinant inbred lines (RIL) developed from two crosses between one high-oil maize inbred and two normal dent maize inbreds were evaluated for grain starch content and its correlation with oil content under four environments. Single-trait QTL for starch content in single-population and joint-population analysis, and multiple-trait QTL for both starch and oil content were detected, and compared with the result obtained in the two related F(2∶3) populations. Totally, 20 single-population QTL for grain starch content were detected. No QTL was simultaneously detected across all ten cases. QTL at bins 5.03 and 9.03 were all detected in both populations and in 4 and 5 cases, respectively. Only 2 of the 16 joint-population QTL had significant effects in both populations. Three single-population QTL and 8 joint-population QTL at bins 1.03, 1.04-1.05, 3.05, 8.04-8.05, 9.03, and 9.05 could be considered as fine-mapped. Common QTL across F(2∶3) and RIL generations were observed at bins 5.04, 8.04 and 8.05 in population 1 (Pop.1), and at bin 5.03 in population 2 (Pop.2). QTL at bins 3.02-3.03, 3.05, 8.04-8.05 and 9.03 should be focused in high-starch maize breeding. In multiple-trait QTL analysis, 17 starch-oil QTL were detected, 10 in Pop.1 and 7 in Pop.2. And 22 single-trait QTL failed to show significance in multiple-trait analysis, 13 QTL for starch content and 9 QTL for oil content. However, QTL at bins 1.03, 6.03-6.04 and 8.03-8.04 might increase grain starch content and/or grain oil content without reduction in another trait. Further research should be conducted to validate the effect of these QTL in the simultaneous improvement of grain starch and oil content in maize.

Quantitative trait loci analysis of osteochondrosis traits in the elbow joint of pigs

Osteochondrosis is a growth disorder in the cartilage of young animals and is characterised by lesions found in the cartilage and bone. This study identified quantitative trait loci (QTLs) associated with six osteochondrosis lesion traits in the elbow joint of finishing pigs. The traits were: thickening of the cartilage, lesion in the subchondral bone, irregular cartilage surface, fissure under the cartilage, an irregular sagittal central groove and depression of the proximal edge of the radius. The study comprised 7172 finishing pigs from crossing 12 Duroc boars with 600 crossbred Landrace × Large White sows and included 462 single nucleotide polymorphism markers. The results showed 18 QTLs exceeding the 5% genome-wide threshold. The QTLs associated with lesions in the medial part of the condylus humeri (assumed to be the four main osteochondrosis traits) were, in most cases, at common locations, whereas the QTLs associated with depression of the proximal edge of the radius in general were on the same chromosomes but at separate locations. The detected QTLs explain a large part of the genetic variation, which is promising for incorporating osteochondrosis into a breeding programme using marker-assisted selection.

Environmental stress increases selection against and dominance of deleterious mutations in inbred families of the Pacific oyster Crassostrea gigas

The deleterious effects of inbreeding are well documented and of major concern in conservation biology. Stressful environments have generally been shown to increase inbreeding depression; however, little is known about the underlying genetic mechanisms of the inbreeding-by-stress interaction and to what extent the fitness of individual deleterious mutations is altered under stress. Using microsatellite marker segregation data and quantitative trait locus (QTL) mapping methods, I performed a genome scan for deleterious mutations affecting viability (viability or vQTL) in two inbred families of the Pacific oyster Crassostrea gigas, reared in a stressful, nutrient-poor diet and a favourable, nutrient-rich diet, which had significant effects on growth and survival. Twice as many vQTL were detected in the stressful diet compared with the favourable diet, resulting primarily from substantially greater mortality of homozygous genotypes. At vQTL, estimates of selection (s) and dominance (h) were greater in the stressful environment (= 0.86 vs. 0.54 and = 0.35 vs. 0.18, in stressful and nonstressful diets, respectively). There was no evidence of interaction between vQTL. Individual vQTL differed across diets in selection only, or in both selection and dominance, and some vQTL were not affected by diet. These results suggest that stress-associated increases in selection against individual deleterious alleles underlie greater inbreeding depression with stress. Furthermore, the finding that inbreeding-by-environment interaction appears, to some extent, to be locus specific, helps to explain previous observations of lineage-specific expression of inbreeding depression and environment-specific purging, which have important implications for conservation and evolutionary biology.

A statistical procedure to map high-order epistasis for complex traits

Genetic interactions or epistasis have been thought to play a pivotal role in shaping the formation, development and evolution of life. Previous work focused on lower-order interactions between a pair of genes, but it is obviously inadequate to explain a complex network of genetic interactions and pathways. We review and assess a statistical model for characterizing high-order epistasis among more than two genes or quantitative trait loci (QTLs) that control a complex trait. The model includes a series of start-of-the-art standard procedures for estimating and testing the nature and magnitude of QTL interactions. Results from simulation studies and real data analysis warrant the statistical properties of the model and its usefulness in practice. High-order epistatic mapping will provide a routine procedure for charting a detailed picture of the genetic regulation mechanisms underlying the phenotypic variation of complex traits.

A two-step approach to map quantitative trait loci for meat quality in connected porcine F(2) crosses considering main and epistatic effects

The aim of this study was to map QTL for meat quality traits in three connected porcine F(2) crosses comprising around 1000 individuals. The three crosses were derived from the founder breeds Chinese Meishan, European Wild Boar and Pietrain. The animals were genotyped genomewide for approximately 250 genetic markers, mostly microsatellites. They were phenotyped for seven meat quality traits (pH at 45 min and 24 h after slaughter, conductivity at 45 min and 24 h after slaughter, meat colour, drip loss and rigour). QTL mapping was conducted using a two-step procedure. In the first step, the QTL were mapped using a multi-QTL multi-allele model that was tailored to analyse multiple connected F(2) crosses. It considered additive, dominance and imprinting effects. The major gene RYR1:g.1843C>T affecting the meat quality on SSC6 was included as a cofactor in the model. The mapped QTL were tested for pairwise epistatic effects in the second step. All possible epistatic effects between additive, dominant and imprinting effects were considered, leading to nine orthogonal forms of epistasis. Numerous QTL were found. The most interesting chromosome was SSC6. Not all genetic variance of meat quality was explained by RYR1:g.1843C>T. A small confidence interval was obtained, which facilitated the identification of candidate genes underlying the QTL. Epistasis was significant for the pairwise QTL on SSC12 and SSC14 for pH24 and for the QTL on SSC2 and SSC5 for rigour. Some evidence for additional pairwise epistatic effects was found, although not significant. Imprinting was involved in epistasis.

Quantitative trait locus analysis of stage-specific inbreeding depression in the Pacific oyster Crassostrea gigas

Inbreeding depression and genetic load have been widely observed, but their genetic basis and effects on fitness during the life cycle remain poorly understood, especially for marine animals with high fecundity and high, early mortality (type-III survivorship). A high load of recessive mutations was previously inferred for the Pacific oyster Crassostrea gigas, from massive distortions of zygotic, marker segregation ratios in F(2) families. However, the number, genomic location, and stage-specific onset of mutations affecting viability have not been thoroughly investigated. Here, we again report massive distortions of microsatellite-marker segregation ratios in two F(2) hybrid families, but we now locate the causative deleterious mutations, using a quantitative trait locus (QTL) interval-mapping model, and we characterize their mode of gene action. We find 14-15 viability QTL (vQTL) in the two families. Genotypic frequencies at vQTL generally suggest selection against recessive or partially recessive alleles, supporting the dominance theory of inbreeding depression. No epistasis was detected among vQTL, so unlinked vQTL presumably have independent effects on survival. For the first time, we track segregation ratios of vQTL-linked markers through the life cycle, to determine their stage-specific expression. Almost all vQTL are absent in the earliest life stages examined, confirming zygotic viability selection; vQTL are predominantly expressed before the juvenile stage (90%), mostly at metamorphosis (50%). We estimate that, altogether, selection on vQTL caused 96% mortality in these families, accounting for nearly all of the actual mortality. Thus, genetic load causes substantial mortality in inbred Pacific oysters, particularly during metamorphosis, a critical developmental transition warranting further investigation.

Mapping quantitative trait loci for metabolic and cytological fatness traits of connected F2 crosses in pigs

In the present study 3 connected F(2) crosses were used to map QTL for classical fat traits as well as fat-related metabolic and cytological traits in pigs. The founder breeds were Chinese Meishan, European Wild Boar, and Pietrain with to some extent the same founder animals in the different crosses. The different selection history of the breeds for fatness traits as well as the connectedness of the crosses led to a high statistical power. The total number of F(2) animals varied between 694 and 966, depending on the trait. The animals were genotyped for around 250 genetic markers, mostly microsatellites. The statistical model was a multi-allele, multi-QTL model that accounted for imprinting. The model was previously introduced from plant breeding experiments. The traits investigated were backfat depth and fat area as well as relative number of fat cells with different sizes and 2 metabolic traits (i.e., soluble protein content as an indicator for the level of metabolic turnover and NADP-malate dehydrogenase as an indicator for enzyme activity). The results revealed in total 37 significant QTL on chromosomes 1, 2, 4, 5, 6, 7, 8, 9, 14, 17, and 18, with often an overlap of confidence intervals of several traits. These confidence intervals were in some cases remarkably small, which is due to the high statistical power of the design. In total, 18 QTL showed significant imprinting effects. The small and overlapping confidence intervals for the classical fatness traits as well as for the cytological and metabolic traits enabled positional and functional candidate gene identification for several mapped QTL.

Large scale genome-wide association and LDLA mapping study identifies QTLs for boar taint and related sex steroids

Background

Boar taint is observed in a high proportion of uncastrated male pigs and is characterized by an unpleasant odor/flavor in cooked meat, primarily caused by elevated levels of androstenone and skatole. Androstenone is a steroid produced in the testis in parallel with biosynthesis of other sex steroids like testosterone and estrogens. This represents a challenge when performing selection against androstenone in breeding programs, without simultaneously decreasing levels of other steroids. The aim of this study was to use high-density genome wide association (GWA) in combination with linkage disequilibrium-linkage analysis (LDLA) to identify quantitative trait loci (QTL) associated with boar taint compounds and related sex steroids in commercial Landrace (n = 1,251) and Duroc (n = 918) breeds.

Results

Altogether, 14 genome wide significant (GWS) QTL regions for androstenone in subcutaneous fat were obtained from the LDLA study in Landrace and 14 GWS QTL regions in Duroc. LDLA analysis revealed that 7 of these QTL regions, located on SSC 1, 2, 3, 7 and 15, were obtained in both breeds. All 14 GWS androstenone QTLs in Landrace are also affecting the estrogens at chromosome wise significance (CWS) or GWS levels, while in Duroc, 3 of the 14 QTLs affect androstenone without affecting any of the estrogens. For skatole, 10 and 4 QTLs were GWS in the LDLA analysis for Landrace and Duroc respectively, with 4 of these detected in both breeds. The GWS QTLs for skatole obtained by LDLA are located at SSC 1, 5, 6, 7, 10, 11, 13 and 14.

Conclusion

This is the first report applying the Porcine 60 K SNP array for simultaneous analysis of boar taint compounds and related sex hormones, using both GWA and LDLA approaches. Several QTLs are involved in regulation of androstenone and skatole, and most of the QTLs for androstenone are also affecting the levels of estrogens. Seven QTLs for androstenone were detected in one breed and confirmed in the other, i.e. in an independent sample, although the majority of QTLs are breed specific. Most QTLs for skatole do not negatively affect other sex hormones and should be easier to implement into the breeding scheme.

Revealing genetic relationships between compounds affecting boar taint and reproduction in pigs

Boar taint is characterized by an unpleasant taste or odor in intact male pigs and is primarily attributed to increased concentrations of androstenone and skatole and to a lesser extent by increased indole. The boar taint compounds skatole and indole are produced by gut bacteria, metabolized in the liver, and stored in the fat tissue. Androstenone, on the other hand, is synthesized in the testis along with testosterone and estrogens, which are known to be important factors affecting fertility. The main goal of this study was to investigate the relationship between genetic factors involved in the primary boar taint compounds in an attempt to discover ways to reduce boar taint without decreasing fertility-related compounds. Heritabilities and genetic correlations between traits were estimated for compounds related to boar taint (androstenone, skatole, indole) and reproduction (testosterone, 17β-estradiol, and estrone sulfate). Heritabilities in the range of 0.47 to 0.67 were detected for androstenone concentrations in both fat and plasma, whereas those for skatole and indole were slightly less (0.27 to 0.41). The genetic correlations between androstenone in plasma and fat were extremely high (0.91 to 0.98) in Duroc and Landrace. In addition, genetic correlations between androstenone (both plasma and fat) and the other sex steroids (estrone sulfate, 17β-estradiol, and testosterone) were very high, in the range of 0.80 to 0.95. Furthermore, a genome-wide association study (GWA) and a combined linkage disequilibrium and linkage analysis (LDLA) were conducted on 1,533 purebred Landrace and 1,027 purebred Duroc to find genome regions involved in genetic control of the boar taint compounds androstenone, skatole, and indole, and sex hormones related to fertility traits. Up to 3,297 informative SNP markers were included for both breeds, including SNP from several boar taint candidate genes. From the GWA study, we found that altogether 27 regions were significant at a genome-wide level (P < 0.05) and an additional 7 regions were significant at a chromosomal level. From the LDLA study, 7 regions were significant on a genome-wide level and an additional 7 regions were significant at a chromosomal level. The most convincing associations were obtained in 6 regions affecting skatole and indole in fat on chromosomes 1, 2, 3, 7, 13, and 14, 1 region on chromosome 6 affecting androstenone in plasma only, and 5 regions on chromosomes 3, 4, 13, and 15 affecting androstenone, testosterone, and estrogens.

Joint QTL analysis of three connected F2-crosses in pigs

Background

Numerous QTL mapping resource populations are available in livestock species. Usually they are analysed separately, although the same founder breeds are often used. The aim of the present study was to show the strength of analysing F₂-crosses jointly in pig breeding when the founder breeds of several F₂-crosses are the same.

Methods

Three porcine F₂-crosses were generated from three founder breeds (i.e. Meishan, Pietrain and wild boar). The crosses were analysed jointly, using a flexible genetic model that estimated an additive QTL effect for each founder breed allele and a dominant QTL effect for each combination of alleles derived from different founder breeds. The following traits were analysed: daily gain, back fat and carcass weight. Substantial phenotypic variation was observed within and between crosses. Multiple QTL, multiple QTL alleles and imprinting effects were considered. The results were compared to those obtained when each cross was analysed separately.

Results

For daily gain, back fat and carcass weight, 13, 15 and 16 QTL were found, respectively. For back fat, daily gain and carcass weight, respectively three, four, and five loci showed significant imprinting effects. The number of QTL mapped was much higher than when each design was analysed individually. Additionally, the test statistic plot along the chromosomes was much sharper leading to smaller QTL confidence intervals. In many cases, three QTL alleles were observed.

Conclusions

The present study showed the strength of analysing three connected F₂-crosses jointly. In this experiment, statistical power was high because of the reduced number of estimated parameters and the large number of individuals. The applied model was flexible and was computationally fast.

A general method for controlling the genome-wide type I error rate in linkage and association mapping experiments in plants

Control of the genome-wide type I error rate (GWER) is an important issue in association mapping and linkage mapping experiments. For the latter, different approaches, such as permutation procedures or Bonferroni correction, were proposed. The permutation test, however, cannot account for population structure present in most association mapping populations. This can lead to false positive associations. The Bonferroni correction is applicable, but usually on the conservative side, because correlation of tests cannot be exploited. Therefore, a new approach is proposed, which controls the genome-wide error rate, while accounting for population structure. This approach is based on a simulation procedure that is equally applicable in a linkage and an association-mapping context. Using the parameter settings of three real data sets, it is shown that the procedure provides control of the GWER and the generalized genome-wide type I error rate (GWER(k)).

Generalized linear model for interval mapping of quantitative trait loci

We developed a generalized linear model of QTL mapping for discrete traits in line crossing experiments. Parameter estimation was achieved using two different algorithms, a mixture model-based EM (expectation-maximization) algorithm and a GEE (generalized estimating equation) algorithm under a heterogeneous residual variance model. The methods were developed using ordinal data, binary data, binomial data and Poisson data as examples. Applications of the methods to simulated as well as real data are presented. The two different algorithms were compared in the data analyses. In most situations, the two algorithms were indistinguishable, but when large QTL are located in large marker intervals, the mixture model-based EM algorithm can fail to converge to the correct solutions. Both algorithms were coded in C++ and interfaced with SAS as a user-defined SAS procedure called PROC QTL.

Quantitative trait loci analysis of swine meat quality traits

A QTL study was performed in large half-sib families to characterize the genetic background of variation in pork quality traits as well as to examine the possibilities of including QTL in a marker-assisted selection scheme. The quality traits included ultimate pH in LM and the semimembranosus, drip loss, and the Minolta color measurements L*, a*, and b* representing meat lightness, redness, and yellowness, respectively. The families consist of 3,883 progenies of 12 Duroc boars that were evaluated to identify the QTL. The linkage map consists of 462 SNP markers on 18 porcine autosomes. Quantitative trait loci were mapped using a linear mixed model with fixed factors (sire, sex, herd, month, sow age) and random factors (polygenic effect, QTL effects, and litter). Chromosome-wide and genome-wide significance thresholds were determined by Peipho's approach, and 95% Bayes credibility intervals were estimated from a posterior distribution of the QTL position. In total, 31 QTL for the 6 meat quality traits were found to be significant at the 5% chromosome-wide level, among which 11 QTL were significant at the 5% genome-wide level and 5 of these were significant at the 0.1% genome-wide level. Segregation of the identified QTL in different families was also investigated. Most of the identified QTL segregated in 1 or 2 families. For the QTL affecting ultimate pH in LM and semimembranosus and L* and b* value on SSC6, the positions of the QTL and the shapes of the likelihood curves were almost the same. In addition, a strong correlation of the estimated effects of these QTL was found between the 4 traits, indicating that the same genes control these traits. A similar pattern was seen on SSC15 for the QTL affecting ultimate pH in the 2 muscles and drip loss. The results from this study will be helpful for fine mapping and identifying genes affecting meat quality traits, and tightly linked markers may be incorporated into marker-assisted selection programs.

Mapping quantitative trait Loci using distorted markers

Quantitative trait locus (QTL) mapping is usually performed using markers that follow a Mendelian segregation ratio. We developed a new method of QTL mapping that can use markers with segregation distortion (non-Mendelian markers). An EM (expectation-maximization) algorithm is used to estimate QTL and SDL (segregation distortion loci) parameters. The joint analysis of QTL and SDL is particularly useful for selective genotyping. Application of the joint analysis is demonstrated using a real life data from a wheat QTL mapping experiment.

QTL mapping in intercross and backcross populations

In the past two decades, various statistical approaches have been developed to identify quantitative trait locus with experimental organisms. In this chapter, we introduce several commonly used QTL mapping methods for intercross and backcross populations. Important issues related to QTL mapping, such as threshold and confidence interval calculations are also discussed. We list and describe five public domain QTL software packages commonly used by biologists.

QTL for body weight, morphometric traits and stress response in European sea bass Dicentrarchus labrax

Natural mating and mass spawning in the European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) complicate genetic studies and the implementation of selective breeding schemes. We utilized a two-step experimental design for detecting QTL in mass-spawning species: 2122 offspring from natural mating between 57 parents (22 males, 34 females and one missing) phenotyped for body weight, eight morphometric traits and cortisol levels, had been previously assigned to parents based on genotypes of 31 DNA microsatellite markers. Five large full-sib families (five sires and two dams) were selected from the offspring (570 animals), which were genotyped with 67 additional markers. A new genetic map was compiled, specific to our population, but based on the previously published map. QTL mapping was performed with two methods: half-sib regression analysis (paternal and maternal) and variance component analysis accounting for all family relationships. Two significant QTL were found for body weight on linkage group 4 and 6, six significant QTL for morphometric traits on linkage groups 1B, 4, 6, 7, 15 and 23 and three suggestive QTL for stress response on linkage groups 3, 14 and 23. The QTL explained between 8% and 38% of phenotypic variance. The results are the first step towards identifying genes involved in economically important traits like body weight and stress response in European sea bass.

Fine mapping quantitative trait loci under selective phenotyping strategies based on linkage and linkage disequilibrium criteria

In fine mapping of a large-scale experimental population where collection of phenotypes are very expensive, difficult to record or time-demanding, selective phenotyping could be used to phenotype the most informative individuals. Linkage analyses based sampling criteria (LAC) and linkage disequilibrium-based sampling criteria (LDC) for selecting individuals to phenotype are compared to random phenotyping in a quantitative trait loci (QTL) verification experiment using stochastic simulation. Several strategies based on LAC and LDC for selecting the most informative 30%, 40% or 50% of individuals for phenotyping to extract maximum power and precision in a QTL fine mapping experiment were developed and assessed. Linkage analyses for the mapping was performed for individuals sampled on LAC within families and combined linkage disequilibrium and linkage analyses was performed for individuals sampled across the whole population based on LDC. The results showed that selecting individuals with similar haplotypes to the paternal haplotypes (minimum recombination criterion) using LAC compared to random phenotyping gave at least the same power to detect a QTL but decreased the accuracy of the QTL position. However, in order to estimate unbiased QTL parameters based on LAC in a large half-sib family, prior information on QTL position was required. The LDC improved the accuracy to estimate the QTL position but not significantly compared to random phenotyping with the same sample size. When applying LDC (all phenotyping levels), the estimated QTL effect were closer to the true value in comparison to LAC. The results showed that the LDC were better than the LAC to select individuals for phenotyping and contributed to detection of the QTL.

A study on the mapping of quantitative trait loci in advanced populations derived from two inbred lines

In genetic and biological studies, the F2 population is one of the most popular and commonly used experimental populations mainly because it can be readily produced and its genome structure possesses several niceties that allow for productive investigation. These niceties include the equivalence between the proportion of recombinants and recombination rates, the capability of providing a complete set of three genotypes for every locus and an analytically attractive first-order Markovian property. Recently, there has been growing interest in using the progeny populations from F2 (advanced populations) because their genomes can be managed to meet specific purposes or can be used to enhance investigative studies. These advanced populations include recombinant inbred populations, advanced intercrossed populations, intermated recombinant inbred populations and immortalized F2 populations. Due to an increased number of meiosis cycles, the genomes of these advanced populations no longer possess the Markovian property and are relatively more complicated and different from the F2 genomes. Although issues related to quantitative trait locus (QTL) mapping using advanced populations have been well documented, still these advanced populations are often investigated in a manner similar to the way F2 populations are studied using a first-order Markovian assumption. Therefore, more efforts are needed to address the complexities of these advanced populations in more details. In this article, we attempt to tackle these issues by first modifying current methods developed under this Markovian assumption to propose an ad hoc method (the Markovian method) and explore its possible problems. We then consider the specific genome structures present in the advanced populations without invoking this assumption to propose a more adequate method (the non-Markovian method) for QTL mapping. Further, some QTL mapping properties related to the confounding problems that result from ignoring epistasis and to mapping closely linked QTL are derived and investigated across the different populations. Simulations show that the non-Markovian method outperforms the Markovian method, especially in the advanced populations subject to selfing. The results presented here may give some clues to the use of advanced populations for more powerful and precise QTL mapping.