Search

Phenotype diversity within cultivated Capsicum chinense is particularly evident for fruit shape and size. We used this diversity in C. chinense to further unravel the genetic mechanisms underlying fruit shape variation in pepper and related Solanaceous species. We identified candidate genes for C. chinense fruit shape, explored their contribution to population structure, and characterized their potential function in pepper fruit shape. Using genotyping by sequencing, we identified 43,081 single nucleotide polymorphisms (SNPs) from diverse collections of C. chinense. Principal component, neighbor-joining tree, and population structure analyses resolved 3 phylogenetically robust clusters associated with fruit shapes. Genome-wide association study (GWAS) was used to identify associated genomic regions with various fruit shape traits obtained from image analysis with Tomato Analyzer software. In our GWAS, we selected 12 SNPs associated with locule number trait and 8 SNP markers associated with other fruit shape traits such as perimeter, area, obovoid, ellipsoid and morphometrics (5y, 6y and 7y). The SNPs in CLAVATA1, WD-40, Auxin receptor, AAA type ATPase family protein, and RNA polymerase III genes were the major markers identified for fruit locule number from our GWAS results. Furthermore, we found SNPs in tetratricopeptide-repeat thioredoxin-like 3, enhancer of ABA co-receptor 1, subunit of exocyst complex 8 and pleiotropic drug resistance proteins associated with various fruit shape traits. CLAVATA1, WD-40 and Auxin receptor genes are known genes that affect tomato fruit shape. In this study, we used Arabidopsis thaliana T-DNA insertion knockout mutants and expression profiles for functional characterization of newly identified genes and to understand their role in fruit shape.

Background Watermelon seeds are a powerhouse of value-added traits such as proteins, free amino acids, vitamins, and essential minerals, offering a paleo-friendly dietary option. Despite the availability of substantial genetic variation, there is no sufficient information on the natural variation in seed-bound amino acids or proteins across the watermelon germplasm. This study aimed to analyze the natural variation in watermelon seed amino acids and total protein and explore underpinning genetic loci by genome-wide association study (GWAS). Methods The study evaluated the distribution of seed-bound free amino acids and total protein in 211 watermelon accessions of Citrullus spp, including 154 of Citrullus lanatus, 54 of Citrullus mucosospermus (egusi) and three of Citrullus amarus. We used the GWAS approach to associate seed phenotypes with 11,456 single nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing (GBS). Results Our results demonstrate a significant natural variation in different free amino acids and total protein content across accessions and geographic regions. The accessions with high protein content and proportion of essential amino acids warrant its use for value-added benefits in the food and feed industries via biofortification. The GWAS analysis identified 188 SNPs coinciding with 167 candidate genes associated with watermelon seed-bound amino acids and total protein. Clustering of SNPs associated with individual amino acids found by principal component analysis was independent of the speciation or cultivar groups and was not selected during the domestication of sweet watermelon. The identified candidate genes were involved in metabolic pathways associated with amino acid metabolism, such as Argininosuccinate synthase, explaining 7% of the variation in arginine content, which validate their functional relevance and potential for marker-assisted analysis selection. This study provides a platform for exploring potential gene loci involved in seed-bound amino acids metabolism, useful in genetic analysis and development of watermelon varieties with superior seed nutritional values.

MicroRNAs (miRNAs) are regulators of the post-transcription stage of gene activity documented to play central roles in flower and fruit development in model plant species. However, little is known about their roles and differences in domesticated and wild Capsicum species. In this study, we used high-throughput sequencing to analyze the miRNA content at three developmental stages (flower, small fruit, and middle fruit) from two cultivated (C. baccatum and C. annuum) and two wild (C. chacoense and C. eximium) pepper species. This analysis revealed 22 known and 27 novel miRNAs differentially expressed across species and tissues. A number of stage- and species-specific miRNAs were identified, and Gene Ontology terms were assigned to 138 genes targeted by the miRNAs. Most Gene Ontology terms were for the categories “genetic information processing”, “signaling and cellular processes”, “amino acid metabolism”, and “carbohydrate metabolism”. Enriched KEGG analysis revealed the pathways amino acids, sugar and nucleotide metabolism, starch and sucrose metabolism, and fructose-mannose metabolism among the principal ones regulated by miRNAs during pepper fruit ripening. We predicted miRNA–target gene interactions regulating flowering time and fruit development, including miR156/157 with SPL genes, miR159 with GaMYB proteins, miR160 with ARF genes, miR172 with AP2-like transcription factors, and miR408 with CLAVATA1 gene across the different Capsicum species. In addition, novel miRNAs play an important role in regulating interactions potentially controlling plant pathogen defense and fruit quality via fructokinase, alpha-L-arabinofuranosidase, and aromatic and neutral amino acid transporter. Overall, the small RNA-sequencing results from this study represent valuable information that provides a solid foundation for uncovering the miRNA-mediated mechanisms of flower and fruit development between domesticated and wild Capsicum species.

Acid mine drainage (AMD) is a huge environmental problem in mountain-top mining regions worldwide, including the Appalachian Mountains in the United States. This study applied a genome-wide association study (GWAS) to uncover genomic loci in Arabidopsis associated with tolerance to AMD toxicity. We characterized five major root phenotypes&mdash, cumulative root length, average root diameter, root surface area, root volume, and primary root length&mdash, in 180 Arabidopsis accessions in response to AMD-supplemented growth medium. GWAS of natural variation in the panel revealed genes associated with tolerance to an acidic environment. Most of these genes were transcription factors, anion/cation transporters, metal transporters, and unknown proteins. Two T-DNA insertion mutants, At1g63005 (miR399b) and At2g05635 (DEAD helicase RAD3), showed enhanced acidity tolerance. Our GWAS and the reverse genetic approach revealed genes involved in conferring tolerance to coal AMD. Our results indicated that proton resistance in hydroponic conditions could be an important index to improve plant growth in acidic soil, at least in acid-sensitive plant species.

The ankyrin (ANK) repeat protein family is largely distributed across plants and has been found to participate in multiple processes such as plant growth and development, hormone response, response to biotic and abiotic stresses. It is considered as one of the major markers of capsaicin content in pepper fruits. In this study, we performed a genome-wide identification and expression analysis of genes encoding ANK proteins in three Capsicum species: Capsicum baccatum, Capsicum annuum and Capsicum chinense. We identified a total of 87, 85 and 96 ANK genes in C. baccatum, C. annuum and C. chinense genomes, respectively. Next, we performed a comprehensive bioinformatics analysis of the Capsicum ANK gene family including gene chromosomal localization, Cis-elements, conserved motif identification, intron/exon structural patterns and gene ontology classification as well as profile expression. Phylogenetic and domain organization analysis grouped the Capsicum ANK gene family into ten subfamilies distributed across all 12 pepper chromosomes at different densities. Analysis of the expression of ANK genes in leaf and pepper fruits suggested that the ANKs have specific expression patterns at various developmental stages in placenta tissue. Our results provide valuable information for further studies of the evolution, classification and putative functions of ANK genes in pepper.

The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.

Elevated CO2 along with drought is a serious global threat to crop productivity. Therefore, understanding the molecular mechanisms plants use to protect these stresses is the key for plant growth and development. In this study, we mimicked natural stress conditions under a controlled Soil-Plant-Atmosphere-Research (SPAR) system and provided the evidence for how miRNAs regulate target genes under elevated CO2 and drought conditions. Significant physiological and biomass data supported the effective utilization of source-sink (leaf to root) under elevated CO2. Additionally, elevated CO2 partially rescued the effect of drought on total biomass. We identified both known and novel miRNAs differentially expressed during drought, CO2, and combined stress, along with putative targets. A total of 32 conserved miRNAs belonged to 23 miRNA families, and 25 novel miRNAs were identified by deep sequencing. Using the existing sweet potato genome database and stringent analyses, a total of 42 and 22 potential target genes were predicted for the conserved and novel miRNAs, respectively. These target genes are involved in drought response, hormone signaling, photosynthesis, carbon fixation, sucrose and starch metabolism, etc. Gene ontology and KEGG ontology functional enrichment revealed that these miRNAs might target transcription factors (MYB, TCP, NAC), hormone signaling regulators (ARF, AP2/ERF), cold and drought factors (corA), carbon metabolism (ATP synthase, fructose-1,6-bisphosphate), and photosynthesis (photosystem I and II complex units). Our study is the first report identifying targets of miRNAs under elevated CO2 levels and could support the molecular mechanisms under elevated CO2 in sweet potato and other crops in the future.

Watermelon is a good source of citrulline, a non-protein amino acid. Citrulline has several therapeutic and clinical implications as it produces nitric oxide via arginine. In plants, citrulline plays a pivotal role in nitrogen transport and osmoprotection. The purpose of this study was to identify single nucleotide polymorphism (SNP) markers associated with citrulline metabolism using a genome-wide association study (GWAS) and understand the role of citrulline in watermelon domestication. A watermelon collection consisting of 187 wild, landraces, and cultivated accessions was used to estimate citrulline content. An association analysis involved a total of 12,125 SNPs with a minor allele frequency (MAF)&thinsp, &gt, 0.05 in understanding the population structure and phylogeny in light of citrulline accumulation. Wild egusi types and landraces contained low to medium citrulline content, whereas cultivars had higher content, which suggests that obtaining higher content of citrulline is a domesticated trait. GWAS analysis identified candidate genes (ferrochelatase and acetolactate synthase) showing a significant association of SNPs with citrulline content. Haplotype networking indicated positive selection from wild to domesticated watermelon. To our knowledge, this is the first study showing genetic regulation of citrulline variation in plants by using a GWAS strategy. These results provide new insights into the citrulline metabolism in plants and the possibility of incorporating high citrulline as a trait in watermelon breeding programs.

Summary Years of selection for desirable fruit quality traits in dessert watermelon (Citrullus lanatus) has resulted in a narrow genetic base in modern cultivars. Development of novel genomic and genetic resources offers great potential to expand genetic diversity and improve important traits in watermelon. Here, we report a high‐quality genome sequence of watermelon cultivar ‘Charleston Gray’, a principal American dessert watermelon, to complement the existing reference genome from ‘97103’, an East Asian cultivar. Comparative analyses between genomes of ‘Charleston Gray’ and ‘97103’ revealed genomic variants that may underlie phenotypic differences between the two cultivars. We then genotyped 1365 watermelon plant introduction (PI) lines maintained at the U.S. National Plant Germplasm System using genotyping‐by‐sequencing (GBS). These PI lines were collected throughout the world and belong to three Citrullus species, C. lanatus, C. mucosospermus and C. amarus. Approximately 25 000 high‐quality single nucleotide polymorphisms (SNPs) were derived from the GBS data using the ‘Charleston Gray’ genome as the reference. Population genomic analyses using these SNPs discovered a close relationship between C. lanatus and C. mucosospermus and identified four major groups in these two species correlated to their geographic locations. Citrullus amarus was found to have a distinct genetic makeup compared to C. lanatus and C. mucosospermus. The SNPs also enabled identification of genomic regions associated with important fruit quality and disease resistance traits through genome‐wide association studies. The high‐quality ‘Charleston Gray’ genome and the genotyping data of this large collection of watermelon accessions provide valuable resources for facilitating watermelon research, breeding and improvement.

ATP-binding cassette (ABC) transporter genes act as transporters for different molecules across biological membranes and are involved in a diverse range of biological processes. In this study, we performed a genome-wide identification and expression analysis of genes encoding ABC transporter proteins in three Capsicum species, i.e., Capsicum annuum, Capsicum baccatum and Capsicum chinense. Capsicum is a valuable horticultural crop worldwide as an important constituent of many foods while containing several medicinal compounds including capsaicin and dihydrocapsaicin. Our results identified the presence of a total of 200, 185 and 187 ABC transporter genes in C. annuum, C. baccatum and C. chinense genomes, respectively. Capsaicin and dihydrocapsaicin content were determined in green pepper fruits (16 dpa). Additionally, we conducted different bioinformatics analyses including ABC genes classification, gene chromosomal location, Cis elements, conserved motifs identification and gene ontology classification, as well as profile expression of selected genes. Based on phylogenetic analysis and domain organization, the Capsicum ABC gene family was grouped into eight subfamilies. Among them, members within the ABCG, ABCB and ABCC subfamilies were the most abundant, while ABCD and ABCE subfamilies were less abundant throughout all species. ABC members within the same subfamily showed similar motif composition. Furthermore, common cis-elements involved in the transcriptional regulation were also identified in the promoter regions of all Capsicum ABC genes. Gene expression data from RNAseq and reverse transcription-semi-quantitative PCR analysis revealed development-specific stage expression profiles in placenta tissues. It suggests that ABC transporters, specifically the ABCC and ABCG subfamilies, may be playing important roles in the transport of secondary metabolites such as capsaicin and dihydrocapsaicin to the placenta vacuoles, effecting on their content in pepper fruits. Our results provide a more comprehensive understanding of ABC transporter gene family in different Capsicum species while allowing the identification of important candidate genes related to capsaicin content for subsequent functional validation.

The plant microbiome is a key determinant of plant health and productivity, and changes in the plant microbiome can alter the tolerance to biotic and abiotic stresses and the quality of end produce. Little is known about the microbial diversity and its effect on carbohydrate metabolism in ripe fruits. In this study, we aimed to understand the diversity and function of microorganisms in relation to carbohydrate metabolism of ripe watermelon fruits. We used 16S metagenomics and RNAseq metatranscriptomics for analysis of red (PI459074, Congo, and SDRose) and yellow fruit-flesh cultivars (PI227202, PI435990, and JBush) of geographically and metabolically diverse watermelon cultivars. Metagenomics data showed that Proteobacteria were abundant in SDRose and PI227202, whereas Cyanobacteria were most abundant in Congo and PI4559074. In the case of metatranscriptome data, Proteobacteria was the most abundant in all cultivars. High expression of genes linked to infectious diseases and the expression of peptidoglycan hydrolases associated to pathogenicity of eukaryotic hosts was observed in SDRose, which could have resulted in low microbial diversity in this cultivar. The presence of GH28, associated with polygalacturonase activity in JBush and SDRose could be related to cell wall modifications including de-esterification and depolymerization, and consequent loss of galacturonic acid and neutral sugars. Moreover, based on the KEGG annotation of the expressed genes, nine α-galactosidase genes involved in key processes of galactosyl oligosaccharide metabolism, such as raffinose family were identified and galactose metabolism pathway was reconstructed. Results of this study underline the links between the host and fruit-associated microbiome in carbohydrate metabolism of the ripe fruits. The cultivar difference in watermelon reflects the quantum and diversity of the microbiome, which would benefit watermelon and other plant breeders aiming at the holobiont concept to incorporate associated microbiomes in breeding programs.

Our genetic diversity study uses microsatellites of known map position to estimate genome level population structure and linkage disequilibrium, and to identify genomic regions that have undergone selection during watermelon domestication and improvement. Thirty regions that showed evidence of selective sweep were scanned for the presence of candidate genes using the watermelon genome browser (www.icugi.org). We localized selective sweeps in intergenic regions, close to the promoters, and within the exons and introns of various genes. This study provided an evidence of convergent evolution for the presence of diverse ecotypes with special reference to American and European ecotypes. Our search for location of linked markers in the whole-genome draft sequence revealed that BVWS00358, a GA repeat microsatellite, is the GAGA type transcription factor located in the 5′ untranslated regions of a structure and insertion element that expresses a Cys 2 His 2 Zinc finger motif, with presumed biological processes related to chitin response and transcriptional regulation. In addition, BVWS01708, an ATT repeat microsatellite, located in the promoter of a DTW domain-containing protein (Cla002761); and 2 other simple sequence repeats that association mapping link to fruit length and rind thickness.

The genome size for watermelon (Citrullus lanatus), melon (Cucumis melo), cucumber (Cuucmis sativus) and pumpkin (Cucurbita pepo) is 425, 454, 367 and 502 Mbp, respectively, and considered medium size as compared with most other crops. Whole-genome duplication is common in angiosperm plants. Research has revealed a paleohexaploidy (ϒ) event in the common ancestor of eudicots after the divergence of monocotyledons and dicotyledons. While analysis of published whole-genome sequences of cucumber, melon and watermelon showed traces of these ancient duplication events, there was not evidence of more recent whole-genome duplications in these species. Analysis of the syntenic relationships among watermelon, cucumber, melon and grape has identified 3543 orthologous relationships covering the watermelon, melon, and cucumber genomes. Comparison of melon and cucumber genomes synteny to detect shorter regions of rearrangements confirmed previously reported ancestral fusions of five melon chromosome pairs in cucumber, and several inter- and intra-chromosome rearrangements between the two species. Sequenced genomes of cucumber, melon and watermelon revealed a comparable range of genes from 23 k to 27 k protein coding genes with an average of four to six exons per gene. Current progress in gene mapping technologies such as optical mapping, which produces maps of large individual DNA molecules, can improve cucurbit comparative genomics to detect large structural variations and DNA rearrangements across the species.

We used genotyping by sequencing to identify a set of 10,480 single nucleotide polymorphism (SNP) markers for constructing a high-resolution genetic map of 1096 cM for watermelon. We assessed the genome-wide variation in recombination rate (GWRR) across the map and found an association between GWRR and genome-wide nucleotide diversity. Collinearity between the map and the genome-wide reference sequence for watermelon was studied to identify inconsistency and chromosome rearrangements. We assessed genome-wide nucleotide diversity, linkage disequilibrium (LD), and selective sweep for wild, semi-wild, and domesticated accessions of Citrullus lanatus var. lanatus to track signals of domestication. Principal component analysis combined with chromosome-wide phylogenetic study based on 1563 SNPs obtained after LD pruning with minor allele frequency of 0.05 resolved the differences between semi-wild and wild accessions as well as relationships among worldwide sweet watermelon. Population structure analysis revealed predominant ancestries for wild, semi-wild, and domesticated watermelons as well as admixture of various ancestries that were important for domestication. Sliding window analysis of Tajima’s D across various chromosomes was used to resolve selective sweep. LD decay was estimated for various chromosomes. We identified a strong selective sweep on chromosome 3 consisting of important genes that might have had a role in sweet watermelon domestication.

Watermelon (Citrullus lanatus var. lanatus) is one of the most important vegetable crops in the world. Molecular markers have become the tools of choice for resolving watermelon taxonomic relationships and evolution. Increased numbers of single nucleotide polymorphism (SNP) markers together with simple sequence repeat (SSR) markers would be useful for phylogenetic analyses of germplasm accessions and for linkage mapping for marker-assisted breeding with quantitative trait loci and single genes. We aimed to construct a genetic map based on SNPs (generated by Illumina Veracode multiplex assays for genotyping) and SSR markers and evaluate relationships inferred from SNP genotypes between 130 watermelon accessions collected throughout the world. We incorporated 282 markers (232 SNPs and 50 SSRs) into the linkage map. The genetic map consisted of 11 linkage groups spanning 924.72 cM with an average distance of 3.28 cM between markers. Because all of the SNP-containing sequences were assembled with the whole-genome sequence draft for watermelon, chromosome numbers could be readily assigned for all the linkage groups. We found that 134 SNPs were polymorphic in 130 watermelon accessions chosen for diversity studies. The current 384-plex SNP set is a powerful tool for characterizing genetic relatedness and for developing medium-resolution genetic maps.

Principal component analysis (PCA) with 36,621 polymorphic genome-anchored single nucleotide polymorphisms (SNPs) identified collectively for Capsicum annuum and Capsicum baccatum was used to characterize population structure and species domestication of these two important incompatible cultivated pepper species. Estimated mean nucleotide diversity (π) and Tajima's D across various chromosomes revealed biased distribution toward negative values on all chromosomes (except for chromosome 4) in cultivated C. baccatum, indicating a population bottleneck during domestication of C. baccatum. In contrast, C. annuum chromosomes showed positive π and Tajima's D on all chromosomes except chromosome 8, which may be because of domestication at multiple sites contributing to wider genetic diversity. For C. baccatum, 13,129 SNPs were available, with minor allele frequency (MAF) ≥0.05; PCA of the SNPs revealed 283 C. baccatum accessions grouped into 3 distinct clusters, for strong population structure. The fixation index (FST) between domesticated C. annuum and C. baccatum was 0.78, which indicates genome-wide divergence. We conducted extensive linkage disequilibrium (LD) analysis of C. baccatum var. pendulum cultivars on all adjacent SNP pairs within a chromosome to identify regions of high and low LD interspersed with a genome-wide average LD block size of 99.1 kb. We characterized 1742 haplotypes containing 4420 SNPs (range 9–2 SNPs per haplotype). Genome-wide association study (GWAS) of peduncle length, a trait that differentiates wild and domesticated C. baccatum types, revealed 36 significantly associated genome-wide SNPs. Population structure, identity by state (IBS) and LD patterns across the genome will be of potential use for future GWAS of economically important traits in C. baccatum peppers.

Use of 10,129 singleton SNPs of known genomic location in tetraploid cotton provided unique opportunities to characterize genome-wide diversity among 440 Gossypium hirsutum and 219 G. barbadense cultivars and landrace accessions of widespread origin. Using the SNPs distributed genome-wide, we examined genetic diversity, haplotype distribution and linkage disequilibrium patterns in the G. hirsutum and G. barbadense genomes to clarify population demographic history. Diversity and identity-by-state analyses have revealed little sharing of alleles between the two cultivated allotetraploid genomes, with a few exceptions that indicated sporadic gene flow. We found a high number of new alleles, representing increased nucleotide diversity, on chromosomes 1 and 2 in cultivated G. hirsutum as compared with low nucleotide diversity on these chromosomes in landrace G. hirsutum. In contrast, G. barbadense chromosomes showed negative Tajima’s D on several chromosomes for both cultivated and landrace types, which indicate that speciation of G. barbadense itself, might have occurred with relatively narrow genetic diversity. The presence of conserved linkage disequilibrium (LD) blocks and haplotypes between G. hirsutum and G. barbadense provides strong evidence for comparable patterns of evolution in their domestication processes. Our study illustrates the potential use of population genetic techniques to identify genomic regions for domestication.

This genetic diversity study aimed to estimate the population structure and explore the use of association mapping strategies to identify linked markers for bacterial resistance, growth and fruit quality in pomegranate collections from India. In total, 88 accessions including 37 cultivated types were investigated. A total of 112 alleles were amplified by use of 44 publicly available microsatellites for estimating molecular genetic diversity and population structure. Neighbor-joining analysis, model-based population structure and principal component analysis corroborated the genetic relationships among wild-type and cultivated pomegranate collections from India. Our study placed all 88 germplasm into four clusters. We identified a cultivated clade of pomegranates in close proximity to Daru types of wild-type pomegranates that grow naturally near the foothills of the Himalayas. Admixture analysis sorted various lineages of cultivated pomegranates to their respective ancestral forms. We identified four linked markers for fruit weight, titratable acidity and bacterial blight severity. PGCT001 was found associated with both fruit weight and bacterial blight, and the association with fruit weight during both seasons analyzed was significant after Bonferroni correction. This research demonstrates effectiveness of microsatellites to resolve population structure among the wild and cultivar collection of pomegranates and future use for association mapping studies.

Background A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise FST distributions among worldwide cultivated watermelons to track signals of domestication. Results We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise FST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model. Conclusions Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.

Pepper (Capsicum annuum L.) is an economically important crop with added nutritional value. Production of capsaicin is an important quantitative trait with high environmental variance, so the development of markers regulating capsaicinoid accumulation is important for pepper breeding programs. In this study, we performed association mapping at the gene level to identify single nucleotide polymorphisms (SNPs) associated with capsaicin pathway metabolites in a diverse Capsicum annuum collection during two seasons. The genes Pun1, CCR, KAS and HCT were sequenced and matched with the whole-genome sequence draft of pepper to identify SNP locations and for further characterization. The identified SNPs for each gene underwent candidate gene association mapping. Association mapping results revealed Pun1 as a key regulator of major metabolites in the capsaicin pathway mainly affecting capsaicinoids and precursors for acyl moieties of capsaicinoids. Six different SNPs in the promoter sequence of Pun1 were found associated with capsaicin in plants from both seasons. Our results support that CCR is an important control point for the flux of p-coumaric acid to specific biosynthesis pathways. KAS was found to regulate the major precursors for acyl moieties of capsaicinoids and may play a key role in capsaicinoid production. Candidate gene association mapping of Pun1 suggested that the accumulation of capsaicinoids depends on the expression of Pun1, as revealed by the most important associated SNPs found in the promoter region of Pun1.

Knowledge of population structure and linkage disequilibrium among the worldwide collections of peppers currently classified as hot, mild, sweet and ornamental types is indispensable for applying association mapping and genomic selection to improve pepper. The current study aimed to resolve the genetic diversity and relatedness of Capsicum annuum germplasm by use of simple sequence repeat (SSR) loci across all chromosomes in samples collected in 2011 and 2012. The physical distance covered by the entire set of SSRs used was 2,265.9 Mb from the 3.48-Gb hot-pepper genome size. The model-based program STRUCTURE was used to infer five clusters, which was further confirmed by classical molecular-genetic diversity analysis. Mean heterozygosity of various loci was estimated to be 0.15. Linkage disequilibrium (LD) was used to identify 17 LD blocks across various chromosomes with sizes from 0.154 Kb to 126.28 Mb. CAMS-142 of chromosome 1 was significantly associated with both capsaicin (CA) and dihydrocapsaicin (DCA) levels. Further, CAMS-142 was located in an LD block of 98.18 Mb. CAMS-142 amplified bands of 244, 268, 283 and 326 bp. Alleles 268 and 283 bp had positive effects on both CA and DCA levels, with an average R 2 of 12.15 % (CA) and 12.3 % (DCA). Eight markers from seven different chromosomes were significantly associated with fruit weight, contributing an average effect of 15 %. CAMS-199, HpmsE082 and CAMS-190 are the three major quantitative trait loci located on chromosomes 8, 9, and 10, respectively, and were associated with fruit weight in samples from both years of the study. This research demonstrates the effectiveness of using genome-wide SSR-based markers to assess features of LD and genetic diversity within C. annuum.