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.

Ancestral Chromosomes for Family Peronosporaceae Inferred from a Telomere-to-Telomere Genome Assembly of Peronospora effusa.

Downy mildew disease of spinach, caused by the oomycete Peronospora effusa, causes major losses to spinach production. In this study, the 17 chromosomes of P. effusa were assembled telomere-to-telomere, using Pacific Biosciences high-fidelity reads. Of these, 16 chromosomes are complete and gapless; chromosome 15 contains one gap bridging the nucleolus organizer region. This is the first telomere-to-telomere genome assembly for an oomycete. Putative centromeric regions were identified on all chromosomes. This new assembly enables a reevaluation of the genomic composition of Peronospora spp.; the assembly was almost double the size and contained more repeat sequences than previously reported for any Peronospora species. Genome fragments consistently underrepresented in six previously reported assemblies of P. effusa typically encoded repeats. Some genes annotated as encoding effectors were organized into multigene clusters on several chromosomes. Putative effectors were annotated on 16 of the 17 chromosomes. The intergenic distances between annotated genes were consistent with compartmentalization of the genome into gene-dense and gene-sparse regions. Genes encoding putative effectors were enriched in gene-sparse regions. The near-gapless assembly revealed apparent horizontal gene transfer from Ascomycete fungi. Gene order was highly conserved between P. effusa and the genetically oriented assembly of the oomycete Bremia lactucae; high levels of synteny were also detected with Phytophthora sojae. Extensive synteny between phylogenetically distant species suggests that many other oomycete species may have similar chromosome organization. Therefore, this assembly provides the foundation for genomic analyses of diverse oomycetes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Pseudomonas syringae effector HopZ3 suppresses the bacterial AvrPto1-tomato PTO immune complex via acetylation.

The plant pathogen Pseudomonas syringae secretes multiple effectors that modulate plant defenses. Some effectors trigger defenses due to specific recognition by plant immune complexes, whereas others can suppress the resulting immune responses. The HopZ3 effector of P. syringae pv. syringae B728a (PsyB728a) is an acetyltransferase that modifies not only components of plant immune complexes, but also the Psy effectors that activate these complexes. In Arabidopsis, HopZ3 acetylates the host RPM1 complex and the Psy effectors AvrRpm1 and AvrB3. This study focuses on the role of HopZ3 during tomato infection. In Psy-resistant tomato, the main immune complex includes PRF and PTO, a RIPK-family kinase that recognizes the AvrPto effector. HopZ3 acts as a virulence factor on tomato by suppressing AvrPto1Psy-triggered immunity. HopZ3 acetylates AvrPto1Psy and the host proteins PTO, SlRIPK and SlRIN4s. Biochemical reconstruction and site-directed mutagenesis experiments suggest that acetylation acts in multiple ways to suppress immune signaling in tomato. First, acetylation disrupts the critical AvrPto1Psy-PTO interaction needed to initiate the immune response. Unmodified residues at the binding interface of both proteins and at other residues needed for binding are acetylated. Second, acetylation occurs at residues important for AvrPto1Psy function but not for binding to PTO. Finally, acetylation reduces specific phosphorylations needed for promoting the immune-inducing activity of HopZ3's targets such as AvrPto1Psy and PTO. In some cases, acetylation competes with phosphorylation. HopZ3-mediated acetylation suppresses the kinase activity of SlRIPK and the phosphorylation of its SlRIN4 substrate previously implicated in PTO-signaling. Thus, HopZ3 disrupts the functions of multiple immune components and the effectors that trigger them, leading to increased susceptibility to infection. Finally, mass spectrometry used to map specific acetylated residues confirmed HopZ3's unusual capacity to modify histidine in addition to serine, threonine and lysine residues.

Fast Fluorescence Titration Quantification of Plasmid DNA with DNA Attractive Magnetic Nanoparticles.

Fluorescence titration using magnetic nanoparticles (FTMN) was performed as a rapid, inexpensive, and simple method for quantifying the amount of fluorophore-intercalated plasmid DNA on these DNA attractive nanoparticles. Binding of the propidium iodide (PI)-intercalated DNA (PI/DNA) to polyethylenimine (PEI)-coated monodisperse iron oxide magnetic nanoparticles (PEI-MNs) was confirmed with transmission electron microscopy after the two species were mixed in water for less than a minute. The amount of DNA on PEI-MNs in aqueous solution, however, could not be easily determined using direct fluorescence measurements due to strong scattering by aggregated MNs, especially at high nanoparticle concentrations. Instead, fluorescence measurements were taken immediately after the solution of PI/DNA and PEI-MN mixtures was treated with a magnet to pull the PEI-MNs out of the solution. The detected fluorescence signal of the remaining free PI/DNA in the solution decreased as the concentration of PEI-MNs in the pre-treated solutions increased, resulting in a titration curve, which was used to determine the amount of DNA on MNs, the dissociation constant, and binding energy after the concentration of PEI-MNs was calibrated with microwave-plasma atomic emission spectroscopy. Quantitative polymerase chain reaction was used to understand the binding of DNA to MNs and to measure the amount of free PI/DNA in solution, and the results were similar to those obtained with the FTMN method.

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.

Isolating an active and inactive CACTA transposon from lettuce color mutants and characterizing their family.

Dietary flavonoids play an important role in human nutrition and health. Flavonoid biosynthesis genes have recently been identified in lettuce (Lactuca sativa); however, few mutants have been characterized. We now report the causative mutations in Green Super Lettuce (GSL), a natural light green mutant derived from red cultivar NAR; and GSL-Dark Green (GSL-DG), an olive-green natural derivative of GSL. GSL harbors CACTA 1 (LsC1), a 3.9-kb active nonautonomous CACTA superfamily transposon inserted in the 5' untranslated region of anthocyanidin synthase (ANS), a gene coding for a key enzyme in anthocyanin biosynthesis. Both terminal inverted repeats (TIRs) of this transposon were intact, enabling somatic excision of the mobile element, which led to the restoration of ANS expression and the accumulation of red anthocyanins in sectors on otherwise green leaves. GSL-DG harbors CACTA 2 (LsC2), a 1.1-kb truncated copy of LsC1 that lacks one of the TIRs, rendering the transposon inactive. RNA-sequencing and reverse transcription quantitative PCR of NAR, GSL, and GSL-DG indicated the relative expression level of ANS was strongly influenced by the transposon insertions. Analysis of flavonoid content indicated leaf cyanidin levels correlated positively with ANS expression. Bioinformatic analysis of the cv Salinas lettuce reference genome led to the discovery and characterization of an LsC1 transposon family with a putative transposon copy number greater than 1,700. Homologs of tnpA and tnpD, the genes encoding two proteins necessary for activation of transposition of CACTA elements, were also identified in the lettuce genome.

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.

High-resolution genetic dissection of the major QTL for tipburn resistance in lettuce, Lactuca sativa.

Tipburn is an important physiological disorder of lettuce, Lactuca sativa L., related to calcium deficiency that can result in leaf necrosis and unmarketable crops. The major quantitative trait locus (QTL), qTPB5.2, can account for up to 70% of the phenotypic variance for tipburn incidence in the field. This QTL was genetically dissected to identify candidate genes for tipburn by creating lines with recombination events within the QTL and assessing their resistance to tipburn. By comparing lines with contrasting haplotypes, the genetic region was narrowed down to ∼877 Kb that was associated with a reduction of tipburn by ∼60%. Analysis of the lettuce reference genome sequence revealed 12 genes in this region, one of which is a calcium transporter with a single nucleotide polymorphism in an exon between haplotypes with contrasting phenotypes. RNA-seq analysis of recombinants revealed two genes that were differentially expressed between contrasting haplotypes consistent with the tipburn phenotype. One encodes a Teosinte branched1/Cycloidea/Proliferating Cell factor transcription factor; however, differential expression of the calcium transporter was detected. The phenotypic data indicated that there is a second region outside of the ∼877 Kb region but within the QTL, at which a haplotype from the susceptible parent decreased tipburn by 10-20%. A recombinant line was identified with beneficial haplotypes in each region from both parents that showed greater tipburn resistance than the resistant parent; this line could be used as the foundation for breeding cultivars with more resistance than is currently available.

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.

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.

Identification of Major Quantitative Trait Loci Controlling Field Resistance to Downy Mildew in Cultivated Lettuce (Lactuca sativa).

Lettuce downy mildew, caused by Bremia lactucae Regel, is the most economically important foliar disease of lettuce (Lactuca sativa L.). The deployment of resistant cultivars carrying dominant resistance genes (Dm genes) plays a crucial role in integrated downy mildew disease management; however, high variability in pathogen populations leads to the defeat of plant resistance conferred by Dm genes. Some lettuce cultivars exhibit field resistance that is only manifested in adult plants. Two populations of recombinant inbred lines (RILs), originating from crosses between the field resistant cultivars Grand Rapids and Iceberg and susceptible cultivars Salinas and PI491224, were evaluated for downy mildew resistance under field conditions. In all, 160 RILs from the Iceberg × PI491224 and 88 RILs from the Grand Rapids × Salinas populations were genotyped using genotyping by sequencing, which generated 906 and 746 high-quality markers, respectively, that were used for quantitative trait locus (QTL) analysis. We found a QTL in chromosome 4 that is present in both Grand Rapids × Salinas and Iceberg × PI491224 populations that has a major effect on field resistance. We also found two additional significant QTLs in chromosomes 2 and 5 in the Iceberg × PI491224 RIL population. Marker-assisted gene pyramiding of multiple Dm genes in combination with QTLs for field resistance provide the opportunity to develop cultivars with more durable resistance to B. lactucae.

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.

Phytopathogen Effectors Use Multiple Mechanisms to Manipulate Plant Autophagy.

Autophagy is a central part of immunity and hence is a key target of pathogens. However, the precise molecular mechanisms by which plant pathogens manipulate autophagy remain elusive. We identify a network of 88 interactions between 184 effectors from bacterial, fungal, oomycete, and nematode pathogens with 25 Arabidopsis autophagy (ATG) proteins. Notably, Pseudomonas syringae pv tomato (Pto) bacterial effectors HrpZ1, HopF3, and AvrPtoB employ distinct molecular strategies to modulate autophagy. Calcium-dependent HrpZ1 oligomerization targets ATG4b-mediated cleavage of ATG8 to enhance autophagy, while HopF3 also targets ATG8 but suppresses autophagy, with both effectors promoting infection. AvrPtoB affects ATG1 kinase phosphorylation and enhances bacterial virulence. Since pathogens inject limited numbers of effectors into hosts, our findings establish autophagy as a key target during infection. Additionally, as autophagy is enhanced and inhibited by these effectors, autophagy likely has different functions throughout infection and, thus, must be temporally and precisely regulated for successful infection.

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.

Linked-read sequencing of gametes allows efficient genome-wide analysis of meiotic recombination.

Meiotic crossovers (COs) ensure proper chromosome segregation and redistribute the genetic variation that is transmitted to the next generation. Large populations and the demand for genome-wide, fine-scale resolution challenge existing methods for CO identification. Taking advantage of linked-read sequencing, we develop a highly efficient method for genome-wide identification of COs at kilobase resolution in pooled recombinants. We first test this method using a pool of Arabidopsis F2 recombinants, and recapitulate results obtained from the same plants using individual whole-genome sequencing. By applying this method to a pool of pollen DNA from an F1 plant, we establish a highly accurate CO landscape without generating or sequencing a single recombinant plant. The simplicity of this approach enables the simultaneous generation and analysis of multiple CO landscapes, accelerating the pace at which mechanisms for the regulation of recombination can be elucidated through efficient comparisons of genotypic and environmental effects on recombination.

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.