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).

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.

The genetic basis of water-use efficiency and yield in lettuce

BACKGROUND

Water supply limits agricultural productivity of many crops including lettuce. Identifying cultivars within crop species that can maintain productivity with reduced water supply is a significant challenge, but central to developing resilient crops for future water-limited climates. We investigated traits known to be related to water-use efficiency (WUE) and yield in lettuce, a globally important leafy salad crop, in a recombinant inbred line (RIL) lettuce mapping population, produced from a cross between the cultivated Lactuca sativa L. cv. Salinas and its wild progenitor L. serriola L.

RESULTS

Wild and cultivated lettuce differed in their WUE and we observed transgressive segregation in yield and water-use traits in the RILs. Quantitative trait loci (QTL) analysis identified genomic regions controlling these traits under well-watered and droughted conditions. QTL were detected for carbon isotope discrimination, transpiration, stomatal conductance, leaf temperature and yield, controlling 4-23 % of the phenotypic variation. A QTL hotspot was identified on chromosome 8 that controlled carbon isotope discrimination, stomatal conductance and yield under drought. Several promising candidate genes in this region were associated with WUE, including aquaporins, late embryogenesis abundant proteins, an abscisic acid-responsive element binding protein and glutathione S-transferases involved in redox homeostasis following drought stress were also identified.

CONCLUSIONS

For the first time, we have characterised the genetic basis of WUE of lettuce, a commercially important and water demanding crop. We have identified promising candidate genomic regions determining WUE and yield under well-watered and water-limiting conditions, providing important pre-breeding data for future lettuce selection and breeding where water productivity will be a key target.

Genetics of Partial Resistance Against Verticillium dahliae Race 2 in Wild and Cultivated Lettuce

Lettuce (Lactuca sativa) is one of the most economically important vegetables in the United States, with approximately 50% of the domestic production concentrated in the Salinas Valley of California. Verticillium wilt, caused by races 1 and 2 of the fungal pathogen Verticillium dahliae, poses a major threat to lettuce production in this area. Although resistance governed by a single dominant gene against race 1 has previously been identified and is currently being incorporated into commercial cultivars, identification of resistance against race 2 has been challenging and no lines with complete resistance have been identified. In this study, we screened germplasm for resistance and investigated the genetics of partial resistance against race 2 using three mapping populations derived from crosses involving L. sativa × L. sativa and L. serriola × L. sativa. The inheritance of resistance in Lactuca species against race 2 is complex but a common quantitative trait locus (QTL) on linkage group 6, designated qVERT6.1 (quantitative Verticillium dahliae resistance on LG 6, first QTL), was detected in multiple populations. Additional race 2 resistance QTLs located in several linkage groups were detected in individual populations and environments. Because resistance in lettuce against race 2 is polygenic with a large genotype by environment interaction, breeding programs to incorporate these resistance genes should be aware of this complexity as they implement strategies to control race 2.

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.

The genetics of resistance to lettuce drop (Sclerotinia spp.) in lettuce in a recombinant inbred line population from Reine des Glaces × Eruption

KEY MESSAGE

Two QTLs for resistance to lettuce drop, qLDR1.1 and qLDR5.1, were identified. Associated SNPs will be useful in breeding for lettuce drop and provide the foundation for future molecular analysis. Lettuce drop, caused by Sclerotinia minor and S. sclerotiorum, is an economically important disease of lettuce. The association of resistance to lettuce drop with the commercially undesirable trait of fast bolting has hindered the integration of host resistance in control of this disease. Eruption is a slow-bolting cultivar that exhibits a high level of resistance to lettuce drop. Eruption also is completely resistant to Verticillium wilt caused by race 1 of Verticillium dahliae. A recombinant inbred line population from the cross Reine des Glaces × Eruption was genotyped by sequencing and evaluated for lettuce drop and bolting in separate fields infested with either S. minor or V. dahliae. Two quantitative trait loci (QTLs) for lettuce drop resistance were consistently detected in at least two experiments, and two other QTLs were identified in another experiment; the alleles for resistance at all four QTLs originated from Eruption. A QTL for lettuce drop resistance on linkage group (LG) 5, qLDR5.1, was consistently detected in all experiments and explained 11 to 25% of phenotypic variation. On LG1, qLDR1.1 was detected in two experiments explaining 9 to 12% of the phenotypic variation. Three out of four resistance QTLs are distinct from QTLs for bolting; qLDR5.1 is pleiotropic or closely linked with a QTL for early bolting; however, the rate of bolting shows only a small effect on the variance in resistance observed at this locus. The SNP markers linked with these QTLs will be useful in breeding for resistance through marker-assisted selection.

Genetic architecture of tipburn resistance in lettuce

KEY MESSAGE

Two major QTLs for tipburn were identified in LGs 1 and 5 contributing to resistance in cv. Salinas. The findings suggest pleiotropic effects between leaf crinkliness/savoy and tipburn. Tipburn is a physiological disorder in lettuce that is thought to be caused by a localized deficiency of calcium in leaf tissues. To elucidate the genetic architecture of resistance to tipburn in lettuce, seven recombinant inbred line populations were analyzed in multiple environments and years to identify quantitative trait loci (QTLs) for tipburn. Core height, head firmness, head closure, leaf crinkliness, plant fresh weight, and leaf savoy were also analyzed to investigate whether QTLs for these morphological traits collocated with QTLs for tipburn, which would be indicative of pleiotropic effects. Twenty-three major, intermediate, and minor unique QTLs for tipburn were identified in one or more populations scattered throughout the genome. Two major QTLs for tipburn incidence were identified in linkage groups (LGs) 1 and 5, which determined up to 45 and 66% of the phenotypic variance. The major QTL in LG 1 collocated with the head firmness QTL. The major QTL in LG 5 collocated with the QTL for core height, leaf crinkliness, and head firmness. Further research is needed to determine whether these associations are due to pleiotropic effects of the same gene or if the genes determining these traits are tightly linked. The beneficial alleles at the QTLs in LGs 1 and 5 are present in Lactuca sativa cv. Salinas, the genotype sequenced for the reference genome assembly. Therefore, these QTLs are good targets to identify genes causing tipburn as well as regions for marker-assisted selection to improve resistance to tipburn in lettuce.

Genome-wide functional analyses of plant coiled-coil NLR-type pathogen receptors reveal essential roles of their N-terminal domain in oligomerization, networking, and immunity

The ability to induce a defense response after pathogen attack is a critical feature of the immune system of any organism. Nucleotide-binding leucine-rich repeat receptors (NLRs) are key players in this process and perceive the occurrence of nonself-activities or foreign molecules. In plants, coevolution with a variety of pests and pathogens has resulted in repertoires of several hundred diverse NLRs in single individuals and many more in populations as a whole. However, the mechanism by which defense signaling is triggered by these NLRs in plants is poorly understood. Here, we show that upon pathogen perception, NLRs use their N-terminal domains to transactivate other receptors. Their N-terminal domains homo- and heterodimerize, suggesting that plant NLRs oligomerize upon activation, similar to the vertebrate NLRs; however, consistent with their large number in plants, the complexes are highly heterometric. Also, in contrast to metazoan NLRs, the N-terminus, rather than their centrally located nucleotide-binding (NB) domain, can mediate initial partner selection. The highly redundant network of NLR interactions in plants is proposed to provide resilience to perturbation by pathogens.

Elucidating the genetic basis of antioxidant status in lettuce (Lactuca sativa)

A diet rich in phytonutrients from fruit and vegetables has been acknowledged to afford protection against a range of human diseases, but many of the most popular vegetables are low in phytonutrients. Wild relatives of crops may contain allelic variation for genes determining the concentrations of these beneficial phytonutrients, and therefore understanding the genetic basis of this variation is important for breeding efforts to enhance nutritional quality. In this study, lettuce recombinant inbred lines, generated from a cross between wild and cultivated lettuce (Lactuca serriola and Lactuca sativa, respectively), were analysed for antioxidant (AO) potential and important phytonutrients including carotenoids, chlorophyll and phenolic compounds. When grown in two environments, 96 quantitative trait loci (QTL) were identified for these nutritional traits: 4 for AO potential, 2 for carotenoid content, 3 for total chlorophyll content and 87 for individual phenolic compounds (two per compound on average). Most often, the L. serriola alleles conferred an increase in total AOs and metabolites. Candidate genes underlying these QTL were identified by BLASTn searches; in several cases, these had functions suggesting involvement in phytonutrient biosynthetic pathways. Analysis of a QTL on linkage group 3, which accounted for >30% of the variation in AO potential, revealed several candidate genes encoding multiple MYB transcription factors which regulate flavonoid biosynthesis and flavanone 3-hydroxylase, an enzyme involved in the biosynthesis of the flavonoids quercetin and kaempferol, which are known to have powerful AO activity. Follow-up quantitative RT-PCR of these candidates revealed that 5 out of 10 genes investigated were significantly differentially expressed between the wild and cultivated parents, providing further evidence of their potential involvement in determining the contrasting phenotypes. These results offer exciting opportunities to improve the nutritional content and health benefits of lettuce through marker-assisted breeding.

Quantitative trait loci associated with tipburn, heat stress-induced physiological disorders, and maturity traits in crisphead lettuce

Crisphead lettuce (Lactuca sativa L.) crops exhibit several economically important, physiological disorders when grown in high temperature conditions. These include tipburn, rib discoloration, premature bolting, ribbiness, and internal rib cracking. We evaluated seven physiological disorders and three agronomic traits segregating in a recombinant inbred line (RIL) population consisting of 152 F7 RILs derived from an intra-specific cross between two crisphead cultivars, L. sativa cv. Emperor x L. sativa cv. El Dorado; evaluations were carried out at each of two parental maturities in one planting and at one intermediate maturity in a second planting in each of 2 years for a total of six evaluations. A genetic map was developed using 449 polymorphic SNP markers; it comprises 807 cM in 20 linkage groups that covered 51 % of the nine lettuce chromosomes. Composite interval mapping revealed a total of 36 significant QTLs for eight out of the ten traits evaluated. Significant QTLs were distributed in 11 linkage groups on seven of the chromosomes and accounted for up to 83 % of the phenotypic variation observed. The three largest QTLs for rib discoloration, which accounted individually for 7-21 % of the variation, were clustered with stem length, two with ribbiness and one with head firmness. Three major clusters of QTLs revealed pleiotropic effects or tight linkage between tipburn incidence and severity, head type, stem length, head firmness and ribbiness. One QTL, qTPB5.2, was detected in multiple trials and described 38-70 % of the variation in tipburn incidence. qTPB5.2 is, therefore, a useful candidate gene for breeding for tipburn resistance using marker-assisted selection.

An Ultra-High-Density, Transcript-Based, Genetic Map of Lettuce

We have generated an ultra-high-density genetic map for lettuce, an economically important member of the Compositae, consisting of 12,842 unigenes (13,943 markers) mapped in 3696 genetic bins distributed over nine chromosomal linkage groups. Genomic DNA was hybridized to a custom Affymetrix oligonucleotide array containing 6.4 million features representing 35,628 unigenes of Lactuca spp. Segregation of single-position polymorphisms was analyzed using 213 F7:8 recombinant inbred lines that had been generated by crossing cultivated Lactuca sativa cv. Salinas and L. serriola acc. US96UC23, the wild progenitor species of L. sativa The high level of replication of each allele in the recombinant inbred lines was exploited to identify single-position polymorphisms that were assigned to parental haplotypes. Marker information has been made available using GBrowse to facilitate access to the map. This map has been anchored to the previously published integrated map of lettuce providing candidate genes for multiple phenotypes. The high density of markers achieved in this ultradense map allowed syntenic studies between lettuce and Vitis vinifera as well as other plant species.

The genomic architecture of disease resistance in lettuce

Genbank and The Compositae Genome Project database, containing over 42,000 lettuce unigenes from Lactuca sativa cv. Salinas and L. serriola accession UC96US23 were mined to identify 702 candidate genes involved in pathogen recognition (RGCs), resistance signal transduction, defense responses, and disease susceptibility. In addition, to identify sequences representing additional sub-families of nucleotide binding site (NBS)-leucine-rich repeat encoding genes; the major classes of resistance genes (R-genes), NBS-encoding sequences were amplified by PCR using degenerate oligonucleotides designed to NBS sub-families specific to the subclass Asteridae, which includes the Compositae family. These products were cloned and sequenced resulting in 18 novel NBS sequences from cv. Salinas and 15 novel NBS sequences from UC96US23. Using a variety of marker technologies, 294 of the 735 candidate disease resistance genes were mapped in our primary mapping population, which consisted of 119 F7 recombinant inbred lines derived from an interspecific cross between cv. Salinas and UC96US23. Using markers shared across multiple genetic maps, 36 resistance phenotypic loci, including two new loci for resistance to downy mildew and two quantitative trait loci for resistance to anthracnose were positioned onto the reference map to provide a global view of the genomic architecture of disease resistance in lettuce and to identify candidate genes for resistance phenotypes. The majority but not all of the resistance phenotypes were genetically associated with RGCs.

The genomic architecture of disease resistance in lettuce

Genbank and The Compositae Genome Project database, containing over 42,000 lettuce unigenes from Lactuca sativa cv. Salinas and L. serriola accession UC96US23 were mined to identify 702 candidate genes involved in pathogen recognition (RGCs), resistance signal transduction, defense responses, and disease susceptibility. In addition, to identify sequences representing additional sub-families of nucleotide binding site (NBS)-leucine-rich repeat encoding genes; the major classes of resistance genes (R-genes), NBS-encoding sequences were amplified by PCR using degenerate oligonucleotides designed to NBS sub-families specific to the subclass Asteridae, which includes the Compositae family. These products were cloned and sequenced resulting in 18 novel NBS sequences from cv. Salinas and 15 novel NBS sequences from UC96US23. Using a variety of marker technologies, 294 of the 735 candidate disease resistance genes were mapped in our primary mapping population, which consisted of 119 F7 recombinant inbred lines derived from an interspecific cross between cv. Salinas and UC96US23. Using markers shared across multiple genetic maps, 36 resistance phenotypic loci, including two new loci for resistance to downy mildew and two quantitative trait loci for resistance to anthracnose were positioned onto the reference map to provide a global view of the genomic architecture of disease resistance in lettuce and to identify candidate genes for resistance phenotypes. The majority but not all of the resistance phenotypes were genetically associated with RGCs.

A high-density, integrated genetic linkage map of lettuce (Lactuca spp.)

An integrated map for lettuce comprising of 2,744 markers was developed from seven intra- and inter-specific mapping populations. A total of 560 markers that segregated in two or more populations were used to align the individual maps. 2,073 AFLP, 152 RFLP, 130 SSR, and 360 RAPD as well as 29 other markers were assigned to nine chromosomal linkage groups that spanned a total of 1,505 cM and ranged from 136 to 238 cM. The maximum interval between markers in the integrated map is 43 cM and the mean interval is 0.7 cM. The majority of markers segregated close to Mendelian expectations in the intra-specific crosses. In the two L. saligna x L. sativa inter-specific crosses, a total of 155 and 116 markers in 13 regions exhibited significant segregation distortion. Data visualization tools were developed to curate, display and query the data. The integrated map provides a framework for mapping ESTs in one core mapping population relative to phenotypes that segregate in other populations. It also provides large numbers of markers for marker assisted selection, candidate gene identification, and studies of genome evolution in the Compositae.

Inter- and intra-genomic homology of the Brassica genomes: implications for their origin and evolution

In order to determine the homologous regions shared by the cultivated Brassica genomes, linkage maps of the diploid cultivated B. rapa (A genome, n = 10), B. nigra (B genome, n = 8) and B. oleracea (C genome, n = 9), were compared. We found intergenomic conserved regions but with extensitve reordering among the genomes. Eighteen linkage groups from all three species could be associated on the basis of homologous segments based on at least three common markers. Intragenomic homologous conservation was also observed for some of the chromosomes of the A, B and C genomes. A possible chromosome phylogenetic pathway based on an ancestral genome of at least five, and no more than seven chromosomes, was drawn from the chromosomal inter-relationships observed. These results demonstrate that extensive duplication and rearrangement have been involved in the formation of the Brassica genomes from a smaller ancestral genome.

Sequence comparison of two codominant RAPD markers in Brassica nigra: deletions, substitutions and microsatellites

Two RAPD fragments segregating codominantly were investigated in a F2 population of Brassica nigra. Southern hybridization of these DNA fragments to genomic B. nigra DNA digested with several endonucleases revealed similar restriction profiles. Sequencing of the two fragments disclosed 93% homology. The differences were due mainly to an internal 41 nucleotide deletion in one of the fragments. Minor deletions of one to three bases, including a microsatellite of CTT motif were also observed. In addition, base substitutions, mostly transitions were detected. These relatively small differences suggested that the two RAPD products were indeed different versions of the same sequence. The larger fragment of 1154 bp was denominated A05(1) whereas the shorter one, denominated A05(2), had 1116 bp.

Structure and organization of the B genome based on a linkage map in Brassica nigra

We constructed a genetic map on Brassica nigra based on a segregating population of 83 F2 individuals. Three different types of molecular markers were used to build the map including isozymes, restriction fragment length polymorphisms (RFLP), and random amplified polymorphic DNA (RAPD). The final map contained 124 markers distributed in 11 linkage groups. The map covered a total distance of 677 cM with the markers distributed within a mean distance of 5.5cM. Of the sequences found in the B. nigra map, 40% were duplicated and organized into three different types of arrangements. They were either scattered throughout the genome, organized in tandem, or organized in blocks of duplicated loci conserved in more than 1 linkage group.