PBGworks T80
PBGworks T80
PBGworks T202,197
Author:
Matthew Robbins, The Ohio State University
The inbred backcross (IBC) population was proposed by Wehrhahn and Allard (1965) as a way of identifying genes or quantitative trait loci (QTL) that contribute to a quantitatively inherited trait. This is accomplished by developing a population that collectively contains most of the genome of a donor parent, divided among each individual line in the population. The majority of the genome of each line is from the recurrent parent, with a small portion from the donor parent. IBC breeding has also been employed for the introgression of exotic germplasm to improve quantitative traits in crop plants. This method has been utilized in bean (Bliss, 1981; Sullivan and Bliss, 1983), oilseed rape (Butruille et al., 1999), rice (Lin et al., 1998), cucumber (Robbins et al., 2008) and tomato (Hartman and St Clair, 1999; Doganlar et al., 2002; Kabelka et al., 2002; Kabelka et al., 2004; Yang et al., 2005, Robbins et al., 2009) for classical breeding and QTL studies.
The first stage of generating an IBC population (Fig. 1, steps 1–3) is similar to generating a backcross breeding population. One distinction is that many individuals are backcrossed to the recurrent parent to generate an IBC population. The second stage (Fig. 1, step 4) is similar to single-seed descent to generate recombinant inbred lines (RILs).
Figure 1. Schematic illustrating the development of an inbred backcross (IBC) population. Figure credit: Matthew Robbins, The Ohio State University.
An important consideration in creating an IBC population is the number of backcrossing generations. More backcrossing ensures that the IBC lines will be more like the recurrent parent, since the percentage of the genome from donor parent is reduced by half with each generation of backcrossing (see article on backcrossing). However, the probability of recovering the genes from the donor parent is reduced by half each generation due to the backcrossing process. The probability of recovering the gene(s) from the donor parent is (1/2)k+1 for a single gene and (1/2)2k+2 for two unlinked genes.
Development of this lesson was supported in part by the National Institute of Food and Agriculture (NIFA) Solanaceae Coordinated Agricultural Project, agreement 2009-85606-05673, administered by Michigan State University. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the United States Department of Agriculture.
PBGworks 647
Authors:
David M. Francis, The Ohio State University
Heather L. Merk, The Ohio State University
Matthew Robbins, The Ohio State University
This page is a continuation of the Overview of Analysis of Variance page and is intended to help plant breeders consider the notions of fixed and random effects and the impacts these can have on ANOVA in the context of plant breeding. Briefly, ANOVA is a statistical test that takes the total variation and assigns it to known causes, leaving a residual portion allocated to uncontrolled or unexplained variation, called the experimental error. By measuring variability as sums of squares deviating from the mean sum of squares for all observations, the variation assigned to different controlled causes will be additive. It is therefore important to completely define the statistical model. Otherwise, the experimental error may be unnecessarily inflated (McIntosh, 1983).
In the Overview of Analysis of Variance page, we considered the following linear model:
Y = m + f(treatment) + error
where
Intuitively, we may think about the treatments as being under our control and as “fixed.” Usually we are interested in comparing the dependent variable among factors/levels of the fixed effect. For example, we may want to evaluate whether differences in yield (dependent variable) between field locations for some elite cultivars we’ve been developing. To conduct this experiment, we would select the cultivars we want to evaluate and find suitable locations for our trial. We could think of the cultivars and locations as being fixed; we purposely chose to study different cultivars and locations. In this case, we are only interested in the performance of the elite cultivars we’re testing in the specific locations we’re testing.
Random effects, in contrast to fixed effects, are typically used to account for variance in the dependent variable. Also, unlike fixed effects, we aren’t looking to compare one level of the random effect to another. In our example, we could also consider location as a random effect. In the case of random effects, levels are chosen randomly from an infinite population and we want to make inferences that can extend beyond the sample. If this were the case, the cultivars would still be fixed effects, but location would be random. If we felt our locations were representative of all possible locations, we could use the different locations to help us make an evaluation of how well cultivars perform across locations as a whole, not just at the locations we’ve tested. The classification of effects as fixed or random determines the appropriate F-test.
McIntosh (1983) provides a set of reference tables for use during experimental design and analysis. These tables are intended for field experiments conducted over two or more locations or years. Some of the tables are replicated below.
Sources of variation | df | Mean squares | Expected mean squares1 | ||
---|---|---|---|---|---|
RL-RT | RL-FT | FL-FT | |||
Locations (l) | l-1 | M1 | σ2e + rσ2TL + tσ2R(L) + rtσ2L | σ2e + tσ2R(L) + rtσ2L | σ2e + tσ2R(L) + rtσ2L |
Blocks(Location) (r) | l(r-1) | M2 | σ2e + tσ2R(L) | σ2e + tσ2R(L) | σ2e + tσ2R(L) |
Treatment (t) | t-1 | M3 | σ2e + rσ2TL + rlσ2T | σ2e + rσ2TL + rlσ2T | σ2e + rlσ2T |
Location x treatment | (l-1)(t-1) | M4 | σ2e + rσ2TL | σ2e + rσ2TL | σ2e + rσ2TL |
Pooled error | l(r-1)(t-1) | M5 | σ2e | σ2e | σ2e |
1 R = random, F = fixed, L = location, T = treatment
Sources of variation | Mean squares | Expected mean squares1 | ||
---|---|---|---|---|
RL-RT | RL-FT | FL-FT | ||
Locations (l) | M1 | (M1+M5)/(M2+M4) | M1/M2 | M1/M2 |
Blocks(Location) (r) | M2 | |||
Treatment (t) | M3 | M3/M4 | M3/M4 | M3/M5 |
Location x treatment | M4 | M4/M5 | M4/M5 | M4/M5 |
Pooled error | M5 |
1 R = random, F = fixed, L = location, T = treatment
In a genetic/breeding experiment, treatments would likely be genotypes or varieties.
When designing experiments, plant breeders must consider the question they want to answer. Consequently, plant breeders must consider what type of statistical analyses are appropriate to answer the desired question. With regards to ANOVA, two important points should be considered in this context.
Many statistics textbooks provide a good discussion of theory and applications of ANOVA. Two examples are listed below.
Development of this page was supported in part by the National Institute of Food and Agriculture (NIFA) Solanaceae Coordinated Agricultural Project, agreement 2009-85606-05673, administered by Michigan State University. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the United States Department of Agriculture.
PBGworks 865
Author:
Audrey Sebolt, Michigan State University
Michigan State University’s tart cherry breeding and genetics program, led by Dr. Amy Iezzoni (Fig. 1), is the only one of its kind in the United States. After being hired in 1981, Dr. Iezzoni set out to determine industry needs for tart cherry and to access available germplasm to develop improved tart cherry varieties.
Figure 1. Dr. Amy Iezzoni, lead of Michigan State University’s tart cherry breeding and genetics program. Photo credit: Michigan State University Tart Cherry Breeding and Genetics Program.
Virtually all tart cherries are processed. Processed tart cherry products include jams, dried cherries, individually quick frozen cherries, cherry juice, and—making up the largest portion—pie filling. The tart cherry industry is a monoculture, consisting essentially of one cultivar: ‘Montmorency’, a 400-year-old cultivar from France (Fig. 2).
Figure 2. Montmorency tart cherry branch with fruit. Photo credit: Michigan State University Tart Cherry Breeding and Genetics Program.
First, most of the cherry germplasm and excellent varieties that would have provided alternatives to ‘Montmorency’ evolved or were bred in Eastern Europe. Prior to the cold war, they were essentially unavailable to the U.S. Second, ‘Montmorency’ is extremely productive. The trees flourish in the sandy soils and harsh winters of Western Michigan, which produces 75% of the nation’s tart cherries. ‘Montmorency’ requires very little horticultural management and can withstand trunk damage inflicted by mechanical harvesting. Fruit produced from this cultivar are generally uniform in size and have clear flesh and bright red skin, characteristics which have become the standard for ‘American cherry pie’. There are limitations to ‘Montmorency’; the fruit can be soft and the trees are highly susceptible to cherry leaf spot (Blumeriella jaapii) (Fig. 3), which is a major financial cost to the ~$39 million tart cherry industry.
Figure 3. Leaves susceptible (left) and resistant (right) to cherry leaf spot. Photo credit: Michigan State University Tart Cherry Breeding and Genetics Program.
When Dr. Iezzoni joined Michigan State, only a small collection of ‘Montmorency’ sports and varieties from Western Europe were available. Over a 15-year period, Dr. Iezzoni collected cherry accessions from Eastern Europe, the center of diversity for tart cherry, to expand her germplasm base. That effort led to the establishment of the world’s largest tart cherry germplasm collection, located at Michigan State University’s Clarksville Horticultural Research Station.
To incorporate germplasm from Eastern Europe, Dr. Iezzoni overcame genetically controlled self-incompatibility (Fig. 4). Today, Dr. Iezzoni’s tart cherry breeding program focuses on increased firmness, pit size and shape, late bloom time, disease resistance, processing savings due to less use of colorants and/or sugar, freestone or “airfree” (Fig. 5), and high yield.
Figure 4. Cherry pollination. Photo credit: Michigan State University Tart Cherry Breeding and Genetics Program.
Figure 5. Freestone or “air free” tart cherry. Photo credit: Cameron Peace, Washington State University.
For more information, visit Dr. Iezzoni’s website.
Development of this page was supported in part by the Michigan Cherry Committee and the USDA’s National Institute of Food and Agriculture (NIFA). Project title: RosBREED: Enabling marker-assisted breeding in Rosaceae is provided by the Specialty Crops Research Initiative Competitive Grant 2009-51181-05808. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s)and do not necessarily reflect the view of the United States Department of Agriculture.
PBGworks 937
Author:
Matthew Robbins, The Ohio State University
The polymerase chain reaction (PCR) is a procedure that mimics the cellular process of DNA replication using the machinery of heat-resistant bacteria in a cyclic manner, resulting in several million copies of a specific DNA sequence that can then be visualized through electrophoresis and staining with a dye. PCR is commonly used in plant genetics and molecular breeding to copy a specific DNA fragment from the genome of an individual as a step in the process of molecular marker assisted selection. The use of PCR to copy a specific portion of a genome is analogous to photocopying a specific page of a book. Table 1 illustrates this analogy by comparing the component required to copy DNA by PCR to those needed to photocopy a page of a book.
Photocopier items | PCR components |
---|---|
The book | The entire genome (called the DNA template) |
The page | A portion of the genome (fragment) we are interested in |
A bookmark | Primers that “mark” the specific fragment |
The copy machine |
The enzyme that copies DNA |
Paper and toner |
The four bases that make up DNA |
In the same way that a bookmark identifies the specific page to photocopy out of a book, PCR primers identify the specific fragment to be copied from the entire genome. In order to copy a page, the photocopier uses the paper and toner to make the copy. Similarly, the polymerase requires nucleotides to produce a replicate of the original DNA fragment.
To understand in more detail how these components function in PCR, the Plant and Soil Sciences eLibrary at the University of Nebraska-Lincoln has an informative lesson on PCR including an animation of the process:
Photo credit: Plant and Soil Science eLibrary
Another animation on PCR can be found at the Dolan DNA Learning Center, part of The Cold Spring Harbor Laboratory.
Photo credit: The Dolan DNA Learning Center
The Genetics Science Learning Center at the University of Utah also has an animation on PCR.
Photo credit: The Genetics Science Learning Center
When using PCR for genotyping, the amplified DNA fragments can be analyzed several different ways. DNA amplified by PCR can be:
For some PCR related entertainment, we recommend “The PCR Song“. With lyrics such as “PCR, when you need to find out who’s your Daddy; PCR, when you need to solve a crime…” this video produced by BioRad features characterizations of famous and not-so-famous folk singers. If you like the musical theme, the “GTCA Song” song rocks to the tune of YMCA while reviewing the biochemistry of PCR.
Development of this lesson was supported in part by the National Institute of Food and Agriculture (NIFA) Solanaceae Coordinated Agricultural Project, agreement 2009-85606-05673, administered by Michigan State University. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the United States Department of Agriculture.
Mention of specific companies is not intended for promotion purposes.
PBGworks 653
Authors:
Heather Merk, The Ohio State University; Deana Namuth-Covert, University of Nebraska-Lincoln; Matthew Robbins, The Ohio State University
At the end of this lesson you should:
The purpose of this article is to provide an example of how to genotype individual tomato plants with a molecular DNA marker. There are several different molecular marker systems available to assist plant breeding programs. For the purposes of this lesson, the marker chosen as an example is a cleaved amplified polymorphism (CAP) marker, a type of marker that is often visualized by gel electrophoresis. Briefly, a CAP marker exploits differences in DNA sequences between two polymerase chain reaction (PCR) products based on the presence or absence of restriction enzyme cutting sites found within that segment of DNA. To genotype a CAP marker, the segment of DNA is amplified using PCR then cut with a restriction enzyme (referred to as digestion, or restriction enzyme digestion), which only cuts at a specific DNA sequence. After digestion, the DNA is separated on agarose gel. CAP markers are designed so that the restriction enzyme will cut the DNA of one genotype, but not another.
Although different breeding program schemes can be used, in this particular case, the individual plants are from an F2 population that is segregating for the marker. In all breeding programs, the specific marker being used must be segregating among the plant population being used in order to be useful.
CAP markers are generally visualized using gel electrophoresis. When scoring any molecular DNA marker using gel electrophoresis, keep the following considerations in mind:
All these considerations will make it easier to score a marker from a gel photo. Next we will follow a specific CAP marker example in a tomato breeding program.
The gel photo below (Fig.1) is a CAP marker, CosOH57, genotyped in 30 individuals that were part of a larger F2 population developed from the parents OH88119 and 06.8068. The population was developed as part of a breeding project to incorporate bacterial spot resistance into elite germplasm. In order to score the gel, the bands are evaluated based on the considerations listed above:
Figure 1. Example gel photo of CAP marker CosOH57. The gel includes a DNA ladder, the parental genotypes (OH88119 and 6.8068), and 30 F2 individuals. Photo credit: Matthew Robbins, The Ohio State University.
Knowing the information outlined above, the gel can be scored. Most computer programs that use marker data in subsequent analyses have a specified data format. For segregating populations, many programs code the data in relation to the parents. For example, Joinmap and MapMaker, two programs that are commonly used for mapping, code genotypes from an F2 population as follows:
Code | Genotype |
---|---|
A | homozygous for parent 1 allele |
B | homozygous for parent 2 allele |
H | heterozygous |
C | not genotype A (dominant B allele, so could be a genotype like parent 2 or heterozygous) |
D | not genotype B (dominant A allele, so could be a genotype like parent 1 or heterozygous) |
“.” | genotype unknown (missing data) |
Keep in mind the following when scoring the genotypes:
Using the genotypic codes, each individual tomato plant is scored (Fig. 1). In the example we are following, CosOH57 is a codominant marker, so the 30 F2 individuals are coded as “A” when only the 216 bp band is present, “B” when a plant has both the 145 and 71 bp bands present, or “H” when all three bands are showing for an individual tomato plant.
Genotypic scores can also be coded by the molecular weight of the fragment. This is useful when genotyping a set of individuals without common parents, and especially if multiple alleles of the marker are present. In this simpler CosOH57 example, using the molecular weight scoring method, parent 1 would be scored as “216” and parent 2 could be scored as either “145” or “71.”
Once the molecular marker is scored, it is useful to organize the data in a spreadsheet or table format. This allows data from other markers genotyped in the same population to be combined in preparation for mapping or other analyses. The individual genotypes for CosOH57 have to be reorganized into a table with markers as rows and individual plant genotypes as columns (Table 2). It is important that “F2 Plant #1” is always the same plant, no matter the particular marker being genotyped. This is a common format for mapping software. The rows for Marker2 and Marker3 indicate that genotypic data can be added for additional markers. Although parental genotypes are not included in mapping analysis, it is useful to keep them with the data for reference.
Marker | OH88119 | 6.8068 | F2 Plant 1 | F2 Plant 2 | F2 Plant 3 | … |
---|---|---|---|---|---|---|
CosOH57 | A | B | A | A | H | |
Marker2 | ||||||
Marker3 | ||||||
… |
Data summaries are also useful to check whether the data collected seems reasonable based on what you expect for a particular population, or if something else may be going on, such as the marker being linked to a trait we are selecting for or forces such as natural selection are distorting the expected segregation pattern. In our example, we may want to verify that the CosOH57 marker genotypes segregate as expected—1:2:1—using a chi-square goodness-of-fit test (note: For a refresher on how to use chi-square, you may want to take a look at the chi-square lesson). The data for the gel photo above, not including the parents, is summarized in Table 3. The observed column is determined simply by counting the number of individual plants with each genotype. The expected number of each genotype is calculated by multiplying the expected frequency of the genotype by the total number of plants being genotyped:
Expected = Expected Frequency x Total
The expected frequency is determined based on the segregation ratio of 1:2:1 for our F2 population, which is 0.25: 0.5 :0.25. Thus, the expected frequency of the “A” genotype for CosOH57 is:
Expected “A” Genotype = Expected Frequency of “A” Genotype x Total Number of F2 Plants Being Genotyped
or
Expected “A” Genotype = 0.25 x 30 = 7.5
The expected frequencies and number of each genotype are also presented in Table 3.
Genotype | Observed | Expected frequency | Expected |
---|---|---|---|
A | 13 | 0.25 | 7.5 |
B | 7 | 0.25 | 7.5 |
H | 10 | 0.5 | 15 |
Total | 30 | 1 | 30 |
When the observed and expected numbers are used in a chi-squared goodness-of-fit test, the calculated p value is 0.057. Since this p value is a little greater than 0.05, a common level to declare significance, there is some evidence that CosOH57 may segregate as expected. Closer inspection of the data indicates that the actual observed frequency of genotype “A” may be higher than expected, while the H genotype may be lower than expected. Additional caution should be exercised because the relatively small number of F2 individuals make it difficult to interpret this chi-square test. Ideally, statisticians recommend genotyping an F2 population using at least 50 individuals.
In this tutorial we learned how to genotype a CAP marker that was scored in an F2 population. The principles we used apply to any other molecular marker that we may genotype, particularly molecular markers genotyped on a gel. These general principles also apply to other plant breeding schemes. We also learned how to organize data so that we can use it for genetic mapping. Finally, we learned how to perform a chi-square analysis as an additional test to help us determine the reliability of a specific marker in our breeding population.
For additional practice scoring an agarose gel:
Development of this page was supported in part by the National Institute of Food and Agriculture (NIFA) Solanaceae Coordinated Agricultural Project, agreement 2009-85606-05673, administered by Michigan State University. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the United States Department of Agriculture.
PBGworks 659
Author:
Matthew Robbins, The Ohio State University
Gel electrophoresis is commonly used in plant breeding and genomics for genotyping with molecular markers, but there are several other applications as well (see below). For example, specific DNA fragments used as markers and isolated from individual plants are amplified by the polymerase chain reaction (PCR) and the resulting DNA fragments are subsequently loaded on a gel. The gel is a solid, gelatin-like substance used to separate DNA fragments based on size. The gel is placed in a conductive salt buffer to which an electrical field is applied. As the negatively-charged DNA fragments migrate toward the positive pole, the gel acts as a size filter, with smaller fragments migrating faster than larger fragments.
In addition, this video illustrates the basics of DNA extraction and gel electrophoresis in tomato:
The Plant and Soil Sciences eLibrary at the University of Nebraska-Lincoln has an informative lesson on gel electrophoresis, including an animation of the process:
Photo credit: Plant and Soil Sciences eLibrary
Another animation on gel electrophoresis can be found at the Dolan DNA Learning Center, part of The Cold Spring Harbor Laboratory:
Photo credit: The Dolan DNA Learning Center
The Genetics Science Learning Center at the University of Utah also has an animation on gel electrophoresis:
Photo credit: The Genetics Science Learning Center
DNA can be separated by electrophoresis to:
Development of this lesson was supported in part by the National Institute of Food and Agriculture (NIFA) Solanaceae Coordinated Agricultural Project, agreement 2009-85606-05673, administered by Michigan State University. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the United States Department of Agriculture.
PBGworks 654
Authors:
David M. Francis, The Ohio State University; Heather L. Merk, The Ohio State University; Matthew Robbins, The Ohio State University
The analysis of variance (ANOVA) is a statistical tool that has two common applications in a plant breeding context. First, ANOVA can be used to test for differences between treatments in an experiment. Common examples of treatments are genotype, location, and variety. Second, ANOVA can be used to aid in estimates of heritability by partitioning variances. This module focuses on simple ANOVA models to evaluate differences between treatments.
Like other statistical tests, ANOVA assumes that certain assumptions are met. One of the principal assumptions of ANOVA is that the samples come from normally distributed populations, each with the same variance. In addition, it is assumed that the residuals come from a normally distributed population with equal variances (σ2). The Kruskal–Wallis test is an alternative to ANOVA when the above assumptions cannot be met.
ANOVA is a tool that can be used to test for differences among treatment means when the independent variable is categorical (e.g., genotypes could be AA, Aa, aa) and the dependent variable is continuous (e.g., yield measured in tons/acre). How does this work?
In ANOVA, the total variance of all samples is calculated. Portions of the total variance can be attributed to known causes (e.g., genotype). This leaves a residual portion of the variance that is uncontrolled or unexplained and is referred to as experimental error. Then the between-treatment variation (e.g., AA genotype variation vs. Aa genotype variation vs. aa genotype variation) is compared to the within-treatment variation (experimental error) (e.g., variation within the aa genotype) to assess whether differences in mean value between treatments are due to the treatment effects or chance.
In the simplest case, linear equations can be developed to describe the relationship between a trait and treatment. The question can then be asked, “which linear equation best fits the data for each treatment?” These linear equations take the following form:
Y = µ + f(treatment) + error
where
In this module we provide two examples of ANOVA and sample data sets to assess differences in treatment effect. In the first example, four methods of soybean transformation are evaluated to determine whether transformation method affects expression of a stress-response gene. In the second example, two molecular markers are evaluated to determine whether genotype of each molecular marker results in differences in disease severity in a BC1 population.
ANOVA is a statistical tool that has applications to experiments in which we want to assess whether there is a difference in a continuous variable between treatment groups. In a plant breeding context, this page demonstrated the utility of ANOVA in gene expression studies and molecular marker analysis.
Many statistics textbooks provide a good discussion of theory and applications of ANOVA. A few examples are listed below.
The following videos provide detailed instructions for calculating components of ANOVA tables (ANOVA1 and 2) and hypothesis testing (ANOVA3).
Development of this page was supported in part by the National Institute of Food and Agriculture (NIFA) Solanaceae Coordinated Agricultural Project, agreement 2009-85606-05673, administered by Michigan State University. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the United States Department of Agriculture.
PBGworks 650
Author:
David M. Francis, The Ohio State University
If you have problems viewing this video connect with our YouTube channel or see the YouTube troubleshooting guide.
Development of this page was supported in part by the National Institute of Food and Agriculture (NIFA) Solanaceae Coordinated Agricultural Project, agreement 2009-85606-05673, administered by Michigan State University. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the United States Department of Agriculture.
PBGworks 620