Skip to main content
Biology LibreTexts

8: Genome Construction

  • Page ID
    107316
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)
    Photo of a microscope on a table in front of a line of maize samples in bowls.
    Fig. 1 Operations research is a tool to address crop breeding objectives. Photo by Iowa State University.

    One of the main challenges in plant breeding is the development of the best marker-assisted breeding method for complex traits. At the present, marker-based approaches are limited in their ability to detect and quantify marker-trait relationships, in particular for traits that are under the influence of gene x gene and gene x environment interactions. Also, as you have learned in previous lessons of this course, QTL estimates are biased by population size and a limited set of environments, making QTL estimates less suitable for crop improvement. For this reason, simulation modeling is an emerging important tool to choose among proposed breeding methods because experimental evaluation of breeding methods is time and resource-limited. Another challenge is management of multiple breeding objectives for several complex traits, making it more likely that an operations research approach called multi-objective optimization will gain favor in crop breeding. Thus, this lesson will introduce operations research as a tool to address multiple crop breeding objectives.

    Learning Objectives
    • Summarize and state the concept of genetic gain
    • Introduce the concept of multi-objective optimization
    • Introduce the concept of operations research in plant breeding
     
    Plants in test tubes.
    Fig. 2 Evaluating materials requires expert training, patience, equipment and time. Photo by Iowa State University.

    Recapitulation of the Concept of Genetic Gain

    Definition

    Genetic gain (∆G) is defined as the predicted change in the mean value of a trait within a population as a result of selection. The ∆G equation (Fig. 3) allows the comparison of predicted effectiveness of particular breeding methods and helps breeders decide how resources should be allocated for achieving various breeding objectives.

    Simple distribution graph with equation for genetic gain. Described in caption.
    Fig. 3 The genetic gain equation and its components. The curve illustrates the distribution showing frequency of individuals in a breeding population (y axis) that display various phenotypic values (x axis) of individuals in a breeding population. For the above curve, the mean phenotypic value of the original population is denoted µ0, and the mean phenotypic value for the selected individuals is denoted µS. Genetic components (σ2) and phenotypic distribution (σP) are indicated. Adapted from Moose and Mumm, 2008.

    Commercialization Challenges

    Figure 4 illustrates a generic plant breeding program involving mating, evaluation, selection, and testing of breeding materials resulting in commercialization of a cultivar. Such a program faces the challenges of time to commercialization of a cultivar, and resources allocated to obtain such cultivar from thousands of individuals.

    Process for a breeding program. Parents are selected, image of distribution graph, create genetic variability, evaluate, test, strip trials, and commercialization.
    Fig. 4 An example of a breeding program.

    Crop Yield Progress

    Despite such challenges, from the 1940s, the yields of corn and soybean in the United States have continued to rise (Fig. 5) mainly due to improvement in crop genetics and agronomic practices.

    Corn and soybean yield graphs. Corn yields similar results in both states; soybeans are more reliable in high yield states.
    Fig. 5 Average corn and soybean yields (in kg/ha) for the U.S. from 1866 to 2015 in high-yield potential states (Iowa, Illinois, and Indiana) and low-yield potential states (Kentucky, Missouri, and Tennessee). Data from National Agricultural Statistics Service.

    Global Food Demand Trends

    Layered bar chart of food production from 1964-66, 1997-99, and 2023, estimated. General growth across all areas.
    Fig. 6 Global progress in food consumption. Adapted from FAO, 2002.

    Despite the upward trend in crop yields in the US and other parts of the world, rising human and animal populations will pose a greater demand for more to be produced per unit of land. The growing global demand for food (Fig. 6) raises the question of whether it is possible to double the current level of production in the next 20 years (Fig. 7). Undoubtedly, to reach 300 bushels/acre of corn by 2030 will require cutting-edge approaches in genomics and breeding. But the problem will be the cost of reaching such a high level of yield with limited time and resources. Thus, integration of new approaches, for example, Genomic Selection, transgenics, and operations research, may be necessary. The next lesson sections entail application of operations research tools in plant breeding as a novel approach to increase ∆G.

    Need for Advancement

    Historically plant breeding has been a form of art: to create new varieties. Thus, ∆G has depended on management of resources, to produce new varieties; while optimization has been ignored. Nonetheless, plant breeding has the potential to become an engineering discipline, relying on operations research, which will be necessary for average yields to double by 2030 (Fig. 7).

    Line chart of bushels per acre over timefor US corn yields, 1985 to 2014. If current trends continue, the goal of 300 bushels per acre won't be hit until 2092.
    Fig. 7 New plant breeding tools will be needed to produce 300 bushels (Bu) per acre by 2030. Data from National Agricultural Statistics Service.

    Multi-Objective Optimization

    Area graph showing genetic diversity by number of successions. A steep rise is shown in the graph.
    Fig. 8 Some multi-objective optimization problems involve maximizing one variable while minimizing another.

    Introduction

    Multi-objective optimization (MO) is an operations research approach used in various fields, including engineering, finance, biomedicine and management. Optimization involves application of more than one objective processes for evaluation that can take into account multiple criteria that need to be considered for making a decision. Therefore, as information on plant genomes continue to emerge, it is now possible to apply the MO approach for large scale plant breeding (Xu et al., 2011). For example, a plant breeding goal may have two objectives, 1) selection and fixation of desirable genes at a set of loci controlling a trait of interest, and 2) keeping genetic variability at the remaining loci to retain adaptability. The challenge of applying MO to solve these competing objectives is identification of optimal solutions to the problems (Chinchuluun and Pardalos, 2007). Such solutions are called Pareto optimal solutions, and they are a measure of MO optimization efficiency.

    Pareto optimal solutions

    We will not dwell on the mathematics used to derive Pareto optimal solutions in this lesson. But it is important to know that there usually exist multiple Pareto optimal solutions for MO problems, and searching for all Pareto optimal solutions can be expensive and time consuming (Chinchuluun and Pardalos, 2007). Nonetheless, recent advances in computational research suggest that it is possible to obtain Pareto optimal solutions for plant breeding problems within reasonable computation time (Xu et al., 2011). Such solutions will be useful tools to help plant breeders make informed decisions in the world of large amounts of genomics data for multiple breeding objectives for complex traits. 2011). Such solutions will be useful tools to help plant breeders make informed decisions in the world of large amounts of genomics data for multiple breeding objectives for complex traits.

    Decorative image.
    Fig. 9 Computer systems can now yield Pareto optimal solutions in acceptable time frames.

    Operations Research in Plant Breeding

    Operations Research involves the application of mathematical models to provide optimal solutions to a problem. An OR approach (Fig. 10) consists four components, 1) Problem, 2), Model, 3) Algorithm, and 4) Solver.

    The words Problem, model, algorithm, and solver with arrows pointing between each item.
    Fig. 10 A multi-objective optimization plant breeding problem requires the use of optimization models, algorithms, and computer technologies. Complexity of the problem, robustness of model, algorithm used, and computer solvers influence the cost of solving the problem.

    Step 1: Defining the Problem

    There is an original population of individuals (Fig. 11). Each individual has a pair of chromosomes, and each chromosome has a number of genes. Some genes are undesirable, while the desirable ones have different variants. The desirable genes will be assigned a value of 1 and undesirable 0. What is the best way to assemble all variants of desirable genes into a target population?

    Integer Programming from an original population to a target population. Rows of numbers and letters are represented.
    Fig. 11 Operations research can help assess the possibility of stacking genes into multiple backgrounds.

    Step 2: Developing a Model

    A model has four key elements – data, decisions, objective, and constraints.

    1. Data
    2. Decisions – A decision would have to be made about number of data and recombination points, and number of chromosomes in the target population.
    3. Objective – The objective is to maximize probability of getting the target population.
    4. Constraints – Constraints can be, for example, the number of chromosomes in target population without undesirable alleles, but such that all desirable variants are retained. Also, the maximum number of recombination events could be another constraint.

    Step 3: Designing a Suitable Algorithm

    The problem in this example belongs to a class of so-called non-deterministic polynomial-time hard (NP-hard) problems (Xu et al. 2011). Importantly, if an algorithm solves one NP-hard problem, it can be used to solve all other NP-hard problems.

    Decorative image.
    Fig. 12 Often a computer program will be used as an algorithm for solving NP-hard problems, such this one for the Traveling Salesman Problem.

    Step 4: Solving the Problem

    Computation time spent to solve the problem in Figure 11 was 0.03 seconds (W.D. Beavis, personal communication).

    Integer programming image with a photo of someone holding a stopwatch overlaid on it.

    Genome Construction vs. Genomic Selection

    The hypothesis is that genome construction is better than genomic selection (Fig. 13). The hypothesis is developed from the premise that a target genotype can be defined. However, if the target genotype has to be determined using experimental methods, then GS will be more effective because experimental methods are underpowered and biased.

    Two lines graphs for genome construction and selection. Genome construction sees genetic gain over generations in a step-pattern, whereas selection sees genetic gain in a curve over time.
    Fig. 13 Comparison of genome construction and genomic selection approaches

    References

    Chinchuluun, A. P. M. Pardalos. 2007. A survey of recent developments in multiobjective optimization. Ann. Oper. Res. 154: 29-50.

    Egli, D. B. 2008. Comparison of corn and soybean yields in the United States: Historical trends and future prospects. Agron J. 100: S-79-S-88.

    Expert Meeting on “How to Feed the World in 2050,” FAO, Rome, 24-26 June 2009.

    FAO, 2002. World agriculture: towards 2015/2030. United Nations, 2002.

    Moose, S. P., and R. H. Mumm. 2008. Molecular plant breeding as the foundation for 21st century crop improvement. Plant Physiol. 147: 969-977.

    Xu, P., L. Wang, and W. D. Beavis. 2011. An optimization approach to gene stacking. Europ. J. Oper. Res. 214: 168-178.

    How to cite this module: Lübberstedt, T., W. Beavis, and W. Suza. (2023). Genome Construction. In W. P. Suza, & K. R. Lamkey (Eds.), Molecular Plant Breeding. Iowa State University Digital Press.


    8: Genome Construction is shared under a CC BY-NC 4.0 license and was authored, remixed, and/or curated by Walter Suza, Thomas Lübberstedt, William Beavis, & William Beavis.