Salinity is a rising concern in many lettuce-growing regions. Lettuce (Lactuca sativa L.) is sensitive to salinity, which reduces plant biomass, and causes leaf burn and early senescence. We sought to identify physiological traits important in salt tolerance that allows lettuce adaptation to high salinity while maintaining its productivity. Based on previous salinity tolerance studies, one sensitive and one tolerant genotype each was selected from crisphead, butterhead, and romaine, as well as leaf types of cultivated lettuce and its wild relative, L. serriola L. Physiological parameters were measured four weeks after transplanting two-day old seedlings into 350 mL volume pots filled with sand, hydrated with Hoagland nutrient solution and grown in a growth chamber. Salinity treatment consisted of gradually increasing concentrations of NaCl and CaCl₂ from 0 mM/0 mM at the time of transplanting, to 30 mM/15 mM at the beginning of week three, and maintaining it until harvest. Across the 10 genotypes, leaf area and fresh weight decreased 0–64% and 16–67%, respectively, under salinity compared to the control. Salinity stress increased the chlorophyll index by 4–26% in the cultivated genotypes, while decreasing it by 5–14% in the two wild accessions. Tolerant lines less affected by elevated salinity were characterized by high values of the chlorophyll fluorescence parameters Fv/Fm and instantaneous photosystem II quantum yield (QY), and lower leaf transpiration.

The C-repeat binding factor (CBF)/dehydration-responsive element binding (DREB1) proteins play a prominent role in freezing tolerance and are highly conserved in higher plants. Here we performed a genome-wide search of the CBF/DREB1 gene family in lettuce (Lactuca sativa L.) and identified 14 members of the family with one member gene containing a non-sense mutation within the AP2 DNA-binding domain. A comprehensive phylogenetic analysis of the CBF/DREB1 family members in 20 plant species from the Asterid or Rosid clade provided evidence that tandem duplication played an important role in the expansion of the CBF/DREB1 family. Expression analysis showed that twelve of the lettuce CBF genes were responsive to low temperature (4 °C), and that three and six of them could also be responsive to salt and heat stresses, respectively. Unlike Arabidopsis thaliana whose members of the CBF/DREB1 family respond only to a particular stress, lettuce CBFs provide wider protection from combinations of abiotic stresses. A global transcriptome analysis revealed distinctive temporal expression patterns among the cold-regulated genes in lettuce plants exposed to low temperature. Genes induced throughout the cold treatment are enriched in functions associated with protection from UV and high-light intensity and the genes suppressed after 7 days of cold exposure are enriched in photosynthesis-associated functions. These results provide insight into the molecular evolutionary properties of the CBF/DREB1 gene family in lettuce and a reference for genetic improvement of the lettuce response to cold acclimation.

White rust, caused by Albugo occidentalis, is one of the major yield-limiting diseases of spinach (Spinacia oleracea) in some major commercial production areas, particularly in southern Texas in the United States. The use of host resistance is the most economical and environment-friendly approach to managing white rust in spinach production. The objectives of this study were to conduct a genome-wide associating study (GWAS), to identify single nucleotide polymorphism (SNP) markers associated with white rust resistance in spinach, and to perform genomic prediction (GP) to estimate the prediction accuracy (PA). A GWAS panel of 346 USDA (US Dept. of Agriculture) germplasm accessions was phenotyped for white rust resistance under field conditions and GWAS was performed using 13 235 whole-genome resequencing (WGR) generated SNPs. Nine SNPs, chr2_53 049 132, chr3_58 479 501, chr3_95 114 909, chr4_9 176 069, chr4_17 807 168, chr4_83 938 338, chr4_87 601 768, chr6_1 877 096, and chr6_31 287 118, located on chromosomes 2, 3, 4, and 6 were associated with white rust resistance in this GWAS panel. Four scenarios were tested for PA using Pearson’s correlation coefficient (r) between the genomic estimation breeding value (GEBV) and the observed values: (1) different ratios between the training set and testing set (fold), (2) different GP models, (3) different SNP numbers in three different SNP sets, and (4) the use of GWAS-derived significant SNP markers. The results indicated that a 2- to 10-fold difference in the various GP models had similar, although not identical, averaged r values in each SNP set; using GWAS-derived significant SNP markers would increase PA with a high r-value up to 0.84. The SNP markers and the high PA can provide valuable information for breeders to improve spinach by marker-assisted selection (MAS) and genomic selection (GS).

Downy mildew, caused by the oomycete Peronospora effusa, is the most economically important disease on spinach. Fourteen new races of P. effusa have been identified in the last three decades. The frequent emergence of new races of P. effusa continually overcome the genetic resistance to the pathogen. The objectives of this research were to more clearly map the downy mildew resistance locus RPF1 in spinach, to identify single nucleotide polymorphism (SNP) markers associated with the resistance, and to refine the candidate genes responsible for the resistance. Progeny from populations generated from crosses of cultivars resistant (due to RPF1) to race 13 of P. effusa (Swan, T-Bird, Squirrel, and Tonga) with race 13 susceptible cultivars (Whale and Polka) were inoculated and the downy mildew disease response determined. Association analysis was performed in TASSEL, GAPIT, PLINK, and GENESIS programs using SNP markers identified from genotyping by sequencing (GBS). Association analysis mapped the race 13 resistance loci (RPF1) to positions 0.39, 0.69, 0.94-0.98, and 1.2 Mb of chromosome 3. The associated SNPs were within 1–7 kb of the disease resistance genes Spo12784, Spo12719, Spo12905, and Spo12821, and 11–18 Kb from Spo12903. This study extended our understanding of the genetic basis of downy mildew resistance in spinach and provided the most promising candidate genes Spo12784 and Spo12903 near the RPF1 locus, to pursue functional validation. The SNP markers may be used to select for the resistant lines to improve genetic resistance against the downy mildew pathogen and in developing durably resistant cultivars.

Genetic diversity is an important resource in crop breeding to improve cultivars with desirable traits. Selective breeding can lead to a reduction of genetic diversity. However, our understanding on this subject remains limited in lettuce ( Lactuca sativa L.). Genotyping‐by‐sequencing (GBS) can provide a reduced version of the genome as a cost‐effective method to identify genetic variants across the genome. We genotyped a diverse set of 441 lettuce accessions using the GBS method. Phylogenetic and population genetic analyses indicated substantial genetic divergence among four horticultural types of lettuce: butterhead, crisphead, leaf, and romaine. Genetic‐diversity estimates between and within the four types indicated that the crisphead type was the most differentiated from other types, whereas its population was the most homogenous with the slowest linkage disequilibrium (LD) decay among the four types. These results suggested that crisphead lettuces had relatively less genetic variation across the genome as well as low gene flow from other types. We identified putative selective sweep regions that showed low genetic variation in the crisphead type. Genome‐wide association study (GWAS) and quantitative trait loci (QTL) analyses provided evidence that these genomic regions were, in part, associated with delayed bolting, implicating the positive selection of delayed bolting in reducing variation. Our findings enhance the current understanding of genetic diversity and the impacts of selective breeding on patterning genetic variation in lettuce.

Deep understanding of genetic architecture of water-stress tolerance is critical for efficient and optimal development of water-stress tolerant cultivars, which is the most economical and environmentally sound approach to maintain lettuce production with limited irrigation. Lettuce (Lactuca sativa L.) production in areas with limited precipitation relies heavily on the use of ground water for irrigation. Lettuce plants are highly susceptible to water-stress, which also affects their nutrient uptake efficiency. Water stressed plants show reduced growth, lower biomass, and early bolting and flowering resulting in bitter flavors. Traditional phenotyping methods to evaluate water-stress are labor intensive, time-consuming and prone to errors. High throughput phenotyping platforms using kinetic chlorophyll fluorescence and hyperspectral imaging can effectively attain physiological traits related to photosynthesis and secondary metabolites that can enhance breeding efficiency for water-stress tolerance. Kinetic chlorophyll fluorescence and hyperspectral imaging along with traditional horticultural traits identified genomic loci affected by water-stress. Supervised machine learning models were evaluated for their accuracy to distinguish water-stressed plants and to identify the most important water-stress related parameters in lettuce. Random Forest (RF) had classification accuracy of 89.7% using kinetic chlorophyll fluorescence parameters and Neural Network (NN) had classification accuracy of 89.8% using hyperspectral imaging derived vegetation indices. The top ten chlorophyll fluorescence parameters and vegetation indices selected by sequential forward selection by RF and NN were genetically mapped using a L. sativa × L. serriola interspecific recombinant inbred line (RIL) population. A total of 25 quantitative trait loci (QTL) segregating for water-stress related horticultural traits, 26 QTL for the chlorophyll fluorescence traits and 34 QTL for spectral vegetation indices (VI) were identified. The percent phenotypic variation (PV) explained by the horticultural QTL ranged from 6.41 to 19.5%, PV explained by chlorophyll fluorescence QTL ranged from 6.93 to 13.26% while the PV explained by the VI QTL ranged from 7.2 to 17.19%. Eight QTL clusters harboring co-localized QTL for horticultural traits, chlorophyll fluorescence parameters and VI were identified on six lettuce chromosomes. Molecular markers linked to the mapped QTL clusters can be targeted for marker-assisted selection to develop water-stress tolerant lettuce.