Combined detection and introgression of QTL in outbred populations

Background

Detecting a QTL is only the first step in genetic improvement programs. When a QTL with desirable characteristics is found, e.g. in a wild or unimproved population, it may be interesting to introgress the detected QTL into the commercial population. One approach to shorten the time needed for introgression is to combine both QTL identification and introgression, into a single step. This combines the strengths of fine mapping and backcrossing and paves the way for introgression of desirable but unknown QTL into recipient animal and plant lines.

Methods

The method consisting in combining QTL mapping and gene introgression has been extended from inbred to outbred populations in which QTL allele frequencies vary both in recipient and donor lines in different scenarios and for which polygenic effects are included in order to model background genes. The effectiveness of the combined QTL detection and introgression procedure was evaluated by simulation through four backcross generations.

Results

The allele substitution effect is underestimated when the favourable QTL allele is not fixed in the donor line. This underestimation is proportional to the frequency differences of the favourable QTL allele between the lines. In most scenarios, the estimates of the QTL location are unbiased and accurate. The retained donor chromosome segment and linkage drag are similar to expected values from other published studies.

Conclusions

In general, our results show that it is possible to combine QTL detection and introgression even in outbred species. Separating QTL mapping and introgression processes is often thought to be longer and more costly. However, using a combined process saves at least one generation. With respect to the linkage drag and obligatory drag, the results of the combined detection and introgression scheme are very similar to those of traditional introgression schemes.     

The effects of selective genotyping on estimates of proportion of recombination between linked quantitative trait loci

Selective genotyping is the marker assay of only the more extreme phenotypes for a quantitative trait and is intended to increase the efficiency of quantitative trait loci (QTL) mapping. We show that selective genotyping can bias estimates of the recombination frequency between linked QTLs — upwardly when QTLs are in repulsion phase, and downwardly when QTLs are in coupling phase. We examined these biases under simple models involving two QTLs segregating in a backcross or F₂ population, using both analytical models and computer simulations. We found that bias is a function of the proportion selected, the magnitude of QTL effects, distance between QTLs and the dominance of QTLs. Selective genotyping thus may decrease the power of mapping multiple linked QTLs and bias the construction of a marker map. We suggest a large proportion than previously suggested (50%) or the entire population be genotyped if linked QTLs of large effects (explain > 10% phenotypic variance) are evident. New models need to be developed to explicitly incorporate selection into QTL map construction.     

Complementary genetic and genomic approaches help characterize the linkage group I seed protein QTL in soybean

Background

Cotton fibers (produced by Gossypium species) are the premier natural fibers for textile production. The two tetraploid species, G. barbadense (Gb) and G. hirsutum (Gh), differ significantly in their fiber properties, the former having much longer, finer and stronger fibers that are highly prized. A better understanding of the genetics and underlying biological causes of these differences will aid further improvement of cotton quality through breeding and biotechnology. We evaluated an inter-specific Gh × Gb recombinant inbred line (RIL) population for fiber characteristics in 11 independent experiments under field and glasshouse conditions. Sites were located on 4 continents and 5 countries and some locations were analyzed over multiple years.

Results

The RIL population displayed a large variability for all major fiber traits. QTL analyses were performed on a per-site basis by composite interval mapping. Among the 651 putative QTLs (LOD > 2), 167 had a LOD exceeding permutation based thresholds. Coincidence in QTL location across data sets was assessed for the fiber trait categories strength, elongation, length, length uniformity, fineness/maturity, and color. A meta-analysis of more than a thousand putative QTLs was conducted with MetaQTL software to integrate QTL data from the RIL and 3 backcross populations (from the same parents) and to compare them with the literature. Although the global level of congruence across experiments and populations was generally moderate, the QTL clustering was possible for 30 trait x chromosome combinations (5 traits in 19 different chromosomes) where an effective co-localization of unidirectional (similar sign of additivity) QTLs from at least 5 different data sets was observed. Most consistent meta-clusters were identified for fiber color on chromosomes c6, c8 and c25, fineness on c15, and fiber length on c3.

Conclusions

Meta-analysis provided a reliable means of integrating phenotypic and genetic mapping data across multiple populations and environments for complex fiber traits. The consistent chromosomal regions contributing to fiber quality traits constitute good candidates for the further dissection of the genetic and genomic factors underlying important fiber characteristics, and for marker-assisted selection.     

Mapping of QTL conferring resistance to sharp eyespot (Rhizoctonia cerealis) in bread wheat at the adult plant growth stage

Key message

Seven sharp eyespot resistance QTL were detected consistently across five environments and delimited to seven DNA marker intervals, respectively, six of which were independent of plant height and heading time.

Abstract

Sharp eyespot, caused mainly by the soil-borne fungus Rhizoctonia cerealis, is one of the important diseases of bread wheat (Triticum aestivum L.). This disease has escalated into a major threat to wheat production in some regions of the world. Wheat resistance to sharp eyespot can be a potential means to reduce the needs for application of fungicides and agricultural inputs. In the present study, the winter wheat lines, Luke and AQ24788-83, both of which possess quantitative resistance to sharp eyespot, were crossed and a population consisting 241 recombinant-inbred lines (RILs) was constructed. These RILs were assessed for sharp eyespot resistance by conducting five field and greenhouse trials during the period from 2008 to 2012, and they were genotyped with 549 simple-sequence repeat DNA markers. Seven quantitative trait loci (QTL) were detected consistently across the five trial environments to be associated with the sharp eyespot resistance. They were mapped on chromosomes 1A, 2B, 3B, 4A, 5D, 6B, and 7B. Four of these QTL are unequivocally novel, while it is possible that the other three might also be novel. Plant height and heading date of the 241 RILs were recorded in the four field trials. All of the seven disease resistance QTL were independent of plant height and heading time except one that was significantly associated with plant heading time. This association might be attributed genetically to a single QTL, or to different but closely linked QTL. In the case of single QTL, pleiotropism might be involved or the sharp eyespot resistance might be conferred in a physical instead of physiological nature.     

Testing for homogeneity of gametic disequilibrium across strata

Background

Variance component (VC) models are commonly used for Quantitative Trait Loci (QTL) mapping in outbred populations. Here, the QTL effect is given as a random effect and a critical part of the model is the relationship between the phenotypic values and the random effect. In the traditional VC model, each individual has a unique QTL effect and the relationship between these random effects is given as a covariance structure (known as the identity-by-descent (IBD) matrix).

Results

We present an alternative notation of the variance component model, where the elements of the random effect are independent base generation allele effects and sampling term effects. The relationship between the phenotypic vales and the random effect is given by an incidence matrix, which results in a novel, but statistically equivalent, version of the traditional VC model. A general algorithm to estimate this incidence matrix is presented. Since the model is given in terms of base generation allele effects and sampling term effects, these effects can be estimated separately using best linear unbiased prediction (BLUP). From simulated data, we showed that biallelic QTL effects could be accurately clustered using the BLUP obtained from our model notation when markers are fully informative, and that the accuracy increased with the size of the QTL effect. We also developed a measure indicating whether a base generation marker homozygote is a QTL heterozygote or not, by comparing the variances of the sampling term BLUP and the base generation allele BLUP. A ratio greater than one gives strong support for a QTL heterozygote.

Conclusion

We developed a simple presentation of the VC QTL model for identification of base generation allele effects in QTL linkage analysis. The base generation allele effects and sampling term effects were separated in our model notation. This clarifies the assumptions of the model and should also enhance the development of genome scan methods.     

Whole-genome QTL analysis for MAGIC

Key message

An efficient whole genome method of QTL analysis is presented for Multi-parent advanced generation integrated crosses.

Abstract

Multi-parent advanced generation inter-cross (MAGIC) populations have been developed for mice and several plant species and are useful for the genetic dissection of complex traits. The analysis of quantitative trait loci (QTL) in these populations presents some additional challenges compared with traditional mapping approaches. In particular, pedigree and marker information need to be integrated and founder genetic data needs to be incorporated into the analysis. Here, we present a method for QTL analysis that utilizes the probability of inheriting founder alleles across the whole genome simultaneously, either for intervals or markers. The probabilities can be found using three-point or Hidden Markov Model (HMM) methods. This whole-genome approach is evaluated in a simulation study and it is shown to be a powerful method of analysis. The HMM probabilities lead to low rates of false positives and low bias of estimated QTL effect sizes. An implementation of the approach is available as an R package. In addition, we illustrate the approach using a bread wheat MAGIC population.     

Threatened edible insects in Hidalgo, Mexico and some measures to preserve them

Background

The Caatinga (dry land vegetation) is one of the most characteristic vegetation types in northeastern Brazil. It occupies a large percentage of the semi-arid region there, and generally supports two major types of economic activity: seasonal agriculture and the harvesting of plant products. However, very little information is available concerning the interaction of people with the plants of the Caatinga.

Methods

A study was undertaken with the participation of 31 adults from a rural community in the state of Pernambuco, Brazil, in order to analyze the patterns of use of medicinal plant resources, and to test a number of hypotheses concerning their use and local knowledge about them. The sources of medicinal plants used by the local community, the differences in oral information concerning the use of plants with their effective uses, and the role of exotic plants in local folk medicine practices were examined.

Results

Forty-eight plant species were cited as having medicinal uses, of which 56.25% are native to the Caatinga region. The patterns of harvesting and the importance of these trees and shrubs as medicinal plants seem to be compatible with a hypothesis based on the seasonal availability of plant resources. There is no direct correlation between known medicinal plants and those used by the local population, which agrees with observations made in different tropical regions. However, this observation was not interpreted in terms of the idea of "erosion" of knowledge (commonly used to explain this lack of correlation), but rather to propose two new concepts: "mass knowledge" and "stock knowledge".

Conclusion

Native plants are a very significant component of locally used medicinal plants, although exotic plants are important for treating specific health problems – which leads the proposal of a hypothesis of diversification.     

Identification of agronomically important QTL in tetraploid potato cultivars using a marker–trait association analysis

Key message

Nineteen tuber quality traits in potato were phenotyped in 205 cultivars and 299 breeder clones. Association analysis using 3364 AFLP loci and 653 SSR-alleles identified QTL for these traits.

Abstract

Two association mapping panels were analysed for marker–trait associations to identify quantitative trait loci (QTL). The first panel comprised 205 historical and contemporary tetraploid potato cultivars that were phenotyped in field trials at two locations with two replicates (the academic panel). The second panel consisted of 299 potato cultivars and included recent breeds obtained from five Dutch potato breeding companies and reference cultivars (the industrial panel). Phenotypic data for the second panel were collected during subsequent clonal selection generations at the individual breeding companies. QTL were identified for 19 agro-morphological and quality traits. Two association mapping models were used: a baseline model without, and a more advanced model with correction for population structure and genetic relatedness. Correction for population structure and genetic relatedness was performed with a kinship matrix estimated from marker information. The detected QTL partly not only confirmed previous studies, e.g. for tuber shape and frying colour, but also new QTL were found like for after baking darkening and enzymatic browning. Pleiotropic effects could be discerned for several QTL.     

QTL mapping of adult-plant resistance to leaf rust in a RIL population derived from a cross of wheat cultivars Shanghai 3/Catbird and Naxos

Key message

Six QTL for adult plant resistance to leaf rust, including two QTL effective against additional diseases, were identified in a RIL population derived from a cross between Shanghai 3/Catbird and Naxos.

Abstract

Leaf rust is an important wheat disease and utilization of adult-plant resistance (APR) may be the best approach to achieve long-term protection from the disease. The CIMMYT spring wheat line Shanghai 3/Catbird (SHA3/CBRD) showed a high level of APR to Chinese Puccinia triticina pathotypes in the field. To identify APR genes in this line, a mapping population of 164 recombinant inbred lines (RILs) was developed from a cross of this line and Naxos, a moderately susceptible German cultivar. The RILs were evaluated for final disease severity (FDS) at Baoding, Hebei province, and Zhoukou, Henan province, in the 2010–2011 and 2011–2012 cropping seasons. QTL analysis detected one major QTL derived from SHA3/CBRD on chromosome 2BS explaining from 15 to 37 % of the phenotypic variance across environments. In addition one minor resistance QTL on chromosome 1AL from SHA3/CBRD and four minor QTL from Naxos on chromosomes 2DL, 5B, 7BS, and 7DS were also detected. SHA3/CBRD also possessed seedling resistance gene Lr26, and Naxos contained Lr1 based on gene postulation following tests with an array of P. triticina pathotypes and molecular marker assays. These seedling resistance and APR genes and their closely linked molecular markers are potentially useful for improving leaf rust resistance in wheat breeding programs.     

A multi-marker test based on family data in genome-wide association study

Background

Requirements for successful implementation of multivariate animal threshold models including phenotypic and genotypic information are not known yet. Here simulated horse data were used to investigate the properties of multivariate estimators of genetic parameters for categorical, continuous and molecular genetic data in the context of important radiological health traits using mixed linear-threshold animal models via Gibbs sampling. The simulated pedigree comprised 7 generations and 40000 animals per generation. Additive genetic values, residuals and fixed effects for one continuous trait and liabilities of four binary traits were simulated, resembling situations encountered in the Warmblood horse. Quantitative trait locus (QTL) effects and genetic marker information were simulated for one of the liabilities. Different scenarios with respect to recombination rate between genetic markers and QTL and polymorphism information content of genetic markers were studied. For each scenario ten replicates were sampled from the simulated population, and within each replicate six different datasets differing in number and distribution of animals with trait records and availability of genetic marker information were generated. (Co)Variance components were estimated using a Bayesian mixed linear-threshold animal model via Gibbs sampling. Residual variances were fixed to zero and a proper prior was used for the genetic covariance matrix.

Results

Effective sample sizes (ESS) and biases of genetic parameters differed significantly between datasets. Bias of heritability estimates was -6% to +6% for the continuous trait, -6% to +10% for the binary traits of moderate heritability, and -21% to +25% for the binary traits of low heritability. Additive genetic correlations were mostly underestimated between the continuous trait and binary traits of low heritability, under- or overestimated between the continuous trait and binary traits of moderate heritability, and overestimated between two binary traits. Use of trait information on two subsequent generations of animals increased ESS and reduced bias of parameter estimates more than mere increase of the number of informative animals from one generation. Consideration of genotype information as a fixed effect in the model resulted in overestimation of polygenic heritability of the QTL trait, but increased accuracy of estimated additive genetic correlations of the QTL trait.

Conclusion

Combined use of phenotype and genotype information on parents and offspring will help to identify agonistic and antagonistic genetic correlations between traits of interests, facilitating design of effective multiple trait selection schemes.     

A proteomics study of the response of North Ronaldsay sheep to copper challenge

Background

The accuracy of predicting disease occurrence using epidemic models relies on an understanding of the system or population under investigation. At the time of the Foot and Mouth disease (FMD) outbreak of 2001, there were limited reports in the literature as to the cattle population structure in Britain. In this paper we examine the temporal patterns of cattle births, deaths, imports and movements occurring within Britain, reported to the Department for the Environment, Food and Rural Affairs (DEFRA) through the British Cattle Movement service (BCMS) during the period 1^st January 2002 to 28^th February 2005.

Results

In Britain, the number of reported cattle births exhibit strong seasonality characterised by a large spring peak followed by a smaller autumn peak. Other event types also exhibit strong seasonal trends; both the reported number of cattle slaughtered and "on-farm" cattle deaths increase during the final part of the year. After allowing for seasonal components by smoothing the data, we illustrate that there is very little remaining non-seasonal trend in the number of cattle births, "on-farm" deaths, slaughterhouse deaths, on- and off-movements. However after allowing for seasonal fluctuations the number of cattle imports has been decreasing since 2002. Reporting of movements, births and deaths was more frequent on certain days of the week. For instance, greater numbers of cattle were slaughtered on Tuesdays, Wednesdays and Thursdays. Evidence for digit preference was found in the reporting of births and "on-farm" deaths with particular bias towards over reporting on the 1^st, 10^th and 20^th of each month.

Conclusion

This study provides insight into the population and movement dynamics of the British cattle population. Although the population is in constant flux, seasonal and long term trends can be identified in the number of reported births, deaths and movements of cattle. Incorporating this temporal variation in epidemic disease modelling may result in more accurate model predictions and may usefully inform future surveillance strategies.     

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Background

In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.

Methods

An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.

Results

Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.

Conclusion

This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.     

An Empirical Method for Establishing Positional Confidence Intervals Tailored for Composite Interval Mapping of QTL

Background

Improved genetic resolution and availability of sequenced genomes have made positional cloning of moderate-effect QTL realistic in several systems, emphasizing the need for precise and accurate derivation of positional confidence intervals (CIs) for QTL. Support interval (SI) methods based on the shape of the QTL likelihood curve have proven adequate for standard interval mapping, but have not been shown to be appropriate for use with composite interval mapping (CIM), which is one of the most commonly used QTL mapping methods.

Results

Based on a non-parametric confidence interval (NPCI) method designed for use with the Haley-Knott regression method for mapping QTL, a CIM-specific method (CIM-NPCI) was developed to appropriately account for the selection of background markers during analysis of bootstrap-resampled data sets. Coverage probabilities and interval widths resulting from use of the NPCI, SI, and CIM-NPCI methods were compared in a series of simulations analyzed via CIM, wherein four genetic effects were simulated in chromosomal regions with distinct marker densities while heritability was fixed at 0.6 for a population of 200 isolines. CIM-NPCIs consistently capture the simulated QTL across these conditions while slightly narrower SIs and NPCIs fail at unacceptably high rates, especially in genomic regions where marker density is high, which is increasingly common for real studies. The effects of a known CIM bias toward locating QTL peaks at markers were also investigated for each marker density case. Evaluation of sub-simulations that varied according to the positions of simulated effects relative to the nearest markers showed that the CIM-NPCI method overcomes this bias, offering an explanation for the improved coverage probabilities when marker densities are high.

Conclusions

Extensive simulation studies herein demonstrate that the QTL confidence interval methods typically used to positionally evaluate CIM results can be dramatically improved by accounting for the procedural complexity of CIM via an empirical approach, CIM-NPCI. Confidence intervals are a critical measure of QTL utility, but have received inadequate treatment due to a perception that QTL mapping is not sufficiently precise for procedural improvements to matter. Technological advances will continue to challenge this assumption, creating even more need for the current improvement to be refined.     

Quantitative trait loci that control the oil content variation of rapeseed (Brassica napus L.)

Key message

This report describes an integrative analysis of seed-oil-content quantitative trait loci (QTL) in Brassica napus , using a high-density genetic map to align QTL among different populations.

Abstract

Rapeseed (Brassica napus) is an important source of edible oil and sustainable energy. Given the challenge involved in using only a few genes to substantially increase the oil content of rapeseed without affecting the fatty acid composition, exploitation of a greater number of genetic loci that regulate the oil content variation among rapeseed germplasm is of fundamental importance. In this study, we investigated variation in the seed-oil content among two related genetic populations of Brassica napus, the TN double-haploid population and its derivative reconstructed-F₂ population. Each population was grown in multiple experiments under different environmental conditions. Mapping of quantitative trait loci (QTL) identified 41 QTL in the TN populations. Furthermore, of the 20 pairs of epistatic interaction loci detected, approximately one-third were located within the QTL intervals. The use of common markers on different genetic maps and the TN genetic map as a reference enabled us to project QTL from an additional three genetic populations onto the TN genetic map. In summary, we used the TN genetic map of the B. napus genome to identify 46 distinct QTL regions that control seed-oil content on 16 of the 19 linkage groups of B. napus. Of these, 18 were each detected in multiple populations. The present results are of value for ongoing efforts to breed rapeseed with high oil content, and alignment of the QTL makes an important contribution to the development of an integrative system for genetic studies of rapeseed.     

Genetic diversity and association mapping of seed vigor in rice (Oryza sativa L.)

Key message

Twenty-seven QTLs were identified for rice seed vigor, in which 16 were novel QTLs. Fifteen elite parental combinations were designed for improving seed vigor in rice.

Abstract

Seed vigor is closely related to direct seeding in rice (Oryza sativa L.). Previous quantitative trait locus (QTL) studies for seed vigor were mainly derived from bi-parental segregating populations and no report from natural populations. In this study, association mapping for seed vigor was performed on a selected sample of 540 rice cultivars (419 from China and 121 from Vietnam). Population structure was estimated on the basis of 262 simple sequence repeat (SSR) markers. Seed vigor was evaluated by root length (RL), shoot length (SL) and shoot dry weight in 2011 and 2012. Abundant phenotypic and genetic diversities were found in the studied population. The population was divided into seven subpopulations, and the levels of linkage disequilibrium (LD) ranged from 10 to 80 cM. We identified 27 marker–trait associations involving 18 SSR markers for three traits. According to phenotypic effects for alleles of the detected QTLs, elite alleles were mined. These elite alleles could be used to design parental combinations and the expected results would be obtained by pyramiding or substituting the elite alleles per QTL (apart from possible epistatic effects). Our results demonstrate that association mapping can complement and enhance previous QTL information for marker-assisted selection and breeding by design.     

Variation at the vernalisation genes Vrn-H1 and Vrn-H2 determines growth and yield stability in barley (Hordeum vulgare) grown under dryland conditions in Syria

Key message

Spring growth in barley controlled by natural variation at Vrn-H1 and Vrn-H2 improved yield stability in marginal Syrian environments.

Abstract

The objective of the present study was to identify QTL influencing agronomic performance in rain-fed Mediterranean environments in a recombinant inbred line (RIL) population, ARKE derived from the Syrian barley landrace, Arta and the Australian feed cultivar, Keel. The population was field tested for agronomic performance at two locations in Syria for 4 years with two sowing dates, in autumn and winter. Genotypic variability in yield of the RIL population was mainly affected by year-to-year variation presumably caused by inter-annual differences in rainfall distribution. The spring growth habit and early flowering inherited from the Australian cultivar Keel increased plant height and biomass and improved yield stability in Syrian environments. QTL for yield and biomass coincided with the map location of flowering time genes, in particular the vernalisation genes Vrn-H1 and Vrn-H2. In marginal environments with terminal drought, the Vrn-H1 allele inherited from Keel improved final biomass and yield. Under changing climate conditions, such as shorter winters, reduced rainfall, and early summer drought, spring barley might thus outperform the traditional vernalisation-sensitive Syrian landraces. We present the ARKE population as a valuable genetic resource to further elucidate the genetics of drought adaptation of barley in the field.     

On the efficiency of composite interval mapping in an outbred population

Composite interval mapping (CIM) has been successfully applied to the detection of QTL in experimental animals and plants. However, practical analyses based on CIM have been reported mainly for populations derived from cross between inbred lines. There are few studies on QTL analyses with CIM in outbred populations. To evaluate the applicability of CIM to outbred populations is prerequisite for the fine mapping of QTL in industrial animals such as pig and chicken. Some markers are usually not fully informative in outbred populations. In application of CIM to outbred populations, the influence of inclusion of such uninformative markers used as covariates on the efficiency of CIM should be investigated. In this paper a least-squares method for CIM was formalized in an F(2) population derived by crossing two outbred lines. The efficiencies of CIM were evaluated for outbred populations in comparison with simple interval mapping (SIM) for several cases of marker informativeness using simulations. By incorporating markers linked to a tested position as well as those unlinked, CIM showed a higher efficiency to separate two linked QTL over SIM. The efficiency of dissection was enhanced as the marker informativeness was increased. The power of CIM to detect an isolated QTL was improved by excluding markers linked to a tested position from covariates and higher than SIM regardless of marker informativeness. In conclusion, CIM is a useful procedure for the analysis of QTL in outbred populations even under low marker informativeness.     

Identification, distribution and molecular evolution of the pacifastin gene family in Metazoa

Background

This study was motivated by the observation of unusual mitochondrial haplotype distributions and associated physiological differences between populations of the killifish Fundulus heteroclitus distributed along the Atlantic coast of North America. A distinct "northern" haplotype is fixed in all populations north of New Jersey, and does not appear south of New Jersey except in extreme upper-estuary fresh water habitats, and northern individuals are known to be more tolerant of hyposmotic conditions than southern individuals. Complete mitochondrial genomes were sequenced from individuals from northern coastal, southern coastal, and fresh water populations (and from out-groups). Comparative genomics approaches were used to test multiple evolutionary hypotheses proposed to explain among-population genome variation including directional selection and hybridization.

Results

Structure and organization of the Fundulus mitochondrial genome is typical of animals, yet subtle differences in substitution patterns exist among populations. No signals of directional selection or hybridization were detected. Mitochondrial genes evolve at variable rates, but all genes exhibit very low dN/dS ratios across all lineages, and the southern population harbors more synonymous polymorphism than other populations.

Conclusion

Evolution of mitochondrial genomes within Fundulus is primarily governed by interaction between strong purifying selection and demographic influences, including larger historical population size in the south. Though directional selection and hybridization hypotheses were not supported, adaptive processes may indirectly contribute to partitioning of variation between populations.     

Persistence of accuracy of genome-wide breeding values over generations when including a polygenic effect

Background

When estimating marker effects in genomic selection, estimates of marker effects may simply act as a proxy for pedigree, i.e. their effect may partially be attributed to their association with superior parents and not be linked to any causative QTL. Hence, these markers mainly explain polygenic effects rather than QTL effects. However, if a polygenic effect is included in a Bayesian model, it is expected that the estimated effect of these markers will be more persistent over generations without having to re-estimate the marker effects every generation and will result in increased accuracy and reduced bias.

Methods

Genomic selection using the Bayesian method, ''BayesB'' was evaluated for different marker densities when a polygenic effect is included (GWpEBV) and not included (GWEBV) in the model. Linkage disequilibrium and a mutation drift balance were obtained by simulating a population with a Ne of 100 over 1,000 generations.

Results

Accuracy of selection was slightly higher for the model including a polygenic effect than for the model not including a polygenic effect whatever the marker density. The accuracy decreased in later generations, and this reduction was stronger for lower marker densities. However, no significant difference in accuracy was observed between the two models. The linear regression of TBV on GWEBV and GWpEBV was used as a measure of bias. The regression coefficient was more stable over generations when a polygenic effect was included in the model, and was always between 0.98 and 1.00 for the highest marker density. The regression coefficient decreased more quickly with decreasing marker density.

Conclusions

Including a polygenic effect had no impact on the selection accuracy, but showed reduced bias, which is especially important when estimates of genome-wide markers are used to estimate breeding values over more than one generation.