Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection.

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

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.

PolyHaplotyper: haplotyping in polyploids based on bi-allelic marker dosage data.

BACKGROUND: For genetic analyses, multi-allelic markers have an advantage over bi-allelic markers like SNPs (single nucleotide polymorphisms) in that they carry more information about the genetic constitution of individuals. This is especially the case in polyploids, where individuals carry more than two alleles at each locus. Haploblocks are multi-allelic markers that can be derived by phasing sets of closely-linked SNP markers. Phased haploblocks, similarly to other multi-allelic markers, will therefore be advantageous in genetic tasks like linkage mapping, QTL mapping and genome-wide association studies. RESULTS: We present a new method to reconstruct haplotypes from SNP dosages derived from genotyping arrays, which is applicable to polyploids. This method is implemented in the software package PolyHaplotyper. In contrast to existing packages for polyploids it makes use of full-sib families among the samples to guide the haplotyping process. We show that in this situation it is much more accurate than other available software, using experimental hexaploid data and simulated tetraploid data. CONCLUSIONS: Our method and the software package PolyHaplotyper in which it is implemented extend the available tools for haplotyping in polyploids. They perform especially well in situations where one or more full-sib families are present.

Genetic dissection of morphological variation in rosette leaves and leafy heads in cabbage (Brassica oleracea var. capitata).

KEY MESSAGE: Correlations between morphological traits of cabbage rosette leaves and heads were found. Genome-wide association studies of these traits identified 50 robust quantitative trait loci in multiple years. Half of these loci affect both organs. Cabbage (Brassica oleracea var. capitata) is an economically important vegetable crop cultivated worldwide. Cabbage plants go through four vegetative stages: seedling, rosette, folding and heading. Rosette leaves are the largest leaves of cabbage plants and provide most of the energy needed to produce the leafy head. To understand the relationship and the genetic basis of leaf development and leafy head formation, 308 cabbage accessions were scored for rosette leaf and head traits in three-year field trials. Significant correlations were found between morphological traits of rosette leaves and heads, namely leaf area with the head area, height and width, and leaf width with the head area and head height, when heads were harvested at a fixed number of days after sowing. Fifty robust quantitative trait loci (QTLs) for rosette leaf and head traits distributed over all nine chromosomes were identified with genome-wide association studies. All these 50 loci were identified in multiple years and generally affect multiple traits. Twenty-five of the QTL were associated with both rosette leaf and leafy head traits. We discuss thirteen candidate genes identified in these QTL that are expressed in heading leaves, with an annotation related to auxin and other phytohormones, leaf development, and leaf polarity that likely play a role in leafy head development or rosette leaf expansion.

Multiallelic models for QTL mapping in diverse polyploid populations.

Quantitative trait locus (QTL) analysis allows to identify regions responsible for a trait and to associate alleles with their effect on phenotypes. When using biallelic markers to find these QTL regions, two alleles per QTL are modelled. This assumption might be close to reality in specific biparental crosses but is unrealistic in situations where broader genetic diversity is studied. Diversity panels used in genome-wide association studies or multi-parental populations can easily harbour multiple QTL alleles at each locus, more so in the case of polyploids that carry more than two alleles per individual. In such situations a multiallelic model would be closer to reality, allowing for different genetic effects for each potential allele in the population. To obtain such multiallelic markers we propose the usage of haplotypes, concatenations of nearby SNPs. We developed "mpQTL" an R package that can perform a QTL analysis at any ploidy level under biallelic and multiallelic models, depending on the marker type given. We tested the effect of genetic diversity on the power and accuracy difference between bi-allelic and multiallelic models using a set of simulated multiparental autotetraploid, outbreeding populations. Multiallelic models had higher detection power and were more precise than biallelic, SNP-based models, particularly when genetic diversity was higher. This confirms that moving to multi-allelic QTL models can lead to improved detection and characterization of QTLs. KEY MESSAGE: QTL detection in populations with more than two functional QTL alleles (which is likely in multiparental and/or polyploid populations) is more powerful when using multiallelic models, rather than biallelic models.

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.

Detecting quantitative trait loci and exploring chromosomal pairing in autopolyploids using polyqtlR.

MOTIVATION: The investigation of quantitative trait loci (QTL) is an essential component in our understanding of how organisms vary phenotypically. However, many important crop species are polyploid (carrying more than two copies of each chromosome), requiring specialized tools for such analyses. Moreover, deciphering meiotic processes at higher ploidy levels is not straightforward, but is necessary to understand the reproductive dynamics of these species, or uncover potential barriers to their genetic improvement. RESULTS: Here, we present polyqtlR, a novel software tool to facilitate such analyses in (auto)polyploid crops. It performs QTL interval mapping in F1 populations of outcrossing polyploids of any ploidy level using identity-by-descent probabilities. The allelic composition of discovered QTL can be explored, enabling favourable alleles to be identified and tracked in the population. Visualization tools within the package facilitate this process, and options to include genetic co-factors and experimental factors are included. Detailed information on polyploid meiosis including prediction of multivalent pairing structures, detection of preferential chromosomal pairing and location of double reduction events can be performed. AVAILABILITYAND IMPLEMENTATION: polyqtlR is freely available from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polyqtlR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Using probabilistic genotypes in linkage analysis of polyploids.

KEY MESSAGE: In polyploids, linkage mapping is carried out using genotyping with discrete dosage scores. Here, we use probabilistic genotypes and we validate it for the construction of polyploid linkage maps. Marker genotypes are generally called as discrete values: homozygous versus heterozygous in the case of diploids, or an integer allele dosage in the case of polyploids. Software for linkage map construction and/or QTL analysis usually relies on such discrete genotypes. However, it may not always be possible, or desirable, to assign definite values to genotype observations in the presence of uncertainty in the genotype calling. Here, we present an approach that uses probabilistic marker dosages for linkage map construction in polyploids. We compare our method to an approach based on discrete dosages, using simulated SNP array and sequence reads data with varying levels of data quality. We validate our approach using experimental data from a potato (Solanum tuberosum L.) SNP array applied to an F1 mapping population. In comparison to the approach based on discrete dosages, we mapped an additional 562 markers. All but three of these were mapped to the expected chromosome and marker position. For the remaining three markers, no physical position was known. The use of dosage probabilities is of particular relevance for map construction in polyploids using sequencing data, as these often result in a higher level of uncertainty regarding allele dosage.

Fine mapping of a thrips resistance QTL in Capsicum and the role of diterpene glycosides in the underlying mechanism.

KEY MESSAGE: A major thrips resistance QTL in Capsicum was fine-mapped to a region of 0.4 Mbp, and a multidisciplinary approach has been used to study putative underlying mechanisms. Resistance to thrips is an important trait for pepper growers. These insects can cause extensive damage to fruits, flowers and leaves on field and greenhouse grown plants worldwide. Two independent studies in Capsicum identified diterpene glycosides as metabolites that are correlated with thrips resistance. In this study, we fine-mapped a previously defined thrips resistance QTL on chromosome 6, to a region of 0.4 Mbp harbouring 15 genes. Two of these 15 candidate genes showed differences in gene expression upon thrips induction, when comparing plants carrying the resistance allele in homozygous state to plants with the susceptibility allele in homozygous state for the QTL region. Three genes, including the two genes that showed difference in gene expression, contained a SNP that was predicted to lead to changes in protein structure. Therefore, these three genes, i.e. an acid phosphatase 1 (APS1), an organic cation/carnitine transporter 7 (OCT7) and an uncharacterized locus LOC107874801, are the most likely candidates for playing a role in thrips resistance and are a first step in elucidating the genetic basis of thrips resistance in Capsicum. In addition, we show that the diterpene glycoside profiles did not differ between plants with the resistance and susceptibility allele for the chromosome 6 QTL, suggesting that these compounds do not play a role in the resistance conferred by the genes located in the major thrips resistance QTL studied.

The ability to manipulate ROS metabolism in pepper may affect aphid virulence.

Myzus persicae has severe economic impact on pepper (Capsicum) cultivation. Previously, we identified two populations of M. persicae, NL and SW, that were avirulent and virulent, respectively on C. baccatum accession PB2013071. The transcriptomics approach used in the current study, which is the first study to explore the pepper-aphid interaction at the whole genome gene expression level, revealed genes whose expression is differentially regulated in pepper accession PB2013071 upon infestation with these M. persicae populations. The NL population induced ROS production genes, while the SW population induced ROS scavenging genes and repressed ROS production genes. We also found that the SW population can induce the removal of ROS which accumulated in response to preinfestion with the NL population, and that preinfestation with the SW population significantly improved the performance of the NL population. This paper supports the hypothesis that M. persicae can overcome the resistance in accession PB2013071 probably because of its ability to manipulate plant defense response especially the ROS metabolism and such ability may benefit avirulent conspecific aphids.

Aphid populations showing differential levels of virulence on Capsicum accessions.

The green peach aphid, Myzus persicae, is one of the most threatening pests in pepper cultivation and growers would benefit from resistant varieties. Previously, we identified two Capsicum accessions as susceptible and three as resistant to M. persicae using an aphid population originating from the Netherlands (NL). Later on we identified an aphid population originating from a different geographical region (Switserland, SW) that was virulent on all tested Capsicum accessions. The objective of the current work is to describe in detail different aspects of the interaction between two aphid populations and two selected Capsicum accessions (one that was susceptible [PB2013046] and one that was resistant [PB2013071] to population NL), including biochemical processes involved. Electrical penetration graph (EPG) recordings showed similar feeding activities for both aphid populations on PB2013046. On accession PB2013071 the aphid population SW was able to devote significantly more time to phloem ingestion than population NL. We also studied plant defense response and found that plants of accession PB2013046 could not induce an accumulation of reactive oxygen species and callose formation after infestation with either aphid population. However, plants of PB2013071 induced a stronger defense response after infestation by population NL than after infestation by population SW. Based on these results, population SW of M. persicae seems to have overcome the resistance of PB2013071 that prevented feeding of aphids from NL population. The potential mechanism by which SW population overcomes the resistance is discussed.

Aphid resistance in Capsicum maps to a locus containing LRR-RLK gene analogues.

KEY MESSAGE: A QTL for aphid resistance on pepper chromosome 2 was identified and validated. This QTL affects aphid survival and reproduction, and was fine mapped to a locus containing LRR-RLK analogues. Myzus persicae is one of the most threatening insect pests that adversely affects pepper (Capsicum) cultivation. Resistance to aphids was previously identified in Capsicum baccatum. This study aimed at elucidating the genetics of aphid resistance in C. baccatum. A QTL analysis was carried out for M. persicae resistance in an F2 population derived from an intraspecific cross between a highly resistant plant and a susceptible plant. Survival and reproduction were used as resistance parameters. Interval mapping detected two QTLs affecting aphid survival (Rmpas-1) and reproduction (Rmprp-1), respectively, both localized in the same area and sharing the same top marker on chromosome 2. Use of this marker as co-factor in multiple-QTL mapping analysis revealed a second, minor QTL (Rmprp-2) only affecting aphid reproduction, on chromosome 4. Fine mapping confirmed the effects of Rmpas-1 and Rmprp-1 and narrowed the major QTL Rmprp-1 down to a genomic region of 96 kb which is predicted to encode four analogues of resistance genes of the receptor-like kinase family containing a leucine-rich repeat domain (LRR-RLKs). This work provides not only initial information for breeding aphid-resistant pepper varieties, but also forms the basis for future molecular analysis of gene(s) involved in aphid resistance.

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

New genotyping technologies, offering the possibility of high genetic resolution at low cost, have helped fuel a surge in interest in the genetic analysis of polyploid species. Nevertheless, autopolyploid species present extra challenges not encountered in diploids and allopolyploids, such as polysomic inheritance or double reduction. Here we investigate the power and precision of quantitative trait locus (QTL) analysis in outcrossing autopolyploids, comparing the results of a model that assumes random bivalent chromosomal pairing during meiosis to one that also allows for multivalents and double reduction. Through a series of simulation studies we found that marginal gains in QTL detection power are achieved using the double reduction model when multivalent pairing occurs. However, when exploring the effect of variable genotypic information across parental homologs, we found that both QTL detection power and precision require high and uniform genotypic information contents. This effect far outweighed considerations regarding bivalent or multivalent pairing (and double reduction) during meiosis. We propose that autopolyploid QTL studies be accompanied by both marker coverage information and per-homolog genotypic information coefficients (GIC). Application of these methods to an autotetraploid potato (Solanum tuberosum L.) mapping population confirmed our ability to locate and dissect QTL in highly heterozygous outcrossing autotetraploid populations.

FitTetra 2.0 - improved genotype calling for tetraploids with multiple population and parental data support.

BACKGROUND: Genetic studies in tetraploids are lagging behind in comparison with studies of diploids as the complex genetics of tetraploids require much more elaborated computational methodologies. Recent advancements in development of molecular techniques and computational tools facilitate new methods for automated, high-throughput genotype calling in tetraploid species. We report on the upgrade of the widely-used fitTetra software aiming to improve its accuracy, which to date is hampered by technical artefacts in the data. RESULTS: Our upgrade of the fitTetra package is designed for a more accurate modelling of complex collections of samples. The package fits a mixture model where some parameters of the model are estimated separately for each sub-collection. When a full-sib family is analyzed, we use parental genotypes to predict the expected segregation in terms of allele dosages in the offspring. More accurate modelling and use of parental data increases the accuracy of dosage calling. We tested the package on data obtained with an Affymetrix Axiom 60 k array and compared its performance with the original version and the recently published ClusterCall tool, showing that at least 20% more SNPs could be called with our updated. CONCLUSION: Our updated software package shows clearly improved performance in genotype calling accuracy. Estimation of mixing proportions of the underlying dosage distributions is separated for full-sib families (where mixture proportions can be estimated from the parental dosages and inheritance model) and unstructured populations (where they are based on the assumption of Hardy-Weinberg equilibrium). Additionally, as the distributions of signal ratios of the dosage classes can be assumed to be the same for all populations, including parental data for some subpopulations helps to improve fitting other populations as well. The R package fitTetra 2.0 is freely available under the GNU Public License as Additional file with this article.

QTL mapping of insect resistance components of Solanum galapagense.

KEY MESSAGE: QTLs for insect resistance parameters, trichome type IV development, and more than 200 non-volatile metabolites, including 76 acyl sugars, all co-locate at the end of Chromosome 2 of Solanum galapagense. Host plant resistance is gaining importance as more and more insecticides are being banned due to environmental concerns. In tomato, resistance towards insects is found in wild relatives and has been attributed to the presence of glandular trichomes and their specific phytochemical composition. In this paper, we describe the results from a large-scale QTL mapping of data from whitefly resistance tests, trichome phenotyping and a comprehensive metabolomics analysis in a recombinant inbred line population derived from a cross between the cultivated Solanum lycopersicum and the wild relative S. galapagense, which is resistant to a range of pest insects. One major QTL (Wf-1) was found to govern the resistance against two different whitefly species. This QTL co-localizes with QTLs for the presence of trichomes type IV and V, as well as all 76 acyl sugars detected and about 150 other non-volatile phytochemicals, including methyl esters of the flavonols myricetin and quercetin. Based on these results, we hypothesize that Wf-1 is regulating the formation of glandular trichome type IV on the leaf epidermis, enabling the production and accumulation of bioactive metabolites in this type of trichomes.

Combining QTL mapping with transcriptome and metabolome profiling reveals a possible role for ABA signaling in resistance against the cabbage whitefly in cabbage.

Whiteflies are among the world's most significant agricultural pests and chemical insecticides are extensively used to reduce crop damage to acceptable levels. However, nearly all insecticides pose a threat to the environment and alternative control methods, such as breeding of crop varieties that are inherently insect-resistant, are needed. Previously, a strong source of plant-age dependent resistance to the cabbage whitefly (Aleyrodes proletella) has been identified in the modern white cabbage (Brassica oleracea var. capitata) variety Rivera. However, nothing is known about the molecular mechanisms or the genes involved in this resistance. In the present study, a multidisciplinary approach combining transcriptome and metabolome profiling with genetic mapping was used to identify the molecular players of whitefly resistance in cabbage. Transcriptome profiles of young (susceptible) and older (resistant) Rivera plants were analyzed using RNA sequencing. While many genes involved in general processes were differentially expressed between both ages, several defense-related processes were overrepresented in the transcriptome profile of older plants. Hormone measurements revealed that jasmonic acid (JA) levels decreased upon whitefly infestation at both plant ages. Interestingly, abscisic acid (ABA) levels showed contrasting effects in response to whitefly infestation: ABA levels were reduced in young plants but induced in older plants upon whitefly feeding. Auxin levels were significantly lower in older plants compared with young plants, independent of whitefly presence, while glucosinolate levels were higher. Additionally, whitefly performance was monitored in an F2 population derived from a cross between Rivera and the susceptible white cabbage variety Christmas Drumhead. Significant QTL intervals were mapped on chromosome 2 and 9 for oviposition rate and whitefly adult survival, respectively. Several genes that were higher expressed in older plants and located in the identified QTL intervals were orthologous to Arabidopsis genes that have been related to ABA signaling, suggesting a role for ABA in the regulation of resistance towards whiteflies. Our results show that combining different omics approaches is a useful strategy to identify candidate genes underlying insect resistance.

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.

Reduced phloem uptake of Myzus persicae on an aphid resistant pepper accession.

BACKGROUND: The green peach aphid (GPA), Myzus persicae, is economically one of the most threatening pests in pepper cultivation, which not only causes direct damage but also transmits many viruses. Breeding aphid resistant pepper varieties is a promising and environmentally friendly method to control aphid populations in the field and in the greenhouse. Until now, no strong sources of resistance against the GPA have been identified. Therefore the main aims of this study were to identify pepper materials with a good level of resistance to GPA and to elucidate possible resistance mechanisms. RESULTS: We screened 74 pepper accessions from different geographical areas for resistance to M. persicae. After four rounds of evaluation we identified one Capsicum baccatum accession (PB2013071) as highly resistant to M. persicae, while the accessions PB2013062 and PB2012022 showed intermediate resistance. The resistance of PB2013071 resulted in a severely reduced uptake of phloem compared to the susceptible accession, as determined by Electrical Penetration Graph (EPG) studies. Feeding of M. persicae induced the expression of callose synthase genes and resulted in callose deposition in the sieve elements in resistant, but not in susceptible plants. CONCLUSIONS: Three aphid resistant pepper accessions were identified, which will be important for breeding aphid resistant pepper varieties in the future. The most resistant accession PB2013071 showed phloem-based resistance against aphid infestation.