Sensitivity of methods for estimating breeding values using genetic markers to the number of QTL and distribution of QTL variance

The objective of this simulation study was to compare the effect of the number of QTL and distribution of QTL variance on the accuracy of breeding values estimated with genomewide markers (MEBV). Three distinct methods were used to calculate MEBV: a Bayesian Method (BM), Least Angle Regression (LARS) and Partial Least Square Regression (PLSR). The accuracy of MEBV calculated with BM and LARS decreased when the number of simulated QTL increased. The accuracy decreased more when QTL had different variance values than when all QTL had an equal variance. The accuracy of MEBV calculated with PLSR was affected neither by the number of QTL nor by the distribution of QTL variance. Additional simulations and analyses showed that these conclusions were not affected by the number of individuals in the training population, by the number of markers and by the heritability of the trait. Results of this study show that the effect of the number of QTL and distribution of QTL variance on the accuracy of MEBV depends on the method that is used to calculate MEBV.     

Application of non-parametric bootstrap methods to estimate confidence intervals for QTL location in a beef cattle QTL experimental population

Empirical confidence intervals (CIs) for the estimated quantitative trait locus (QTL) location from selective and non-selective non-parametric bootstrap resampling methods were compared for a genome scan involving an Angus x Brahman reciprocal fullsib backcross population. Genetic maps, based on 357 microsatellite markers, were constructed for 29 chromosomes using CRI-MAP V2.4. Twelve growth, carcass composition and beef quality traits (n = 527-602) were analysed to detect QTLs utilizing (composite) interval mapping approaches. CIs were investigated for 28 likelihood ratio test statistic (LRT) profiles for the one QTL per chromosome model. The CIs from the non-selective bootstrap method were largest (87 7 cM average or 79-2% coverage of test chromosomes). The Selective II procedure produced the smallest CI size (42.3 cM average). However, CI sizes from the Selective II procedure were more variable than those produced by the two LOD drop method. CI ranges from the Selective II procedure were also asymmetrical (relative to the most likely QTL position) due to the bias caused by the tendency for the estimated QTL position to be at a marker position in the bootstrap samples and due to monotonicity and asymmetry of the LRT curve in the original sample.     

Transformation of QTL genotypic effects to allelic effects

The genotypic and allelic effect models are equivalent in terms of QTL detection in a simple additive model, but the QTL allelic model has the advantage of providing direct information for marker-assisted selection. However, the allelic matrix is four times as large as the genotypic IBD matrix, causing computational problems, especially in genome scans examining multiple positions. Transformation from genotypic to allelic effects, after estimating the genotypic effects with a smaller IBD matrix, can solve this problem. Although the validity of transformation from genotypic to allelic effects has been disputed, this work proves that transformation can successfully yield unique allelic effects when genotypic and allelic IBD matrixes exist.     

Estimation,variance and optimal sampling of gene diversity

An extension of Nei's analysis of diversity in a subdivided population is proposed for a haploid locus. The differentiation G _STbecomes a natural extension of Wright's F _STand generalizes Weir and Cockerham's parameter of co-ancestry by relaxing the assumption of identical correlation for all the alleles. Inter- and intrapopulation variances of the estimated diversities and differentiation are derived. Finally, the optimal sampling strategy for measuring G _STwhen a fixed number of individuals can be analysed is considered. It is shown that, at a given locus, there is a unique sample size per population which yields the smallest variance of G _ST,regardless of the number of populations studied. These theoretical developments are illustrated with an analysis of chloroplast DNA diversity in a forest tree. The results emphasize the necessity of sampling many populations, rather than many individuals per population, for an accurate measurement of the subdivision of gene diversity at a single locus.     

Phosphate number and acyl chain length determine the subcellular location and lateral mobility of phosphoinositides

Phosphoinositides are critical regulators of ion channel and transporter activity. There are multiple isomers of biologically active phosphoinositides in the plasma membrane and the different lipid species are non-randomly distributed. However, the mechanism by which cells impose selectivity and directionality on lipid movements and so generate a non-random lipid distribution remains unclear. In the present study we investigated which structural elements of phosphoinositides are responsible for their subcellular location and movement. We incubated phosphatidylinositol (PI), phosphatidylinositol 4-monophosphate (PI(4)P) and phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) with short or long acyl chains in CHO and HEK cells. We show that phosphate number and acyl chain length determine cellular location and translocation movement. In CHO cells, PI(4,5)P2 with a long acyl chain was released into the cytosol easily because of a low partition coefficient whereas long chain PI was released more slowly because of a high partition coefficient. In HEK cells, the cellular location and translocation movement of PI were similar to those of PI in CHO cells, whereas those of PI(4,5)P2 were different; some mechanism restricted the translocation movement of PI(4,5)P2, and this is in good agreement with the extremely low lateral diffusion of PI(4,5)P2. In contrast to the dependence on the number of phosphates of the phospholipid head group of long acyl chain analogs, short acyl chain phospholipids easily undergo translocation movement regardless of cell type and number of phosphates in the lipid headgroup.     

Evolution of variance in offspring number: The effects of population size and migration

It was shown by Gillespie [1974. Am. Nat. 108, 145–151], that if two genotypes produce the same average number of offspring on but have a different variance associated within each generation, the genotype with a lower variance will have a higher effective fitness. Specifically, the effective fitness is {ei65-1}, where w is the mean fitness, {ei65-2} is the variance in offspring number, and N is the total population size. The model also predicts that if a strategy has a higher arithmetic mean fitness and a higher variance than the competitor, the outcome of selection will depend on the population size (with larger population sizes favoring the highvariance, high-mean genotype). This suggests that for metapopulation sizes favoring the high-variance, high-mean genotype). This suggests that for metapopulations with large numbers of (relatively) small demes, a strategy with lower variance and lower mean may be favored if the migration rate is low while higher migration rates (consistent with a larger effective population size) favor the opposite strategy. Individual-based simulation confirms that this is indeed the case for an island model of migration, though the effect of migration differs greatly depending on whether migration precedes or follows selection. It is noted in the appendix that while Gillespie [1974. Am. Nat. 108, 145–151] does seem to be heuristically accurate, it is not clear that the definition of effective fitness follows from his derivation.     

A strategy to determine HLA class II restriction broadly covering the DR, DP, and DQ allelic variants most commonly expressed in the general population

Classic ways to determine MHC restriction involve inhibition with locus-specific antibodies and antigen presentation assays with panels of cell lines matched or mismatched at the various loci of interest. However, these determinations are often complicated by T cell epitope degeneracy and promiscuity. We describe a selection of 46 HLA DR, DQ, and DP specificities that provide worldwide population (phenotypic) coverage of almost 90 % at each locus, and account for over 66 % of all genes at each locus. This panel afforded coverage of at least four HLA class II alleles in over 95 % of the individuals in four study populations of diverse ethnicity from the USA and South Africa. Next, a panel of single HLA class II-transfected cell lines, corresponding to these 46 allelic variants was assembled, consisting of lines previously developed and 15 novel lines generated for the present study. The novel lines were validated by assessing their HLA class II expression by FACS analysis, the in vitro peptide binding activity of HLA molecules purified from the cell lines, and their antigen presenting capacity to T cell lines of known restriction. We also show that these HLA class II-transfected cell lines can be used to rapidly and unambiguously determine HLA restriction of epitopes recognized by an individual donor in a single experiment. This panel of lines will enable high throughput determination of HLA restriction, enabling better characterization of HLA class II-restricted T cell responses and facilitating the development of HLA tetrameric staining reagents.     

Genome‐wide association QTL mapping for teat number in a purebred population of Duroc pigs

Because of increasing litter size in Western pig breeds, additional teats are desirable to increase the capacity for nursing offspring. We applied genome‐wide SNP markers to detect QTL regions that affect teat number in a Duroc population. We phenotyped 1024 animals for total teat number. A total of 36 588 SNPs on autosomes were used in the analysis. The estimated heritability for teat number was 0.34 ± 0.05 on the basis of a genomic relationship matrix constructed from all SNP markers. Using a BayesC method, we identified a total of 18 QTL regions that affected teat number in Duroc pigs; 9 of the 18 regions were newly detected.     

Using haplotypes for the prediction of allelic identity to fine-map QTL: characterization and properties

Background

Numerous methods have been developed over the last decade to predict allelic identity at unobserved loci between pairs of chromosome segments along the genome. These loci are often unobserved positions tested for the presence of quantitative trait loci (QTL). The main objective of this study was to understand from a theoretical standpoint the relation between linkage disequilibrium (LD) and allelic identity prediction when using haplotypes for fine mapping of QTL. In addition, six allelic identity predictors (AIP) were also compared in this study to determine which one performed best in theory and application.

Results

A criterion based on a simple measure of matrix distance was used to study the relation between LD and allelic identity prediction when using haplotypes. The consistency of this criterion with the accuracy of QTL localization, another criterion commonly used to compare AIP, was evaluated on a set of real chromosomes. For this set of chromosomes, the criterion was consistent with the mapping accuracy of a simulated QTL with either low or high effect. As measured by the matrix distance, the best AIP for QTL mapping were those that best captured LD between a tested position and a QTL. Moreover the matrix distance between a tested position and a QTL was shown to decrease for some AIP when LD increased. However, the matrix distance for AIP with continuous predictions in the [0,1] interval was algebraically proven to decrease less rapidly up to a lower bound with increasing LD in the simplest situations, than the discrete predictor based on identity by state between haplotypes (IBS _hap), for which there was no lower bound. The expected LD between haplotypes at a tested position and alleles at a QTL is a quantity that increases naturally when the tested position gets closer to the QTL. This behavior was demonstrated with pig and unrelated human chromosomes.

Conclusions

When the density of markers is high, and therefore LD between adjacent loci can be assumed to be high, the discrete predictor IBS _hap is recommended since it predicts allele identity correctly when taking LD into account.     

QTL for root angle and number in a population developed from bread wheats (Triticum aestivum) with contrasting adaptation to water-limited environments

Root architecture traits in wheat are important in deep soil moisture acquisition and may be used to improve adaptation to water-limited environments. The genetic architecture of two root traits, seminal root angle and seminal root number, were investigated using a doubled haploid population derived from SeriM82 and Hartog. Multiple novel quantitative trait loci (QTL) were identified, each one having a modest effect. For seminal root angle, four QTL (?log₁₀(P) >3) were identified on 2A, 3D, 6A and 6B, and two suggestive QTL (?log₁₀(P) >2) on 5D and 6B. For root number, two QTL were identified on 4A and 6A with four suggestive QTL on 1B, 3A, 3B and 4A. QTL for root angle and root number did not co-locate. Transgressive segregation was found for both traits. Known major height and phenology loci appear to have little effect on root angle and number. Presence or absence of the T1BL.1RS translocation did not significantly influence root angle. Broad sense heritability (h ²) was estimated as 50 % for root angle and 31 % for root number. Root angle QTL were found to be segregating between wheat cultivars adapted to the target production region indicating potential to select for root angle in breeding programs.     

Copy number alterations and allelic ratio in relation to recurrence of rectal cancer

Background

In rectal cancer, total mesorectal excision surgery combined with preoperative (chemo)radiotherapy reduces local recurrence rates but does not improve overall patient survival, a result that may be due to the harmful side effects and/or co-morbidity of preoperative treatment. New biomarkers are needed to facilitate identification of rectal cancer patients at high risk for local recurrent disease. This would allow for preoperative (chemo)radiotherapy to be restricted to high-risk patients, thereby reducing overtreatment and allowing personalized treatment protocols. We analyzed genome-wide DNA copy number (CN) and allelic alterations in 112 tumors from preoperatively untreated rectal cancer patients. Sixty-six patients with local and/or distant recurrent disease were compared to matched controls without recurrence. Results were validated in a second cohort of tumors from 95 matched rectal cancer patients. Additionally, we performed a meta-analysis that included 42 studies reporting on CN alterations in colorectal cancer and compared results to our own data.

Results

The genomic profiles in our study were comparable to other rectal cancer studies. Results of the meta-analysis supported the hypothesis that colon cancer and rectal cancer may be distinct disease entities. In our discovery patient study cohort, allelic retention of chromosome 7 was significantly associated with local recurrent disease. Data from the validation cohort were supportive, albeit not statistically significant, of this finding.

Conclusions

We showed that retention of heterozygosity on chromosome 7 may be associated with local recurrence in rectal cancer. Further research is warranted to elucidate the mechanisms and effect of retention of chromosome 7 on the development of local recurrent disease in rectal cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1550-0) contains supplementary material, which is available to authorized users.     

Conditional QTL underlying resistance to late blight in a diploid potato population

A large number of quantitative trait loci (QTL) for resistance to late blight of potato have been reported with a "conventional" method in which each phenotypic trait reflects the cumulative genetic effects for the duration of the disease process. However, as genes controlling response to disease may have unique contributions with specific temporal features, it is important to consider the phenotype as dynamic. Here, using the net genetic effects evidenced at consecutive time points during disease development, we report the first conditional mapping of QTL underlying late blight resistance in potato under five environments in Peru. Six conditional QTL were mapped, one each on chromosome 2, 7 and 12 and three on chromosome 9. These QTL represent distinct contributions to the phenotypic variation at different stages of disease development. By comparison, when conventional mapping was conducted, only one QTL was detected on chromosome 9. This QTL was the same as one of the conditional QTL. The results imply that conditional QTL reflect genes that function at particular stages during the host-pathogen interaction. The dynamics revealed by conditional QTL mapping could contribute to the understanding of the molecular mechanism of late blight resistance and these QTL could be used to target genes for marker development or manipulation to improve resistance.     

QTL affecting stress response to crowding in a rainbow trout broodstock population

Background

Genomic analyses have the potential to impact selective breeding programs by identifying markers that serve as proxies for traits which are expensive or difficult to measure. Also, identifying genes affecting traits of interest enhances our understanding of their underlying biochemical pathways. To this end we conducted genome scans of seven rainbow trout families from a single broodstock population to identify quantitative trait loci (QTL) having an effect on stress response to crowding as measured by plasma cortisol concentration. Our goal was to estimate the number of major genes having large effects on this trait in our broodstock population through the identification of QTL.

Results

A genome scan including 380 microsatellite markers representing 29 chromosomes resulted in the de novo construction of genetic maps which were in good agreement with the NCCCWA genetic map. Unique sets of QTL were detected for two traits which were defined after observing a low correlation between repeated measurements of plasma cortisol concentration in response to stress. A highly significant QTL was detected in three independent analyses on Omy16, many additional suggestive and significant QTL were also identified. With linkage-based methods of QTL analysis such as half-sib regression interval mapping and a variance component method, we determined that the significant and suggestive QTL explain about 40-43% and 13-27% of the phenotypic trait variation, respectively.

Conclusions

The cortisol response to crowding stress is a complex trait controlled in a sub-sample of our broodstock population by multiple QTL on at least 8 chromosomes. These QTL are largely different from others previously identified for a similar trait, documenting that population specific genetic variants independently affect cortisol response in ways that may result in different impacts on growth. Also, mapping QTL for multiple traits associated with stress response detected trait specific QTL which indicate the significance of the first plasma cortisol measurement in defining the trait. Fine mapping these QTL can lead towards the identification of genes affecting stress response and may influence approaches to selection for this economically important stress response trait.

     

The mean and variance of environmental temperature interact to determine physiological tolerance and fitness

Global climate change poses one of the greatest threats to biodiversity. Most analyses of the potential biological impacts have focused on changes in mean temperature, but changes in thermal variance will also impact organisms and populations. We assessed the combined effects of the mean and variance of temperature on thermal tolerances, organismal survival, and population growth in Drosophila melanogaster. Because the performance of ectotherms relates nonlinearly to temperature, we predicted that responses to thermal variation (±0° or ±5°C) would depend on the mean temperature (17° or 24°C). Consistent with our prediction, thermal variation enhanced the rate of population growth (r(max)) at a low mean temperature but depressed this rate at a high mean temperature. The interactive effect on fitness occurred despite the fact that flies improved their heat and cold tolerances through acclimation to thermal conditions. Flies exposed to a high mean and a high variance of temperature recovered from heat coma faster and survived heat exposure better than did flies that developed at other conditions. Relatively high survival following heat exposure was associated with low survival following cold exposure. Recovery from chill coma was affected primarily by the mean temperature; flies acclimated to a low mean temperature recovered much faster than did flies acclimated to a high mean temperature. To develop more realistic predictions about the biological impacts of climate change, one must consider the interactions between the mean environmental temperature and the variance of environmental temperature.     

mRNA location and translation rate determine protein targeting to dual destinations

1. Download : Download high-res image (143KB)
2. Download : Download full-size image