Genome-wide SNP discovery and population structure analysis in pepper (Capsicum annuum) using genotyping by sequencing

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BACKGROUND

Phytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, 'PWB-185', 'PWB-186', and '6347', at the 4-8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.

RESULTS

Most of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.

CONCLUSIONS

Results will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

Revitalizing agriculture: next-generation genotyping and -omics technologies enabling molecular prediction of resilient traits in the Solanaceae family

This review highlights -omics research in Solanaceae family, with a particular focus on resilient traits. Extensive research has enriched our understanding of Solanaceae genomics and genetics, with historical varietal development mainly focusing on disease resistance and cultivar improvement but shifting the emphasis towards unveiling resilience mechanisms in genebank-preserved germplasm is nowadays crucial. Collecting such information, might help researchers and breeders developing new experimental design, providing an overview of the state of the art of the most advanced approaches for the identification of the genetic elements laying behind resilience. Building this starting point, we aim at providing a useful tool for tackling the global agricultural resilience goals in these crops.

Ridge regression and deep learning models for genome-wide selection of complex traits in New Mexican Chile peppers

BACKGROUND

Genomewide prediction estimates the genomic breeding values of selection candidates which can be utilized for population improvement and cultivar development. Ridge regression and deep learning-based selection models were implemented for yield and agronomic traits of 204 chile pepper genotypes evaluated in multi-environment trials in New Mexico, USA.

RESULTS

Accuracy of prediction differed across different models under ten-fold cross-validations, where high prediction accuracy was observed for highly heritable traits such as plant height and plant width. No model was superior across traits using 14,922 SNP markers for genomewide selection. Bayesian ridge regression had the highest average accuracy for first pod date (0.77) and total yield per plant (0.33). Multilayer perceptron (MLP) was the most superior for flowering time (0.76) and plant height (0.73), whereas the genomic BLUP model had the highest accuracy for plant width (0.62). Using a subset of 7,690 SNP loci resulting from grouping markers based on linkage disequilibrium coefficients resulted in improved accuracy for first pod date, ten pod weight, and total yield per plant, even under a relatively small training population size for MLP and random forest models. Genomic and ridge regression BLUP models were sufficient for optimal prediction accuracies for small training population size. Combining phenotypic selection and genomewide selection resulted in improved selection response for yield-related traits, indicating that integrated approaches can result in improved gains achieved through selection.

CONCLUSIONS

Accuracy values for ridge regression and deep learning prediction models demonstrate the potential of implementing genomewide selection for genetic improvement in chile pepper breeding programs. Ultimately, a large training data is relevant for improved genomic selection accuracy for the deep learning models.

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.

Narrow genetic diversity in germplasm from the Guinean and Sudano-Guinean zones in Benin indicates the need to broaden the genetic base of sweet fig banana (Musa acuminata cv Sotoumon)

Sweet fig (M. acuminata cv. Sotoumon) is an economically important dessert banana in Benin, with high nutritional, medicinal, and cultural values. Nevertheless, its productivity and yield are threatened by biotic and abiotic stresses. Relevant knowledge of the genetic diversity of this economically important crop is essential for germplasm conservation and the development of breeding programs. However, very little is known about the genetic makeup of this cultivar in Benin. To advance the understanding of genetic diversity in sweet fig banana germplasm, a Genotype-By-Sequencing (GBS) was performed on a panel of 273 accessions collected in different phytogeographical zones of Benin. GBS generated 8,457 quality SNPs, of which 1992 were used for analysis after filtering. The results revealed a low diversity in the studied germplasm (He = 0.0162). Genetic differentiation was overall very low in the collection as suggested by the negative differentiation index (Fstg = -0.003). The Analysis of Molecular Variance (AMOVA) indicated that the variation between accessions within populations accounted for 83.8% of the total variation observed (P < 0.001). The analysis of population structure and neighbor-joining tree partitioned the germplasm into three clusters out of which a predominant major one contained 98.1% of all accessions. These findings demonstrate that current sweet fig banana genotypes shared a common genetic background, which made them vulnerable to biotic and abiotic stress. Therefore, broadening the genetic base of the crop while maintaining its quality attributes and improving yield performance is of paramount importance. Moreover, the large genetic group constitutes an asset for future genomic selection studies in the crop and can guide the profiling of its conservation strategies.

Association mapping in bambara groundnut [Vigna subterranea (L.) Verdc.] reveals loci associated with agro-morphological traits

BACKGROUND

Genome-wide association studies (GWAS) are important for the acceleration of crop improvement through knowledge of marker-trait association (MTA). This report used DArT SNP markers to successfully perform GWAS on agro-morphological traits using 270 bambara groundnut [Vigna subterranea (L.) Verdc.] landraces sourced from diverse origins. The study aimed to identify marker traits association for nine agronomic traits using GWAS and their candidate genes. The experiment was conducted at two different locations laid out in alpha lattice design. The cowpea [Vigna unguiculata (L.) Walp.] reference genome (i.e. legume genome most closely related to bambara groundnut) assisted in the identification of candidate genes.

RESULTS

The analyses showed that linkage disequilibrium was found to decay rapidly with an average genetic distance of 148 kb. The broadsense heritability was relatively high and ranged from 48.39% (terminal leaf length) to 79.39% (number of pods per plant). The GWAS identified a total of 27 significant marker-trait associations (MTAs) for the nine studied traits explaining 5.27% to 24.86% of phenotypic variations. Among studied traits, the highest number of MTAs was obtained from seed coat colour (6) followed by days to flowering (5), while the least is days to maturity (1), explaining 5.76% to 11.03%, 14.5% to 19.49%, and 11.66% phenotypic variations, respectively. Also, a total of 17 candidate genes were identified, varying in number for different traits; seed coat colour (6), days to flowering (3), terminal leaf length (2), terminal leaf width (2), number of seed per pod (2), pod width (1) and days to maturity (1).

CONCLUSION

These results revealed the prospect of GWAS in identification of SNP variations associated with agronomic traits in bambara groundnut. Also, its present new opportunity to explore GWAS and marker assisted strategies in breeding of bambara groundnut for acceleration of the crop improvement.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of wild tea plant in Guizhou Plateau

BACKGROUND

Tea, the second largest consumer beverage in the world after water, is widely cultivated in tropical and subtropical areas. However, the effect of environmental factors on the distribution of wild tea plants is unclear.

RESULTS

A total of 159 wild tea plants were collected from different altitudes and geological types of the Guizhou Plateau. Using the genotyping-by-sequencing method, 98,241 high-quality single nucleotide polymorphisms were identified. Genetic diversity, population structure analysis, principal component analysis, phylogenetic analysis, and linkage disequilibrium were performed. The genetic diversity of the wild tea plant population from the Silicate Rock Classes of Camellia gymnogyna was higher than that from the Carbonate Rock Classes of Camellia tachangensis. In addition, the genetic diversity of wild tea plants from the second altitude gradient was significantly higher than that of wild tea plants from the third and first altitude gradients. Two inferred pure groups (GP01 and GP02) and one inferred admixture group (GP03) were identified by population structure analysis and were verified by principal component and phylogenetic analyses. The highest differentiation coefficients were determined for GP01 vs. GP02, while the lowest differentiation coefficients were determined for GP01 vs. GP03.

CONCLUSIONS

This study revealed the genetic diversity and geographical distribution characteristics of wild tea plants in the Guizhou Plateau. There are significant differences in genetic diversity and evolutionary direction between Camellia tachangensis with Carbonate Rock Classes at the first altitude gradient and Camellia gymnogyna with Silicate Rock Classes at the third altitude gradient. Geological environment, soil mineral element content, soil pH, and altitude markedly contributed to the genetic differentiation between Camellia tachangensis and Camellia gymnogyna.

Genome-wide association study of leaf-related traits in tea plant in Guizhou based on genotyping-by-sequencing

BACKGROUND

Studying the genetic characteristics of tea plant (Camellia spp.) leaf traits is essential for improving yield and quality through breeding and selection. Guizhou Plateau, an important part of the original center of tea plants, has rich genetic resources. However, few studies have explored the associations between tea plant leaf traits and single nucleotide polymorphism (SNP) markers in Guizhou.

RESULTS

In this study, we used the genotyping-by-sequencing (GBS) method to identify 100,829 SNP markers from 338 accessions of tea germplasm in Guizhou Plateau, a region with rich genetic resources. We assessed population structure based on high-quality SNPs, constructed phylogenetic relationships, and performed genome-wide association studies (GWASs). Four inferred pure groups (G-I, G-II, G-III, and G-IV) and one inferred admixture group (G-V), were identified by a population structure analysis, and verified by principal component analyses and phylogenetic analyses. Through GWAS, we identified six candidate genes associated with four leaf traits, including mature leaf size, texture, color and shape. Specifically, two candidate genes, located on chromosomes 1 and 9, were significantly associated with mature leaf size, while two genes, located on chromosomes 8 and 11, were significantly associated with mature leaf texture. Additionally, two candidate genes, located on chromosomes 1 and 2 were identified as being associated with mature leaf color and mature leaf shape, respectively. We verified the expression level of two candidate genes was verified using reverse transcription quantitative polymerase chain reaction (RT-qPCR) and designed a derived cleaved amplified polymorphism (dCAPS) marker that co-segregated with mature leaf size, which could be used for marker-assisted selection (MAS) breeding in Camellia sinensis.

CONCLUSIONS

In the present study, by using GWAS approaches with the 338 tea accessions population in Guizhou, we revealed a list of SNPs markers and candidate genes that were significantly associated with four leaf traits. This work provides theoretical and practical basis for the genetic breeding of related traits in tea plant leaves.

Genetic Diversity and Population Structure of Chinese Chestnut (Castanea mollissima Blume) Cultivars Revealed by GBS Resequencing

Chinese chestnut (Castanea mollissima Bl.) is one of the earliest domesticated and cultivated fruit trees, and it is widely distributed in China. Because of the high quality of its nuts and its high resistance to abiotic and biotic stresses, Chinese chestnut could be used to improve edible chestnut varieties worldwide. However, the unclear domestication history and highly complex genetic background of Chinese chestnut have prevented the efficiency of breeding efforts. To explore the genetic diversity and structure of Chinese chestnut populations and generate new insights that could aid chestnut breeding, heterozygosity statistics, molecular variance analysis, ADMIXTURE analysis, principal component analysis, and phylogenetic analysis were conducted to analyze single nucleotide polymorphism data from 185 Chinese chestnut landraces from five geographical regions in China via genotyping by sequencing. Results showed that the genetic diversity level of the five populations from different regions was relatively high, with an observed heterozygosity of 0.2796-0.3427. The genetic diversity level of the population in the mid-western regions was the highest, while the population north of the Yellow River was the lowest. Molecular variance analysis showed that the variation among different populations was only 2.07%, while the intra-group variation reached 97.93%. The Chinese chestnut samples could be divided into two groups: a northern and southern population, separated by the Yellow River; however, some samples from the southern population were genetically closer to samples from the northern population. We speculate that this might be related to the migration of humans during the Han dynasty due to the frequent wars that took place during this period, which might have led to the introduction of chestnut to southern regions. Some samples from Shandong Province and Beijing City were outliers that did not cluster with their respective groups, and this might be caused by the special geographical, political, and economic significance of these two regions. The findings of our study showed the complex genetic relationships among Chinese chestnut landraces and the high genetic diversity of these resources.

Genetic diversity and marker trait association for phytophthora resistance in chilli

BACKGROUND

Chilli is an important commercial crop with positive returns tendency. Phytophthora root rot causes drastic damage to chilli plant. Dearth of detecting marker trait associations is a major hinderance in practicing marker assisted selection in chilli breeding.

METHODS AND RESULTS

Herein, 110 chilli accessions were assessed for 15 agronomic traits under control and disease infected conditions for two crop seasons (2018-2019). The SSR genotyping revealed high values of major allele frequency (MAF = 0.70), genetic diversity (GD = 0.39) and Polymorphic Information Content (PIC = 0.31). Principal coordinate analysis and population structure analysis showed distribution of diverse genotypes in all groups by dividing 110 genotypes in three populations and nine sub-populations. The UPGMA based Archaeopteryx tree was in concordance with population structure analysis. Linkage disequilibrium analysis evaluated that LD decays within 3-10 bp. Marker trait association (MTA) revealed the associations of 35 SSRs with 14 morphological traits. The significant MTA for marker CAeMS073 with relative leaf damage (RLD, 0.183 R²) under control and treated conditions was consistently observed in both models. The markers, CAMS173 and CAMS194 were found to be strongly associated with RLD and Disease Index (DI), respectively. The absence of MTA was detected for height of first branch.

CONCLUSION

The MTAs reported in this study can facilitate marker assisted breeding for developing chilli germplasm resistant against Phytophthora capsici.

Chile Pepper (Capsicum) Breeding and Improvement in the "Multi-Omics" Era

Chile pepper (Capsicum spp.) is a major culinary, medicinal, and economic crop in most areas of the world. For more than hundreds of years, chile peppers have "defined" the state of New Mexico, USA. The official state question, "Red or Green?" refers to the preference for either red or the green stage of chile pepper, respectively, reflects the value of these important commodities. The presence of major diseases, low yields, decreased acreages, and costs associated with manual labor limit production in all growing regions of the world. The New Mexico State University (NMSU) Chile Pepper Breeding Program continues to serve as a key player in the development of improved chile pepper varieties for growers and in discoveries that assist plant breeders worldwide. Among the traits of interest for genetic improvement include yield, disease resistance, flavor, and mechanical harvestability. While progress has been made, the use of conventional breeding approaches has yet to fully address producer and consumer demand for these traits in available cultivars. Recent developments in "multi-omics," that is, the simultaneous application of multiple omics approaches to study biological systems, have allowed the genetic dissection of important phenotypes. Given the current needs and production constraints, and the availability of multi-omics tools, it would be relevant to examine the application of these approaches in chile pepper breeding and improvement. In this review, we summarize the major developments in chile pepper breeding and present novel tools that can be implemented to facilitate genetic improvement. In the future, chile pepper improvement is anticipated to be more data and multi-omics driven as more advanced genetics, breeding, and phenotyping tools are developed.

Genome-wide diversity analysis to infer population structure and linkage disequilibrium among Colombian coconut germplasm

Genetic diversity and relatedness of accessions for coconut growing in Colombia was unknown until this study. Here we develop single nucleotide polymorphisms (SNPs) along the coconut genome based on Genotyping by Sequencing (GBS) with the goal of analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diverse coconut panel consisting of 112 coconut accessions from the Atlantic and Pacific coasts of Colombia. A comprehensive catalog of approximately 40,000 SNPs with a minor allele frequency (MAF) of > 0.05 is presented. A total of 40,614 SNPs were found but only 19,414 anchored to chromosomes. Of these, 10,338 and 4606 were exclusive to the Atlantic and Pacific gene pools, respectively, and 3432 SNPs could differentiate both gene pools. A filtered subset of unlinked and anchored SNPs (1271) showed a population structure at K = 4, separating accessions from the Pacific and Atlantic coasts that can also be distinguished by palm height, as found in previous studies. The Pacific groups had a slow LD decay, low Fixation Index (Fst) and low nucleotide diversity (π), while the Atlantic group had slightly higher genetic diversity and faster LD decay. Genome-wide diversity analyses are of importance to promote germplasm conservation and breeding programs aimed at developing new cultivars better adapted to the region.

Genetic Diversity and Population Structure of Doum Palm (Hyphaene compressa) Using Genotyping by Sequencing

Doum palm (Hyphaene compressa) is a perennial economic plant primarily growing in Kenya’s Arid and Semi-Arid Lands (ASALs). It is heavily relied upon for food, animal feed, construction materials and medicine, making it an ideal plant for resource sustainability. However, the limited information on its genetic resources has hindered its breeding and conservation studies. This study used the genotyping by sequencing approach to identify Single Nucleotide Polymorphisms. These SNPs were further used to assess the genetic diversity and population structure of 96 H. compressa accessions from Coastal, Northern and Eastern ASAL regions of Kenya using two approaches; reference-based and de novo-based assemblies. STRUCTURE analysis grouped the sampled accessions into two genetic clusters (Cluster 1 and Cluster 2). Cluster 1 included accessions from the Northern region, whereas Cluster 2 included all accessions from Eastern and Coastal regions. Accessions from Kwale (Coastal) had mixed ancestry from both Cluster 1 and Cluster 2. These STRUCTURE findings were further supported by principal components analysis, discriminant analysis of principal components and phylogenetic analysis. Analysis of molecular variance indicated greater genetic variation within populations (92.7%) than among populations (7.3%). An overall F_ST of 0.074 was observed, signifying moderate genetic differentiation among populations. The results of this study will provide information useful in breeding, marker-assisted selection and conservation management of H. compressa.

Genetic diversity and population structure of rice (Oryza sativa L.) landraces from Kerala, India analyzed through genotyping-by-sequencing

Researchers stand at the vanguard of advancement and application of next-generation sequencing technology for developing dominant strategies for the sustainable management of genetically diverse crops. We attempt to fill the existing research lacuna in the molecular characterization of potent rice landraces in Kerala. Genotyping-by-sequencing (GBS) was performed on 96 Kerala rice accessions to identify single-nucleotide polymorphisms (SNPs), to examine the genetic diversity, population structure, and to delineate linkage disequilibrium (LD) pattern. GBS identified 5856 high-quality SNPs. The structure analysis indicated three subpopulations with the highest probability for population clustering with significant genetic differentiation, confirmed by principal component analysis. The genome-wide LD decay distance was 772 kb, at which the r² dropped to half its maximum value. The analysis of genetic properties of the identified SNP panel with an average polymorphism information content (PIC) value of 0.22 and a minor allele frequency (MAF) > 0.1 unveiled their efficacy in genome-wide association studies (GWAS). High F_ST (0.266) and low Nm (0.692) portray a strong genetic differentiation among the rice landraces, complementing the genetic structuring observed in the studied population. Slow LD decay in the rice landraces reflects their self-pollinating behavior and the indirect selection of desired traits by domestication. Moreover, the high LD entails only a minimum number of SNP markers for detecting marker-trait association. The diverse germplasm utilized in this study can be further utilized to disclose genetic variants associated with phenotypic traits and define signatures of selection via GWAS and selective sweep, respectively.

Harnessing Crop Wild Diversity for Climate Change Adaptation

Warming and drought are reducing global crop production with a potential to substantially worsen global malnutrition. As with the green revolution in the last century, plant genetics may offer concrete opportunities to increase yield and crop adaptability. However, the rate at which the threat is happening requires powering new strategies in order to meet the global food demand. In this review, we highlight major recent 'big data' developments from both empirical and theoretical genomics that may speed up the identification, conservation, and breeding of exotic and elite crop varieties with the potential to feed humans. We first emphasize the major bottlenecks to capture and utilize novel sources of variation in abiotic stress (i.e., heat and drought) tolerance. We argue that adaptation of crop wild relatives to dry environments could be informative on how plant phenotypes may react to a drier climate because natural selection has already tested more options than humans ever will. Because isolated pockets of cryptic diversity may still persist in remote semi-arid regions, we encourage new habitat-based population-guided collections for genebanks. We continue discussing how to systematically study abiotic stress tolerance in these crop collections of wild and landraces using geo-referencing and extensive environmental data. By uncovering the genes that underlie the tolerance adaptive trait, natural variation has the potential to be introgressed into elite cultivars. However, unlocking adaptive genetic variation hidden in related wild species and early landraces remains a major challenge for complex traits that, as abiotic stress tolerance, are polygenic (i.e., regulated by many low-effect genes). Therefore, we finish prospecting modern analytical approaches that will serve to overcome this issue. Concretely, genomic prediction, machine learning, and multi-trait gene editing, all offer innovative alternatives to speed up more accurate pre- and breeding efforts toward the increase in crop adaptability and yield, while matching future global food demands in the face of increased heat and drought. In order for these 'big data' approaches to succeed, we advocate for a trans-disciplinary approach with open-source data and long-term funding. The recent developments and perspectives discussed throughout this review ultimately aim to contribute to increased crop adaptability and yield in the face of heat waves and drought events.

Single nucleotide polymorphisms reveal genetic diversity in New Mexican chile peppers (Capsicum spp.)

BACKGROUND

Chile peppers (Capsicum spp.) are among the most important horticultural crops in the world due to their number of uses. They are considered a major cultural and economic crop in the state of New Mexico in the United States. Evaluating genetic diversity in current New Mexican germplasm would facilitate genetic improvement for different traits. This study assessed genetic diversity, population structure, and linkage disequilibrium (LD) among 165 chile pepper genotypes using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS).

RESULTS

A GBS approach identified 66,750 high-quality SNP markers with known map positions distributed across the 12 chromosomes of Capsicum. Principal components analysis revealed four distinct clusters based on species. Neighbor-joining phylogenetic analysis among New Mexico State University (NMSU) chile pepper cultivars showed two main clusters, where the C. annuum genotypes grouped together based on fruit or pod type. A Bayesian clustering approach for the Capsicum population inferred K = 2 as the optimal number of clusters, where the C. chinense and C. frutescens grouped in a single cluster. Analysis of molecular variance revealed majority of variation to be between the Capsicum species (76.08 %). Extensive LD decay (~ 5.59 Mb) across the whole Capsicum population was observed, demonstrating that a lower number of markers would be required for implementing genome wide association studies for different traits in New Mexican type chile peppers. Tajima's D values demonstrated positive selection, population bottleneck, and balancing selection for the New Mexico Capsicum population. Genetic diversity for the New Mexican chile peppers was relatively low, indicating the need to introduce new alleles in the breeding program to broaden the genetic base of current germplasm.

CONCLUSIONS

Genetic diversity among New Mexican chile peppers was evaluated using GBS-derived SNP markers and genetic relatedness on the species level was observed. Introducing novel alleles from other breeding programs or from wild species could help increase diversity in current germplasm. We present valuable information for future association mapping and genomic selection for different traits for New Mexican chile peppers for genetic improvement through marker-assisted breeding.

Genotyping by Sequencing-Based Discovery of SNP Markers and Construction of Linkage Map from F5 Population of Pepper with Contrasting Powdery Mildew Resistance Trait

Powdery mildew (PM) is a common fungal disease infecting pepper plants worldwide. Molecular breeding of pepper cultivars with powdery mildew resistance is desirable for the economic improvement of pepper cultivation. In the present study, 188 F5 population derived from AR1 (PM resistant) and TF68 (PM sensitive) parents were subjected to high-throughput genotyping by sequencing (GBS) for the identification of single nucleotide polymorphism (SNP) markers. Further, the identified SNP markers were utilized for the construction of genetic linkage map and QTL analysis. Overall read mapping percentage of 87.29% was achieved in this study with the total length of mapped region ranging from 2,956,730 to 25,537,525 bp. A total of 41,111 polymorphic SNPs were identified, and a final of 1,841 SNPs were filtered for the construction of a linkage map. A total of 12 linkage groups were constructed corresponding to each chromosome with 1,308 SNP markers with the map length of 2506.8 cM. Further, two QTLs such as Pm-2.1 and Pm-5.1 were identified in chromosomes 2 and 5, respectively, for the PM resistance. Overall, the outcomes of the present endeavor can be utilized for the marker-assisted selection of pepper with powdery mildew-resistant trait.

Genotyping-by-Sequencing Technology in Plant Taxonomy and Phylogeny

Genotyping-by-sequencing (GBS) is a method to discover and genotype simultaneous genome-wide high-throughput single nucleotide polymorphisms (SNPs). GBS is based on reducing genome complexity with restriction enzymes. Here we describe a method developed by Elshire et al. for constructing simplified GBS libraries and recent bioinformatic approaches developed to analyze the large volume of polymorphism data generated by this method. GBS approach is suitable for population studies, taxonomic and phylogenic studies, germplasm characterization, and breeding and trait mapping for a wide range of organisms, including plants with complex genomes.

Development of a Panel of Genotyping-in-Thousands by Sequencing in Capsicum

Genotyping by sequencing (GBS) enables genotyping of multiple loci at low cost. However, the single nucleotide polymorphisms (SNPs) revealed by GBS tend to be randomly distributed between individuals, limiting their direct comparisons without applying the various filter options to obtain a comparable dataset of SNPs. Here, we developed a panel of a multiplex targeted sequencing method, genotyping-in-thousands by sequencing (GT-seq), to genotype SNPs in Capsicum spp. Previously developed Fluidigm^® SNP markers were converted to GT-seq markers and combined with new GT-seq markers developed using SNP information obtained through GBS. We then optimized multiplex PCR conditions: we obtained the highest genotyping rate when the first PCR consisted of 25 cycles. In addition, we determined that 101 primer pairs performed best when amplifying target sequences of 79 bp. We minimized interference of multiplex PCR by primer dimer formation using the PrimerPooler program. Using our GT-seq pipeline on Illumina Miseq and Nextseq platforms, we genotyped up to 1,500 (Miseq) and 1,300 (Nextseq) samples for the optimum panel size of 100 loci. To allow the genotyping of Capsicum species, we designed 332 informative GT-seq markers from Fluidigm SNP markers and GBS-derived SNPs. This study illustrates the first application of GT-seq in crop plants. The GT-seq marker set developed here will be a useful tool for molecular breeding of peppers in the future.

Large scale SNP unearthing and genetic architecture analysis in sea-captured and cultured populations of Cynoglossus semilaevis

Knowledge on population structure and genetic diversity is a focal point for association mapping studies and genomic selection. Genotyping by sequencing (GBS) represents an innovative method for large scale SNP detection and genotyping of genetic resources. Here we used the GBS approach for the genome-wide identification of SNPs in a collection of Cynoglossus semilaevis and for the assessment of the level of genetic diversity in C. semilaevis genotypes. GBS analysis generated a total of 55.12 Gb high-quality sequence data, with an average of 0.63 Gb per sample. The total number of SNP markers was 563, 109. In order to explore the genetic diversity of C. semilaevis and to select a minimal core set representing most of the total genetic variation with minimum redundancy, C. semilaevis sequences were analyzed using high quality SNPs. Based on hierarchical clustering, it was possible to divide the collection into 2 clusters. The marine fishing populations were clustered and clearly separated from the cultured populations, and the cultured populations from Hebei was also distinct from the other two local populations. These analyses showed that genotypes were clustered based on species-related features. Differential significant SNPs were also captured and validated by GBS and SNaPshot, with linkage disequilibrium and haplotype analysis, seven SNPs have been confirmed to have obvious differentiation in two populations, which may be used as the characteristic evaluation sites of sea-captured and cultured Cynoglossus semilaevis populations. And SNP markers and information on population structure developed in this study will undoubtedly support genome-wide association mapping studies and marker-assisted selection programs. These differential SNPs could be also employed as the characteristic evaluation sites of sea-captured and cultured Cynoglossus semilaevis populations in future.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Resolving complicated relationships of the Panax bipinnatifidus complex in southwestern China by RAD-seq data

The P. binpinnatifidus complex included most of the Panax species distributed in Sino-Himalaya regions except for P. pseudoginseng, P. stipuleanatus and P. notoginseng. However, the delimitation and identification of these taxa within the species complex are very difficult due to the existence of morphological intermediates, and their evolutionary relationships remain unresolved despite several studies have been carried out based on traditional DNA markers. The taxonomic uncertainty hinders the identification, conservation and exploration of these wild populations of Panax. To study this species complex, we employed ddRAD-seq data of these taxa from 18 different localities of southwestern China, using two RAD analysis pipelines, STACKS and pyRAD. Based on the results of phylogenetic analysis, the species complex was divided into four clades with high supports, which largely agreed with morphologically described species. Two clades, corresponding to P. vietnamensis and P. zingiberensis, respectively, were sister groups, indicating that these two species had a closer genetic relationship; the third clade was consisted of samples with bamboo-like rhizomes named as P. wangianus clade, and the fourth one with moniliform rhizomes was named as P. bipinnatifidus clade. The population genetic structure analysis and D-statistics test showed the localized admixture among these species, which indicated that introgression had occurred among the related lineages continuously distributed in southeastern Yunnan and adjacent regions.

Heterogeneity in the distribution of 159 drug-response related SNPs in world populations and their genetic relatedness

Interethnic variability in drug response arises from genetic differences associated with drug metabolism, action and transport. These genetic variations can affect drug efficacy as well as cause adverse drug reactions (ADRs). We retrieved drug-response related single nucleotide polymorphism (SNP) associated data from databases and analyzed to elucidate population specific distribution of 159 drug-response related SNPs in twenty six populations belonging to five super-populations (African, Admixed Americans, East Asian, European and South Asian). Significant interpopulation differences exist in the minor (variant) allele frequencies (MAFs), linkage disequilibrium (LD) and haplotype distributions among these populations. 65 of the drug-response related alleles, which are considered as minor (variant) in global population, are present as the major alleles (frequency ≥0.5) in at least one or more populations. Populations that belong to the same super-population have similar distribution pattern for majority of the variant alleles. These drug response related variant allele frequencies and their pairwise LD measure (r²) can clearly distinguish the populations in a way that correspond to the known evolutionary history of human and current geographic distributions, while D' cannot. The data presented here may aid in identifying drugs that are more appropriate and/or require pharmacogenetic testing in these populations. Our findings emphasize on the importance of distinct, ethnicity-specific clinical guidelines, especially for the African populations, to avoid ADRs and ensure effective drug treatment.

Genome-wide association mapping of QTLs implied in potato virus Y population sizes in pepper: evidence for widespread resistance QTL pyramiding

In this study, we looked for genetic factors in the pepper (Capsicum annuum) germplasm that control the number of potato virus Y (PVY) particles entering the plant (i.e. effective population size at inoculation) and the PVY accumulation at the systemic level (i.e. census population size). Using genotyping-by-sequencing (GBS) in a core collection of 256 pepper accessions, we obtained 10 307 single nucleotide polymorphisms (SNPs) covering the whole genome. Genome-wide association studies (GWAS) detected seven SNPs significantly associated with the virus population size at inoculation and/or systemic level on chromosomes 4, 6, 9 and 12. Two SNPs on chromosome 4 associated with both PVY population sizes map closely to the major resistance gene pvr2 encoding the eukaryotic initiation factor 4E. No obvious candidates for resistance were identified in the confidence intervals for the other chromosomes. SNPs detected on chromosomes 6 and 12 colocalized with resistance quantitative trait loci (QTLs) previously identified with a biparental population. These results show the efficiency of GBS and GWAS in C. annuum, indicate highly consistent results between GWAS and classical QTL mapping, and suggest that resistance QTLs identified with a biparental population are representative of a much larger collection of pepper accessions. Moreover, the resistance alleles at these different loci were more frequently combined than expected by chance in the core collection, indicating widespread pyramiding of resistance QTLs and widespread combination of resistance QTLs and major effect genes. Such pyramiding may increase resistance efficiency and/or durability.

Genome-Wide SNP discovery and genomic characterization in avocado (Persea americana Mill.)

Modern crop breeding is based on the use of genetically and phenotypically diverse plant material and, consequently, a proper understanding of population structure and genetic diversity is essential for the effective development of breeding programs. An example is avocado, a woody perennial fruit crop native to Mesoamerica with an increasing popularity worldwide. Despite its commercial success, there are important gaps in the molecular tools available to support on-going avocado breeding programs. In order to fill this gap, in this study, an avocado 'Hass' draft assembly was developed and used as reference to study 71 avocado accessions which represent the three traditionally recognized avocado horticultural races or subspecies (Mexican, Guatemalan and West Indian). An average of 5.72 M reads per individual and a total of 7,108 single nucleotide polymorphism (SNP) markers were produced for the 71 accessions analyzed. These molecular markers were used in a study of genetic diversity and population structure. The results broadly separate the accessions studied according to their botanical race in four main groups: Mexican, Guatemalan, West Indian and an additional group of Guatemalan × Mexican hybrids. The high number of SNP markers developed in this study will be a useful genomic resource for the avocado community.

Target sequencing reveals genetic diversity, population structure, core-SNP markers, and fruit shape-associated loci in pepper varieties

BACKGROUND

The widely cultivated pepper (Capsicum spp.) is one of the most diverse vegetables; however, little research has focused on characterizing the genetic diversity and relatedness of commercial varieties grown in China. In this study, a panel of 92 perfect single-nucleotide polymorphisms (SNPs) was identified using re-sequencing data from 35 different C. annuum lines. Based on this panel, a Target SNP-seq genotyping method was designed, which combined multiplex amplification of perfect SNPs with Illumina sequencing, to detect polymorphisms across 271 commercial pepper varieties.

RESULTS

The perfect SNPs panel had a high discriminating capacity due to the average value of polymorphism information content, observed heterozygosity, expected heterozygosity, and minor allele frequency, which were 0.31, 0.28, 0.4, and 0.31, respectively. Notably, the studied pepper varieties were morphologically categorized based on fruit shape as blocky-, long horn-, short horn-, and linear-fruited. The long horn-fruited population exhibited the most genetic diversity followed by the short horn-, linear-, and blocky-fruited populations. A set of 35 core SNPs were then used as kompetitive allele-specific PCR (KASPar) markers, another robust genotyping technique for variety identification. Analysis of genetic relatedness using principal component analysis and phylogenetic tree construction indicated that the four fruit shape populations clustered separately with limited overlaps. Based on STRUCTURE clustering, it was possible to divide the varieties into five subpopulations, which correlated with fruit shape. Further, the subpopulations were statistically different according to a randomization test and F_st statistics. Nine loci, located on chromosomes 1, 2, 3, 4, 6, and 12, were identified to be significantly associated with the fruit shape index (p < 0.0001).

CONCLUSIONS

Target SNP-seq developed in this study appears as an efficient power tool to detect the genetic diversity, population relatedness and molecular breeding in pepper. Moreover, this study demonstrates that the genetic structure of Chinese pepper varieties is significantly influenced by breeding programs focused on fruit shape.

Genome Wide Association Study and Genomic Selection of Amino Acid Concentrations in Soybean Seeds

Soybean is a major source of protein for human consumption and animal feed. Releasing new cultivars with high nutritional value is one of the major goals in soybean breeding. To achieve this goal, genome-wide association studies of seed amino acid contents were conducted based on 249 soybean accessions from China, US, Japan, and South Korea. The accessions were evaluated for 15 amino acids and genotyped by sequencing. Significant genetic variation was observed for amino acids among the accessions. Among the 231 single nucleotide polymorphisms (SNPs) significantly associated with variations in amino acid contents, fifteen SNPs localized near 14 candidate genes involving in amino acid metabolism. The amino acids were classified into two groups with five in one group and seven amino acids in the other. Correlation coefficients among the amino acids within each group were high and positive, but the correlation coefficients of amino acids between the two groups were negative. Twenty-five SNP markers associated with multiple amino acids can be used to simultaneously improve multi-amino acid concentration in soybean. Genomic selection analysis of amino acid concentration showed that selection efficiency of amino acids based on the markers significantly associated with all 15 amino acids was higher than that based on random markers or markers only associated with individual amino acid. The identified markers could facilitate selection of soybean varieties with improved seed quality.

Genetic diversity and structure of Capsicum annuum as revealed by start codon targeted and directed amplified minisatellite DNA markers

BACKGROUND

Identification of high resolving DNA-based markers is of paramount importance to unlock the potential of genetic diversity and selection of unique accessions of Capsicum annuum L., within Cross River and Ebonyi States of Nigeria, for breeding and conservation. Therefore, we comparatively explored the effectiveness of start codon targeted (SCoT) and directed amplified minisatellite DNA (DAMD) markers for diversity analysis of the accessions. Fifteen accessions were collected for DNA extraction and amplifications with the markers.

RESULTS

Dendrograms from SCoT and DAMD categorized the accessions into five and three genetic groups, respectively, while the principal component analysis identified five genetic clusters, each from the markers. The average values of allele, gene diversity and polymorphic information content detected with SCoT and DAMD demonstrate that the two markers were effective and efficient, especially, SCoT in genetic diversity study of the accessions of pepper. Number of polymorphic loci (NPL) and percentage polymorphic loci (PPL) from SCoT (NPL = 64, PPL = 80.00-95.73%) and DAMD (NPL = 56, PPL = 53.33-86.67%) were high, but higher in SCoT markers. Other effective genetic parameters (effective number of alleles, Nei's genetic diversity and Shannon's information indices) identified with the two marker systems elucidated the allelic richness, rich genetic diversity within the populations and informative nature of the markers, especially SCoT. The intraspecific genetic diversity, interspecific genetic diversity, and coefficient of differentiation obtained with SCoT and DAMD further exposed the genetic structure with more genetic divergence within than among the populations of the accessions. Estimate of gene flow from the SCoT markers was 3.8375 and 0.6.2042 for the DAMD markers. The estimate of gene flow values from the markers indicated extensiveness with SCoT (Nm = 3.8375) and extremely extensive with DAMD (Nm = 6.2042) among the populations.

CONCLUSION

This study shows that SCoT markers may be more useful and informative than DAMD in measuring genetic diversity and differentiation of the accessions of the genus Capsicum. Genetic parameters obtained with SCoT showed that the accessions from Cross River were more genetically diverse than the ones from Ebonyi State. Therefore, SCoT may be a preferred marker in evaluating genetic diversity for improvement and conservation of this spicy crop, C. capsicum.

GWAS Discovery Of Candidate Genes for Yield-Related Traits in Peanut and Support from Earlier QTL Mapping Studies

Peanut (Arachis hypogaea L.) is one of the most important oil crops worldwide, and its yet increasing market demand may be met by genetic improvement of yield related traits, which may be facilitated by a good understanding of the underlying genetic base of these traits. Here, we have carried out a genome-wide association study (GWAS) with the aim to identify genomic regions and the candidate genes within these regions that may be involved in determining the phenotypic variation at seven yield-related traits in peanut. For the GWAS analyses, 195 peanut accessions were phenotyped and/or genotyped; the latter was done using a genotyping-by-sequencing approach, which produced a total of 13,435 high-quality single nucleotide polymorphisms (SNPs). Analyses of these SNPs show that the analyzed peanut accessions can be approximately grouped into two big groups that, to some extent, agree with the botanical classification of peanut at the subspecies level. By taking this genetic structure as well as the relationships between the analyzed accessions into consideration, our GWAS analyses have identified 93 non-overlapping peak SNPs that are significantly associated with four of the studied traits. Gene annotation of the genome regions surrounding these peak SNPs have found a total of 311 unique candidate genes. Among the 93 yield-related-trait-associated SNP peaks, 12 are found to be co-localized with the quantitative trait loci (QTLs) that were identified by earlier related QTL mapping studies, and these 12 SNP peaks are only related to three traits and are almost all located on chromosomes Arahy.05 and Arahy.16. Gene annotation of these 12 co-localized SNP peaks have found 36 candidates genes, and a close examination of these candidate genes found one very interesting gene (arahy.RI9HIF), the rice homolog of which produces a protein that has been shown to improve rice yield when over-expressed. Further tests of the arahy.RI9HIF gene, as well as other candidate genes especially those within the more confident co-localized genomic regions, may hold the potential for significantly improving peanut yield.

Genotyping by Sequencing of Cultivated Lentil (Lens culinaris Medik.) Highlights Population Structure in the Mediterranean Gene Pool Associated With Geographic Patterns and Phenotypic Variables

Cultivated lentil (Lens culinaris Medik.) is one of the oldest domesticated crops and one of the most important grain legumes worldwide. The Mediterranean Basin holds large part of lentil biodiversity; however, no genetic structure was defined within the Mediterranean gene pool. In this study, we used high-throughput genotyping by sequencing to resolve the genetic structure of the Mediterranean ex situ lentil collection held at the Italian National Research Council. Sequencing of a 188-plex genotyping-by-sequencing library and bioinformatics treatment of data yielded 6,693 single nucleotide polymorphisms. Analysis of nonredundant genotypes with nonparametric and parametric methods highlighted the occurrence of five highly differentiated genetic clusters. Clustering could be related to geographic patterns and phenotypic traits, indicating that post-domestication routes introducing cultivation in Mediterranean countries and selection were major forces shaping lentil population structure. The estimation of the fixation index F_ST at individual single nucleotide polymorphism loci allowed the identification of distinctive alleles across clusters, suggesting the possibility to set up molecular keys for the assignment of lentil germplasm to specific genetic groups. Finally, significant associations between markers and phenotypic data were identified. Overall, the results of this study are of major importance for lentil conservation genetics and breeding and provide insights on the lentil evolutionary history.

Genetic Characterization of Apulian Olive Germplasm as Potential Source in New Breeding Programs

The olive is a fruit tree species with a century-old history of cultivation in theMediterranean basin. In Apulia (Southern Italy), the olive is of main social, cultural and economicimportance, and represents a hallmark of the rural landscape. However, olive cultivation in thisregion is threatened by the recent spread of the olive quick decline syndrome (OQDS) disease, thusthere is an urgent need to explore biodiversity and search for genetic sources of resistance. Herein,a genetic variation in Apulian olive germplasm was explored, as a first step to identify genotypeswith enhanced bio-agronomic traits, including resistance to OQDS. A preselected set of nuclearmicrosatellite markers allowed the acquisition of genotypic profiles, and to define geneticrelationships between Apulian germplasm and widespread cultivars. The analysis highlighted thebroad genetic variation in Apulian accessions and the presence of different unique genetic profiles.The results of this study lay a foundation for the organization of new breeding programs for olivegenetic improvement.

Genomic diversity and novel genome-wide association with fruit morphology in Capsicum, from 746k polymorphic sites

Capsicum is one of the major vegetable crops grown worldwide. Current subdivision in clades and species is based on morphological traits and coarse sets of genetic markers. Broad variability of fruits has been driven by breeding programs and has been mainly studied by linkage analysis. We discovered 746k variable sites by sequencing 1.8% of the genome in a collection of 373 accessions belonging to 11 Capsicum species from 51 countries. We describe genomic variation at population-level, confirm major subdivision in clades and species, and show that the known major subdivision of C. annuum separates large and bulky fruits from small ones. In C. annuum, we identify four novel loci associated with phenotypes determining the fruit shape, including a non-synonymous mutation in the gene Longifolia 1-like (CA03g16080). Our collection covers all the economically important species of Capsicum widely used in breeding programs and represent the widest and largest study so far in terms of the number of species and number of genetic variants analyzed. We identified a large set of markers that can be used for population genetic studies and genetic association analyses. Our results provide a comprehensive and precise perspective on genomic variability in Capsicum at population-level and suggest that future fine genetic association studies will yield useful results for breeding.

Genome-Wide Association Study of Resistance to Xanthomonas gardneri in the USDA Pepper (Capsicum) Collection

Bacterial spot Xanthomonas species cause significant disease outbreaks on tomato and pepper in tropical and subtropical regions throughout the world. Host resistance has been one of the key components of integrated disease management approaches to mitigate plant pathogens. Although a number of resistance genes have been identified in pepper against bacterial spot xanthomonads, emergence of bacterial spot pathogen variants capable of overcoming these sources and changing pathogen distributions reinforce the importance of identifying novel candidates to incorporate into breeding programs. We conducted a genome-wide association study (GWAS) on a diverse U.S. Department of Agriculture collection of pepper germplasm including different species of Capsicum to identify novel sources of disease resistance against a highly virulent X. gardneri strain isolated from a recent outbreak. GWAS identified highly significant single nucleotide polymorphisms associated with defoliation in response to infection with X. gardneri. Functionally relevant candidate genes encoded products involved in disease resistance/susceptibility, hormone signaling, and basal resistance against multiple pathogens in various host-pathogen systems. The X. gardneri-resistant genotypes and quantitative trait loci identified in this study provide alleles that could be used for a resistance gene pyramiding effort against different species of bacterial spot xanthomonads in pepper.

Translation initiation codon (ATG) or SCoT markers-based polymorphism study within and across various Capsicum accessions: insight from their amplification, cross-transferability and genetic diversity

Being an economical and nutritional crop, Capsicum appeases people's peppery taste and is found to bewidely distributed all over the world having vast diversity. In the present study, genetic polymorphism, cross transferability (CT) and genetic diversity were examined among the 54 different accessions of Capsicum species including 49 of Capsicum annuum, three of C. baccatum and two of C. frutescens, using a set of 36 start codon targeted (SCoT) primers. Of the total, 35 SCoT markers showed successful amplification profile among chilli germplasms and an average primer polymorphism was reported as 81.52% which ranged from 50% (SCoT-6) to 100% (SCoT-11). A total of 365 amplicons were obtained with an average of 10.43 bands per primer and the length of the bands ranged from 150 bp to 1.2 kb. Further, polymorphic information content value of SCoT markers ranged from 0.42 (for SCoT-25) to 0.86 (SCoT-27) with an average of 0.78. The average value of CT of SCoT markers was 44.08% ranged from 14.25% to 57.26% among different chilli accessions. A dendrogram was constructed and established genetic relationship among 54 capsicum species, with the help of translation initiation codon polymorphisms or SCoT primer amplification. This study suggests the effectiveness of SCoT marker system for characterizing and assessing genetic diversity of Capsicum germplasm, which can be used for evolutionary studies and to identify agronomically important traits.

Genetic diversity and population structure of Ethiopian Capsicum germplasms

We established a collection of 142 Capsicum genotypes from different geographical areas of Ethiopia with the aim of capturing genetic diversity. Morphological traits and high-resolution melting analysis distinguished one Capsicum baccatum, nine Capsicum frutescens and 132 Capsicum annuum accessions in the collection. Measurement of plant growth parameters revealed variation between germplasms in parameters including plant height, stem thickness, internode length, number of side branches, fruit width, and fruit length. Broad-sense heritability was maximum for fruit weight, followed by length and width of leaves. We used genotyping by sequencing (GBS) to identify single-nucleotide polymorphisms (SNPs) in the panel of 142 Capsicum germplasms and found 2,831,791 genome-wide SNP markers. Among these, we selected 53,284 high-quality SNPs and used them to estimate the level of genetic diversity, population structure, and phylogenetic relationships. From model-based ancestry analysis, the phylogenetic tree and principal-coordinate analysis (PCoA), we identified two distinct genetic populations: one comprising 132 C. annuum accessions and the other comprising the nine C. frutescens accessions. GWAS analysis detected 509 SNP markers that were significantly associated with fruit-, stem- and leaf-related traits. This is the first comprehensive report of the analysis of genetic variation in Ethiopian Capsicum species involving a large number of accessions. The results will help breeders utilize the germplasm collection to improve existing commercial cultivars.

Development and Application of InDel Markers for Capsicum spp. Based on Whole-Genome Re-Sequencing

Genome-wide identification of Insertion/Deletion polymorphisms (InDels) in Capsicum spp. was performed through comparing whole-genome re-sequencing data from two Capsicum accessions, C. annuum cv. G29 and C. frutescens cv. PBC688, with the reference genome sequence of C. annuum cv. CM334. In total, we identified 1,664,770 InDels between CM334 and PBC688, 533,523 between CM334 and G29, and 1,651,856 between PBC688 and G29. From these InDels, 1605 markers of 3-49 bp in length difference between PBC688 and G29 were selected for experimental validation: 1262 (78.6%) showed polymorphisms, 90 (5.6%) failed to amplify, and 298 (18.6%) were monomorphic. For further validation of these InDels, 288 markers were screened across five accessions representing five domesticated species. Of these assayed markers, 194 (67.4%) were polymorphic, 87 (30.2%) monomorphic and 7 (2.4%) failed. We developed three interspecific InDels, which associated with three genes and showed specific amplification in five domesticated species and clearly differentiated the interspecific hybrids. Thus, our novel PCR-based InDel markers provide high application value in germplasm classification, genetic research and marker-assisted breeding in Capsicum species.

Capsaicinoids: Pungency beyond Capsicum

Capsaicinoids are metabolites responsible for the appealing pungency of Capsicum (chili pepper) species. The completion of the Capsicum annuum genome has sparked new interest into the development of biotechnological applications involving the manipulation of pungency levels. Pungent dishes are already part of the traditional cuisine in many countries, and numerous health benefits and industrial applications are associated to capsaicinoids. This raises the question of how to successfully produce more capsaicinoids, whose biosynthesis is strongly influenced by genotype-environment interactions in fruits of Capsicum. In this Opinion article we propose that activating the capsaicinoid biosynthetic pathway in a more amenable species such as tomato could be the next step in the fascinating story of pungent crops.

Genetic diversity, population structure, and relationships in a collection of pepper (Capsicum spp.) landraces from the Spanish centre of diversity revealed by genotyping-by-sequencing (GBS)

Pepper (Capsicum spp.) is one of the most important vegetable crops; however, pepper genomic studies lag behind those of other important Solanaceae. Here we present the results of a high-throughput genotyping-by-sequencing (GBS) study of a collection of 190 Capsicum spp. accessions, including 183 of five cultivated species (C. annuum, C. chinense, C. frutescens, C. baccatum, and C. pubescens) and seven of the wild form C. annuum var. glabriusculum. Sequencing generated 6,766,231 high-quality read tags, of which 40.7% were successfully aligned to the reference genome. SNP calling yielded 4083 highly informative segregating SNPs. Genetic diversity and relationships of a subset of 148 accessions, of which a complete passport information was available, was studied using principal components analysis (PCA), discriminant analysis of principal components (DAPC), and phylogeny approaches. C. annuum, C. baccatum, and C. chinense were successfully separated by all methods. Our population was divided into seven clusters by DAPC, where C. frutescens accessions were clustered together with C. chinense. C. annuum var. glabriusculum accessions were spread into two distinct genetic pools, while European accessions were admixed and closely related. Separation of accessions was mainly associated to differences in fruit characteristics and origin. Phylogeny studies showed a close relation between Spanish and Mexican accessions, supporting the hypothesis that the first arose from a main genetic flow from the latter. Tajima's D statistic values were consistent with positive selection in the C. annuum clusters, possibly related to domestication or selection towards traits of interest. This work provides comprehensive and relevant information on the origin and relationships of Spanish landraces and for future association mapping studies in pepper.

Biotechnological and Digital Revolution for Climate-Smart Plant Breeding

Climate change, associated with global warming, extreme weather events, and increasing incidence of weeds, pests and pathogens, is strongly influencing major cropping systems. In this challenging scenario, miscellaneous strategies are needed to expedite the rate of genetic gains with the purpose of developing novel varieties. Large plant breeding populations, efficient high-throughput technologies, big data management tools, and downstream biotechnology and molecular techniques are the pillars on which next generation breeding is based. In this review, we describe the toolbox the breeder has to face the challenges imposed by climate change, remark on the key role bioinformatics plays in the analysis and interpretation of big “omics” data, and acknowledge all the benefits that have been introduced into breeding strategies with the biotechnological and digital revolution.

Large Scale Phenotyping Provides Insight into the Diversity of Vegetative and Reproductive Organs in a Wide Collection of Wild and Domesticated Peppers (Capsicum spp.)

In the past years, the diversity of Capsicum has been mainly investigated through genetics and genomics approaches, fewer efforts have been made in the field of plant phenomics. Assessment of crop traits with high-throughput methodologies could enhance the knowledge of the plant phenome, giving at the same time a key contribution to the understanding of the function of many genes. In this study, a wide germplasm collection of 307 accessions retrieved from 48 world regions, and belonging to nine Capsicum species was characterized for 54 plant, leaf, flower and fruit traits. Conventional descriptors and semi-automated tools based on image analysis and colour coordinate detection were used. Significant differences were found among accessions, between species and between sweet and spicy cultivated types, revealing a large diversity. The results highlighted how the domestication process and the continued selection have increased the variability of fruit shape and colour. Hierarchical clustering based on conventional and fruit morphological descriptors reflected the separation of species on the basis of their phylogenetic relationships. These observations suggested that the flow between distinct gene pools could have contributed to determine the similarity of the species on the basis of morphological plant and fruit parameters. The approach used represents the first high-throughput phenotyping effort in Capsicum spp. aimed at broadening the knowledge of the diversity of domesticated and wild peppers. The data could help to select best the candidates for breeding and provide new insight into the understanding of the genetic base of the fruit shape of pepper.

GBS-derived SNP catalogue unveiled wide genetic variability and geographical relationships of Italian olive cultivars

Information on the distribution of genetic variation is essential to preserve olive germplasm from erosion and to recover alleles lost through selective breeding. In addition, knowledge on population structure and genotype-phenotype associations is crucial to support modern olive breeding programs that must respond to new environmental conditions imposed by climate change and novel biotic/abiotic stressors. To further our understanding of genetic variation in the olive, we performed genotype-by-sequencing on a panel of 94 Italian olive cultivars. A reference-based and a reference-independent SNP calling pipeline generated 22,088 and 8,088 high-quality SNPs, respectively. Both datasets were used to model population structure via parametric and non parametric clustering. Although the two pipelines yielded a 3-fold difference in the number of SNPs, both described wide genetic variability among our study panel and allowed individuals to be grouped based on fruit weight and the geographical area of cultivation. Multidimensional scaling analysis on identity-by-state allele-sharing values as well as inference of population mixtures from genome-wide allele frequency data corroborated the clustering pattern we observed. These findings allowed us to formulate hypotheses about geographical relationships of Italian olive cultivars and to confirm known and uncover novel cases of synonymy.

Molecular Markers Associated to Two Non-allelic Genic Male Sterility Genes in Peppers (Capsicum annuum L.)

Male sterility is of high importance in hybrid seed production of hot and sweet peppers. Genic (or nuclear) male sterility (GMS) is a simply inherited (usually monogenic recessive) and highly stable trait. However, one major disadvantage of using GMS is 1:1 segregation of male sterile to male fertile plants in every subsequent generation. Molecular markers tightly linked to genic male sterility (ms) genes would facilitate an efficient and rapid transfer of ms genes into different genetic backgrounds through marker-assisted backcrossing. The two non-allelic genic male sterility genes ms3 and ms w in hot and sweet pepper backgrounds, respectively, are monogenic recessive. Genotyping by sequencing (GBS) in an F2 population segregating for ms3 gene in hot pepper and in an F6 inbred near-isogenic line (NIL) population segregating for ms w gene in sweet pepper yielded 9,713 and 7,453 single nucleotide polymorphism markers, respectively. Four candidate SNPs co-segregating with ms3 gene and one co-segregating with ms w gene were identified by bulk segregant analysis and physically mapped to chromosomes 1 and 5, respectively. In hot pepper, two markers [HPGMS2 (CAPS) and HPGMS3 (dCAPS)] located 3.83 cM away from the ms3 gene and in sweet pepper the dCAPS marker SPGMS1 co-segregated (completely linked) with the ms w gene were developed. These markers will increase the efficacy of the male sterility genes for pepper breeding, as they can be useful in developing the genic male sterile lines in parental inbred lines of commercial hybrids through marker-assisted backcrossing, hybrid seed production, and genetic purity testing of hybrid seeds.

Population genomics and morphometric assignment of western honey bees (Apis mellifera L.) in the Republic of South Africa

BACKGROUNDS

Apis mellifera scutellata and A.m. capensis (the Cape honey bee) are western honey bee subspecies indigenous to the Republic of South Africa (RSA). Both bees are important for biological and economic reasons. First, A.m. scutellata is the invasive "African honey bee" of the Americas and exhibits a number of traits that beekeepers consider undesirable. They swarm excessively, are prone to absconding (vacating the nest entirely), usurp other honey bee colonies, and exhibit heightened defensiveness. Second, Cape honey bees are socially parasitic bees; the workers can reproduce thelytokously. Both bees are indistinguishable visually. Therefore, we employed Genotyping-by-Sequencing (GBS), wing geometry and standard morphometric approaches to assess the genetic diversity and population structure of these bees to search for diagnostic markers that can be employed to distinguish between the two subspecies.

RESULTS

Apis mellifera scutellata possessed the highest mean number of polymorphic SNPs (among 2449 informative SNPs) with minor allele frequencies > 0.05 (Np = 88%). The RSA honey bees generated a high level of expected heterozygosity (Hexp = 0.24). The mean genetic differentiation (FST; 6.5%) among the RSA honey bees revealed that approximately 93% of the genetic variation was accounted for within individuals of these subspecies. Two genetically distinct clusters (K = 2) corresponding to both subspecies were detected by Model-based Bayesian clustering and supported by Principal Coordinates Analysis (PCoA) inferences. Selected highly divergent loci (n = 83) further reinforced a distinctive clustering of two subspecies across geographical origins, accounting for approximately 83% of the total variation in the PCoA plot. The significant correlation of allele frequencies at divergent loci with environmental variables suggested that these populations are adapted to local conditions. Only 17 of 48 wing geometry and standard morphometric parameters were useful for clustering A.m. capensis, A.m. scutellata, and hybrid individuals.

CONCLUSIONS

We produced a minimal set of 83 SNP loci and 17 wing geometry and standard morphometric parameters useful for identifying the two RSA honey bee subspecies by genotype and phenotype. We found that genes involved in neurology/behavior and development/growth are the most prominent heritable traits evolved in the functional evolution of honey bee populations in RSA. These findings provide a starting point for understanding the functional basis of morphological differentiations and ecological adaptations of the two honey bee subspecies in RSA.

Occurrence of variable levels of health-promoting fruit compounds in horn-shaped Italian sweet pepper varieties assessed by a comprehensive approach

BACKGROUND

Cultivated pepper is a rich source of diverse bioactive compounds with potential properties related to the prevention of major degenerative and proliferative diseases. In the present study, 15 sweet pepper varieties, highly appreciated on the Italian market, were assessed for variation in the content of ascorbic acid, isoprenoids and flavonoids, as well as morpho-agronomic performances and molecular diversity.

RESULTS

The collection under study showed a wide variability for all traits considered. Traditional cultivars, although less productive, were characterized by a high content of β-carotene and ascorbic acid, reaching maximum concentrations of 230.5 mg kg^-1 fresh weight (fw) and 2750 mg kg^-1 fw, respectively. Strong correlations were demonstrated between neoxanthin and luteolins. Fruit weight was positively correlated with α-tocopherol content and negatively correlated with luteolins content. The genotyping by sequencing platform allowed the identification of 1833 single nucleotide polymorphism, which better defined the relationships among cultivars, based on provenance and improvement rate.

CONCLUSIONS

The present study provides an overview of the variability in the expression of fruit nutritional traits in a collection of horn-shaped pepper cultivars, integrating agronomic and molecular data. The impact for breeding and consumers is discussed. © 2017 Society of Chemical Industry.

Population structure, genetic diversity and downy mildew resistance among Ocimum species germplasm

BACKGROUND

The basil (Ocimum spp.) genus maintains a rich diversity of phenotypes and aromatic volatiles through natural and artificial outcrossing. Characterization of population structure and genetic diversity among a representative sample of this genus is severely lacking. Absence of such information has slowed breeding efforts and the development of sweet basil (Ocimum basilicum L.) with resistance to the worldwide downy mildew epidemic, caused by the obligate oomycete Peronospora belbahrii. In an effort to improve classification of relationships 20 EST-SSR markers with species-level transferability were developed and used to resolve relationships among a diverse panel of 180 Ocimum spp. accessions with varying response to downy mildew.

RESULTS

Results obtained from nested Bayesian model-based clustering, analysis of molecular variance and unweighted pair group method using arithmetic average (UPGMA) analyses were synergized to provide an updated phylogeny of the Ocimum genus. Three (major) and seven (sub) population (cluster) models were identified and well-supported (P < 0.001) by PhiPT (Φ_PT) values of 0.433 and 0.344, respectively. Allelic frequency among clusters supported previously developed hypotheses of allopolyploid genome structure. Evidence of cryptic population structure was demonstrated for the k1 O. basilicum cluster suggesting prevalence of gene flow. UPGMA analysis provided best resolution for the 36-accession, DM resistant k3 cluster with consistently strong bootstrap support. Although the k3 cluster is a rich source of DM resistance introgression of resistance into the commercially important k1 accessions is impeded by reproductive barriers as demonstrated by multiple sterile F1 hybrids. The k2 cluster located between k1 and k3, represents a source of transferrable tolerance evidenced by fertile backcross progeny. The 90-accession k1 cluster was largely susceptible to downy mildew with accession 'MRI' representing the only source of DM resistance.

CONCLUSIONS

High levels of genetic diversity support the observed phenotypic diversity among Ocimum spp. accessions. EST-SSRs provided a robust evaluation of molecular diversity and can be used for additional studies to increase resolution of genetic relationships in the Ocimum genus. Elucidation of population structure and genetic relationships among Ocimum spp. germplasm provide the foundation for improved DM resistance breeding strategies and more rapid response to future disease outbreaks.

Resolving relationships and phylogeographic history of the Nyssa sylvatica complex using data from RAD-seq and species distribution modeling

Nyssa sylvatica complex consists of several woody taxa occurring in eastern North America. These taxa were recognized as two or three species including three or four varieties by different authors. Due to high morphological similarities and complexity of morphological variation, classification and delineation of taxa in the group have been difficult and controversial. Here we employ data from RAD-seq to elucidate the genetic structure and phylogenetic relationships within the group. Using the genetic evidence, we evaluate previous classifications and delineate species. We also employ Species Distribution Modeling (SDM) to evaluate impacts of climatic changes on the ranges of the taxa and to gain insights into the relevant refugia in eastern North America. Results from Molecular Variance Analysis (AMOVA), STRUCTURE, phylogenetic analyses using Maximum likelihood, Bayesian Inference, and Splittree methods of RAD-seq data strongly support a two-clade pattern, largely separating samples of N. sylvatica from those of N. biflora-N. ursina mix. Divergence time analysis with BEAST suggests the two clades diverged in the mid Miocene. The ancestor of the present trees of N. sylvatica was suggested to be in the Pliocene and that of N. biflora-N. ursina mix in the end of the Miocene. Results from SDM predicted a smaller range in the southern part of the species present range of each clade during the Last Glacial Maximum (LGM). A northward expansion of the ranges during interglacial period and a northward shift of the ranges in the future under a model of global warming were also predicted. Our results support the recognition of two species in the complex, N. sylvatica and N. biflora, following the phylogenetic species concept. We found no genetic evidence supporting recognitions of intraspecific taxa. However, we propose subsp. ursina and subsp. biflora within N. biflora due to their distinction in habits, distributions, and habitats. Our results further support movements of trees in eastern North America in response to climatic changes. Finally, our study demonstrates that RAD-seq data and a combination of population genomics and SDM are valuable in resolving relationship and biogeographic history of closely related species that are taxonomically difficult.