Quantifying the Power and Precision of QTL Analysis in Autopolyploids Under Bivalent and Multivalent Genetic Models

GenAPoPop 1.0: A user-friendly software to analyse genetic diversity and structure from partially clonal and selfed autopolyploid organisms

Autopolyploidy is quite common in most clades of eukaryotes. The emergence of sequence-based genotyping methods with individual and marker tags now enables confident allele dosage, overcoming the main obstacle to the democratization of the population genetic approaches when studying ecology and evolution of autopolyploid populations and species. Reproductive modes, including clonality, selfing and allogamy, have deep consequences on the ecology and evolution of population and species. Analysing genetic diversity and its dynamics over generations is one efficient way to infer the relative importance of clonality, selfing and allogamy in populations. GenAPoPop is a user-friendly solution to compute the specific corpus of population genetic indices, including indices about genotypic diversity, needed to analyse partially clonal, selfed and allogamous polysomic populations genotyped with confident allele dosage. It also easily provides the posterior probabilities of quantitative reproductive modes in autopolyploid populations genotyped at two-time steps and a graphical representation of the minimum spanning trees of the genetic distances between polyploid individuals, facilitating the interpretation of the genetic coancestry between individuals in hierarchically structured populations. GenAPoPop complements the previously existing solutions, including SPAGEDI and POLYGENE, to use genotypings to study the ecology and evolution of autopolyploid populations. It was specially developed with a simple graphical interface and workflow, and comes with a simulator to facilitate practical courses and teaching of population genetics for autopolyploid populations.

Development of a 135K SNP genotyping array for Actinidia arguta and its applications for genetic mapping and QTL analysis in kiwifruit

Kiwifruit (Actinidia spp) is a woody, perennial and deciduous vine. In this genus, there are multiple ploidy levels but the main cultivated cultivars are polyploid. Despite the availability of many genomic resources in kiwifruit, SNP genotyping is still a challenge given these different levels of polyploidy. Recent advances in SNP array technologies have offered a high-throughput genotyping platform for genome-wide DNA polymorphisms. In this study, we developed a high-density SNP genotyping array to facilitate genetic studies and breeding applications in kiwifruit. SNP discovery was performed by genome-wide DNA sequencing of 40 kiwifruit genotypes. The identified SNPs were stringently filtered for sequence quality, predicted conversion performance and distribution over the available Actinidia chinensis genome. A total of 134 729 unique SNPs were put on the array. The array was evaluated by genotyping 400 kiwifruit individuals. We performed a multidimensional scaling analysis to assess the diversity of kiwifruit germplasm, showing that the array was effective to distinguish kiwifruit accessions. Using a tetraploid F1 population, we constructed an integrated linkage map covering 3060.9 cM across 29 linkage groups and performed QTL analysis for the sex locus that has been identified on Linkage Group 3 (LG3) in Actinidia arguta. Finally, our dataset presented evidence of tetrasomic inheritance with partial preferential pairing in A. arguta. In conclusion, we developed and evaluated a 135K SNP genotyping array for kiwifruit. It has the advantage of a comprehensive design that can be an effective tool in genetic studies and breeding applications in this high-value crop.

QTL mapping in outbred tetraploid (and diploid) diallel populations

Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multiallelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and nonadditive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with three tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher-order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.

Genomics-based discrimination of 2n gamete formation mechanisms in polyploids: a case study in nonaploid Diospyros kaki 'Akiou'

Unreduced gametes (2n gametes), possessing double the haploid genome, whatever ploidy that happens to be, are a common source of ploidy variation in plant populations. First and second division restitution (FDR and SDR) are the dominant mechanisms of 2n gamete production; all else being equal, FDR gametes have a higher degree of heterozygosity, thus they are advantageous in breeding. The discrimination of these mechanisms from the consequence of hybridization is challenging, especially in higher polyploids, and usually requires information on centromere location. In this study, we propose a genotyping-based strategy to uncover the mechanisms of 2n gamete formation in progeny that has a higher ploidy than its parents. Simulation of 2n gamete production revealed that FDR and SDR pathways can be discriminated based on allele transmission patterns alone without information on centromere location. We applied this strategy to study the formation mechanism of a nonaploid Diospyros kaki ‘Akiou', which was bred via hybridization between D. kaki hexaploid cultivars. The result demonstrated that ‘Akiou' was derived from the fertilization of a normal female gamete by a 2n male gamete and that this 2n gamete was produced through FDR. Consequently, the distinct duplex transmission pattern in the FDR gamete enabled us to infer the genomic characteristics of polyploid persimmon. The method could be tested only for the plant being polypoid, which allows for the ability to discriminate causes of 2n gamete formation using allele dosage in progeny, and will be useful in future studies of polyploid genomics.

A tetrasomic inheritance model and likelihood-based method for mapping quantitative trait loci in autotetraploid species

Dissecting the genetic architecture of quantitative traits in autotetraploid species is a methodologically challenging task, but a pivotally important goal for breeding globally important food crops, including potato and blueberry, and ornamental species such as rose. Mapping quantitative trait loci (QTLs) is now a routine practice in diploid species but is far less advanced in autotetraploids, largely due to a lack of analytical methods that account for the complexities of tetrasomic inheritance. We present a novel likelihood-based method for QTL mapping in outbred segregating populations of autotetraploid species. The method accounts properly for sophisticated features of gene segregation and recombination in an autotetraploid meiosis. It may model and analyse molecular marker data with or without allele dosage information, such as that from microarray or sequencing experiments. The method developed outperforms existing bivalent-based methods, which may fail to model and analyse the full spectrum of experimental data, in the statistical power of QTL detection, and accuracy of QTL location, as demonstrated by an intensive simulation study and analysis of data sets collected from a segregating population of potato (Solanum tuberosum). The study enables QTL mapping analysis to be conducted in autotetraploid species under a rigorous tetrasomic inheritance model.

QTL Mapping for Resistance to Cankers Induced by Pseudomonas syringae pv. actinidiae (Psa) in a Tetraploid Actinidia chinensis Kiwifruit Population

Polyploidy is a key driver of significant evolutionary changes in plant species. The genus Actinidia (kiwifruit) exhibits multiple ploidy levels, which contribute to novel fruit traits, high yields and resistance to the canker-causing dieback disease incited by Pseudomonas syringae pv. actinidiae (Psa) biovar 3. However, the genetic mechanism for resistance to Psa observed in polyploid kiwifruit is not yet known. In this study we performed detailed genetic analysis of a tetraploid Actinidia chinensis var. chinensis population derived from a cross between a female parent that exhibits weak tolerance to Psa and a highly Psa-resistant male parent. We used the capture-sequencing approach across the whole kiwifruit genome and generated the first ultra-dense maps in a tetraploid kiwifruit population. We located quantitative trait loci (QTLs) for Psa resistance on these maps. Our approach to QTL mapping is based on the use of identity-by-descent trait mapping, which allowed us to relate the contribution of specific alleles from their respective homologues in the male and female parent, to the control of Psa resistance in the progeny. We identified genes in the diploid reference genome whose function is suggested to be involved in plant defense, which underly the QTLs, including receptor-like kinases. Our study is the first to cast light on the genetics of a polyploid kiwifruit and suggest a plausible mechanism for Psa resistance in this species.

Diversity in Metabolites and Fruit Quality Traits in Blueberry Enables Ploidy and Species Differentiation and Establishes a Strategy for Future Genetic Studies

Blueberry is well recognized as a rich source of health promoting phytochemicals such as flavonoids and phenolic acids. Multiple studies in blueberry and other crops indicated that flavonoids and phenolic acids function as bioactive compounds in the human body promoting multiple health effects. Despite their importance, information is limited about the levels of variation in bioactive compounds within and between ploidy level and species, and their association with fruit quality traits. Such information is crucial to define a strategy to study the genetic mechanisms controlling these traits and to select for these traits in blueberry breeding programs. Here we evaluated 33 health related phytochemicals belonging to four major groups of flavonoids and phenolic acids across 128 blueberry accessions over two years together with fruit quality traits, including fruit weight, titratable acidity, total soluble acids and pH. Highly significant variation between accessions, years, and accession by year interaction were identified for most of the traits. Cluster analysis grouped phytochemicals by their functional structure (e.g., anthocyanins, flavanols, flavonols, and phenolic acids). Multivariate analysis of the traits resulted in separation of diploid, tetraploid and hexaploid accessions. Broad sense heritability of the traits estimated in 100 tetraploid accessions, ranged from 20 to 90%, with most traits revealing moderate to high broad sense heritability (H² > 40%), suggesting that strong genetic factors control these traits. Fruit size can be estimated as a proxy of fruit weight or volume and vice versa, and it was negatively correlated with content of most of phytochemicals evaluated here. However, size-independent variation for anthocyanin content and profile (e.g., acylated vs. non-acylated anthocyanin) exists in the tetraploid accessions and can be explored to identify other factors such as genes related to the biosynthetic pathway that control this trait. This result also suggests that metabolite concentrations and fruit size, to a certain degree can be improved simultaneously in breeding programs. Overall, the results of this study provide a framework to uncover the genetic basis of bioactive compounds and fruit quality traits and will be useful to advance blueberry-breeding programs focusing on integrating these traits.

High-Resolution Linkage Map With Allele Dosage Allows the Identification of Regions Governing Complex Traits and Apospory in Guinea Grass (Megathyrsus maximus)

Forage grasses are mainly used in animal feed to fatten cattle and dairy herds, and guinea grass (Megathyrsus maximus) is considered one of the most productive of the tropical forage crops that reproduce by seeds. Due to the recent process of domestication, this species has several genomic complexities, such as autotetraploidy and aposporous apomixis. Consequently, approaches that relate phenotypic and genotypic data are incipient. In this context, we built a linkage map with allele dosage and generated novel information of the genetic architecture of traits that are important for the breeding of M. maximus. From a full-sib progeny, a linkage map containing 858 single nucleotide polymorphism (SNP) markers with allele dosage information expected for an autotetraploid was obtained. The high genetic variability of the progeny allowed us to map 10 quantitative trait loci (QTLs) related to agronomic traits, such as regrowth capacity and total dry matter, and 36 QTLs related to nutritional quality, which were distributed among all homology groups (HGs). Various overlapping regions associated with the quantitative traits suggested QTL hotspots. In addition, we were able to map one locus that controls apospory (apo-locus) in HG II. A total of 55 different gene families involved in cellular metabolism and plant growth were identified from markers adjacent to the QTLs and APOSPORY locus using the Panicum virgatum genome as a reference in comparisons with the genomes of Arabidopsis thaliana and Oryza sativa. Our results provide a better understanding of the genetic basis of reproduction by apomixis and traits important for breeding programs that considerably influence animal productivity as well as the quality of meat and milk.

Multi-environment QTL analysis of plant and flower morphological traits in tetraploid rose

KEY MESSAGE

Rose morphological traits such as prickles or petal number are influenced by a few key QTL which were detected across different growing environments-necessary for genomics-assisted selection in non-target environments. Rose, one of the world's most-loved and commercially important ornamental plants, is predominantly tetraploid, possessing four rather than two copies of each chromosome. This condition complicates genetic analysis, and so the majority of previous genetic studies in rose have been performed at the diploid level. However, there may be advantages to performing genetic analyses at the tetraploid level, not least because this is the ploidy level of most breeding germplasm. Here, we apply recently developed methods for quantitative trait loci (QTL) detection in a segregating tetraploid rose population (F1 = 151) to unravel the genetic control of a number of key morphological traits. These traits were measured both in the Netherlands and Kenya. Since ornamental plant breeding and selection are increasingly being performed at locations other than the production sites, environment-neutral QTL are required to maximise the effectiveness of breeding programmes. We detected a number of robust, multi-environment QTL for such traits as stem and petiole prickles, petal number and stem length that were localised on the recently developed high-density SNP linkage map for rose. Our work explores the complex genetic architecture of these important morphological traits at the tetraploid level, while helping to advance the methods for marker-trait exploration in polyploid species.