A next-generation marker genotyping platform (AmpSeq) in heterozygous crops: a case study for marker-assisted selection in grapevine

Integrating targeted genetic markers to genotyping-by-sequencing for an ultimate genotyping tool

New selection methods, using trait-specific markers (marker-assisted selection (MAS)) and/or genome-wide markers (genomic selection (GS)), are becoming increasingly widespread in breeding programs. This new era requires innovative and cost-efficient solutions for genotyping. Reduction in sequencing cost has enhanced the use of high-throughput low-cost genotyping methods such as genotyping-by-sequencing (GBS) for genome-wide single-nucleotide polymorphism (SNP) profiling in large breeding populations. However, the major weakness of GBS methodologies is their inability to genotype targeted markers. Conversely, targeted methods, such as amplicon sequencing (AmpSeq), often face cost constraints, hindering genome-wide genotyping across a large cohort. Although GBS and AmpSeq data can be generated from the same sample, an efficient method to achieve this is lacking. In this study, we present the Genome-wide & Targeted Amplicon (GTA) genotyping platform, an innovative way to integrate multiplex targeted amplicons into the GBS library preparation to provide an all-in-one cost-effective genotyping solution to breeders and research communities. Custom primers were designed to target 23 and 36 high-value markers associated with key agronomical traits in soybean and barley, respectively. The resulting multiplex amplicons were compatible with the GBS library preparation enabling both GBS and targeted genotyping data to be produced efficiently and cost-effectively. To facilitate data analysis, we have introduced Fast-GBS.v3, a user-friendly bioinformatic pipeline that generates comprehensive outputs from data obtained following sequencing of GTA libraries. This high-throughput low-cost approach will greatly facilitate the application of DNA markers as it provides required markers for both MAS and GS in a single assay.

Rpv10.2: A Haplotype Variant of Locus Rpv10 Enables New Combinations for Pyramiding Downy Mildew Resistance Traits in Grapevine

In viticulture, pathogens like the oomycete Plasmopara viticola, the causal agent of downy mildew, can cause severe yield loss and require extensive application of plant protection chemicals. Breeders are generating pathogen-resistant varieties exploiting American and Asian wild Vitis germplasm as sources of resistance. Several loci mediating resistance to P. viticola have been identified in the past but may be overcome by specifically adapted strains of the pathogen. Aiming to find and characterize novel loci, a cross population with Vitis amurensis ancestry was investigated searching for resistance-correlated quantitative trait loci (QTL). As a prerequisite, a genetic map was generated by analyzing the 244 F1 individuals derived from a cross of the downy mildew susceptible Vitis vinifera cultivar 'Tigvoasa' and the resistant V. amurensis pBC1 breeding line We 90-06-12. This genetic map is based on the information from 627 molecular markers including 56 simple sequence repeats and 571 rhAmpSeq markers. A phenotypic characterization of the progeny showed a clear segregation of the resistance traits in the F1 population after an experimental inoculation of leaf discs with downy mildew. Combining genetic and phenotypic data, an analysis for QTL revealed a major locus on linkage Group 9 that correlates strongly with the resistance to downy mildew. The locus was mapped to a region of about 80 kb on the PN40024 (12x.V2) grapevine reference genome. This genomic region co-localizes with the formerly identified locus Rpv10 from the grapevine cultivar 'Solaris'. As we found different allele sizes of the locus-linked SSR markers than those characterizing the known Rpv10 locus and differences in the sequence of a candidate gene, it was regarded as a haplotype variant and named Rpv10.2.

Multiple routes to fungicide resistance: Interaction of Cyp51 gene sequences, copy number and expression

We examined the molecular basis of triazole resistance in Blumeria graminis f. sp. tritici (wheat mildew, Bgt), a model organism among powdery mildews. Four genetic models for responses to triazole fungicides were identified among US and UK isolates, involving multiple genetic mechanisms. Firstly, only two amino acid substitutions in CYP51B lanosterol demethylase, the target of triazoles, were associated with resistance, Y136F and S509T (homologous to Y137F and S524T in the reference fungus Zymoseptoria tritici). As sequence variation did not explain the wide range of resistance, we also investigated Cyp51B copy number and expression, the latter using both reverse transcription-quantitative PCR and RNA-seq. The second model for resistance involved higher copy number and expression in isolates with a resistance allele; thirdly, however, moderate resistance was associated with higher copy number of wild-type Cyp51B in some US isolates. A fourth mechanism was heteroallelism with multiple alleles of Cyp51B. UK isolates, with significantly higher mean resistance than their US counterparts, had higher mean copy number, a high frequency of the S509T substitution, which was absent from the United States, and in the most resistant isolates, heteroallelism involving both sensitivity residues Y136+S509 and resistance residues F136+T509. Some US isolates were heteroallelic for Y136+S509 and F136+S509, but this was not associated with higher resistance. The obligate biotrophy of Bgt may constrain the tertiary structure and thus the sequence of CYP51B, so other variation that increases resistance may have a selective advantage. We describe a process by which heteroallelism may be adaptive when Bgt is intermittently exposed to triazoles.

Phenotyping, genetics, and "-omics" approaches to unravel and introgress enhanced resistance against apple scab (Venturia inaequalis) in apple cultivars (Malus × domestica)

Apple scab disease, caused by the fungus Venturia inaequalis, endangers commercial apple production globally. It is predominantly managed by frequent fungicide sprays that can harm the environment and promote the development of fungicide-resistant strains. Cultivation of scab-resistant cultivars harboring diverse qualitative Rvi resistance loci and quantitative trait loci associated with scab resistance could reduce the chemical footprint. A comprehensive understanding of the host-pathogen interaction is, however, needed to efficiently breed cultivars with enhanced resistance against a variety of pathogenic strains. Breeding efforts should not only encompass pyramiding of Rvi loci and their corresponding resistance alleles that directly or indirectly recognize pathogen effectors, but should also integrate genes that contribute to effective downstream defense mechanisms. This review provides an overview of the phenotypic and genetic aspects of apple scab resistance, and currently known corresponding defense mechanisms. Implementation of recent "-omics" approaches has provided insights into the complex network of physiological, molecular, and signaling processes that occur before and upon scab infection, thereby revealing the importance of both constitutive and induced defense mechanisms. Based on the current knowledge, we outline advances toward more efficient introgression of enhanced scab resistance into novel apple cultivars by conventional breeding or genetic modification techniques. However, additional studies integrating different "-omics" approaches combined with functional studies will be necessary to unravel effective defense mechanisms as well as key regulatory genes underpinning scab resistance in apple. This crucial information will set the stage for successful knowledge-based breeding for enhanced scab resistance.

Rapid and easy construction of a simplified amplicon sequencing (simplified AmpSeq) library for marker-assisted selection

Marker-assisted selection (MAS) is fundamental for plant breeding programs, as it can identify desirable seedlings at a young stage and reduce the cost, time and space needed for plant maintenance, especially for perennial crops. To facilitate the process of genotyping, which is time consuming and laborious, we developed a simplified amplicon sequencing (simplified AmpSeq) library construction method for next-generation sequencing that can be applied to MAS in breeding programs. The method is based on one-step PCR with a mixture of two primer sets: the first consisting of tailed target primers, the second of primers that contain flow-cell binding sites, indexes and tail sequences complementary to those in the first set. To demonstrate the process of MAS using s implified AmpSeq, we created databases of genotypes for important traits by using cultivar collections including triploid cultivars and segregating seedlings of Japanese pear (Pyrus pyrifolia Nakai), Japanese chestnut (Castanea crenata Sieb. et Zucc.) and apple (Malus domestica Borkh.). Simplified AmpSeq has the advantages of high repeatability, ability to estimate allele number in polyploid species and semi-automatic evaluation using target allele frequencies. Because this method provides high flexibility for designing primer sets and targeting any variant, it will be useful for plant breeding programs.

Fully sequencing the cassava full-length cDNA library reveals unannotated transcript structures and alternative splicing events in regions with a high density of single nucleotide variations, insertions-deletions, and heterozygous sequences

The primary transcript structure provides critical insights into protein diversity, transcriptional modification, and functions. Cassava transcript structures are highly diverse because of alternative splicing (AS) events and high heterozygosity. To precisely determine and characterize transcript structures, fully sequencing cloned transcripts is the most reliable method. However, cassava annotations were mainly determined according to fragmentation-based sequencing analyses (e.g., EST and short-read RNA-seq). In this study, we sequenced the cassava full-length cDNA library, which included rare transcripts. We obtained 8,628 non-redundant fully sequenced transcripts and detected 615 unannotated AS events and 421 unannotated loci. The different protein sequences resulting from the unannotated AS events tended to have diverse functional domains, implying that unannotated AS contributes to the truncation of functional domains. The unannotated loci tended to be derived from orphan genes, implying that the loci may be associated with cassava-specific traits. Unexpectedly, individual cassava transcripts were more likely to have multiple AS events than Arabidopsis transcripts, suggestive of the regulated interactions between cassava splicing-related complexes. We also observed that the unannotated loci and/or AS events were commonly in regions with abundant single nucleotide variations, insertions-deletions, and heterozygous sequences. These findings reflect the utility of completely sequenced FLcDNA clones for overcoming cassava-specific annotation-related problems to elucidate transcript structures. Our work provides researchers with transcript structural details that are useful for annotating highly diverse and unique transcripts and alternative splicing events.

Validation of KASP-SNP markers in cassava germplasm for marker-assisted selection of increased carotenoid content and dry matter content

Provitamin A biofortification and increased dry matter content are important breeding targets in cassava improvement programs worldwide. Biofortified varieties contribute to the alleviation of provitamin A deficiency, a leading cause of preventable blindness common among pre-school children and pregnant women in developing countries particularly Africa. Dry matter content is a major component of dry yield and thus underlies overall variety performance and acceptability by growers, processors, and consumers. Single nucleotide polymorphism (SNP) markers linked to these traits have recently been discovered through several genome-wide association studies but have not been deployed for routine marker-assisted selection (MAS). This is due to the lack of useful information on markers' performances in diverse genetic backgrounds. To overcome this bottleneck, technical and biological validation of the loci associated with increased carotenoid content and dry matter content were carried out using populations independent of the marker discovery population. In the present study, seven previously identified markers for these traits were converted to a robust set of uniplex allele-specific polymerase chain reaction (PCR) assays and validated in two independent pre-breeding and breeding populations. These assays were efficient in discriminating marker genotypic classes and had an average call rate greater than 98%. A high correlation was observed between the predicted and observed carotenoid content as inferred by root yellowness intensity in the breeding (r = 0.92) and pre-breeding (r = 0.95) populations. On the other hand, dry matter content-markers had moderately low predictive accuracy in both populations (r< 0.40) due to the more quantitative nature of the trait. This work confirmed the markers' effectiveness in multiple backgrounds, therefore, further strengthening their value in cassava biofortification to ensure nutritional security as well as dry matter content productivity. Our study provides a framework to guide future marker validation, thus leading to the more routine use of markers in MAS in cassava improvement programs.

Development and validation of an optimized marker set for genomic selection in southern U.S. rice breeding programs

The potential of genomic selection (GS) to increase the efficiency of breeding programs has been clearly demonstrated; however, the implementation of GS in rice (Oryza sativa L.) breeding programs has been limited. In recent years, efforts have begun to work toward implementing GS into the Louisiana State University (LSU) Agricultural Center rice breeding program. One of the first steps for successful GS implementation is to establish a suitable marker set for the target germplasm and a reliable, cost-effective genotyping platform capable of providing informative marker data with an adequate turnaround time. The objective of this study was to develop a marker set for routine GS and demonstrate its effectiveness in southern U.S. rice germplasm. The utility of the resulting marker set, the LSU500, for GS applications was demonstrated using four years of breeding data across 7,607 experimental lines and four elite biparental populations. The predictive ability of GS ranged from 0.13 to 0.78 for key traits across different market classes and yield trials. Comparisons between phenotypic selection and GS within biparental populations demonstrates similar performance of GS compared with phenotypic selection in predicting future performance. The prediction accuracies obtained with the LSU500 marker set demonstrates the utility of this marker set for cost-effective GS applications in southern U.S. rice breeding programs. The LSU500 marker set has been established through the genotyping service provider Agriplex Genomics, and in the future, it will undergo improvements to reduce the cost and increase the accuracy of GS.

Genetic Diversity of Wild and Cultivated Muscadine Grapes (Vitis rotundifolia Michx.)

The muscadine (Vitis rotundifolia syn. Muscadinia rotundifolia) is an American grape species native to the southeastern United States that has been cultivated for centuries. Muscadines are one of three grape species in subgenus Muscadinia with a chromosome number of 2n = 40 (V. rotundifolia, Vitis munsoniana, and Vitis popenoei), making them genetically distinct from the European wine and table grape (Vitis vinifera) and other species in subgenus Euvitis. Crop improvement efforts have been continuous since the late 19th century, yet the germplasm that served as the foundation for early muscadine breeding efforts was sourced from a relatively small portion of their native range, mostly in the coastal plains of North Carolina. This study used the rhAmpSeq Vitis core panel haplotype markers to genotype 194 Muscadinia accessions from five cultivated populations and 15 wild populations collected across their native range. Wild populations from the western half of the native range were generally less genetically differentiated than hypothesized, but were genetically distinct from the material used in both past and present breeding efforts. One population collected from coastal North Carolina grouped closely with V. munsoniana accessions despite being well outside the reported range for that species. Principal coordinate and structure analyses revealed three main groups within the 194 accessions: one for cultivated material, one for wild V. rotundifolia, and one for V. munsoniana and V. popenoei. At K = 5, structure results showed that more recent muscadine cultivars are further differentiated from wild accessions and varieties. These analyses confirmed our hypothesis that muscadine cultivars are genetically differentiated from their wild counterparts. This study also showed that genetic diversity in V. rotundifolia is not equally distributed across its native range and that the limited number of genotypes used in crop improvement efforts has not fully utilized the genetic diversity within the species.

Stable QTL for malate levels in ripe fruit and their transferability across Vitis species

Malate is a major contributor to the sourness of grape berries (Vitis spp.) and their products, such as wine. Excessive malate at maturity, commonly observed in wild Vitis grapes, is detrimental to grape and wine quality and complicates the introgression of valuable disease resistance and cold hardy genes through breeding. This study investigated an interspecific Vitis family that exhibited strong and stable variation in malate at ripeness for five years and tested the separate contribution of accumulation, degradation, and dilution to malate concentration in ripe fruit in the last year of study. Genotyping was performed using transferable rhAmpSeq haplotype markers, based on the Vitis collinear core genome. Three significant QTL for ripe fruit malate on chromosomes 1, 7, and 17, accounted for over two-fold and 6.9 g/L differences, and explained 40.6% of the phenotypic variation. QTL on chromosomes 7 and 17 were stable in all and in three out of five years, respectively. Variation in pre-veraison malate was the major contributor to variation in ripe fruit malate (39%), and based on two and five years of data, respectively, their associated QTL overlapped on chromosome 7, indicating a common genetic basis. However, use of transferable markers on a closely related Vitis family did not yield a common QTL across families. This suggests that diverse physiological mechanisms regulate the levels of this key metabolite in the Vitis genus, a conclusion supported by a review of over a dozen publications from the past decade, showing malate-associated genetic loci on all 19 chromosomes.

Mapping Major Disease Resistance Genes in Soybean by Genome-Wide Association Studies

Soybean is one of the most valuable agricultural crops in the world. Besides, this legume is constantly attacked by a wide range of pathogens (fungi, bacteria, viruses, and nematodes) compromising yield and increasing production costs. One of the major disease management strategies is the genetic resistance provided by single genes and quantitative trait loci (QTL). Identifying the genomic regions underlying the resistance against these pathogens on soybean is one of the first steps performed by molecular breeders. In the past, genetic mapping studies have been widely used to discover these genomic regions. However, over the last decade, advances in next-generation sequencing technologies and their subsequent cost decreasing led to the development of cost-effective approaches to high-throughput genotyping. Thus, genome-wide association studies applying thousands of SNPs in large sets composed of diverse soybean accessions have been successfully done. In this chapter, a comprehensive review of the majority of GWAS for soybean diseases published since this approach was developed is provided. Important diseases caused by Heterodera glycines, Phytophthora sojae, and Sclerotinia sclerotiorum have been the focus of the several GWAS. However, other bacterial and fungi diseases also have been targets of GWAS. As such, this GWAS summary can serve as a guide for future studies of these diseases. The protocol begins by describing several considerations about the pathogens and bringing different procedures of molecular characterization of them. Advice to choose the best isolate/race to maximize the discovery of multiple R genes or to directly map an effective R gene is provided. A summary of protocols, methods, and tools to phenotyping the soybean panel is given to several diseases. We also give details of options of DNA extraction protocols and genotyping methods, and we describe parameters of SNP quality to soybean data. Websites and their online tools to obtain genotypic and phenotypic data for thousands of soybean accessions are highlighted. Finally, we report several tricks and tips in Subheading 4, especially related to composing the soybean panel as well as generating and analyzing the phenotype data. We hope this protocol will be helpful to achieve GWAS success in identifying resistance genes on soybean.

The Genetic Basis of Anthocyanin Acylation in North American Grapes (Vitis spp.)

Hydroxycinnamylated anthocyanins (or simply 'acylated anthocyanins') increase color stability in grape products, such as wine. Several genes that are relevant for anthocyanin acylation in grapes have been previously described; however, control of the degree of acylation in grapes is complicated by the lack of genetic markers quantitatively associated with this trait. To characterize the genetic basis of anthocyanin acylation in grapevine, we analyzed the acylation ratio in two closely related biparental families, Vitis rupestris B38 × 'Horizon' and 'Horizon' × Illinois 547-1, for 2 and 3 years, respectively. The acylation ratio followed a bimodal and skewed distribution in both families, with repeatability estimates larger than 0.84. Quantitative trait locus (QTL) mapping with amplicon-based markers (rhAmpSeq) identified a strong QTL from 'Horizon' on chromosome 3, near 15.85 Mb in both families and across years, explaining up to 85.2% of the phenotypic variance. Multiple candidate genes were identified in the 14.85-17.95 Mb interval, in particular, three copies of a gene encoding an acetyl-CoA-benzylalcohol acetyltransferase-like protein within the two most strongly associated markers. Additional population-specific QTLs were found in chromosomes 9, 10, 15, and 16; however, no candidate genes were described. The rhAmpSeq markers reported here, which were previously shown to be highly transferable among the Vitis genus, could be immediately implemented in current grapevine breeding efforts to control the degree of anthocyanin acylation and improve the quality of grapes and their products.

Multiple independent recombinations led to hermaphroditism in grapevine

We studied the grape sex-determining region (SDR) in 12 Vitis genomes and demonstrated its conservation across 556 genotypes including 193 accessions from 47 world-wide wild grapevine species and 363 accessions of cultivated grapevine. Although the grape SDR is recombination free in all wild species, we found two distinct hermaphrodite (H) haplotypes (H1 and H2) among the cultivated grapevines, both chimeras of male (M) and female (f) haplotypes. The two independent recombinations carry different genetic signatures which long predate the domestication of grapevine, suggesting independent evolutions of this trait in wild European grapevine gene pools prior to human domestication.

Hermaphroditic (perfect) flowers were a key trait in grapevine domestication, enabling a drastic increase in yields due to the efficiency of self-pollination in the domesticated grapevine (Vitis vinifera L. ssp. vinifera). In contrast, all extant wild Vitis species are dioecious, each plant having only male or female flowers. In this study, we identified the male (M) and female (f) haplotypes of the sex-determining region (SDR) in the wild grapevine species V. cinerea and confirmed the boundaries of the SDR. We also demonstrated that the SDR and its boundaries are precisely conserved across the Vitis genus using shotgun resequencing data of 556 wild and domesticated accessions from North America, East Asia, and Europe. A high linkage disequilibrium was found at the SDR in all wild grape species, while different recombination signatures were observed along the hermaphrodite (H) haplotype of 363 cultivated accessions, revealing two distinct H haplotypes, named H1 and H2. To further examine the H2 haplotype, we sequenced the genome of two grapevine cultivars, 'Riesling' and 'Chardonnay'. By reconstructing the first two H2 haplotypes, we estimated the divergence time between H1 and H2 haplotypes at ∼6 million years ago, which predates the domestication of grapevine (∼8,000 y ago). Our findings emphasize the important role of recombination suppression in maintaining dioecy in wild grape species and lend additional support to the hypothesis that at least two independent recombination events led to the reversion to hermaphroditism in grapevine.

Whole genome resequencing and custom genotyping unveil clonal lineages in 'Malbec' grapevines (Vitis vinifera L.)

Grapevine cultivars are clonally propagated to preserve their varietal attributes. However, genetic variations accumulate due to the occurrence of somatic mutations. This process is anthropically influenced through plant transportation, clonal propagation and selection. Malbec is a cultivar that is well-appreciated for the elaboration of red wine. It originated in Southwestern France and was introduced in Argentina during the 1850s. In order to study the clonal genetic diversity of Malbec grapevines, we generated whole-genome resequencing data for four accessions with different clonal propagation records. A stringent variant calling procedure was established to identify reliable polymorphisms among the analyzed accessions. The latter procedure retrieved 941 single nucleotide variants (SNVs). A reduced set of the detected SNVs was corroborated through Sanger sequencing, and employed to custom-design a genotyping experiment. We successfully genotyped 214 Malbec accessions using 41 SNVs, and identified 14 genotypes that clustered in two genetically divergent clonal lineages. These lineages were associated with the time span of clonal propagation of the analyzed accessions in Argentina and Europe. Our results show the usefulness of this approach for the study of the scarce intra-cultivar genetic diversity in grapevines. We also provide evidence on how human actions might have driven the accumulation of different somatic mutations, ultimately shaping the Malbec genetic diversity pattern.

Mining Grapevine Downy Mildew Susceptibility Genes: A Resource for Genomics-Based Breeding and Tailored Gene Editing

Several pathogens continuously threaten viticulture worldwide. Until now, the investigation on resistance loci has been the main trend to understand the interaction between grapevine and the mildew causal agents. Dominantly inherited gene-based resistance has shown to be race-specific in some cases, to confer partial immunity, and to be potentially overcome within a few years since its introgression. Recently, on the footprint of research conducted in Arabidopsis, putative genes associated with downy mildew susceptibility have been discovered also in the grapevine genome. In this work, we deep-sequenced four putative susceptibility genes-namely VvDMR6.1, VvDMR6.2, VvDLO1, VvDLO2-in 190 genetically diverse grapevine genotypes to discover new sources of broad-spectrum and recessively inherited resistance. Identified Single Nucleotide Polymorphisms were screened in a bottleneck analysis from the genetic sequence to their impact on protein structure. Fifty-five genotypes showed at least one impacting mutation in one or more of the scouted genes. Haplotypes were inferred for each gene and two of them at the VvDMR6.2 gene were found significantly more represented in downy mildew resistant genotypes. The current results provide a resource for grapevine and plant genetics and could corroborate genomic-assisted breeding programs as well as tailored gene editing approaches for resistance to biotic stresses.

Population structure of Erysiphe necator on domesticated and wild vines in the Middle East raises questions on the origin of the grapevine powdery mildew pathogen

Plant pathogens usually originate and diversify in geographical regions where hosts and pathogens co-evolve. Erysiphe necator, the causal agent of grape powdery mildew, is a destructive pathogen of grapevines worldwide. Although Eastern US is considered the centre of origin and diversity of E. necator, previous reports on resistant native wild and domesticated Asian grapevines suggest Asia as another possible origin of the pathogen. By using multi-locus sequencing, microsatellites and a novel application of amplicon sequencing (AmpSeq), we show that the population of E. necator in Israel is composed of three genetic groups: Groups A and B that are common worldwide, and a new group IL, which is genetically differentiated from any known group in Europe and Eastern US. Group IL showed distinguished ecological characteristics: it was dominant on wild and traditional vines (95%); its abundance increased along the season; and was more aggressive than A and B isolates on both wild and domesticated vines. The low genetic diversity within group IL suggests that it has invaded Israel from another origin. Therefore, we suggest that the Israeli E. necator population was founded by at least two invasions, of which one could be from a non-East American source, possibly from Asian origin.

Discovery of the REN11 Locus From Vitis aestivalis for Stable Resistance to Grapevine Powdery Mildew in a Family Segregating for Several Unstable and Tissue-Specific Quantitative Resistance Loci

Race-specific resistance loci, whether having qualitative or quantitative effects, present plant-breeding challenges for phenotypic selection and deciding which loci to select or stack with other resistance loci for improved durability. Previously, resistance to grapevine powdery mildew (GPM, caused by Erysiphe necator) was predicted to be conferred by at least three race-specific loci in the mapping family B37-28 × C56-11 segregating for GPM resistance from Vitis aestivalis. In this study, 9 years of vineyard GPM disease severity ratings plus a greenhouse and laboratory assays were genetically mapped, using a rhAmpSeq core genome marker platform with 2,000 local haplotype markers. A new qualitative resistance locus, named REN11, on the chromosome (Chr) 15 was found to be effective in nearly all (11 of 12) vineyard environments on leaves, rachis, berries, and most of the time (7 of 12) stems. REN11 was independently validated in a pseudo-testcross with the grandparent source of resistance, "Tamiami." Five other loci significantly predicted GPM severity on leaves in only one or two environments, which could indicate race-specific resistance or their roles in different timepoints in epidemic progress. Loci on Chr 8 and 9 reproducibly predicted disease severity on stems but not on other tissues and had additive effects with REN11 on the stems. The rhAmpSeq local haplotype sequences published in this study for REN11 and Chr 8 and 9 stem quantitative trait locus (QTL) can be used directly for marker-assisted selection or converted to SNP assays. In screening for REN11 in a diversity panel of 20,651 vines representing the diversity of Vitis, this rhAmpSeq haplotype had a false positive rate of 0.034% or less. The effects of the other foliar resistance loci detected in this study seem too unstable for genetic improvement regardless of quantitative effect size, whether due to race specificity or other environmental variables.

Transcriptome-Derived Amplicon Sequencing Markers Elucidate the U.S. Podosphaera macularis Population Structure Across Feral and Commercial Plantings of Humulus lupulus

Obligately biotrophic plant pathogens pose challenges in population genetic studies due to their genomic complexities and elaborate culturing requirements with limited biomass. Hop powdery mildew (Podosphaera macularis) is an obligately biotrophic ascomycete that threatens sustainable hop production. P. macularis populations of the Pacific Northwest (PNW) United States differ from those of the Midwest and Northeastern United States, lacking one of two mating types needed for sexual recombination and harboring two strains that are differentially aggressive on the cultivar Cascade and able to overcome the Humulus lupulus R-gene R6 (V6), respectively. To develop a high-throughput marker platform for tracking the flow of genotypes across the United States and internationally, we used an existing transcriptome of diverse P. macularis isolates to design a multiplex of 54 amplicon sequencing markers, validated across a panel of 391 U.S. samples and 123 international samples. The results suggest that P. macularis from U.S. commercial hop yards form one population closely related to P. macularis of the United Kingdom, while P. macularis from U.S. feral hop locations grouped with P. macularis of Eastern Europe. Included in this multiplex was a marker that successfully tracked V6-virulence in 65 of 66 samples with a confirmed V6-phenotype. A new qPCR assay for high-throughput genotyping of P. macularis mating type generated the highest resolution distribution map of P. macularis mating type to date. Together, these genotyping strategies enable the high-throughput and inexpensive tracking of pathogen spread among geographical regions from single-colony samples and provide a roadmap to develop markers for other obligate biotrophs.

Candidate gene mapping identifies genomic variations in the fire blight susceptibility genes HIPM and DIPM across the Malus germplasm

Development of apple (Malus domestica) cultivars resistant to fire blight, a devastating bacterial disease caused by Erwinia amylovora, is a priority for apple breeding programs. Towards this goal, the inactivation of members of the HIPM and DIPM gene families with a role in fire blight susceptibility (S genes) can help achieve sustainable tolerance. We have investigated the genomic diversity of HIPM and DIPM genes in Malus germplasm collections and used a candidate gene-based association mapping approach to identify SNPs (single nucleotide polymorphisms) with significant associations to fire blight susceptibility. A total of 87 unique SNP variants were identified in HIPM and DIPM genes across 93 Malus accessions. Thirty SNPs showed significant associations (p < 0.05) with fire blight susceptibility traits, while two of these SNPs showed highly significant (p < 0.001) associations across two different years. This research has provided knowledge about genetic diversity in fire blight S genes in diverse apple accessions and identified candidate HIPM and DIPM alleles that could be used to develop apple cultivars with decreased fire blight susceptibility via marker-assisted breeding or biotechnological approaches.

Novel stable QTLs identification for berry quality traits based on high-density genetic linkage map construction in table grape

BACKGROUND

Aroma, berry firmness and berry shape are three main quality traits in table grape production, and also the important target traits in grapevine breeding. However, the information about their genetic mechanisms is limited, which results in low accuracy and efficiency of quality breeding in grapevine. Mapping and isolation of quantitative trait locus (QTLs) based on the construction of genetic linkage map is a powerful approach to decipher the genetic determinants of complex quantitative traits.

RESULTS

In the present work, a final integrated map consisting of 3411 SLAF markers on 19 linkage groups (LGs) with an average distance of 0.98 cM between adjacent markers was generated using the specific length amplified fragment sequencing (SLAF-seq) technique. A total of 9 significant stable QTLs for Muscat flavor, berry firmness and berry shape were identified on two linkage groups among the hybrids analyzed over three consecutive years from 2016 to 2018. Notably, new stable QTLs for berry firmness and berry shape were found on LG 8 respectively for the first time. Based on biological function and expression profiles of candidate genes in the major QTL regions, 3 genes (VIT_08s0007g00440, VIT_08s0040g02740 and VIT_08s0040g02350) related to berry firmness and 3 genes (VIT_08s0032g01110, VIT_08s0032g01150 and VIT_08s0105g00200) linked to berry shape were highlighted. Overexpression of VIT_08s0032g01110 in transgenic Arabidopsis plants caused the change of pod shape.

CONCLUSIONS

A new high-density genetic map with total 3411 markers was constructed with SLAF-seq technique, and thus enabled the detection of narrow interval QTLs for relevant traits in grapevine. VIT_08s0007g00440, VIT_08s0040g02740 and VIT_08s0040g02350 were found to be related to berry firmness, while VIT_08s0032g01110, VIT_08s0032g01150 and VIT_08s0105g00200 were linked to berry shape.

Biotechnological Approaches: Gene Overexpression, Gene Silencing, and Genome Editing to Control Fungal and Oomycete Diseases in Grapevine

Downy mildew, powdery mildew, and grey mold are some of the phytopathological diseases causing economic losses in agricultural crops, including grapevine, worldwide. In the current scenario of increasing global warming, in which the massive use of agrochemicals should be limited, the management of fungal disease has become a challenge. The knowledge acquired on candidate resistant (R) genes having an active role in plant defense mechanisms has allowed numerous breeding programs to integrate these traits into selected cultivars, even though with some limits in the conservation of the proper qualitative characteristics of the original clones. Given their gene-specific mode of action, biotechnological techniques come to the aid of breeders, allowing them to generate simple and fast modifications in the host, without introducing other undesired genes. The availability of efficient gene transfer procedures in grapevine genotypes provide valid tools that support the application of new breeding techniques (NBTs). The expertise built up over the years has allowed the optimization of these techniques to overexpress genes that directly or indirectly limit fungal and oomycetes pathogens growth or silence plant susceptibility genes. Furthermore, the downregulation of pathogen genes which act as virulence effectors by exploiting the RNA interference mechanism, represents another biotechnological tool that increases plant defense. In this review, we summarize the most recent biotechnological strategies optimized and applied on Vitis species, aimed at reducing their susceptibility to the most harmful fungal and oomycetes diseases. The best strategy for combating pathogenic organisms is to exploit a holistic approach that fully integrates all these available tools.

A new SNP genotyping technology Target SNP-seq and its application in genetic analysis of cucumber varieties

To facilitate the utility of SNP-based genotyping, we developed a new method called target SNP-seq which combines the advantages of multiplex PCR amplification and high throughput sequencing. Compared with KASP, Microarrays, GBS and other SNP genotyping methods, target SNP-seq is flexible both in SNPs and samples, yields high accuracy, especially when genotyping genome wide perfect SNPs with high polymorphism and conserved flanking sequences, and is cost-effective, requiring 3 days and $7 for per DNA sample to genotype hundreds of SNP loci. The present study established a DNA fingerprint of 261 cucumber varieties by target SNP-seq with 163 perfect SNPs from 4,612,350 SNPs based on 182 cucumber resequencing datasets. Four distinct subpopulations were found in 261 Chinese cucumber varieties: the north China type, the south China type, the Europe type, and the Xishuangbanna type. The north China type and Xishuangbanna type harbored lower genetic diversity, indicating greater risk of genetic erosion in these two subpopulations. Furthermore, a core set of 24 SNPs was able to distinguish 99% of the 261 cucumber varieties. 29 core cucumber backbone varieties in China were identified. Therefore, target SNP-seq provides a new way to screen out core SNP loci from the whole genome for DNA fingerprinting of crop varieties. The high efficiency and low cost of target SNP-seq is more competitive than the current SNP genotyping methods, and it has excellent application prospects in genetic research, as well as in promoting plant breeding processes in the near future.

Haplotyping the Vitis collinear core genome with rhAmpSeq improves marker transferability in a diverse genus

Transferable DNA markers are essential for breeding and genetics. Grapevine (Vitis) breeders utilize disease resistance alleles from congeneric species ~20 million years divergent, but existing Vitis marker platforms have cross-species transfer rates as low as 2%. Here, we apply a marker strategy targeting the inferred Vitis core genome. Incorporating seven linked-read de novo assemblies and three existing assemblies, the Vitis collinear core genome is estimated to converge at 39.8 Mb (8.67% of the genome). Adding shotgun genome sequences from 40 accessions enables identification of conserved core PCR primer binding sites flanking polymorphic haplotypes with high information content. From these target regions, we develop 2,000 rhAmpSeq markers as a PCR multiplex and validate the panel in four biparental populations spanning the diversity of the Vitis genus, showing transferability increases to 91.9%. This marker development strategy should be widely applicable for genetic studies in many taxa, particularly those ~20 million years divergent.

Virus-Induced Flowering by Apple Latent Spherical Virus Vector: Effective Use to Accelerate Breeding of Grapevine

Apple latent spherical virus (ALSV) was successfully used in promoting flowering (virus-induced flowering, VIF) in apple and pear seedlings. In this paper, we report the use of ALSV vectors for VIF in seedlings and in vitro cultures of grapevine. After adjusting experimental conditions for biolistic inoculation of virus RNA, ALSV efficiently infected not only progeny seedlings of Vitis spp. ‘Koshu,’ but also in vitro cultures of V. vinifera ‘Neo Muscat’ without inducing viral symptoms. The grapevine seedlings and in vitro cultures inoculated with an ALSV vector expressing the ‘florigen’ gene (Arabidopsis Flowering locus T, AtFT) started to set floral buds 20–30 days after inoculation. This VIF technology was successfully used to promote flowering and produce grapes with viable seeds in in vitro cultures of F₁ hybrids from crosses between V. ficifolia and V. vinifera and made it possible to analyze the quality of fruits within a year after germination. High-temperature (37 °C) treatment of ALSV-infected grapevine disabled virus movement to newly growing tissue to obtain ALSV-free shoots. Thus, the VIF using ALSV vectors can be used to shorten the generation time of grapevine seedlings and accelerate breeding of grapevines with desired traits.

R-Loci Arrangement Versus Downy and Powdery Mildew Resistance Level: A Vitis Hybrid Survey

For the viticulture of the future, it will be an essential prerequisite to manage grapevine diseases with fewer chemical inputs. The development and the deployment of novel mildew resistant varieties are considered one of the most promising strategies towards a sustainable viticulture. In this regard, a collection of 102 accessions derived from crossing Vitis hybrids with V. vinifera varieties was studied. In addition to the true-to-type analysis, an exhaustive genetic characterization was carried out at the 11 reliable mildew resistance (R) loci available in the literature to date. Our findings highlight the pyramiding of R-loci against downy mildew in 15.7% and against powdery mildew in 39.2% of the total accessions. The genetic analysis was coupled with a three-year evaluation of disease symptoms in an untreated field in order to assess the impact of the R-loci arrangement on the disease resistance degree at leaf and bunch level. Overall, our results strongly suggest that R-loci pyramiding does not necessarily mean to increase the overall disease resistance, but it guarantees the presence of further barriers in case of pathogens overcoming the first. Moreover, our survey allows the discovery of new mildew resistance sources useful for novel QTL identifications towards marker-assisted breeding.

The epiphytic microbiota of sour rot-affected grapes differs minimally from that of healthy grapes, indicating causal organisms are already present on healthy berries

Sour rot is a disease complex produced by an interaction between grape berries and various species of yeast and acetic acid bacteria in the presence of Drosophila fruit flies. While yeast and bacteria are consistently found on healthy grape berries worldwide, we explored whether the composition of these epiphytic communities differed depending on the presence or absence of sour rot symptoms. Using high-throughput sequencing, we characterized the microbiome of sour rot-affected grapes from two geographical areas across two years. In 2015 and 2016, both healthy and sour rot-affected berries were collected from commercial and research vineyards in Geneva, NY and commercial vineyards in Tasmania, AUS. In this experiment, all associated organisms grouped together primarily by location, and not by presence/absence of symptoms or cultivar. The predominant difference between asymptomatic and symptomatic samples, regardless of location, was the abundance of Acetobacter species, which were significantly more plentiful in the symptomatic samples. Yeast genera such as Candida, Hanseniaspora, Pichia and Saccharomyces were abundant in both sets of samples, but varied by region. The consistent presence of yeast species and the increased abundance of acetic acid-generating bacteria is consistent with our understanding of their etiological role in sour rot development. In 2016, diseased grapes also were collected from vineyards in Fredonia, NY, and Modesto, CA. Consistent with our comparison study, all associated organisms again grouped together primarily by location. Yeast genera such as Candida, Hanseniaspora, Pichia and Saccharomyces were abundant in both sets of samples, but varied by region. The consistent presence of yeast species and the abundance of acetic acid-generating bacteria in both experiments is consistent with our understanding of their etiological role in sour rot development.

Accelerating crop genetic gains with genomic selection

Genomic prediction based on additive genetic effects can accelerate genetic gain. There are opportunities for further improvement by including non-additive effects that access untapped sources of genetic diversity. Several studies have reported a worrying gap between the projected global future demand for plant-based products and the current annual rates of production increase, indicating that enhancing the rate of genetic gain might be critical for future food security. Therefore, new breeding technologies and strategies are required to significantly boost genetic improvement of future crop cultivars. Genomic selection (GS) has delivered considerable genetic gain in animal breeding and is becoming an essential component of many modern plant breeding programmes as well. In this paper, we review the lessons learned from implementing GS in livestock and the impact of GS on crop breeding, and discuss important features for the success of GS under different breeding scenarios. We highlight major challenges associated with GS including rapid genotyping, phenotyping, genotype-by-environment interaction and non-additivity and give examples for opportunities to overcome these issues. Finally, the potential of combining GS with other modern technologies in order to maximise the rate of crop genetic improvement is discussed, including the potential of increasing prediction accuracy by integration of crop growth models in GS frameworks.

Computational Analysis of AmpSeq Data for Targeted, High-Throughput Genotyping of Amplicons

Amplicon sequencing (AmpSeq) is a practical, intuitive strategy with a semi-automated computational pipeline for analysis of highly multiplexed PCR-derived sequences. This genotyping platform is particularly cost-effective when multiplexing 96 or more samples with a few amplicons up to thousands of amplicons. Amplicons can target from a single nucleotide to the upper limit of the sequencing platform. The flexibility of AmpSeq's wet lab methods make it a tool of broad interest for diverse species, and AmpSeq excels in flexibility, high-throughput, low-cost, accuracy, and semi-automated analysis. Here we provide an open science framework procedure to output data out of an AmpSeq project, with an emphasis on the bioinformatics pipeline to generate SNPs, haplotypes and presence/absence variants in a set of diverse genotypes. Open-access tutorial datasets with actual data and a containerization open source software instance are provided to enable training in each of these genotyping applications. The pipelines presented here should be applicable to the analysis of various target-enriched (e.g., amplicon or sequence capture) Illumina sequence data.

Target SSR-Seq: A Novel SSR Genotyping Technology Associate With Perfect SSRs in Genetic Analysis of Cucumber Varieties

Simple sequence repeats (SSR) - also known as microsatellites - have been used extensively in genetic analysis, fine mapping, quantitative trait locus (QTL) mapping, as well as marker-assisted selection (MAS) breeding and other techniques. Despite a plethora of studies reporting that perfect SSRs with stable motifs and flanking sequences are more efficient for genetic research, the lack of a high throughput technology for SSR genotyping has limited their use as genetic targets in many crops. In this study, we developed a technology called Target SSR-seq that combined the multiplexed amplification of perfect SSRs with high throughput sequencing. This method can genotype plenty of SSR loci in hundreds of samples with highly accurate results, due to the substantial coverage afforded by high throughput sequencing. We also detected 844 perfect SSRs based on 182 resequencing datasets in cucumber, of which 91 SSRs were selected for Target SSR-seq. Finally, 122 SSRs, including 31 SSRs for varieties identification, were used to genotype 382 key cucumber varieties readily available in Chinese markets using our Target SSR-seq method. Libraries of PCR products were constructed and then sequenced on the Illumina HiSeq X Ten platform. Bioinformatics analysis revealed that 111 filtered SSRs were accurately genotyped with an average coverage of 1289× at an extremely low cost; furthermore, 398 alleles were observed in 382 cucumber cultivars. Genetic analysis identified four populations: northern China type, southern China type, European type, and Xishuangbanna type. Moreover, we acquired a set of 16 core SSRs for the identification of 382 cucumber varieties, of which 42 were isolated as backbone cucumber varieties. This study demonstrated that Target SSR-seq is a novel and efficient method for genetic research.

Genetic Mapping of the Incompatibility Locus in Olive and Development of a Linked Sequence-Tagged Site Marker

The genetic control of self-incompatibility (SI) has been recently disclosed in olive. Inter-varietal crossing confirmed the presence of only two incompatibility groups (G1 and G2), suggesting a simple Mendelian inheritance of the trait. A double digest restriction associated DNA (ddRAD) sequencing of a biparental population segregating for incompatibility groups has been performed and high-density linkage maps were constructed in order to map the SI locus and identify gene candidates and linked markers. The progeny consisted of a full-sib family of 229 individuals derived from the cross 'Leccino' (G1) × 'Dolce Agogia' (G2) varieties, segregating 1:1 (G1:G2), in accordance with a diallelic self-incompatibility (DSI) model. A total of 16,743 single nucleotide polymorphisms was identified, 7,006 in the female parent 'Leccino' and 9,737 in the male parent 'Dolce Agogia.' Each parental map consisted of 23 linkage groups and showed an unusual large size (5,680 cM in 'Leccino' and 3,538 cM in 'Dolce Agogia'). Recombination was decreased across all linkage groups in pollen mother cells of 'Dolce Agogia,' the parent with higher heterozygosity, compared to megaspore mother cells of 'Leccino,' in a context of a species that showed exceptionally high recombination rates. A subset of 109 adult plants was assigned to either incompatibility group by a stigma test and the diallelic self-incompatibility (DSI) locus was mapped to an interval of 5.4 cM on linkage group 18. This region spanned a size of approximately 300 Kb in the olive genome assembly. We developed a sequence-tagged site marker in the DSI locus and identified five haplotypes in 57 cultivars with known incompatibility group assignment. A combination of two single-nucleotide polymorphisms (SNPs) was sufficient to predict G1 or G2 phenotypes in olive cultivars, enabling early marker-assisted selection of compatible genotypes and allowing for a rapid screening of inter-compatibility among cultivars in order to guarantee effective fertilization and increase olive production. The construction of high-density linkage maps has led to the development of the first functional marker in olive and provided positional candidate genes in the SI locus.

An aromatic amino acid and associated helix in the C-terminus of the potato leafroll virus minor capsid protein regulate systemic infection and symptom expression

The C-terminal region of the minor structural protein of potato leafroll virus (PLRV), known as the readthrough protein (RTP), is involved in efficient virus movement, tissue tropism and symptom development. Analysis of numerous C-terminal deletions identified a five-amino acid motif that is required for RTP function. A PLRV mutant expressing RTP with these five amino acids deleted (Δ5aa-RTP) was compromised in systemic infection and symptom expression. Although the Δ5aa-RTP mutant was able to move long distance, limited infection foci were observed in systemically infected leaves suggesting that these five amino acids regulate virus phloem loading in the inoculated leaves and/or unloading into the systemically infected tissues. The 5aa deletion did not alter the efficiency of RTP translation, nor impair RTP self-interaction or its interaction with P17, the virus movement protein. However, the deletion did alter the subcellular localization of RTP. When co-expressed with a PLRV infectious clone, a GFP tagged wild-type RTP was localized to discontinuous punctate spots along the cell periphery and was associated with plasmodesmata, although localization was dependent upon the developmental stage of the plant tissue. In contrast, the Δ5aa-RTP-GFP aggregated in the cytoplasm. Structural modeling indicated that the 5aa deletion would be expected to perturb an α-helix motif. Two of 30 plants infected with Δ5aa-RTP developed a wild-type virus infection phenotype ten weeks post-inoculation. Analysis of the virus population in these plants by deep sequencing identified a duplication of sequences adjacent to the deletion that were predicted to restore the α-helix motif. The subcellular distribution of the RTP is regulated by the 5-aa motif which is under strong selection pressure and in turn contributes to the efficient long distance movement of the virus and the induction of systemic symptoms.

HS-SPME-GC-MS Analyses of Volatiles in Plant Populations-Quantitating Compound × Individual Matrix Effects

Headspace solid-phase microextraction (HS-SPME) coupled to gas chromatography⁻mass spectrometry (GC-MS) is widely employed for volatile analyses of plants, including mapping populations used in plant breeding research. Studies often employ a single internal surrogate standard, even when multiple analytes are measured, with the assumption that any relative changes in matrix effects among individuals would be similar for all compounds, i.e., matrix effects do not show Compound × Individual interactions. We tested this assumption using individuals from two plant populations: an interspecific grape (Vitis spp.) mapping population (n = 140) and a tomato (Solanum spp.) recombinant inbred line (RIL) population (n = 148). Individual plants from the two populations were spiked with a cocktail of internal standards (n = 6, 9, respectively) prior to HS-SPME-GC-MS. Variation in the relative responses of internal standards indicated that Compound × Individual interactions exist but were different between the two populations. For the grape population, relative responses among pairs of internal standards varied considerably among individuals, with a maximum of 249% relative standard deviation (RSD) for the pair of [U13C]hexanal and [U13C]hexanol. However, in the tomato population, relative responses of internal standard pairs varied much less, with pairwise RSDs ranging from 8% to 56%. The approach described in this paper could be used to evaluate the suitability of using surrogate standards for HS-SPME-GC-MS studies in other plant populations.

Two dominant loci determine resistance to Phomopsis cane lesions in F1 families of hybrid grapevines

KEY MESSAGE

Rapid characterization of novel NB-LRR-associated resistance to Phomopsis cane spot on grapevine using high-throughput sampling and low-coverage sequencing for genotyping, locus mapping and transcriptome analysis provides insights into genetic resistance to a hemibiotrophic fungus. Phomopsis cane and leaf spot, caused by the hemibiotrophic fungus Diaporthe ampelina (syn = Phomopsis viticola), reduces the productivity in grapevines. Host resistance was studied on three F1 families derived from crosses involving resistant genotypes 'Horizon', Illinois 547-1, Vitis cinerea B9 and V. vinifera 'Chardonnay'. All families had progeny with extremely susceptible phenotypes, developing lesions on both dormant canes and maturing fruit clusters. Segregation of symptoms was observed under natural levels of inoculum in the field, while phenotypes on green shoots were confirmed under controlled inoculations in greenhouse. High-density genetic maps were used to localize novel qualitative resistance loci named Rda1 and Rda2 from V. cinerea B9 and 'Horizon', respectively. Co-linearity between reference genetic and physical maps allowed localization of Rda2 locus between 1.5 and 2.4 Mbp on chromosome 7, and Rda1 locus between 19.3 and 19.6 Mbp of chromosome 15, which spans a cluster of five NB-LRR genes. Further dissection of this locus was obtained by QTL mapping of gene expression values 14 h after inoculation across a subset of the 'Chardonnay' × V. cinerea B9 progeny. This provided evidence for the association between transcript levels of two of these NB-LRR genes with Rda1, with increased NB-LRR expression among susceptible progeny. In resistant parent V. cinerea B9, inoculation with D. ampelina was characterized by up-regulation of SA-associated genes and down-regulation of ethylene pathways, suggesting an R-gene-mediated response. With dominant effects associated with disease-free berries and minimal symptoms on canes, Rda1 and Rda2 are promising loci for grapevine genetic improvement.

Cacao biotechnology: current status and future prospects

Theobroma cacao-The Food of the Gods, provides the raw material for the multibillion dollar chocolate industry and is also the main source of income for about 6 million smallholders around the world. Additionally, cocoa beans have a number of other nonfood uses in the pharmaceutical and cosmetic industries. Specifically, the potential health benefits of cocoa have received increasing attention as it is rich in polyphenols, particularly flavonoids. At present, the demand for cocoa and cocoa-based products in Asia is growing particularly rapidly and chocolate manufacturers are increasing investment in this region. However, in many Asian countries, cocoa production is hampered due to many reasons including technological, political and socio-economic issues. This review provides an overview of the present status of global cocoa production and recent advances in biotechnological applications for cacao improvement, with special emphasis on genetics/genomics, in vitro embryogenesis and genetic transformation. In addition, in order to obtain an insight into the latest innovations in the commercial sector, a survey was conducted on granted patents relating to T. cacao biotechnology.

NGS-Based Genotyping, High-Throughput Phenotyping and Genome-Wide Association Studies Laid the Foundations for Next-Generation Breeding in Horticultural Crops

Demographic trends and changes to climate require a more efficient use of plant genetic resources in breeding programs. Indeed, the release of high-yielding varieties has resulted in crop genetic erosion and loss of diversity. This has produced an increased susceptibility to severe stresses and a reduction of several food quality parameters. Next generation sequencing (NGS) technologies are being increasingly used to explore “gene space” and to provide high-resolution profiling of nucleotide variation within germplasm collections. On the other hand, advances in high-throughput phenotyping are bridging the genotype-to-phenotype gap in crop selection. The combination of allelic and phenotypic data points via genome-wide association studies is facilitating the discovery of genetic loci that are associated with key agronomic traits. In this review, we provide a brief overview on the latest NGS-based and phenotyping technologies and on their role to unlocking the genetic potential of vegetable crops; then, we discuss the paradigm shift that is underway in horticultural crop breeding.

Enhancing genetic gain in the era of molecular breeding

As one of the important concepts in conventional quantitative genetics and breeding, genetic gain can be defined as the amount of increase in performance that is achieved annually through artificial selection. To develop pro ducts that meet the increasing demand of mankind, especially for food and feed, in addition to various industrial uses, breeders are challenged to enhance the potential of genetic gain continuously, at ever higher rates, while they close the gaps that remain between the yield potential in breeders' demonstration trials and the actual yield in farmers' fields. Factors affecting genetic gain include genetic variation available in breeding materials, heritability for traits of interest, selection intensity, and the time required to complete a breeding cycle. Genetic gain can be improved through enhancing the potential and closing the gaps, which has been evolving and complemented with modern breeding techniques and platforms, mainly driven by molecular and genomic tools, combined with improved agronomic practice. Several key strategies are reviewed in this article. Favorable genetic variation can be unlocked and created through molecular and genomic approaches including mutation, gene mapping and discovery, and transgene and genome editing. Estimation of heritability can be improved by refining field experiments through well-controlled and precisely assayed environmental factors or envirotyping, particularly for understanding and controlling spatial heterogeneity at the field level. Selection intensity can be significantly heightened through improvements in the scale and precision of genotyping and phenotyping. The breeding cycle time can be shortened by accelerating breeding procedures through integrated breeding approaches such as marker-assisted selection and doubled haploid development. All the strategies can be integrated with other widely used conventional approaches in breeding programs to enhance genetic gain. More transdisciplinary approaches, team breeding, will be required to address the challenge of maintaining a plentiful and safe food supply for future generations. New opportunities for enhancing genetic gain, a high efficiency breeding pipeline, and broad-sense genetic gain are also discussed prospectively.

Genotyping by Sequencing for SNP-Based Linkage Analysis and Identification of QTLs Linked to Fruit Quality Traits in Japanese Plum (Prunus salicina Lindl.)

Marker-assisted selection (MAS) in stone fruit (Prunus species) breeding is currently difficult to achieve due to the polygenic nature of the most relevant agronomic traits linked to fruit quality. Genotyping by sequencing (GBS), however, provides a large quantity of useful data suitable for fine mapping using Single Nucleotide Polymorphisms (SNPs) from a reference genome. In this study, GBS was used to genotype 272 seedlings of three F1 Japanese plum (Prunus salicina Lindl) progenies derived from crossing "98-99" (as a common female parent) with "Angeleno," "September King," and "September Queen" as male parents. Raw sequences were aligned to the Peach genome v1, and 42,909 filtered SNPs were obtained after sequence alignment. In addition, 153 seedlings from the "98-99" × "Angeleno" cross were used to develop a genetic map for each parent. A total of 981 SNPs were mapped (479 for "98-99" and 502 for "Angeleno"), covering a genetic distance of 688.8 and 647.03 cM, respectively. Fifty five seedlings from this progeny were phenotyped for different fruit quality traits including ripening time, fruit weight, fruit shape, chlorophyll index, skin color, flesh color, over color, firmness, and soluble solids content in the years 2015 and 2016. Linkage-based QTL analysis allowed the identification of genomic regions significantly associated with ripening time (LG4 of both parents and both phenotyping years), fruit skin color (LG3 and LG4 of both parents and both years), chlorophyll degradation index (LG3 of both parents in 2015) and fruit weight (LG7 of both parents in 2016). These results represent a promising situation for GBS in the identification of SNP variants associated to fruit quality traits, potentially applicable in breeding programs through MAS, in a highly heterozygous crop species such as Japanese plum.

Interval mapping for red/green skin color in Asian pears using a modified QTL-seq method

Pears with red skin are attractive to consumers and provide additional health benefits. Identification of the gene(s) responsible for skin coloration can benefit cultivar selection and breeding. The use of QTL-seq, a bulked segregant analysis method, can be problematic when heterozygous parents are involved. The present study modified the QTL-seq method by introducing a |Δ(SNP-index)| parameter to improve the accuracy of mapping the red skin trait in a group of highly heterozygous Asian pears. The analyses were based on mixed DNA pools composed of 28 red-skinned and 27 green-skinned pear lines derived from a cross between the 'Mantianhong' and 'Hongxiangsu' red-skinned cultivars. The 'Dangshansuli' cultivar genome was used as reference for sequence alignment. An average single-nucleotide polymorphism (SNP) index was calculated using a sliding window approach (200-kb windows, 20-kb increments). Nine scaffolds within the candidate QTL interval were in the fifth linkage group from 111.9 to 177.1 cM. There was a significant linkage between the insertions/deletions and simple sequence repeat markers designed from the candidate intervals and the red/green skin (R/G) locus, which was in a 582.5-kb candidate interval that contained 81 predicted protein-coding gene models and was composed of two subintervals at the bottom of the fifth chromosome. The ZFRI 130-16, In2130-12 and In2130-16 markers located near the R/G locus could potentially be used to identify the red skin trait in Asian pear populations. This study provides new insights into the genetics controlling the red skin phenotype in this fruit.

Lessons from a Phenotyping Center Revealed by the Genome-Guided Mapping of Powdery Mildew Resistance Loci

The genomics era brought unprecedented opportunities for genetic analysis of host resistance, but it came with the challenge that accurate and reproducible phenotypes are needed so that genomic results appropriately reflect biology. Phenotyping host resistance by natural infection in the field can produce variable results due to the uncontrolled environment, uneven distribution and genetics of the pathogen, and developmentally regulated resistance among other factors. To address these challenges, we developed highly controlled, standardized methodologies for phenotyping powdery mildew resistance in the context of a phenotyping center, receiving samples of up to 140 grapevine progeny per F₁ family. We applied these methodologies to F₁ families segregating for REN1- or REN2-mediated resistance and validated that some but not all bioassays identified the REN1 or REN2 locus. A point-intercept method (hyphal transects) to quantify colony density objectively at 8 or 9 days postinoculation proved to be the phenotypic response most reproducibly predicted by these resistance loci. Quantitative trait locus (QTL) mapping with genotyping-by-sequencing maps defined the REN1 and REN2 loci at relatively high resolution. In the reference PN40024 genome under each QTL, nucleotide-binding site-leucine-rich repeat candidate resistance genes were identified-one gene for REN1 and two genes for REN2. The methods described here for centralized resistance phenotyping and high-resolution genetic mapping can inform strategies for breeding resistance to powdery mildews and other pathogens on diverse, highly heterozygous hosts.