Genome-Wide Association Analysis of Resistance to Pseudoperonospora cubensis in Citron Watermelon

Genome-Wide Association Mapping and Genomic Prediction of Fusarium Wilt Race 2 Resistance in the USDA Citrullus amarus Collection

Fusarium wilt caused by Fusarium oxysporum f. sp. niveum (Fon) race 2 is a serious disease in watermelon and can reduce yields by 80%. Genome-wide association studies (GWAS) are a valuable tool in dissecting the genetic basis of traits. Citrullus amarus accessions (n = 120) from the USDA germplasm collection were genotyped with whole-genome resequencing, resulting in 2,126,759 single nucleotide polymorphic (SNP) markers that were utilized for GWAS. Three models were used for GWAS with the R package GAPIT. Mixed linear model (MLM) analysis did not identify any significant marker associations. FarmCPU identified four quantitative trait nucleotides (QTN) on three different chromosomes (i.e., chromosomes 1, 5, and 9), and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK) identified one QTN on chromosome 10 as significantly associated with Fon race 2 resistance. FarmCPU identified four QTN that explained 60% of Fon race 2 resistance, and the single QTN from BLINK explained 27%. Relevant candidate genes were found within the linkage disequilibrium (LD) blocks of these significant SNPs, including genes encoding aquaporins, expansins, 2S albumins, and glutathione S-transferases which have been shown to be involved in imparting resistance to Fusarium spp. Genomic predictions (GP) for Fon race 2 resistance using all 2,126,759 SNPs resulted in a mean prediction accuracy of 0.08 with five-fold cross-validation employing genomic best linear unbiased prediction (gBLUP) or ridge-regression best linear unbiased prediction (rrBLUP). Mean prediction accuracy with gBLUP leave-one-out cross-validation was 0.48. Thus, along with identifying genomic regions associated with Fon race 2 resistance among the accessions, this study observed prediction accuracies that were strongly influenced by population size.

Genome-Wide Association Study of Resistance to Pseudomonas syringae in the USDA Collection of Citrullus amarus

Pseudomonas leaf spot (PLS), caused by Pseudomonas syringae pv. syringae, is an emerging disease of watermelon in the United States with the potential to severely reduce yield under humid conditions. The genetic basis of resistance to this disease is not known and no resistant germplasm is available. Because Citrullus amarus is an important reservoir of resistance genes for the cultivated watermelon, C. lanatus, we screened the United States Department of Agriculture plant introduction collection of C. amarus for resistance to PLS. Accessions (n = 117) were phenotyped for their level of resistance to PLS in two separate tests. Accession means of percent leaf area affected ranged from 1.5 to 99.4%. The broad-sense heritability for the trait was 0.51. Whole-genome resequencing generated 2,126,759 single-nucleotide polymorphisms (SNPs) which were used to perform a genome-wide association study (GWAS) aimed at discovering molecular markers for resistance. Three different models-BLINK, FarmCPU, and MLM-were included in the GWAS analyses. BLINK and FarmCPU, which are multilocus models, found eight SNPs, located on chromosomes Ca01, Ca05, Ca06, Ca08, and Ca10, that were significantly associated with resistance to PLS. Two of these SNPs were found by both BLINK and FarmCPU. The MLM model did not detect any significant associations. BLINK and FarmCPU estimated an explained phenotypic variance of 43.6 and 28.5%, respectively, for SNP S6_19327000 and 25.0 and 26.0%, respectively, for SNP S1_33362258, the two most significant SNPs found. In total, 43 candidate genes with known involvement in disease resistance were discovered within the genomic intervals of seven of the eight peak SNPs. Eleven of the candidate genes that were found have been reported to be involved in resistance to P. syringae in other plant species. Two significant SNPs were within resistance genes previously documented to play important roles of plant resistance specific to P. syringae in other pathosystems. The SNPs identified in this study will be instrumental in finding causal genes involved in PLS resistance in watermelon and developing resistant germplasm through breeding.

Genetic control of flowering time and fruit yield in citron watermelon

Flowering time and fruit yield are important traits in watermelon crop improvement. There is limited information on the inheritance and genomic loci underlying flowering time and yield performance, especially in citron watermelon. A total of 125 citron watermelon accessions were evaluated in field trials over two growing seasons for days to male and female flowers, fruit count, fruit weight, and fruit yield. The germplasm was genotyped with more than two million single-nucleotide polymorphism (SNP) markers generated via whole-genome resequencing. Trait mapping was conducted using a genome-wide association study (GWAS). Broad-sense heritability for all traits ranged from moderate to high, indicating that genetic improvement through breeding and selection is feasible. Significant marker-trait associations were uncovered for days to female flower (chromosomes Ca04, Ca05, Ca08, and Ca09), fruit count (on Ca02, Ca03, and Ca05), fruit weight (on Ca02, Ca06, Ca08, Ca10, and Ca11), and fruit yield on chromosomes Ca05, Ca07, and Ca09. The phenotypic variation explained by the significant SNPs ranged from 1.6 to 25.4, highlighting the complex genetic architecture of the evaluated traits. Candidate genes relevant to flowering time and fruit yield component traits were uncovered on chromosomes Ca02, Ca04, Ca05, Ca06, Ca09, and Ca11. These results lay a foundation for marker-assisted trait introgression of flowering time and fruit yield component traits in watermelons.

Distinct Genomic Loci Underlie Quantitative Resistance to Meloidogyne enterolobii Galling and Reproduction in Citrullus amarus

Meloidogyne enterolobii is a virulent species of root-knot nematode that threatens watermelon (Citrullus lanatus) production in the southeastern United States. There are no known sources of root-knot nematode resistance in cultivated C. lanatus. Specific genotypes of a wild watermelon relative, C. amarus, are resistant against M. incognita but the genetics that underly this resistance are still unknown and it is not clear that this same resistance will be effective against M. enterolobii. To identify and characterize new sources of resistance to M. enterolobii, we screened 108 diverse C. amarus lines alongside a susceptible C. lanatus cultivar (Charleston Gray) for resistance against M. enterolobii. Different C. amarus genotypes ranged from resistant to susceptible for the three resistance phenotypes measured; mean percent root system galled ranged from 10 to 73%, mean egg mass counts ranged from 0.3 to 64.5, and mean eggs per gram of root ranged from 326 to 146,160. We used each of these three resistance phenotypes combined with whole-genome resequencing data to conduct a genome-wide association scan that identified significant associations between M. enterolobii resistance and 11 single-nucleotide polymorphisms (SNPs) within the C. amarus genome. Interestingly, SNPs associated with reduced galling and egg masses were located within a single quantitative trait locus (QTL) on chromosome Ca03, while reductions in nematode eggs per gram of root were associated with separate QTL on chromosomes Ca04 and Ca08. The results of this study suggest that multiple genes are involved with M. enterolobii resistance in C. amarus and the SNPs identified will assist with efforts to breed for M. enterolobii resistance in watermelon.

Phenotypic Variability of Root System Architecture Traits for Drought Tolerance among Accessions of Citron Watermelon (Citrullus lanatus var. citroides (L.H. Bailey)

Citron watermelon (Citrullus lanatus var. citroides) is a drought-tolerant cucurbit crop widely grown in sub-Saharan Africa in arid and semi-arid environments. The species is a C3 xerophyte used for multiple purposes, including intercropping with maize, and has a deep taproot system. The deep taproot system plays a key role in the species’ adaptation to dry conditions. Understanding the root system development of this crop could be useful to identify traits for breeding water-use efficient and drought-tolerant varieties. This study compared the root system architecture of citron watermelon accessions under water-stress conditions. Nine selected and drought-tolerant citron watermelon accessions were grown under non-stress (NS) and water stress (WS) conditions using the root rhizotron procedure in a glasshouse. The following root system architecture (RSA) traits were measured: root system width (RSW), root system depth (RSD), convex hull area (CHA), total root length (TRL), root branch count (RBC), total root volume (TRV), leaf area (LA), leaf number (LN), first seminal root length (FSRL), seminal root angle (SRA), root dry mass (RDM), shoot dry mass (SDM), root−shoot mass ratio (RSM), root mass ratio (RMR), shoot mass ratio (SMR) and root tissue density (RTD). The data collected on RSA traits were subjected to an analysis of variance (ANOVA), correlation and principal component analyses. ANOVA revealed a significant (p < 0.05) accession × water stress interaction effect for studied RSA traits. Under WS, RDM exhibited significant and positive correlations with RSM (r = 0.65), RMR (r = 0.66), RSD (r = 0.66), TRL (r = 0.60), RBC (r = 0.72), FSRL (r = 0.73) and LN (r = 0.70). The principal component analysis revealed high loading scores for the following RSA traits: RSW (0.89), RSD (0.97), TRL (0.99), TRV (0.90), TRL (0.99), RMR (0.96) and RDM (0.76). In conclusion, the study has shown that the identified RSA traits could be useful in crop improvement programmes for citron watermelon genotypes with enhanced drought adaptation for improved yield performance under drought-prone environments.