Methods for predicting superior genotypes under multiple environments based on QTL effects

Methods were developed for predicting two kinds of superior genotypes (superior line and superior hybrid) based on quantative trait locus (QTL) effects including epistatic and QTL × environment interaction effects. Formulae were derived for predicting the total genetic effect of any individual with known QTLs genotype derived from the mapping population in a specific environment. Two algorithms, enumeration algorithm and stepwise tuning algorithm, were used to select the best multi-locus combination of all the putative QTLs. Grain weight per plant (GW) in rice was analyzed as a working example to demonstrate the proposed methods. Results showed that the predicted superior lines and superior hybrids had great superiorities over the F₁ hybrid, indicating large breeding potential remained for further improvement on GW. Results also showed that epistatic effects and their interaction with environments largely contributed to the superiorities of the predicted superior lines and superior hybrids. User-friendly software, QTLNetwork, version 1.0, was developed based on the methods in the present paper.     

水、旱条件下水稻茎秆主要抗倒伏性状的QTL分析

以粳型旱稻IRAT109和粳型水稻越富杂交的116个DH株系的群体为材料,利用已构建的分子标记连锁图(包括94个RFLP标记和71个SSR标记),定位了水稻茎秆主要抗倒伏性状的QTL。在水田、旱田栽培条件下,考查了乳熟期DH系及其亲本的茎基粗、茎秆长及茎秆强度等性状。相关分析表明,茎基粗与茎秆长、茎基粗与茎秆强度及茎秆长与茎秆强度间均呈极显著正相关。利用QTLMAPPER进行水、旱田单环境定位分析及水、旱田联合定位分析定位了控制这些性状的QTL。水、旱田单环境定位分析结果表明：3个性状共检测到9个加性QTL和5对上位性QTL;联合定位分析表明：茎基粗、茎秆长共检测到6个加性QTL和6对上位性QTL,其中6个加性QTL和1对上位性QTL在两种方法下都检测到。旱田条件下检测到2个加性及2对上位性QTL(bctla、ct9、c16a-c16c和cs5-cs12)对表型变异的贡献率(简称贡献率)大于30％。这些高贡献率QTL可能对旱田条件下旱稻抗倒伏分子育种有重要意义。     

GSD: a genetic screen database

The systematic assignment of gene function to a sequenced genome is one of the outstanding challenges in the post-genomic era. Large-scale systematic mutagenesis screens are important tools for reaching this goal. Here we describe GSD, a software package that allows storage and integration of data from genetic screens. GSD was initially developed for a large-scale F3 mutagenesis screen for developmental mutants of medaka (Oryzias latipes). The version presented here supports a wide range of different screens (mutagenesis, RNAi, morpholinos, transgenesis and others) using different organisms. Data are stored in a relational database and can be made accessible through web interfaces. Researchers can enter data describing their screened embryos: They can track statistics, submit images and describe the resulting phenotypes using a phenotype classification ontology. We developed a fish phenotype classification ontology of medaka and zebrafish for this software package and made it available to the public. In addition, a list of genetic lines resulting from each screen can be generated. These lines (mutant alleles, transgenic lines) can be described and categorized in the same ways as the screened individuals. Raw data from the screen can be integrated to describe these lines. A query module that searches this list can be used to publish the screen results on the Internet. A test version is available at and the software can be downloaded from this site.     

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.     

Study of strategies for selecting quantitative trait locus mapping procedures by computer simulation

Along with the development and integration of molecular genetics and quantitative genetics, many quantitative trait locus (QTL) mapping studies have been conducted using different mapping populations in various crop species. Existing QTLs can be used for marker-assisted breeding and map-based cloning, whereas the false-positive QTLs are no use. The purpose of this study is to evaluate the suitability of different mapping procedures for data from different genetic models. In this study, four types of recombinant inbred lines (RILs) with different genetic models, viz. additive QTLs (Model I), additive and epistatic QTLs (Model II), additive QTLs and QTL × environment interaction (Model III), additive, epistatic QTLs and QTL × environment interaction (Model IV), were simulated by computer. Six types of QTL mapping procedures, viz. CIM, MIMF, MIMR, ICIM, MQM and NWIM, on four kinds of QTL mapping software, viz. WinQTL Cartographer Version 2.5, IciMapping Version 2.0, MapQTL Version 5.0 and QTLnetwork Version 2.0, were used for screening QTLs of the simulated RILs. The results showed that different mapping procedures have different suitability for different genetic models. CIM and MQM can only screen Model I data. MIMR, MIMF and ICIM can only screen Model I and Model II data. NWIM can screen all four models’ data. It can be concluded that different genetic models’ data have different most suitable mapping procedures. In practical experiments where the genetic model of the data is unknown, a multiple model mapping strategy should be used, that is a full model scanning with complex model procedure followed by verification with other procedures corresponding to the scanning results.     

Marker-based estimation of the coefficient of coancestry in hybrid breeding programmes

Molecular markers allow to estimate the pairwise relatedness between the members of a breeding pool when their selection history is no longer available or has become too complex for a classical pedigree analysis. The field of population genetics has several estimation procedures at its disposal, but when the genotyped individuals are highly selected inbred lines, their application is not warranted as the theoretical assumptions on which these estimators were built, usually linkage equilibrium between marker loci or even Hardy–Weinberg equilibrium, are not met. An alternative approach requires the availability of a genotyped reference set of inbred lines, which allows to correct the observed marker similarities for their inherent upward bias when used as a coancestry measure. However, this approach does not guarantee that the resulting coancestry matrix is at least positive semi-definite (psd), a necessary condition for its use as a covariance matrix. In this paper we present the weighted alikeness in state (WAIS) estimator. This marker-based coancestry estimator is compared to several other commonly applied relatedness estimators under realistic hybrid breeding conditions in a number of simulations. We also fit a linear mixed model to phenotypical data from a commercial maize breeding programme and compare the likelihood of the different variance structures. WAIS is shown to be psd which makes it suitable for modelling the covariance between genetic components in linear mixed models involved in breeding value estimation or association studies. Results indicate that it generally produces a low root mean squared error under different breeding circumstances and provides a fit to the data that is comparable to that of several other marker-based alternatives. Recommendations for each of the examined coancestry measures are provided.     

Industry benefits from recent genetic progress in sheep and beef populations

An analytical model that evaluates the benefits from 10 years of genetic improvement over a 20-year time frame was specified. Estimates of recent genetic trends in recorded traits, industry statistics and published estimates of the economic values of trait changes were used to parameterise the model for the UK sheep and beef industries. Despite rates of genetic change in the relevant performance-recorded breeding populations being substantially less than theoretical predictions, the financial benefits of genetic change were substantial. Over 20 years, the benefits from 10 years of genetic progress at recently achieved rates in recorded hill sheep, sheep crossing sire and sheep terminal sire breeding programmes was estimated to be ￡5.3, ￡1.0 and ￡11.5 million, respectively. If dissemination of genetic material is such that these rates of change are also realised across the entire ram breeding industry, the combined benefits would be ￡110.8 million. For beef cattle, genetic evaluation systems have been operating within all the major breeds for some years with quite widespread use of performance recording, and so genetic trends within the beef breeds were used as predictors of industry genetic change. Benefits from 10 years of genetic progress at recent rates of change, considering a 20-year time frame, in terminal sire beef breeds are expected to be ￡4.9 million. Benefits from genetic progress for growth and carcass characters in dual-purpose beef breeds were ￡18.2 million after subtraction of costs associated with a deterioration in calving traits. These benefits may be further offset by unfavourable associated changes in maternal traits. Additional benefits from identification and use of the best animals available from the breeding sector for commercial matings through performance recording and genetic evaluation could not be quantified. When benefits of genetic improvement were expressed on an annual present value basis and compared with lagged annual investment costs to achieve it, the internal rate of return (IRR) on the combined investment in sheep and beef cattle was 32%. Despite a much higher rate of participation in performance recording, the present value of benefits and the IRR were lower for beef cattle than for sheep. The implications of these results for future national and industry investment in genetic improvement infrastructure were discussed.     

Purification of myelin basic protein from bovine brain

Myelin basic protein (MBP) is a commonly used substrate for in vitro determination of numerous protein kinase activities. Herein we describe a rapid method for isolating relatively large amounts of MBP from bovine brain with a purity greater than that currently available from commercial sources. Lipids were first extracted from the CNS tissue by homogenization in sec-butanol. Washes under neutral and mildly basic conditions were employed to remove neutral and acidic proteins from the defatted residue. MBP was subsequently extracted under acidic conditions and further purified by chromatography on CM Sephadex C-25. Potential contaminating enzyme activities were destroyed by heart treatment. This method typically yields a recovery of 1.0-1.5 mg MBP per gram of starting material with a purity of greater than 95%. The MBP prepared in this manner was suitable for determination of kinase activities by both solution and the "in gel" kinase assay systems.     

Field evaluation of doubled haploid plants in the Apiaceae: dill (Anethum graveolens L.), caraway (Carum carvi L.), and fennel (Foeniculum vulgare Mill.)

The Apiaceae family includes vegetables, as well as herb and spice crops. Compared to major crops, there have been few breeding or genetic improvement programs for any of the Apiaceae, especially the herb and spice species. Haploidy technology can be used to develop uniform, true-breeding lines, as well as to accelerate breeding programs. Field trials of dill (Anethum graveolens L.), caraway (Carum carvi L.), and fennel (Foeniculum vulgare Mill.) doubled haploid (DH) lines were conducted over 2–5 cropping seasons. Several of the DH dill lines had desirable agronomic characteristics such as short uniform stature along with early maturity that would be useful for crop improvement. Seed yields and the essential oil content of the seed harvested from the best performing DH dill lines were either equal to or higher than the parental line. A DH annual caraway line was identified that produced higher seed yields than the industry standard. The main constituents of the essential oil for the DH lines of both dill and caraway were similar to the parental lines. Fennel DH lines exhibited differences in height but were too late in maturity for seed production under prairie conditions. The results indicate that not only were we able to generate DH lines that could be used in a crop improvement program, but we developed DH lines that could be used directly as cultivars as these lines performed better than the industry standard (parental line).     

Antigen processing of myelin basic protein is required prior to recognition by T cells inducing EAE.

The processing and presentation of whole myelin basic protein (MBP) and a 12 amino acid encephalitogenic peptide were investigated using MBP-immune and peptide-immune murine T cell lines. Myelin basic protein is the major component of central nervous system (CNS) white matter capable of inciting an autoimmune response which leads to the disease, experimental allergic encephalomyelitis (EAE), in a number of animal species. MBP-immune T cell lines caused a form of adoptively transferred EAE when injected into naive, syngeneic recipients. It has been found that both whole MBP and peptide required processing in order to induce proliferation of the T cell lines. The proliferative response was greatest when MBP was processed under conditions in which proteolysis was prevented. The demonstration that activation of encephalitogenic MBP immune T cells requires a processed form of MBP may have relevance to the human inflammatory CNS demyelinating condition, multiple sclerosis, for which EAE is the EAE is the prime animal model.     

Advanced optimum contribution selection as a tool to improve regional cattle breeds: a feasibility study for Vorderwald cattle

In the middle of the 20th century, increasing inbreeding rates were identified as a threat to livestock breeding. Consequences include reduced fertility, fitness and phenotypic expression of lethal alleles. An important step in mitigating this inbreeding was the introduction of optimum contribution selection (OCS). OCS facilitates the simultaneous management of genetic gain and inbreeding rates. However, using a standard OCS methodology for regional breeds with historical introgression for upgrading reasons could lead to reinforced selection on introgressed genetic material since those alleles improve the rate of genetic gain and reduce the average kinship in the population. Consequently, regional breeds may become genetically extinct if a standard OCS approach is used. Thus, the advanced OCS (aOCS) approach takes introgressed genetic material into account. The major goals of this study were to (i) gather key information on the feasibility of aOCS under practical conditions of the actual breeding scheme of Vorderwald cattle, (ii) identify superior strategies for implementing the actual scheme and (iii) examine whether historical breeding decisions to increase genetic gain by introgression from commercial breeds could have been avoided by using aOCS. Stochastic simulations were designed in this study to create populations from the historical gene pool by using aOCS. Simultaneously, all practical constraints of a breeding scheme were met. Thus, the simulated populations were comparable with real data. The annual genetic gain was higher in reality (1.56) than in the simulation scenarios (1.12–1.40). The introgressed genetic material increased to 61.3% in reality but was conserved at a final value of 15.3% (±0.78) across simulations. The classical rate of inbreeding and rate of native inbreeding were constrained to 0.092% on an annual basis. This value is equal to an effective population size of 100. The observed values for rates of inbreeding were 0.082–0.087% and 0.087–0.088% for classical and native kinship, respectively. The corresponding figures in reality were 0.067% and 0.184%, respectively. This study suggests that aOCS is feasible for Vorderwald cattle. Strategies for implementation are identified. Finally, we conclude that historical breeding decisions could have been avoided by using aOCS. The genetic gain would have been reduced by at least 12.2%, but the introgressed genetic material, genetic diversity and native genetic diversity would have been more desirable for a breed under conservation.     

Determination of genetic and pathological variance among Tilletia indica isolates and monosporidial lines using PCR based markers and host differentials

Karnal bunt of wheat caused by Tilletia indica is an important international quarantine disease in many countries. In this investigation, genetic and pathological variation among the 10 isolates and 15 monosporidial (Ms) lines belonged to different locations of North-West India was studied. Depending upon the pathogenic potential, most virulent and least aggressive isolate was found from Chaksu (Rajasthan) and Tarau (HP), which scored coefficients of infection 70.98 and 6.22, respectively, on susceptible host genotype HD 2009 under artificially inoculated conditions. Fifteen Ms lines were inoculated in 20 combinations. Most virulent compatible combination was found KB2MsD?×?KB6MsA, which scored co-efficient of infection 74.91%. Out of 32 Inter Simple Sequence Repeats based molecular markers, 28 were polymorphic generating 192 reproducible bands for all the T. indica isolates and Ms lines in this study. A grouping analysis using the unrooted neighbour – joining method was consistent with DARwin software and winboot analysis and combination approach suggested that self-paired Ms lines exhibit narrow genetic diversity. This result will be useful for developing integrated strategies for disease management and breeding programmes for improvement of the varieties.     

Efficient construction of high-density linkage map and its application to QTL analysis in barley

Using a High Efficiency Genome Scanning (HEGS) system and recombinant inbred (RI) lines derived from the cross of Russia 6 and H.E.S. 4, a high-density genetic map was constructed in barley. The resulting 1,595.7-cM map encompassed 1,172 loci distributed on the seven linkage groups comprising 1,134 AFLP, 34 SSR, three STS and vrs1 (kernel row type) loci. Including PCR reactions, gel electrophoresis and data processing, 6 months of work by a single person was sufficient for the whole mapping procedure under a reasonable cost. To make an appraisal of the resolution of genetic analysis for the 95 RI lines based on the constructed linkage map, we measured three agronomic traits: plant height, spike exsertion length and 1,000-kernel weight, and the analyzed quantitative trait loci (QTLs) associated with these traits. The results were compared on the number of detected QTLs and their effects between a high-density map and a skeleton map constructed by selected AFLP and anchor markers. The composite interval mapping on the high-density map detected more QTLs than the other analyses. Closely linked markers with QTLs on the high-density map could be powerful tools for marker-assisted selection in barley breeding programs and further genetic analyses including an advanced backcross analysis or a map-based cloning of QTL. Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by J.S. Heslop-Harrison     

Genetic Analysis of <Emphasis Type="Italic">Porphyra yezoensis</Emphasis> Using Microsatellite Markers

Microsatellites are repetitive genomic elements that show high levels of variation and therefore are useful tools for studying genetic polymorphism and constructing genetic linkage maps of eukaryotic organisms. Porphyra yezoensis Ueda is an economically important seaweed that is being targeted for genetic improvement using marker-assisted breeding. Hence, in an attempt to develop microsatellite markers for P. yezoensis, a microsatellite (or simple sequence repeat)-enriched library was constructed to identify (GA)n and (CA)n motifs. A total of 71 perfect microsatellite clones were identified, of which 30 simple sequence repeat primer pairs were developed. Of these, 24 (80%) amplified polymerase chain reaction products of expected sizes. Twelve primer pairs amplified two to four bands, whereas another 12 primer pairs produced monomorphic banding patterns. Data for 12 loci were analyzed using POPGENE software version 1.32. A total of 29 alleles were produced at 12 loci, with an average of 2.42 alleles (Na) and 1.81 effective alleles (Ne) per locus. These markers were used to analyze the genetic diversity within 11 geographically different lines of P. yezoensis. Overall, these lines were clustered into two divisions with those from close geographic locations clustering together. Further cloning and sequencing of size variant alleles at two microsatellite loci revealed that the variable numbers of motif repeats in different alleles were major sources of polymorphisms.     

A genetic linkage map of microsatellite,gene-specific and morphological markers in diploid <Emphasis Type="Italic">Fragaria</Emphasis>

Knowledge of the genetic diversity of a species is important for the choice of crossing parents in line and hybrid breeding. Our objective was to investigate European winter triticale using simple sequence repeat (SSR) markers and the coancestry coefficient (f) with regard to genetic diversity and grouping of germplasm. Three to five primer pairs for each of the 42 chromosomes were selected to analyse 128 European winter triticale varieties and breeding lines. SSR analysis resulted in the identification of 657 alleles with an average of 6.8 alleles per primer pair. The average polymorphism information content (PIC) for polymorphic markers was 0.54. Correlation between f and genetic similarity (GS) estimates based on Rogers Distance was low (r_f×GS(ABDR)=0.33). The analysis of molecular variance (AMOVA) revealed that 84.7% of the total variation was found within breeding companies, and 15.3% among them. In conclusion, SSR markers from wheat and rye provide a powerful tool for assessing genetic diversity in triticale. Even though no distinct groups within the European winter triticale pool could be detected by principal co-ordinate analysis, this study provides basic information about the genetic relationships for breeding purposes.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by H.C. Becker     

Response of flax genotypes to doubled haploid production

Forty-four flax genotypes with a diverse genetic background were evaluated for anther culture response using a standard anther culture protocol in order to determine the feasibility to initiate a routine haploid production system in applied breeding programs. A strong genotype effect on callus induction and shoot regeneration in anther culture was found in this study. A number of genotypes, including two low cadmium content lines 96-11785 and 96-11826, a high oil content line 96-22109 and a high linolenic acid content line M 4919 were identified as highly responsive. The impact of the findings in this study on flax breeding was discussed. This revised version was published online in June 2006 with corrections to the Cover Date.     

spagedi: a versatile computer program to analyse spatial genetic structure at the individual or population levels

spag e d i version 1.0 is a software primarily designed to characterize the spatial genetic structure of mapped individuals or populations using genotype data of codominant markers. It computes various statistics describing genetic relatedness or differentiation between individuals or populations by pairwise comparisons and tests their significance by appropriate numerical resampling. spag e d i is useful for: (i) detecting isolation by distance within or among populations and estimating gene dispersal parameters; (ii) assessing genetic relatedness between individuals and its actual variance, a parameter of interest for marker based inferences of quantitative inheritance; (iii) assessing genetic differentiation among populations, including the case of haploids or autopolyploids.     

Genome-based prediction of test cross performance in two subsequent breeding cycles

Genome-based prediction of genetic values is expected to overcome shortcomings that limit the application of QTL mapping and marker-assisted selection in plant breeding. Our goal was to study the genome-based prediction of test cross performance with genetic effects that were estimated using genotypes from the preceding breeding cycle. In particular, our objectives were to employ a ridge regression approach that approximates best linear unbiased prediction of genetic effects, compare cross validation with validation using genetic material of the subsequent breeding cycle, and investigate the prospects of genome-based prediction in sugar beet breeding. We focused on the traits sugar content and standard molasses loss (ML) and used a set of 310 sugar beet lines to estimate genetic effects at 384 SNP markers. In cross validation, correlations >0.8 between observed and predicted test cross performance were observed for both traits. However, in validation with 56 lines from the next breeding cycle, a correlation of 0.8 could only be observed for sugar content, for standard ML the correlation reduced to 0.4. We found that ridge regression based on preliminary estimates of the heritability provided a very good approximation of best linear unbiased prediction and was not accompanied with a loss in prediction accuracy. We conclude that prediction accuracy assessed with cross validation within one cycle of a breeding program can not be used as an indicator for the accuracy of predicting lines of the next cycle. Prediction of lines of the next cycle seems promising for traits with high heritabilities.