Source: NORTH CAROLINA STATE UNIV submitted to
A GENETICS-BASED DATA ANALYSIS SYSTEM FOR BREEDERS IN POLYPLOID BREEDING PROGRAMS
Sponsoring Institution
National Institute of Food and Agriculture
Project Status
NEW
Funding Source
Reporting Frequency
Annual
Accession No.
1027948
Grant No.
2022-67013-36269
Project No.
NC09927
Proposal No.
2021-07591
Multistate No.
(N/A)
Program Code
A1141
Project Start Date
Jan 1, 2022
Project End Date
Dec 31, 2024
Grant Year
2022
Project Director
Mollinari, M.
Recipient Organization
NORTH CAROLINA STATE UNIV
(N/A)
RALEIGH,NC 27695
Performing Department
Horticultural Science
Non Technical Summary
Many important agricultural species are polyploids, i.e., have multiple copies of their genomes. They range from staple food crops (potato, sweetpotato) to fruits (strawberry, kiwi, blueberry, banana), ornamental flowers (roses, chrysanthemum), forage crops, turfgrass, and sugar and energy production crops (sugarcane). The transmission of genic material across generations in polyploids is much more intricate and challenging to unravel than diploids, such as maize, rice, and soybeans. Although challenging, understanding inheritance patterns is essential information in breeding programs. With the correct assessment of these patterns, it is possible to associate specific genomic positions to important traits or even find the gene responsible for them and use this information in breeding programs.In the last few years, we have developed a series of computational tools to help breeders and geneticists answer these questions by analyzing genomic data in polyploid species. We developed tools such as VCF2SM and SuperMASSA for processing raw DNA sequences and identifying genetic markers, MAPpoly for constructing genetic maps, and QTLpoly for locating important genes to trait phenotypes also, to perform prediction. Currently, our tools are implemented for limited genetic design: full-sib families, and we are extending to multiple families. This project proposes extending, even more, our previous polyploid genomic tool for general multiple-generation pedigree breeding populations typically present in practical polyploid breeding programs. Moreover, we propose developing a new downstream computational tool, called DecisionPoly, user-friendly and offers clearly illustrated actionable information to assist polyploid breeders in making short- and long-term breeding decisions based on the collected and learned information about their breeding populations different breeding objectives.
Animal Health Component
0%
Research Effort Categories
Basic
50%
Applied
30%
Developmental
20%
Classification

Knowledge Area (KA)Subject of Investigation (SOI)Field of Science (FOS)Percent
2012499108050%
2012410108150%
Goals / Objectives
The main goal of this project is to develop a comprehensive, integrated, open-source, and publicly available pipeline data analysis platform to process genomic data, to infer the complex inheritance patterns from parents to offspring, to map genes that are important for breeding objectives, and to offer breeders actionable information to make short and long-term breeding decisions in practical breeding programs for polyploid species. In its upstream, we aim to develop computational tools to deal with different genomic data, call biallelic and multiallelic markers, combine genetic and genomic information from multiple breeding families, and build complex genetic models from complex population structures. In its downstream, we aim to develop a user-friendly computational tool offering illustrated actionable information to assist breeders in making breeding decisions based on the collected and learned information about their breeding populations for different breeding objectives.For this project, the specific objectives are:To further develop haplotyping algorithms to consider all the relevant information from complex breeding schemes with multiple-generation and partially inter-connected multiple families.To extend the genetic models between genotypes and phenotypes for all scenarios presented in item 1 and build sound and efficient statistical analysis procedures to achieve the purposes of genetic discovery and more accurate prediction.To implement a user-friendly computational tool in Shiny-R language to help breeders make short and long-term decisions in their breeding programs. Breeders will be able to use the results obtained in items 1 and 2 to make educated short-term decisions, such as which individuals to select and mate, as well as long-term decisions which will be attained through interactive breeding exercises (predict outcomes for various breeding decisions); optimize decisions based on breeding objectives; multiple-generation forward simulation exercises (predict results for different breeding strategies).
Project Methods
Haplotype inference in complex pedigrees.To infer the haplotypes in complex pedigrees, we will extend our previous work on constructing genetic maps in full-sib families. Our solution uses multilocus hidden Markov model (HMM) analysis and works for even ploidy levels from diploids up to autooctaploids. Our current model uses input probability-based dosage markers and recovers the multiple polyploid genotypes present in the segregating population, otherwise masked by the biallelic nature of the SNPs. Given the flexibility of the HMM framework, We also expect to extend our model to use multiallelic markers. The concept of HMM in genetic mapping is to use multiple linked markers to estimate the parental linkage phase, the genetic distance between markers, and reconstruct the offspring haplotypes. This procedure can correct many potential genotyping errors.Using the various sources of interconnected evidence (multiple SNPs and individuals), HMM can aggregate multiple SNP information and partially recover the intrinsic high error rate in individual marker dosage callings in polyploid species. Multilocus analysis to construct genetic maps and offspring haplotypes in polyploids is extremely important for information recovery and marker data quality control. We will implement these innovations in our open source publically available R package MAPpoly.Genetic models between genotypes and phenotypes in complex pedigreesTo integrate all the relevant information (genomic markers and trait phenotypes) in a complex pedigree breeding population for joint and informed analysis, we need a cohesive quantitative genetic model that applies to the whole breeding population. A quantitative genetic model can be devised based on the alleles of the population founders, assuming that all the segregating alleles of breeding individuals can be traced back to the founders in probability. The challenge is how to perform the genetic analysis efficiently and informatively.There are two strategies. One is to extend our current QTL random-effect genetic model for a full-sib family (implemented in QTLpoly) to a complex multi-generational breeding population with the alleles and their effects defined in terms of the founders' alleles. Like computing the G matrix for GBLUP (for GS), a corresponding G matrix (or Q matrix) can be computed for each targeted QTL locus. Multiple Q matrices can be built in one model for multiple QTL. The statistical importance of each QTL locus can be evaluated in the variance component. This strategy is similar to the one proposed by Amadeu et al. (2020) for QTL mapping in tetraploid dialleles in their software DiaQTL. However, there needs to be a more efficient way to identify multiple QTL in this case, as a sequential genome search is not computationally efficient. The second strategy may help in this respect.The second strategy is: first to build a large set of the founder's allele effects sampled in the genome at every specified genome interval position; then to put that set of allelic effects in a LASSO analysis (Tibshirani, 1996) to shrink them to a small set; and from this set to identify potential QTL positions for evaluation. Some combinations of the two strategies can help identify both QTL and their relevant allelic effects computationally more manageable way. We can first try to find significant QTL additive effects, and then conditional on those, try to find their essential dominance effects. This analysis can be used for genetic discovery (identification of QTL and additive, dominant, and epistatic allelic effects). It can also be optimized to predict the performance of future generation individuals. This prediction can be used for selection and mating design and optimizing the experimental designs (which can be aided by a forward simulation). Identifying significant dominant allelic effects is particularly important for variety development because they contribute to heterosis or special combining ability. Of course, our computational tool will also have an option to compute the standard genome G matrix and use it to do GBLUP for GS, at least for comparison. We will implement these innovations in our open source publically available R package QTLpoly.Development of a terminal tool for breeders - DecisionPolyTo truly help breeders incorporate genomic information in their breeding programs, we need to make extra efforts to put the downstream tool in the hands of breeders to help them make breeding decisions directly. For this purpose, we ask ourselves this question: What do breeders need? We think they may need a computational tool that can assist them in making the short and long-term breeding decisions based on the collected and learned information about their breeding populations for different breeding objectives.Data analysis of using the upstream computational tools (SuperMASSA, MAPpoly, and QTLpoly) requires quite a bit of proficiency in using R programs and understanding the scientific and technical issues behind the tools. Thus, we will implement DecisonPoly, an easy-to-use interactive application that can access the relevant information from MAPpoly and QTLpoly, display various information and results graphically, provide various options for users to perform further breeding analysis and simulations, and give breeders a direct control on the decision-making process.We will implement this application using the multiplatform language Shiny. Shiny is an R package that makes it easy to build interactive web apps straight from R with direct connections with databases and the other upstream R packages (MAPpoly and QTLpoly). It is highly interactive, can perform analysis, directly communicates with users with graphic results, and tells the story from the data. We can host this standalone app on a webpage or embed it in R Markdown documents or build dashboards.We will develop our tools using the Breading Application Programing Interface (BrAPI) standardized technical specifications (https://brapi.org/) to facilitate communications between different breeding platforms. We also will implement other ways to import and export datasets, such as CSV and Excell files, since these formats are prevalent in many breeding programs. We will also follow the FAIR guidelines to make our analysis, codes, and datasets findable, accessible, interoperable, and reproducible to the scientific community, especially breeders.

Progress 01/01/22 to 12/31/22

Outputs
Target Audience:We aimed to engage breeders and graduate students with an interest or active involvement in polyploid breeding. The diverse plant species studied by these groups spanned blueberries, potatoes, sweetpotatoes, roses, yams, ornamental flowers, and blackberries, among others. Our target audience included participants from both domestic and international personnel and institutions. Changes/Problems: In our initial proposal, we identified the Mwanga Diversity Population (MDP) as the main resource for providing a representative sample of genetic data to be used in developing and testing our models and implementations. However, despite our collaborator's efforts to obtain good-quality DNA from the MDP materials and several adjustments to the genotyping protocol, we were unable to obtain the genotype information of the MDP progeny. As a result, testing the efficiency of the multi-parental model in that scenario has not been possible. This work is ongoing, and we expect to acquire the necessary data in the next reporting cycle. In the meantime, we have established collaborations with other breeding groups, including the International Potato Center (CIP), which provided us with a population that allowed us to construct a consensus genetic map and reconstruct the haplotypes in the progeny. Additionally, we are collaborating with the Texas A&M University Rose Breeding Genetic program, which supplied an interconnected rose population with varying ploidy levels, enabling us to test our procedures. We also decided to change the name of our user-friendly down-stream tool from DecisionPoly to GGSpoly (GGS for Genetic and Genomic Selection) - to be implemented in the next reporting cycles. What opportunities for training and professional development has the project provided?We are working in collaboration with several graduate students on projects directly involved with this project. For the graduate students, we could point out Simon Fraher, who successfully used our tools to identify a single major QTL that explained 70% of the variation in resistance to the nematode Meloidogyne enterolobii and currently developing markers to perform assisted breeding selection at the NCSU sweetpotato breeding program. We are collaborating closely with Gabriel Gesteira, a postdoctoral researcher at NCSU specializing in genotype-phenotype associations, and Cristiani Tanigiuti, a postdoctoral fellow at Texas A&M University, who is involved in multi-population haplotype construction and the development of user-friendly tools. How have the results been disseminated to communities of interest? Formal classroom instruction: We presented our work in guest lectures at NCSU, including Plant Cytogenetics in Plant Breeding, Breeding Asexually Propagated Crops, and Quantitative Genetics Theory and Methods. Workshops and conferences: We participated in the following events to share our knowledge and expertise: a. American Society for Horticultural Science Conference b. 30th Plant and Animal Genome Conference c. Tools for Polyploids workshop, attended by 238 participants d. Advancing Computing Skills in Plant Breeding workshop, organized by the NCSU plant breeding consortium, with an average of 30 attendees Extension and outreach: The talks from the Tools for Polyploids workshop were made available in video format on the Tools for Polyploids webpage at https://www.polyploids.org/2023recordedpresentations. In addition, all source codes for this project are freely available in the git repository https://github.com/mmollina. What do you plan to do during the next reporting period to accomplish the goals?Objective1 and 2. a) Continue to improve our phasing algorithm and explore possibilities to include pedigree with small family sizes and multiple generations. b) Finish the implementation of MAPpoly 2.0 and improve the implementation of MAPpoly-MP and QTLpoly-MP. Objective3. Leveraging our VIEWpoly visualization framework for polyploid genetic analysis (https://cran.r-project.org/package=viewpoly), we plan to commence the development and implementation of GGSpoly in the upcoming reporting cycle. GGSpoly is a computational tool intended to aid breeders in making informed short and long-term breeding decisions based on gathered and processed data about their breeding populations for diverse objectives. We anticipate collaborating closely with breeding groups such as the NCSU Sweetpotato and Potato Breeding Team led by Craig Yencho and the Texas A&M Rose Breeding Group led by Oscar Riera-Lizarazu and David Byrne in order to optimize GGSpoly's effectiveness in real-world breeding situations.

Impacts
What was accomplished under these goals? Objective 1: We successfully implemented mapping and haplotype reconstruction algorithms for polyploid partially inter-connected families, including diploid (2x), tetraploid (4x), and hexaploid (6x) families. Our current implementation allows for the integration of families with varying ploidy levels, enabling the simultaneous reconstruction of genetic maps and offspring haplotypes in single-generation populations given a sufficiently large number of individuals within full-sib populations. The performance of these algorithms was verified using real-world populations, such as a tetraploid potato in a partial diallel 3 x 9 parental configuration and three interconnected hexaploid sweetpotato populations, as well as several in silico simulations. Furthermore, we have optimized the algorithms employed in creating individual genetic maps for individual parents. Our preliminary evaluations demonstrated a substantial decrease in processing time, exhibiting a 10-fold improvement for tetraploid families and a 100-fold improvement for hexaploid families when tested with real data. This refinement is crucial for future implementation. The analysis of markers in single parents will act as an initiation algorithm for merging genetic maps in multiparental families in subsequent phases. Currently, we are implementing these functions in a new version of our mapping software MAPpoly which will incorporate several significant updates designed to enhance its capabilities and efficiency. Objective 2. We also made considerable advancements in connecting the haplotypes obtained in Objective 1 to phenotypes of interconnected populations. Our algorithm is an extension of the algorithm previously developed by Pereira et al. (2020). It relies on a random effect model which is applied in the context of the Multiple Interval Mapping (MIM) procedure. To evaluate and demonstrate the utility of our approach, we applied it to the same three interconnected hexaploid sweetpotato populations described in Objective 1, and also used several in silico simulations. With the posterior haplotype probabilities, we performed multiple QTL mapping and detected significant QTL for beta-carotene, dry matter, and starch content, with consistent allele effects across sub-populations. Objective 3: Nothing to report during this period.

Publications

  • Type: Conference Papers and Presentations Status: Published Year Published: 2023 Citation: Mollinari M. Computational Tools for Genomic Analysis in Polyploids in: 2022 ASHS Annual Conference (presenter) https://ashs.confex.com/ashs/2022/meetingapp.cgi/Paper/38640
  • Type: Conference Papers and Presentations Status: Published Year Published: 2023 Citation: Gesteira G.S., Mollinari M., Pereira G.da S., Olukolu B.A., Oloka B.M., Yencho C., Zeng Z-B. Genetic Mapping in Interconnected Hexaploid Sweetpotato Populations in: International Plant & Animal Genome 30 Conference (presenter) https://plan.core-apps.com/pag_2023/abstract/d8c2bdba-5e59-4f12-ab32-9e01d221a7a0
  • Type: Conference Papers and Presentations Status: Published Year Published: 2023 Citation: Wu S., Sun H., Kitavi M., Hamilton J.P., Gesteira G.S., Mollinari M., Zeng Z-B., Yencho C., Buell R., Fei Z. Advances in the Development of Chromosome-Scale Haplotype-Resolved Genome Assemblies of Hexaploid Sweetpotatoes in: International Plant & Animal Genome 30 Conference https://plan.core-apps.com/pag_2023/abstract/18b6ff6a-ca63-4074-ae87-e4e99fea4f11