Leaf cuticular wax composition of a genetically diverse collection of lettuce (Lactuca sativa L.) cultivars evaluated under field conditions.

Cuticular waxes of plants impart tolerance to many forms of environmental stress and help shed dangerous human pathogens on edible plant parts. Although the chemical composition of waxes on a wide variety of important crops has been described, a detailed wax compositional analysis has yet to be reported for lettuce (Lactuca sativa L.), one of the most widely consumed vegetables. We present herein the leaf wax content and composition of 12 genetically diverse lettuce cultivars sampled across five time points during their vegetative growth phase in the field. Mean total leaf wax amounts across all cultivars varied little over 28 days of vegetative growth, except for a notable decrease in total waxes following a major precipitation event, presumably due to wax degradation from wind and rain. All lettuce cultivars were found to contain a unique wax composition highly enriched in 22- and 24-carbon length 1-alcohols (docosanol and tetracosanol, respectively). In our report, the dominance of these shorter chain length 1-alcohols as wax constituents represents a relatively rare phenotype in plants. The ecological significance of these dominant and relatively short 1-alcohols is still unknown. Although waxes have been a target for improvement of various crops, no such work has been reported for lettuce. This study lays the groundwork for future research that aims to integrate cuticular wax characteristics of field grown plants into the larger context of lettuce breeding and cultivar development.

PhytoOracle: Scalable, modular phenomics data processing pipelines.

As phenomics data volume and dimensionality increase due to advancements in sensor technology, there is an urgent need to develop and implement scalable data processing pipelines. Current phenomics data processing pipelines lack modularity, extensibility, and processing distribution across sensor modalities and phenotyping platforms. To address these challenges, we developed PhytoOracle (PO), a suite of modular, scalable pipelines for processing large volumes of field phenomics RGB, thermal, PSII chlorophyll fluorescence 2D images, and 3D point clouds. PhytoOracle aims to (i) improve data processing efficiency; (ii) provide an extensible, reproducible computing framework; and (iii) enable data fusion of multi-modal phenomics data. PhytoOracle integrates open-source distributed computing frameworks for parallel processing on high-performance computing, cloud, and local computing environments. Each pipeline component is available as a standalone container, providing transferability, extensibility, and reproducibility. The PO pipeline extracts and associates individual plant traits across sensor modalities and collection time points, representing a unique multi-system approach to addressing the genotype-phenotype gap. To date, PO supports lettuce and sorghum phenotypic trait extraction, with a goal of widening the range of supported species in the future. At the maximum number of cores tested in this study (1,024 cores), PO processing times were: 235 minutes for 9,270 RGB images (140.7 GB), 235 minutes for 9,270 thermal images (5.4 GB), and 13 minutes for 39,678 PSII images (86.2 GB). These processing times represent end-to-end processing, from raw data to fully processed numerical phenotypic trait data. Repeatability values of 0.39-0.95 (bounding area), 0.81-0.95 (axis-aligned bounding volume), 0.79-0.94 (oriented bounding volume), 0.83-0.95 (plant height), and 0.81-0.95 (number of points) were observed in Field Scanalyzer data. We also show the ability of PO to process drone data with a repeatability of 0.55-0.95 (bounding area).

Loss-of-function of SAWTOOTH 1 affects leaf dorsiventrality genes to promote leafy heads in lettuce.

The mechanisms underlying leafy heads in vegetables are poorly understood. Here, we cloned a quantitative trait locus (QTL) controlling leafy heads in lettuce (Lactuca sativa). The QTL encodes a transcription factor, SAWTOOTH 1 (LsSAW1), which has a BEL1-like homeodomain and is a homolog of Arabidopsis thaliana. A 1-bp deletion in Lssaw1 contributes to the development of leafy heads. Laser-capture microdissection and RNA-sequencing showed that LsSAW1 regulates leaf dorsiventrality and loss-of-function of Lssaw1 downregulates the expression of many adaxial genes but upregulates abaxial genes. LsSAW1 binds to the promoter region of the adaxial gene ASYMMETRIC LEAVES 1 (LsAS1) to upregulate its expression. Overexpression of LsAS1 compromised the effects of Lssaw1 on heading. LsSAW1 also binds to the promoter region of the abaxial gene YABBY 1 (LsYAB1), but downregulates its expression. Overexpression of LsYAB1 led to bending leaves in LsSAW1 genotypes. LsSAW1 directly interacts with KNOTTED 1 (LsKN1), which is necessary for leafy heads in lettuce. RNA-seq data showed that LsSAW1 and LsKN1 exert antagonistic effects on the expression of thousands of genes. LsSAW1 compromises the ability of LsKN1 to repress LsAS1. Our results suggest that downregulation or loss-of-function of adaxial genes and upregulation of abaxial genes allow for the development of leafy heads.

The upregulated LsKN1 gene transforms pinnately to palmately lobed leaves through auxin, gibberellin, and leaf dorsiventrality pathways in lettuce.

Leaf shape represents a vital agronomic trait for leafy vegetables such as lettuce. Some lettuce cultivars produce lobed leaves, varying from pinnately to palmately lobed, but the genetic mechanisms remain unclear. In this study, we cloned one major quantitative trait locus (QTL) controlling palmately lobed leaves. The candidate gene, LsKN1, encodes a homeobox transcription factor, and has been shown previously to be critical for the development of leafy heads in lettuce. The LsKN1 allele that is upregulated by the insertion of a transposon promotes the development of palmately lobed leaves. We demonstrated that LsKN1 upregulated LsCUC2 and LsCUC3 through different mechanisms, and their upregulation was critical for the development of palmately lobed leaves. LsKN1 binds the promoter of LsPID to promote auxin biosynthesis, which positively contributes to the development of palmately lobed leaves. In contrast, LsKN1 suppresses GA biosynthesis to promote palmately lobed leaves. LsKN1 also binds to the promoter of LsAS1, a dorsiventrality gene, to downregulate its expression. Overexpression of the LsAS1 gene compromised the effects of the LsKN1 gene changing palmately to pinnately lobed leaves. Our study illustrated that the upregulated LsKN1 gene led to palmately lobed leaves in lettuce by integrating several downstream pathways, including auxin, gibberellin, and leaf dorsiventrality pathways.

Quantitative Trait Loci and Candidate Genes Associated with Photoperiod Sensitivity in Lettuce (Lactuca spp.).

KEY MESSAGE: A population of lettuce that segregated for photoperiod sensitivity was planted under long-day and short-day conditions. Genetic mapping revealed two distinct sets of QTLs controlling daylength-independent and photoperiod-sensitive flowering time. The molecular mechanism of flowering time regulation in lettuce is of interest to both geneticists and breeders because of the extensive impact of this trait on agricultural production. Lettuce is a facultative long-day plant which changes in flowering time in response to photoperiod. Variations exist in both flowering time and the degree of photoperiod sensitivity among accessions of wild (Lactuca serriola) and cultivated (L. sativa) lettuce. An F6 population of 236 recombinant inbred lines (RILs) was previously developed from a cross between a late-flowering, photoperiod-sensitive L. serriola accession and an early-flowering, photoperiod-insensitive L. sativa accession. This population was planted under long-day (LD) and short-day (SD) conditions in a total of four field and screenhouse trials; the developmental phenotype was scored weekly in each trial. Using genotyping-by-sequencing (GBS) data of the RILs, quantitative trait loci (QTL) mapping revealed five flowering time QTLs that together explained more than 20% of the variation in flowering time under LD conditions. Using two independent statistical models to extract the photoperiod sensitivity phenotype from the LD and SD flowering time data, we identified an additional five QTLs that together explained more than 30% of the variation in photoperiod sensitivity in the population. Orthology and sequence analysis of genes within the nine QTLs revealed potential functional equivalents in the lettuce genome to the key regulators of flowering time and photoperiodism, FD and CONSTANS, respectively, in Arabidopsis.

A Composite Analysis of Flowering Time Regulation in Lettuce.

Plants undergo profound physiological changes when transitioning from vegetative to reproductive growth. These changes affect crop production, as in the case of leafy vegetables. Lettuce is one of the most valuable leafy vegetable crops in the world. Past genetic studies have identified multiple quantitative trait loci (QTLs) that affect the timing of the floral transition in lettuce. Extensive functional molecular studies in the model organism Arabidopsis provide the opportunity to transfer knowledge to lettuce to explore the mechanisms through which genetic variations translate into changes in flowering time. In this review, we integrated results from past genetic and molecular studies for flowering time in lettuce with orthology and functional inference from Arabidopsis. This summarizes the basis for all known genetic variation underlying the phenotypic diversity of flowering time in lettuce and how the genetics of flowering time in lettuce projects onto the established pathways controlling flowering time in plants. This comprehensive overview reveals patterns across experiments as well as areas in need of further study. Our review also represents a resource for developing cultivars with delayed flowering time.

Drone phenotyping and machine learning enable discovery of loci regulating daily floral opening in lettuce.

Flower opening and closure are traits of reproductive importance in all angiosperms because they determine the success of self- and cross-pollination. The temporal nature of this phenotype rendered it a difficult target for genetic studies. Cultivated and wild lettuce, Lactuca spp., have composite inflorescences that open only once. An L. serriola×L. sativa F6 recombinant inbred line (RIL) population differed markedly for daily floral opening time. This population was used to map the genetic determinants of this trait; the floral opening time of 236 RILs was scored using time-course image series obtained by drone-based phenotyping on two occasions. Floral pixels were identified from the images using a support vector machine with an accuracy >99%. A Bayesian inference method was developed to extract the peak floral opening time for individual genotypes from the time-stamped image data. Two independent quantitative trait loci (QTLs; Daily Floral Opening 2.1 and qDFO8.1) explaining >30% of the phenotypic variation in floral opening time were discovered. Candidate genes with non-synonymous polymorphisms in coding sequences were identified within the QTLs. This study demonstrates the power of combining remote sensing, machine learning, Bayesian statistics, and genome-wide marker data for studying the genetics of recalcitrant phenotypes.

Upregulation of a KN1 homolog by transposon insertion promotes leafy head development in lettuce.

Leafy head is a unique type of plant architecture found in some vegetable crops, with leaves bending inward to form a compact head. The genetic and molecular mechanisms underlying leafy head in vegetables remain poorly understood. We genetically fine-mapped and cloned a major quantitative trait locus controlling heading in lettuce. The candidate gene (LsKN1) is a homolog of knotted 1 (KN1) from Zea mays Complementation and CRISPR/Cas9 knockout experiments confirmed the role of LsKN1 in heading. In heading lettuce, there is a CACTA-like transposon inserted into the first exon of LsKN1 (LsKN1▽). The transposon sequences act as a promoter rather than an enhancer and drive high expression of LsKN1▽. The enhanced expression of LsKN1▽ is necessary but not sufficient for heading in lettuce. Data from ChIP-sequencing, electrophoretic mobility shift assays, and dual luciferase assays indicate that the LsKN1▽ protein binds the promoter of LsAS1 and down-regulates its expression to alter leaf dorsoventrality. This study provides insight into plant leaf development and will be useful for studies on heading in other vegetable crops.

Characterization of four polymorphic genes controlling red leaf colour in lettuce that have undergone disruptive selection since domestication.

Anthocyanins protect plants from biotic and abiotic stressors and provide great health benefits to consumers. In this study, we cloned four genes (Red Lettuce Leaves 1 to 4: RLL1 to RLL4) that contribute to colour variations in lettuce. The RLL1 gene encodes a bHLH transcription factor, and a 5-bp deletion in some cultivars abolishes its function to activate the anthocyanin biosynthesis pathway. The RLL2 gene encodes an R2R3-MYB transcription factor, which was derived from a duplication followed by mutations in its promoter region. The RLL3 gene encodes an R2-MYB transcription factor, which down-regulates anthocyanin biosynthesis through competing with RLL2 for interaction with RLL1; a mis-sense mutation compromises the capacity of RLL3 to bind RLL1. The RLL4 gene encodes a WD-40 transcription factor, homologous to the RUP genes suppressing the UV-B signal transduction pathway in Arabidopsis; a mis-sense mutation in rll4 attenuates its suppressing function, leading to a high concentration of anthocyanins. Sequence analysis of the RLL1-RLL4 genes from wild and cultivated lettuce showed that their function-changing mutations occurred after domestication. The mutations in rll1 disrupt anthocyanin biosynthesis, while the mutations in RLL2, rll3 and rll4 activate anthocyanin biosynthesis, showing disruptive selection for leaf colour during domestication of lettuce. The characterization of multiple polymorphic genes in this study provides the necessary molecular resources for the rational breeding of lettuce cultivars with distinct levels of red pigments and green cultivars with high levels of health-promoting flavonoids.

Correction to: The LsVe1L allele provides a molecular marker for resistance to Verticillium dahliae race 1 in lettuce.

Following publication of the original article [1], the author reported a processing error in Figure 5. This has been corrected in the original article.

The alternative reality of plant mitochondrial DNA: One ring does not rule them all.

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms.

The LsVe1L allele provides a molecular marker for resistance to Verticillium dahliae race 1 in lettuce.

BACKGROUND: Verticillium wilt caused by the fungus Verticillium dahliae race 1 is among the top disease concerns for lettuce in the Salinas and Pajaro Valleys of coastal central California. Resistance of lettuce against V. dahliae race 1 was previously mapped to the single dominant Verticillium resistance 1 (Vr1) locus. Lines of tomato resistant to race 1 are known to contain the closely linked Ve1 and Ve2 genes that encode receptor-like proteins with extracellular leucine-rich repeats; the Ve1 and Ve2 proteins act antagonistically to provide resistance against V. dahliae race 1. The Vr1 locus in lettuce contains a cluster of several genes with sequence similarity to the tomato Ve genes. We used genome sequencing and/or PCR screening along with pathogenicity assays of 152 accessions of lettuce to investigate allelic diversity and its relationship to race 1 resistance in lettuce. RESULTS: This approach identified a total of four Ve genes: LsVe1, LsVe2, LsVe3, and LsVe4. The majority of accessions, however, contained a combination of only three of these LsVe genes clustered on chromosomal linkage group 9 (within ~ 25 kb in the resistant cultivar La Brillante and within ~ 127 kb in the susceptible cultivar Salinas). CONCLUSIONS: A single allele, LsVe1L, was present in all resistant accessions and absent in all susceptible accessions. This allele can be used as a molecular marker for V. dahliae race 1 resistance in lettuce. A PCR assay for rapid detection of race 1 resistance in lettuce was designed based on nucleotide polymorphisms. Application of this assay allows identification of resistant genotypes in early stages of plant development or at seed-level without time- and labor-intensive testing in the field.

RNA sequencing provides insights into the evolution of lettuce and the regulation of flavonoid biosynthesis.

Different horticultural types of lettuce exhibit tremendous morphological variation. However, the molecular basis for domestication and divergence among the different horticultural types of lettuce remains unknown. Here, we report the RNA sequencing of 240 lettuce accessions sampled from the major horticultural types and wild relatives, generating 1.1 million single-nucleotide polymorphisms (SNPs). Demographic modeling indicates that there was a single domestication event for lettuce. We identify a list of regions as putative selective sweeps that occurred during domestication and divergence, respectively. Genome-wide association studies (GWAS) identify 5311 expression quantitative trait loci (eQTL) regulating the expression of 4105 genes, including nine eQTLs regulating genes associated with flavonoid biosynthesis. GWAS for leaf color detects six candidate loci responsible for the variation of anthocyanins in lettuce leaves. Our study provides a comprehensive understanding of the domestication and the accumulation of anthocyanins in lettuce, and will facilitate the breeding of cultivars with improved nutritional value.

Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce.

Lettuce (Lactuca sativa) is a major crop and a member of the large, highly successful Compositae family of flowering plants. Here we present a reference assembly for the species and family. This was generated using whole-genome shotgun Illumina reads plus in vitro proximity ligation data to create large superscaffolds; it was validated genetically and superscaffolds were oriented in genetic bins ordered along nine chromosomal pseudomolecules. We identify several genomic features that may have contributed to the success of the family, including genes encoding Cycloidea-like transcription factors, kinases, enzymes involved in rubber biosynthesis and disease resistance proteins that are expanded in the genome. We characterize 21 novel microRNAs, one of which may trigger phasiRNAs from numerous kinase transcripts. We provide evidence for a whole-genome triplication event specific but basal to the Compositae. We detect 26% of the genome in triplicated regions containing 30% of all genes that are enriched for regulatory sequences and depleted for genes involved in defence.

Consequences of normalizing transcriptomic and genomic libraries of plant genomes using a duplex-specific nuclease and tetramethylammonium chloride.

Several applications of high throughput genome and transcriptome sequencing would benefit from a reduction of the high-copy-number sequences in the libraries being sequenced and analyzed, particularly when applied to species with large genomes. We adapted and analyzed the consequences of a method that utilizes a thermostable duplex-specific nuclease for reducing the high-copy components in transcriptomic and genomic libraries prior to sequencing. This reduces the time, cost, and computational effort of obtaining informative transcriptomic and genomic sequence data for both fully sequenced and non-sequenced genomes. It also reduces contamination from organellar DNA in preparations of nuclear DNA. Hybridization in the presence of 3 M tetramethylammonium chloride (TMAC), which equalizes the rates of hybridization of GC and AT nucleotide pairs, reduced the bias against sequences with high GC content. Consequences of this method on the reduction of high-copy and enrichment of low-copy sequences are reported for Arabidopsis and lettuce.

Analyses of synteny between Arabidopsis thaliana and species in the Asteraceae reveal a complex network of small syntenic segments and major chromosomal rearrangements.

Comparative genomic studies among highly divergent species have been problematic because reduced gene similarities make orthologous gene pairs difficult to identify and because colinearity is expected to be low with greater time since divergence from the last common ancestor. Nevertheless, synteny between divergent taxa in several lineages has been detected over short chromosomal segments. We have examined the level of synteny between the model species Arabidopsis thaliana and species in the Compositae, one of the largest and most diverse plant families. While macrosyntenic patterns covering large segments of the chromosomes are not evident, significant levels of local synteny are detected at a fine scale covering segments of 1-Mb regions of A. thaliana and regions of <5 cM in lettuce and sunflower. These syntenic patches are often not colinear, however, and form a network of regions that have likely evolved by duplications followed by differential gene loss.

Dm3 is one member of a large constitutively expressed family of nucleotide binding site-leucine-rich repeat encoding genes.

The major cluster of resistance genes in lettuce cv. Diana contains approximately 32 nucleotide binding site-leucine-rich repeat encoding genes. Previous molecular dissection of this complex region had identified a large gene, RGC2B, as a candidate for encoding the downy mildew resistance gene, Dm3. This article describes genetic and transgenic complementation data that demonstrated RGC2B is necessary and sufficient to confer resistance with Dm3 specificity. Ethylmethanesulphonate was used to induce mutations to downy mildew susceptibility in cv. Diana (Dm1, Dm3, Dm7, and Dm8). Nineteen families were identified with a complete loss of resistance in one of the four resistance specificities. Sequencing revealed a variety of point mutations in RGC2B in the six dm3 mutants. Losses of resistance were due to single changes in amino acid sequence or a change in an intron splice site. These mutations did not cluster in any particular region of RGC2B. A full-length genomic copy of RGC2B was isolated from a lambdaphage library and introduced into two genotypes of lettuce. Transgenics expressing RGC2B exhibited resistance to all isolates expressing Avr3 from a wide range of geographical origins. In a wildtype Dm3-expressing genotype, many of the RGC2 family members are expressed at low levels throughout the plant.