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.

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.

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.

At least 10 Dm genes conferring resistance to the oomycete downy mildew fungus Bremia lactucae map to the major resistance cluster in lettuce. We investigated the structure of this cluster in the lettuce cultivar Diana, which contains Dm3. A deletion breakpoint map of the chromosomal region flanking Dm3 was saturated with a variety of molecular markers. Several of these markers are components of a family of resistance gene candidates ( RGC2 ) that encode a nucleotide binding site and a leucine-rich repeat region. These motifs are characteristic of plant disease resistance genes. Bacterial artificial chromosome clones were identified by using duplicated restriction fragment length polymorphism markers from the region, including the nucleotide binding site-encoding region of RGC2. Twenty-two distinct members of the RGC2 family were characterized from the bacterial artificial chromosomes; at least two additional family members exist. The RGC2 family is highly divergent; the nucleotide identity was as low as 53% between the most distantly related copies. These RGC2 genes span at least 3.5 Mb. Eighteen members were mapped on the deletion breakpoint map. A comparison between the phylogenetic and physical relationships of these sequences demonstrated that closely related copies are physically separated from one another and indicated that complex rearrangements have shaped this region. Analysis of low-copy genomic sequences detected no genes, including RGC2 , in the Dm3 region, other than sequences related to retrotransposons and transposable elements. The related but divergent family of RGC2 genes may act as a resource for the generation of new resistance phenotypes through infrequent recombination or unequal crossing over. INTRODUCTIONDisease resistance genes frequently occur in tightly linked clusters (Pryor, 1987; Crute and Pink, 1996;Michelmore and Meyers, 1998). Clusters of plant resistance genes were first established by use of classic genetic techniques; detailed molecular analyses are now beginning to unravel the complexity of these loci and the underlying mechanisms determining their structure (Parniske et al., 1997;Song et al., 1997). It is becoming increasingly apparent that such clusters may be both common and complex genomic regions in plants.Clusters of resistance genes have been identified in diverse plant species. More than 30 different resistance specificities to the single fungal pathogen responsible for flax rust disease, Melampsora lini , have been mapped to five linkage groups (Flor, 1971;Islam and Shepherd, 1991). These loci exemplify two possible genetic arrangements that may exist for clusters of resistance genes: the flax L locus contains at least 13 allelic rust resistance specificities, and the more complex M locus exists as a tandem array of at least seven genes (Islam and Shepherd, 1991). In maize, multiple Rp genes, both linked and allelic, have been observed to mediate resistance to the rust fungus Puccinia sorghi ; 16 genetically separable loci were mapped to a single cluster known as th...

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.